Using substitution matrices to estimate probability distributions for biological sequences.
Accurately estimating probabilities from observations is important for probabilistic-based approaches to problems in computational biology. In this paper we present a biologically-motivated method for estimating probability distributions over discrete alphabets from observations using a mixture model of common ancestors. The method is an extension of substitution matrix-based probability estimation methods. In contrast to previous such methods, our method has a simple Bayesian interpretation and has the advantage over Dirichlet mixtures that it is both effective and simple to compute for large alphabets. The method is applied to estimate amino acid probabilities based on observed counts in an alignment and is shown to perform comparably to previous methods. The method is also applied to estimate probability distributions over protein families and improves protein classification accuracy.
Eskin
Eleazar
E
Department of Computer Science, Columbia University, NY 10027, USA. eeskin@cs.columbia.edu
Noble
William Stafford
WS
Singer
Yoram
Y
D000595:Amino Acid Sequence; D019295:Computational Biology; D008433:Mathematics; D008957:Models, Genetic; D011336:Probability; D011506:Proteins / Q000032:analysis / Q000235:genetics
D011506:Proteins