 Software
 Open Access
 Published:
Mocapy++  A toolkit for inference and learning in dynamic Bayesian networks
BMC Bioinformatics volume 11, Article number: 126 (2010)
Abstract
Background
Mocapy++ is a toolkit for parameter learning and inference in dynamic Bayesian networks (DBNs). It supports a wide range of DBN architectures and probability distributions, including distributions from directional statistics (the statistics of angles, directions and orientations).
Results
The program package is freely available under the GNU General Public Licence (GPL) from SourceForge http://sourceforge.net/projects/mocapy. The package contains the source for building the Mocapy++ library, several usage examples and the user manual.
Conclusions
Mocapy++ is especially suitable for constructing probabilistic models of biomolecular structure, due to its support for directional statistics. In particular, it supports the Kent distribution on the sphere and the bivariate von Mises distribution on the torus. These distributions have proven useful to formulate probabilistic models of protein and RNA structure in atomic detail.
Background
A Bayesian network (BN) represents a set of variables and their joint probability distribution using a directed acyclic graph [1, 2]. A dynamic Bayesian network (DBN) is a BN that represents sequences, such as timeseries from speech data or biological sequences [3]. One of the simplest examples of a DBN is the well known hidden Markov model (HMM) [4, 5]. DBNs have been applied with great success to a large number of problems in various fields. In bioinformatics, DBNs are especially relevant because of the sequential nature of biological molecules, and have therefore proven suitable for tackling a large number of problems. Examples are protein homologue detection [6], protein secondary structure prediction [7, 8], gene finding [5], multiple sequence alignment [5] and sampling of protein conformations [9, 10].
Here, we present a general, open source toolkit, called Mocapy++, for inference and learning in BNs and especially DBNs. The main purpose of Mocapy++ is to allow the user to concentrate on the probabilistic model itself, without having to implement customized algorithms. The name Mocapy stands for Markov chain Mo nte Ca rlo and Py thon: the key ingredients in the original implementation of Mocapy (T. Hamelryck, University of Copenhagen, 2004, unpublished). Today, Mocapy has been reimplemented in C++ but the name is kept for historical reasons. Mocapy supports a large range of architectures and probability distributions, and has proven its value in several published applications [9–13]. This article serves as the main single reference for both Mocapy and Mocapy++.
Existing Packages
Kevin Murphy maintains a list of software packages for inference in BNs [14]. Currently, this list contains 54 packages. A small subset of these packages share some key features with Mocapy++ (see Table 1). These packages have an application programming interface (API), perform parameter estimation and are free of charge (at least for academic use). Mocapy++ is mainly intended for use in scientific research, where reproducibility and openness of scientific results are important. Commercial closed source packages are therefore not included in this discussion.
In bioinformatics, models are typically trained using large datasets. Some packages in Table 1 only provide exact inference algorithms that are often not suitable for training models with large datasets. Other packages have no or little support for DBNs, which is important for modelling biomolecular structure. To our knowledge none of the publically available open source toolkits support directional statistics, which has recently become of crucial importance for applications in structural bioinformatics such as modelling protein and RNA structure in 3D detail [9, 10, 12, 15]. Furthermore, Mocapy++ is the only package that uses the stochastic EM [16–18] algorithm for parameter learning (see the Materials and Methods section). These features make Mocapy++ an excellent choice for many tasks in bioinformatics and especially structural bioinformatics.
Implementation
Mocapy++ is implemented as a program library in C++. The library is highly modular and new node types can be added easily. For object serialization and special functions the Boost C++ library [19] is used. All relevant objects are serializable, meaning that Mocapy++ can be suspended and later resumed at any state during training or sampling. The LAPACK library [20] is used for linear algebra routines.
Mocapy++ uses CMake [21] to locate packages and configure the build system and can be used either as a static or shared library. The package includes a Doxygen configuration file for HTML formatted documentation of the source code. An example of a Python interface file for SWIG http://www.swig.org is also included in the package.
Data structures
Most of the internal data is stored in simple Standard Template Library (STL) [22] data structures. However, STL or other public libraries offer little support for multidimensional arrays when the dimension needs to be set at runtime. In Mocapy++ such a multidimensional array is for example needed to store the conditional probability table (CPT) of the discrete nodes. The CPT is a matrix that holds the probabilities of each combination of node and parent values. For example, a discrete node of size 2 with two parents of sizes 3 and 4, respectively, will have a 3 × 4 × 2 matrix as its CPT. Mocapy++ therefore has its own implementation of a multidimensional array, called MDArray. The MDArray class features dynamic allocation of dimensions and provides various slicing operations. The MDArray is also used for storing the training data and other internal data.
Specifying a DBN in Mocapy++
Consider a sequence of observations. Each position in the sequence is characterized by n random variables (called a slice, see Figure 1). Each slice in the sequence can be represented by an ordinary BN, which is duplicated along the sequence as necessary. The sequential dependencies are in turn represented by edges between the consecutive slices. Hence, a DBN in Mocapy++ is defined by three components: a set of nodes that represent all variables for a given slice, the edges between the nodes within a slice (the intra edges) and the edges that connect nodes in two consecutive slices (the inter edges).
Node Types
Mocapy++ supports several node types, each corresponding to a specific probability distribution. The categorical distribution (discrete node), multinomial (for vectors of counts), Gaussian (uni and multivariate), von Mises (uni and bivariate; for data on the circle or the torus, respectively) [23], Kent (5parameter FisherBingham; for data on the sphere) [24] and Poisson distributions are supported. Some node types, such as the bivariate von Mises and Kent nodes, are to our knowledge only available in Mocapy++. The bivariate von Mises and Kent distributions are briefly described here. These distributions belong to the realm of directional statistics, which is concerned with probability distributions on manifolds such as the circle, the sphere or the torus [23, 25].
Kent Distribution
The Kent distribution [9, 24, 26–29], also known as the 5parameter FisherBingham distribution, is a distribution on the 2D sphere (the surface of the 3D ball). It is the 2D member of a larger class of Ndimensional distributions called the FisherBingham distributions. The density function of the Kent distribution is:
where x is a random 3D unit vector that specifies a point on the 2D sphere.
The various parameters can be interpreted as follows:

