An efficient algorithm for accurate computation of the Dirichlet- multinomial log-likelihood function

Bioinformatics

Volumn 30, Issue 11, 2014, Pages 1547-1554

  Yu, Peng ^a   Shaw, Chad A ^b  

^a Texas AandM University   (United States)

^b BAYLOR COLLEGE OF MEDICINE   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ALGORITHM; ARTICLE; HUMAN; METAGENOMICS; METHODOLOGY; MICROFLORA; STATISTICAL MODEL;

ALGORITHMS; HUMANS; LIKELIHOOD FUNCTIONS; METAGENOMICS; MICROBIOTA;

EID: 84901381584     PISSN: 13674803     EISSN: 14602059     Source Type: Journal    
DOI: 10.1093/bioinformatics/btu079     Document Type: Article

References (53)

1. Agresti, A. (2002) Categorical Data Analysis. Wiley Series in Probability and Statistics. 2nd edn. Wiley-Interscience, Hoboken, New Jersey
2. Anders, S. and Huber, W. (2010) Differential expression analysis for sequence count data. Genome Biol., 11, R106
3. Bishop, C.M. (2006) Pattern Recognition and Machine Learning. Information Science and Statistics. Springer
4. Bouguila, N. (2008) Clustering of count data using generalized Dirichlet multinomial distributions. IEEE Trans. Knowl. Data Eng., 20, 462-474
5. Brier, S.S. (1980) Analysis of contingency tables under cluster sampling. Biometrika., 67, 591-596
6. Brown, M. et al. (1993) Using Dirichlet mixture priors to derive hidden Markov models for protein families. In: Proceedings of the First International Conference on Intelligent Systems for Molecular Biology. AAAI Press, pp. 47-55
7. Cameron, C.A. and Trivedi, P.K. (2013) Regression Analysis of Count Data. Econometric Society Monographs. 2nd edn. Cambridge University Press, Cambridge
8. Casella, G. and Berger, R. (2002) Duxbury advanced series in statistics and decision sciences. In: Statistical inference. Thomson Learning, Pacific Grove, CA.
9. Crowder, M.J. (1978) Beta-binomial ANOVA for proportions. Appl. Stat., 27, 34-37
10. Farewell, D.M. and Farewell, V.T. (2013) Dirichlet negative multinomial regression for overdispersed correlated count data. Biostatistics., 14, 395-404
11. Fisher, R.A. (1973) Statistical Methods for Research Workers. 14th edn. Hafner Publishing Company, Edinburgh
12. Freitag, E. and Busam, R. (2009) Complex Analysis. 2nd edn. Springer
13. Gomez-Alvarez, V. et al. (2012) Metagenome analyses of corroded concrete wastewater pipe biofilms reveal a complex microbial system. BMC Microbiol., 12, 122
14. Haseman, J.K. and Kupper, L.L. (1979) Analysis of dichotomous response data from certain toxicological experiments. Biometrics., 35, 281-293. (Pubitemid 9136580
15. Hilbe, J.M. (2011) Negative Binomial Regression. 2nd edn. Cambridge University Press, Cambridge, UK.
16. IEEE Task P754. (2008) IEEE 754-2008, Standard for Floating-Point Arithmetic. IEEE, New York, NY.
17. Kim, B.S. and Margolin, B.H. (1992) Testing goodness of fit of a multinomial model against overdispersed alternatives. Biometrics., 48, 711-719
18. Kleinman, J.C. (1973) Proportions with extraneous variance: Single and independent sample. J. Am. Stat. Assoc., 68, 46-54
19. La Rosa, P.S. et al. (2012) Hypothesis testing and power calculations for taxonomicbased human microbiome data. PLoS One, 7, e52078
20. Lauter, C. and Dinechin, F.D. (2008) Optimizing polynomials for floating-point implementation. In: Proceedings of the 8th Conference on Real Numbers and Computers, Santiago De Compostela, Spain
21. Leckenby, J.D. and Kishi, S. (1984) The Dirichlet multinomial distribution as a magazine exposure model. J. Mark. Res., 21, 100-106
22. Lesnoff, M. and Lancelot, R. (2012) aod: Analysis of Overdispersed Data. R package version 1.3
23. Lu, X. et al. (2013) Son connects the splicing-regulatory network with pluripotency in human embryonic stem cells. Nat. Cell Biol., 15, 1141-52
24. Lutz, C.S. and Moreira, A. (2011) Alternative mRNA polyadenylation in eukaryotes: An effective regulator of gene expression. Wiley Interdiscip. Rev. RNA, 2, 22-31
25. MacKay, D.J.C. and Bauman Peto, L.C. (1994) A hierarchical Dirichlet language model. Nat. Lang. Eng., 1, 1-19
26. Madsen, R.E. et al. (2005)Modeling word burstiness using the Dirichlet distribution. In: Proceedings of the 22nd International Conference on Machine Learning. ICML'05. ACM, New York, NY, pp. 545-552
27. McCullagh, P. and Nelder, J.A. (1989) Generalized Linear Models. Monographs on Statistics and Applied Probability. Chapman and Hall/CRC, New York
28. Mehta, C.R. and Patel, N.R. (1986) Algorithm 643: Fexact: A fortran subroutine for Fisher's exact test on unordered r - c contingency tables. ACM Trans. Math. Softw., 12, 154-161
29. Metzker, M.L. (2010) Sequencing technologies - The next generation. Nat. Rev. Genet., 11, 31-46
30. Mimno, D. and McCallum, A. (2008) Topic Models Conditioned on Arbitrary Features with Dirichlet-Multinomial Regression. UAI, Helsinki, Finland
31. Mosimann, J.E. (1962) On the compound multinomial distribution, the multivariate -distribution, and correlations among proportions. Biometrika., 49, 65-82
32. Mosimann, J.E. (1963) On the compound negative multinomial distribution and correlations among inversely sampled pollen counts. Biometrika., 50, 47-54
33. The NCBI BioProject website. (2010) Human Microbiome Project 16S rRNA Clinical Production Pilot (ID: 48335). http://www.ncbi.nlm.nih.gov/ bioproject?term48335
34. Neerchal, N.K. and Morel, J.G. (2005) An improved method for the computation of maximum likelihood estimates for multinomial overdispersion models. Comput. Stat. Data Anal., 49, 33-43
35. Paul, S.R. et al. (2005) Fisher information matrix of the Dirichlet-multinomial distribution. Biom. J., 47, 230-236
36. Poortema, K. (1999) On modelling overdispersion of counts. Stat. Neerl., 53, 5-20
37. R Core Team. (2013) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria
38. Robinson, M.D. et al. (2010) edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics., 26, 139-40
39. Rosa, P.S.L. et al. (2013) HMP: Hypothesis Testing and Power Calculations for Comparing Metagenomic Samples from HMP. R package version 1.3.1
40. Rowe, C.H. (1931) A proof of the asymptotic series for log -z and log -za. Ann. Math., Second Ser, 32, 10-16
41. Rust, R.T. and Leone, R.P. (1984) The mixed-media Dirichlet multinomial distribution: A model for evaluating television-magazine advertising schedules. J. Mark. Res., 21, 89-99
42. Sjó lander, K. et al. (1996) Dirichlet mixtures: A method for improved detection of weak but significant protein sequence homology. Comput. Appl. Biosci., 12, 327-45
43. Skellam, J.G. (1948) A probability distribution derived from the binomial distribution by regarding the probability of success as variable between the sets of trials. J. R. Stat. Soc. Ser. B Methodol., 10, 257-261
44. Stein, W. et al. (2012) Sage Mathematics Software (Version 5.0.1). The Sage Development Team
45. Tarone, R.E. (1979) Testing the goodness of fit of the binomial distribution. Biometrika., 66, 585-590
46. Tvedebrink, T. (2009) dirmult: Estimation in Dirichlet-Multinomial Distribution. R package version 0.1.2
47. Wan, J. (2012) Global analysis of alternative polyadenylation regulation using highthroughput sequencing. PhD Thesis, University of Iowa
48. Wang, E.T. et al. (2008) Alternative isoform regulation in human tissue transcriptomes. Nature., 456, 470-476
49. Whittaker, E.T. and Watson, G.N. (1927) A Course of Modern Analysis. 4th edn. Reprinted 1990. Cambridge University Press, Cambridge, UK.
50. Winkelmann, R. (2008) Econometric Analysis of Count Data. 5th edn. Springer, Berlin, Germany
51. Yee, T.W. (2010) The VGAM package for categorical data analysis. J. Stat. Softw., 32, 1-34
52. Yee, T.W. (2012) VGAM: Vector Generalized Linear and Additive Models. R package version 0.9-0
53. Yee, T.W. and Wild, C.J. (1996) Vector generalized additive models. J. R. Stat. Soc. B, 58, 481-493