MtreeHMM
Mutagenetic tree hidden Markov models
Current version
0.1.01 Jul 2006
Description
MtreeHMM is a Perl program that implements the Expectation Maximization algorithm described in (Beerenwinkel & Drton, 2007) for estimating the parameters of the mutagenetic tree hidden Markov model from longitudinal clonal sequence data.
Download
MtreeHMM is distributed under the external pageGNU General Public Licensecall_made.
- Downloadmtreehmm-0.1.01.tar.gz (GZ, 15 KB)vertical_align_bottom
(includes source code, toy example, and HIV data set)
Installation and Usage
In order to run MtreeHMM you need, besides the Perl interpreter, the Perl module 'Amoeba.pm' available at external pageCPANcall_made. At the beginning of the PERL script, change the line
use lib "$ENV{HOME}/perllib";
to point to the directory that contains 'Math/Amoeba.pm' (here '~/perllib').
The input data needs to be specified in four different files that share the same filestem:
filestem.model - Mutagenetic tree model, adjacency matrix
filestem.pat - List of binary mutational patterns (clones)
filestem.patient - List of patients
filestem.week - List of time points
The '.model' and '.pat' files are the same as in the external pageMtreemixcall_made package. You may use Mtreemix to estimate the tree structure ('.model' file) from the observed mutational patterns ('.pat' file).
After download and unpacking, type 'mtreehmm.pl -h' for detailed usage instructions. The download package includes a small test dataset named toy, from which the file formats become apparent.
The HIV data set analyzed in (Beerenwinkel & Drton, 2007) is also included.
Authors and Contributors
Citation
- A mutagenetic tree hidden Markov model for longitudinal clonal HIV sequence data
Niko Beerenwinkel and Mathias Drton
external pageBiostatistics 8 (1), 53-71call_made