OnWorks favicon

mcmc_analysis - Online in the Cloud

Run mcmc_analysis in OnWorks free hosting provider over Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator

This is the command mcmc_analysis that can be run in the OnWorks free hosting provider using one of our multiple free online workstations such as Ubuntu Online, Fedora Online, Windows online emulator or MAC OS online emulator



mcmc_analysis - Analysis of the MCMC


mcmc_analysis [OPTIONS] MCMC-data ...


This program reads MCMC output, from one or several files, from programs in the BEEP
package. It is important that the same columns are present in each input file.

Input format:
Output from MCMC iterations on the format
<likelihood> <tab> <iter> <tab> <params>
<likelihood> is the logarithm of the likelihood in float format
<tab> is tab whitespace separating the fields
<iter> is an integer for the ordinal of the iteration
<params> is a list of fields separated by semicolons containing
the parameters of the MCMC.

The MCMC params are typed and given names by the first line in the file,
which is on the format
# L <tab> N <tab> [<name>(<type>)]+
The names ought to be unique, but do not have to be. The <type> is one
and used to infer how to parse and analyze the rest of the lines.

Output format:
A report on the MCMC run, with posterior estimates of the parameters.


The percentage (0 <= x <1), or the number of (x is integer >= 1) of the input to be
discarded as burnin. Default: 0.1.

Plot the parameter named <string>. Output is two columns, the iteration number and
the parameter's value in the iteration.

Ignore the named parameter. See header line in MCMC output for column names. You
can name several columns at once using a comma-delimited format, (e.g. -i
Length,Name). No spaces allowed between column names.

-t Output LaTex for the analysis.

-l Count and report the number of samples in the input file.

-mp NT The maximum number of points to plot.

-sp Explicitly indicate points in plots.

-coda Output file for the CODA package in R

Same as -coda, but includes tree parameters. Each tree is output as an integer ID
(order of visitation in chain, 1,2,...). Make sure to use -i to hide any trees
containing lengths/times/rates.

-P Parallel chains. This implies that several files of (parallel) samples are listed.


The prime-phylo home page: http://prime.sbc.su.se

Use mcmc_analysis online using onworks.net services

Free Servers & Workstations

Download Windows & Linux apps

Linux commands