Skip to main content

Program for estimating amino acid replacement rates

Project description

modelestimator --- Infer sequence evolution rate matrices from a MSA

Example usage

    modelestimator -t 0.001 file1.fa file2.fa file3.fa

Infer a rate matrix (written to stdout) from three alignment files in Fasta format. The output is PAML format by default, and therefore applicable in a number of available phylogenetic softwares.

    modelestimator -b 200 file.fa

Try the experimental bootstrapping feature (200 replicates) on a Fasta multialignment.


modelestimator <options> infiles

<format> should be either FASTA, STOCKHOLM or PHYLIP format.

Output is a rate matrix and residue distribution vector.


  -h, --help            show this help message and exit
  --version             show program's version number and exit
  -a {iqtree,matlab,mrbayes,octave,paml,phyml}, --application {iqtree,matlab,mrbayes,octave,paml,phyml}
                        Choose output format to suit the application you want
                        to use for inference. The 'iqtree', 'paml' and 'phyml'
                        options are identical. The 'matlab' and 'octave'
                        optins are for import into MatLab-compatible programs
                        and are presenting the actual Q matrix rather than the
                        R matrix used by PAML/PhyML, etc. Default: paml
  -f {fasta,clustal,nexus,phylip,stockholm}, --format {fasta,clustal,nexus,phylip,stockholm}
                        Specify sequence format of input files. Default: fasta
  -t T, --threshold T   Stop when consecutive iterations do not change by more
                        than T. Default: 0.001
  -b N, --bootstrap N   Estimate the rate matrix using bootstrapping by
                        computing N resampled replicates of the input
                        multialignment. For each replicate, a rate matrix is
                        computed. The mean matrix the elementwise standard
                        deviation is returned. Only one infile should be given
                        in this mode.
  -B N, --bootstrapped_quality N
                        Estimate the quality of the rate matrix estimate using
                        a bootstrap procedure. The multialignment is resampled
                        N times and a Q matrix is computed for each replicate.
                        Then the difference (matrix norm) between rate matrix
                        estimated without resampling and each bootstrapped Q
                        is computed and the mean difference is returned. Only
                        one infile should be given in this mode. Returns
                        bootstrap norm.

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for modelestimator-v2, version 2.1.3
Filename, size File type Python version Upload date Hashes
Filename, size modelestimator_v2-2.1.3-py3-none-any.whl (23.4 kB) File type Wheel Python version py3 Upload date Hashes View
Filename, size modelestimator-v2-2.1.3.tar.gz (14.4 kB) File type Source Python version None Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page