NAME
Bio::Tools::Run::Phylo::PAML::Yn00 - Wrapper aroud the PAML program yn00
VERSION
version 1.7.3
SYNOPSIS
use Bio::Tools::Run::Phylo::PAML::Yn00;use Bio::AlignIO;my $alignio = Bio::AlignIO->new(-format => 'phylip',-file => 't/data/gf-s85.phylip');my $aln = $alignio->next_aln;my $yn = Bio::Tools::Run::Phylo::PAML::Yn00->new();$yn->alignment($aln);my ($rc,$parser) = $yn->run();while( my $result = $parser->next_result ) {my @otus = $result->get_seqs();my $MLmatrix = $result->get_MLmatrix();# 0 and 1 correspond to the 1st and 2nd entry in the @otus arraymy $dN = $MLmatrix->[0]->[1]->{dN};my $dS = $MLmatrix->[0]->[1]->{dS};my $kaks =$MLmatrix->[0]->[1]->{omega};print "Ka = $dN Ks = $dS Ka/Ks = $kaks\n";}
DESCRIPTION
This is a wrapper around the yn00 (method of Yang and Nielsen, 2000) program of PAML (Phylogenetic Analysis by Maximum Likelihood) package of Ziheng Yang. See http://abacus.gene.ucl.ac.uk/software/paml.html for more information.
This module will generate a proper yn00.ctl file and will run the program in a separate temporary directory to avoid creating temp files all over the place and will cleanup after itself.
Default Values
See the Bio::Tools::Run::Phylo::PAML::Codeml module for documentation of the default values.
ATTRIBUTES
program_name
Title : program_nameUsage : $yn00->program_name()Function: holds the program nameReturns: stringArgs : None
program_dir
Title : program_dirUsage : $yn00->program_dir()Function: returns the program directory, obtained from ENV variable.Returns: stringArgs :
error_string
Title : error_stringUsage : $obj->error_string($newval)Function: Where the output from the last analysus run is stored.Returns : value of error_stringArgs : newvalue (optional)
alignment
Title : alignmentUsage : $codeml->align($aln);Function: Get/Set the L<Bio::Align::AlignI> objectReturns : L<Bio::Align::AlignI> objectArgs : [optional] L<Bio::Align::AlignI>Comment : We could potentially add support for running directly on a filebut we shall keep it simpleSee also: L<Bio::SimpleAlign>
no_param_checks
Title : no_param_checksUsage : $obj->no_param_checks($newval)Function: Boolean flag as to whether or not we shouldtrust the sanity checks for parameter valuesReturns : value of no_param_checksArgs : newvalue (optional)
save_tempfiles
Title : save_tempfilesUsage : $obj->save_tempfiles($newval)Function:Returns : value of save_tempfilesArgs : newvalue (optional)
outfile_name
Title : outfile_nameUsage : my $outfile = $codeml->outfile_name();Function: Get/Set the name of the output file for this run(if you wanted to do something special)Returns : stringArgs : [optional] string to set value to
tempdir
Title : tempdirUsage : my $tmpdir = $self->tempdir();Function: Retrieve a temporary directory name (which is created)Returns : string which is the name of the temporary directoryArgs : none
METHODS
new
Title : newUsage : my $obj = Bio::Tools::Run::Phylo::PAML::Yn00->new();Function: Builds a new Bio::Tools::Run::Phylo::PAML::Yn00 objectReturns : Bio::Tools::Run::Phylo::PAML::Yn00Args : -alignment => the L<Bio::Align::AlignI> object-save_tempfiles => boolean to save the generated tempfiles andNOT cleanup after onesself (default FALSE)
run
Title : runUsage : $yn->run();Function: run the yn00 analysis using the default or updated parametersthe alignment parameter must have been setReturns : 3 values,$rc = 1 for success, 0 for errorshash reference of the Yang calculated Ka/Ks valuesthis is a set of pairwise observations keyed assequencenameA->sequencenameB->datatypehash reference same as the previous one except it for theNei and Gojobori calculated Ka,Ks,omega valuesArgs : none
get_parameters
Title : get_parametersUsage : my %params = $self->get_parameters();Function: returns the list of parameters as a hashReturns : associative array keyed on parameter namesArgs : none
set_parameter
Title : set_parameterUsage : $codeml->set_parameter($param,$val);Function: Sets a codeml parameter, will be validated againstthe valid values as set in the %VALIDVALUES class variable.The checks can be ignored if on turns of param checks like this:$codeml->no_param_checks(1)Returns : boolean if set was success, if verbose is set to -1then no warning will be reportedArgs : $paramname => name of the parameter$value => value to set the parameter toSee also: L<no_param_checks()>
set_default_parameters
Title : set_default_parametersUsage : $codeml->set_default_parameters(0);Function: (Re)set the default parameters from the defaults(the first value in each array in the%VALIDVALUES class variable)Returns : noneArgs : boolean: keep existing parameter values
cleanup
Title : cleanupUsage : $codeml->cleanup();Function: Will cleanup the tempdir directory after a PAML runReturns : noneArgs : none
io
Title : ioUsage : $obj->io($newval)Function: Gets a L<Bio::Root::IO> objectReturns : L<Bio::Root::IO>Args : none
Bio::Tools::Run::Wrapper methods
FEEDBACK
Mailing lists
User feedback is an integral part of the evolution of this and other Bioperl modules. Send your comments and suggestions preferably to the Bioperl mailing list. Your participation is much appreciated.
[email protected] - General discussionhttp://bioperl.org/Support.html - About the mailing lists
Support
Please direct usage questions or support issues to the mailing list: [email protected]
rather than to the module maintainer directly. Many experienced and reponsive experts will be able look at the problem and quickly address it. Please include a thorough description of the problem with code and data examples if at all possible.
Reporting bugs
Report bugs to the Bioperl bug tracking system to help us keep track of the bugs and their resolution. Bug reports can be submitted via the web:
https://github.com/bioperl/bio-tools-phylo-paml/issues
AUTHOR
Jason Stajich <[email protected]>
COPYRIGHT
This software is copyright (c) by Jason Stajich <[email protected]>.
This software is available under the same terms as the perl 5 programming language system itself.