κ: a concentration parameter. The concentration of the density increases with κ.

β: determines the ellipticity of the equal probability contours of the distribution. The ellipticity increases with β. If β = 0, the Kent distribution becomes the von MisesFisher distribution on the 2D sphere.

γ_{1}: the mean direction.

γ_{2}: the main axis of the elliptical equal probability contours.

γ_{3}: the secondary axis of the elliptical equal probability contours.
The normalizing factor C(κ, β) is approximately given by:
The Kent distribution can be fully characterized by 5 independent parameters. The concentration and the shape of the equal probability contours are characterized by the κ and β parameters, respectively. Two angles are sufficient to specify the mean direction on the sphere, and one additional angle fixes the orientation of the elliptical equal probability contours. The latter three angles are in practice specified by the three orthonormal γ vectors, which form a 3 × 3 orthogonal matrix.
The advantage of the Kent distribution over the von MisesFisher distribution on the 2D sphere is that the equal probability contours of the density are not restricted to be circular: they can be elliptical as well. The Kent distribution is equivalent to a Gaussian distribution with unrestricted covariance. Hence, for 2D directional data the Kent distribution is richer than the corresponding von MisesFisher distribution, i.e. it is more suited if the data contains noncircular clusters. The Kent distribution is illustrated in Figure 2. This distribution was used to formulate FB5HMM [9], which is a probabilistic model of the local structure of proteins in terms of the Cα positions.
Bivariate von Mises Distribution
Another distribution from directional statistics is the bivariate von Mises distribution on the torus [23]. This distribution can be used to model bivariate angular data. The density function of the bivariate von Mises (cosine variant) distribution is:
where C(κ_{1}, κ_{2}, κ_{3}) is the normalizing factor and ϕ, ψ are random angles in [0, 2π[. Such an angle pair defines a point on the torus.
The distribution has 5 parameters:

μ and ν are the means for ϕ and ψ respectively.

κ_{1} and κ_{2} are the concentration of ϕ and ψ respectively.

κ_{3} is related to their correlation.
This distribution is illustrated in Figure 3 and described in greater detail in Mardia et al.[23]. The distribution was used by Boomsma et al.[10] to formulate a probabilistic model of the local structure of proteins in atomic detail.
Inference and Learning
Mocapy++ uses a Markov chain Monte Carlo (MCMC) technique called Gibbs sampling [1] to perform inference, i.e. to approximate the probability distribution over the values of the hidden nodes. Sampling methods such as Gibbs sampling are attractive because they allow complicated network architectures and a wide range of probability distributions.
Parameter learning of a DBN with hidden nodes is done using the expectation maximization (EM) method, which provides a maximum likelihood point estimate of the parameters. In the Estep, the values of the hidden nodes are inferred using the current DBN parameters. In the subsequent Mstep, the inferred values of the hidden nodes are used to update the parameters of the DBN using maximum likelihood estimation. The E and Mstep cycle is repeated until convergence. Parameter learning using the EM algorithm requires a method to perform inference over the possible hidden node values. If one uses a stochastic procedure to perform the Estep (as in Mocapy++), a stochastic version of the EM algorithm is obtained. There are two reasons to use a stochastic Estep. First, deterministic inference might be intractable. Second, certain stochastic versions of the EM algorithm are more robust than the classic version of EM [16]. EM algorithms with a stochastic Estep come in two flavors [1, 17]. In Monte Carlo EM (MCEM), a large number of samples is generated in the EM step. In Stochastic EM (SEM) [16–18] only one sample is generated for each hidden node, and a 'completed' dataset is obtained. In contrast to MCEM, SEM has some clear advantages over deterministic EM algorithms: SEM is less dependent on starting conditions, and has a lower tendency to get stuck at saddle points, or insignificant local maxima. Because only one value needs to be sampled for each hidden node in the Estep, SEM can also be considerably faster than MCEM. SEM is especially suited for large datasets, while for small datasets MCEM is a better choice. Mocapy++ supports both forms of EM.
Results and Discussion
Hamelryck et al.[9] sample realistic protein Cαtraces using an HMM with a Kent output node. Boomsma et al.[10] extend this model to full atomic detail using the bivariate von Mises distribution [23]. In both applications, Mocapy was used for parameter estimation and sampling. Zhao et al.[11] used Mocapy for related work. Mocapy has also been used to formulate a probabilistic model of RNA structure [12] (Figure 1) and to infer functional interactions in a biomolecular network [13].
To illustrate the speed of Mocapy++, we use three parameter estimation benchmarks and report the execution time on a standard PC (Intel Core 2 Duo, 2.33 GHz). The first benchmark is an HMM with 50 slices and two discrete nodes in each slice (one hidden node and one output node). All nodes have 5 states. The second benchmark is similar, but with a 4dimensional Gaussian output node and a 10 state hidden node. The third benchmark is more complex and is shown in Figure 4.
Using a training set consisting of 200 sequences, 100 iterations of SEM take 14 seconds for the discrete HMM, 41 seconds for the Gaussian HMM and 195 seconds for the complex BN. The evolution of the loglikelihood during training is shown in Figure 5.
In practice, the most time consuming step in parameter learning is Gibbs sampling of the hidden nodes. The running time for one sweep of Gibbs sampling for a hidden discrete node is O(l × s) where l is the total number of slices in the data and s is the size of the node. The largest model that, to our knowledge, has been successfully trained with Mocapy++ is an extension of TorusDBN [10]. The dataset consisted of 9059 sequences with a total of more than 1.5 million slices. The model has 11897 parameters and one EMiteration takes 860 seconds. The number of SEM iterations needed for likelihood convergence is around 100.
Toolkits for inference and learning in Bayesian networks use many different algorithms and are implemented in a variety of computer languages (Matlab, R, Java,...); comparisons are thus necessarily unfair or even irrelevant. Therefore, we feel it suffices to point out that Mocapy++ has some unique features (such as the support for directional statistics), and that the benchmarks clearly show that its performance is more than satisfactory for real life problems  both with respect to speed and data set size.
Future Directions of Mocapy++
The core of Mocapy++ described here is not expected to change much in future versions of Mocapy++. However, Mocapy++ is an evolving project with room for new features and additions. We therefore encourage people to propose their ideas for improvements and to participate in the development of Mocapy++. Potential directions include:

Additional probability distributions

Structure learning

Graphical user interface

Plugins for reading data in various formats
Conclusions
Mocapy++ has a number of attractive features that are not found together in other toolkits [14]: it is open source, implemented in C++ for optimal speed efficiency and supports directional statistics. This branch of statistics deals with data on unusual manifolds such as the sphere or the torus [25], which is particularly useful to formulate probabilistic models of biomolecular structure in atomic detail [9–12]. Finally, the use of SEM for parameter estimation avoids problems with convergence [16, 17] and allows for the use of large datasets, which are particularly common in bioinformatics. In conclusion, Mocapy++ provides a powerful machine learning tool to tackle a large range of problems in bioinformatics.
Availability and Requirements

Project name: Mocapy++

Project home page: http://sourceforge.net/projects/mocapy

Operating system(s): Linux, Unix, Mac OS X, Windows with Cygwin

Programming language: C++

Other requirements: Boost, CMake and LAPACK, GNU Fortran

License: GNU GPL
References
 1.
Bishop CM: Pattern recognition and machine learning. Springer; 2006.
 2.
Pearl J: Probabilistic reasoning in intelligent systems: networks of plausible inference. Morgan Kaufmann; 1997.
 3.
Ghahramani Z: Learning dynamic Bayesian networks. Lect Notes Comp Sci 1998, 1387: 168–197. full_text
 4.
Rabiner LR: A tutorial on hidden Markov models and selected applications in speech recognition. Proc IEEE 1989, 77(2):257–286. 10.1109/5.18626
 5.
Durbin R, Eddy SR, Krogh A, Mitchison G: Biological sequence analysis. Cambridge University Press; 1999.
 6.
Raval A, Ghahramani Z, Wild DL: A Bayesian network model for protein fold and remote homologue recognition. Bioinformatics 2002, 18(6):788–801. 10.1093/bioinformatics/18.6.788
 7.
Schmidler SC, Liu JS, Brutlag DL: Bayesian segmentation of protein secondary structure. J Comp Biol 2000, 7(1–2):233–248. 10.1089/10665270050081496
 8.
Chu W, Ghahramani Z, Podtelezhnikov A, Wild DL: Bayesian segmental models with multiple sequence alignment profiles for protein secondary structure and contact map prediction. IEEE/ACM Trans Comp Biol Bioinf 2006, 3(2):98–113. 10.1109/TCBB.2006.17
 9.
Hamelryck T, Kent JT, Krogh A: Sampling realistic protein conformations using local structural bias. PLoS Comp Biol 2006., 2(9): 10.1371/journal.pcbi.0020131
 10.
Boomsma W, Mardia KV, Taylor CC, FerkinghoffBorg J, Krogh A, Hamelryck T: A generative, probabilistic model of local protein structure. Proc Natl Acad Sci USA 2008, 105(26):8932–8937. 10.1073/pnas.0801715105
 11.
Zhao F, Li S, Sterner BW, Xu J: Discriminative learning for protein conformation sampling. Prot Struct Func Bioinf 2008, 73: 228–240. 10.1002/prot.22057
 12.
Frellsen J, Moltke I, Thiim M, Mardia KV, FerkinghoffBorg J, Hamelryck T: A probabilistic model of RNA conformational space. PLoS Comp Biol 2009., 5(6): 10.1371/journal.pcbi.1000406
 13.
KwangHyun C, HyungSeok C, SangMok C: Unraveling the functional interaction structure of a biomolecular network through alternate perturbation of initial conditions. J Biochem Biophys Met 2007, 70(4):701–707. 10.1016/j.jbbm.2007.01.008
 14.
Murphy KP: Software for graphical models: A review. Int Soc Bayesian Anal Bull 2007, 14(4):13–15.
 15.
Hamelryck T: Probabilistic models and machine learning in structural bioinformatics. Stat Met Med Res 2009, 18(5):505–526. 10.1177/0962280208099492
 16.
Nielsen SF: The stochastic EM algorithm: estimation and asymptotic results. Bernoulli 2000, 6(3):457–489. 10.2307/3318671
 17.
Gilks WR, Richardson S, Spiegelhalter D: Markov chain Monte Carlo in practice. Chapman & Hall/CRC; 1995.
 18.
Celeux G, Diebolt J: The SEM algorithm: a probabilistic teacher algorithm derived from the EM algorithm for the mixture problem. Comp Stat Quart 1985, 2: 73–92.
 19.
Abrahams D, Gurtovoy A: C++ template metaprogramming: concepts, tools, and techniques from Boost and beyond. AddisonWesley Professional; 2004.
 20.
Angerson E, Bai Z, Dongarra J, Greenbaum A, Mckenney A, Du Croz J, Hammarling S, Demmel J, Bischof C, Sorensen D: LAPACK: A portable linear algebra library for highperformance computers. Proc Supercomp '90 1990, 2–11. full_text
 21.
Wojtczyk M, Knoll A: A cross platform development workflow for C/C++ applications. ICSEA '08 2008, 224–229.
 22.
Musser DR, Saini A: STL: Tutorial and Reference Guide: C++ Programming with the Standard Template Library. AddisonWesley Professional Computing Series; 1996.
 23.
Mardia KV, Taylor CC, Subramaniam G: Protein bioinformatics and mixtures of bivariate von Mises distributions for angular data. Biometrics 2007, 63(2):505–512. 10.1111/j.15410420.2006.00682.x
 24.
Kent JT: The FisherBingham distribution on the sphere. J Roy Stat Soc 1982, 44: 71–80.
 25.
Mardia KV, Jupp PE: Directional statistics. Wiley; 2000.
 26.
Leong P, Carlile S: Methods for spherical data analysis and visualization. J Neurosci Met 1998, 80: 191–200. 10.1016/S01650270(97)00201X
 27.
Peel D, Whiten W, McLachlan G: Fitting mixtures of Kent distributions to aid in joint set identification. J Am Stat Ass 2001, 96: 56–63. 10.1198/016214501750332974
 28.
Kent J, Hamelryck T: Using the FisherBingham distribution in stochastic models for protein structure. In Quantitative Biology, Shape Analysis, and Wavelets. Volume 24. Edited by: Barber S, Baxter P, Mardia K, Walls R. Leeds University Press; 2005:57–60.
 29.
Boomsma W, Kent J, Mardia K, Taylor C, Hamelryck T: Graphical models and directional statistics capture protein structure. In Interdisciplinary Statistics and Bioinformatics. Volume 25. Edited by: Barber S, Baxter P, Mardia K, Walls R. Leeds University Press; 2006:91–94.
 30.
Raiko T, Valpola H, Harva M, Karhunen J: Building blocks for variational Bayesian learning of latent variable models. J Mach Learn Res 2007, 8: 155–201.
 31.
Murphy KP: The Bayes net toolbox for MATLAB. Comp Sci Stat 2001, 33: 2001.
 32.
Lunn DJ, Thomas A, Best N, Spiegelhalter D: WinBUGS  A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comp 2000, 10(4):325–337. 10.1023/A:1008929526011
 33.
Lacave C, Atienza R, Díez FJ: Graphical explanation in Bayesian networks. In ISMDA '00: Proceedings of the First International Symposium on Medical Data Analysis. London, UK: SpringerVerlag; 2000:122–129.
 34.
Druzdzel MJ: SMILE: Structural modeling, inference, and learning engine and GeNIE: A development environment for graphical decisiontheoretic models. AAAI/IAAI; 1999:902–903.
 35.
Bilmes J, Zweig G: The graphical models toolkit: An open source software system for speech and timeseries processing. IEEE ICASSP 2002, 4: 3916–3919.
 36.
Minka T, Winn J, Guiver J, Kannan A: Infer.NET 2.2.Microsoft Research Cambridge; 1999. [http://research.microsoft.com/infernet]
Acknowledgements
The authors thank the colleagues at the Bioinformatics centre who have helped in the development of Mocapy++: Christian Andreetta, Wouter Boomsma, Mikael Borg, Jes Frellsen, Tim Harder and Kasper Stovgaard. We also thank John T. Kent and Kanti Mardia, University of Leeds, UK and Jesper FerkinghoffBorg, Technical University of Denmark for helpful discussions. We acknowledge funding from the Danish Council for Strategic Research (Program Commission on Nanoscience, Biotechnology and IT, Project: simulating proteins on a millisecond time scale, 2106060009).
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
TH designed and implemented Mocapy in Python. MP designed and implemented Mocapy++. MP drafted the manuscript and TH revised the manuscript. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Paluszewski, M., Hamelryck, T. Mocapy++  A toolkit for inference and learning in dynamic Bayesian networks. BMC Bioinformatics 11, 126 (2010). https://doi.org/10.1186/1471210511126
Received:
Accepted:
Published:
Keywords
 Hide Markov Model
 Expectation Maximization
 Hide Node
 Expectation Maximization Algorithm
 Dynamic Bayesian Network