Combining principal component analysis with parameter line-searches to improve the efficacy of Metropolis–Hastings MCMC

Environmental and Ecological Statistics

Volumn 22, Issue 2, 2015, Pages 247-274

  Kennedy, David A ^a,b,c   Dukic, Vanja ^d   Dwyer, Greg ^a  

^a UNIVERSITY OF CHICAGO   (United States)

^b PENNSYLVANIA STATE UNIVERSITY   (United States)

^c NATIONAL INSTITUTES OF HEALTH   (United States)

^d UNIVERSITY OF COLORADO   (United States)

Author keywords

Birth death model; MCMC; Parameter line search; Survival time data; Within host model

Indexed keywords

EID: 84939875610     PISSN: 13528505     EISSN: 15733009     Source Type: Journal    
DOI: 10.1007/s10651-014-0297-0     Document Type: Article

References (78)

1. Alizon S, van Baalen M (2008) Acute or chronic? Within-host models with immune dynamics, infection outcome, and parasite evolution. Am Nat 172:E244–E256
2. Antia R, Levin B, May R (1994) Within-host population-dynamics and the evolution and maintenance of microparasite virulence. Am Nat 144:457–472
3. Armenian H, Lilienfeld A (1983) Incubation period of disease. Epidemiol Rev 5:1–15
4. Ashida M, Brey P (1998) Molecular mechanisms of immune responses in insects. Chapman & Hall, London
5. Baldwin K, Hakim R (1991) Growth and differentiation of the larval midgut epithelium during molting in the moth, Manduca sexta. Tissue Cell 23:411–422
6. Beaumont M, Zhang W, Balding D (2002) Approximate Bayesian computation in population genetics. Genetics 162:2025–2035
7. Bogich T, Shea K (2008) A state-dependent model for the optimal management of an invasive metapopulation. Ecol Appl 18:748–761
8. Bolker B (2008) Ecological models and data in R. Princeton University Press, New Jersey
9. Braun M (1983) Differential equations and their applications, an introduction to applied mathematics, 3rd edn. Springer, New York
10. Brigham C, Power A, Hunter A (2002) Evaluating the internal consistency of recovery plans for federally endangered species. Ecol Appl 12:648–654
11. Brockwell A (2006) Parallel Markov chain Monte Carlo simulation by pre-fetching. J Comput Graph Stat 15:246–261
12. Chakerian J, Holmes S (2012) Computational tools for evaluating phylogenetic and hierarchical clustering trees. J Comput Graph Stat 21:581–599
13. Comon P (1994) Independent component analysis, a new concept. Signal Proces 36:287–314
14. Cory J, Myers J (2003) The ecology and evolution of insect baculoviruses. Annu Rev Ecol Evol Syst 34:239–272
15. Cowles M, Carlin B (1996) Markov chain Monte Carlo convergence diagnostics: a comparative review. J Am Stat Assoc 91:883–904
16. Craiu R, Rosenthal J, Yang C (2009) Learn from thy neighbor: parallel-chain and regional adaptive MCMC. J Am Stat Assoc 104:1454–1466
17. Csillery K, Blum M, Gaggiotti O, Francois O (2010) Approximate Bayesian computation (ABC) in practice. Trends Ecol Evol 25:410–418
18. Doak DF, Morris WF (2010) Demographic compensation and tipping points in climate-induced range shifts. Nature 467:959–962
19. Doob J (1945) Markoff chains: denumerable case. Trans Am Math Soc 58:455–473
20. Dukic V, Lopes H, Polson N (2012) Tracking epidemics with Google Flu trends data and a state-space SEIR model. J Am Stat Assoc 107:1410–1426
21. Feng H, Gould F, Huang Y, Jiang Y, Wu K (2010) Modeling the population dynamics of cotton bollworm Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) over a wide area in northern China. Ecol Model 221:1819–1830
22. Fuller E, Elderd B, Dwyer G (2012) Pathogen persistence in the environment and insect-baculovirus interactions: disease-density thresholds, epidemic burnout and insect outbreaks. Am Nat 179:E70–E96
23. Fuller S, Millet L (2011) Computing performance: Game over or next level? IEEE Comput 44:31–38
24. Geer D (2005) Chip makers turn to multicore processors. IEEE Comput 38:11–13
25. Gelman A, Rubin D (1992) Inference from iterative simulation using multiple sequences. Stat Sci 7:457–472
26. Gilchrist M, Sasaki A (2002) Modeling host-parasite coevolution: a nested approach based on mechanistic models. J Theor Biol 218:289–308
27. Gilks W, Roberts G (1996) Markov chain Monte Carlo in practice, chapter Introducing Markov chain Monte Carlo. Chapman & Hall, London
28. Gillespie D (1977) Exact stochastic simulation of coupled chemical-reactions. J Phys Chem 81:2340–2361
29. Girolami M, Calderhead B (2011) Riemann manifold Langevin and Hamiltonian Monte Carlo methods. J R Stat Soc Ser B 73:123–214
30. Grant A, Restif O, McKinley T, Sheppard M, Maskell D, Mastroeni P (2008) Modelling within-host spatiotemporal dynamics of invasive bacterial disease. PLoS Biol 6:757–770
31. Haario H, Saksman E, Tamminen J (2001) An adaptive Metropolis algorithm. Bernoulli 7:223–242
32. Hartig F, Calabrese JM, Reineking B, Wiegand T, Huth A (2011) Statistical inference for stochastic simulation models—theory and application. Ecol Lett 14:816–827
33. Heidelberger P, Welch P (1983) Simulation run length control in the presence on an initial transient. Oper Res 31:1109–1144
34. Hoover K, Washburn J, Volkman L (2000) Midgut-based resistance of Heliothis virescens to baculovirus infection mediated by phytochemicals in cotton. J Insect Physiol 46:999–1007
35. Hunter-Fujita F, Entwistle P, Evans H, Crook N (1998) Insect viruses and pest management. Wiley, Chichester
36. Ionides E, Breto C, King A (2006) Inference for nonlinear dynamical systems. Proc Natl Sci USA 103:18438–18443
37. Jacob P, Robert C, Smith M (2011) Using parallel computation to improve independent Metropolis-Hastings based estimation. J Comput Graph Stat 20:616–635
38. Jolliffe I (1986) Principal component analysis. Springer, New York
39. Karlin S, Taylor H (1975) A first course in stochastic processes. Academic, New York
40. Kennedy DA, Dukic V, Dwyer G (2014) The mechanisms determining the within-host population dynamics of an insect pathogen. Am Nat 184:407–423
41. Khorsheed E, Hurn M, Jennison C (2011) Mapping electron density in the ionosphere: a principal component MCMC algorithm. Comput Stat Data Anal 55:338–352
42. Kimura M (1983) The neutral theory of molecular evolution. Cambridge University Press, New York
43. King A, Shrestha S, Harvill E, Bjørnstad O (2009) Evolution of acute infections and the invasion-persistence trade-off. Am Nat 173:446–455
44. Kot M (2001) Elements of mathematical ecology. Cambridge University Press, Cambridge
45. Lele S, Dennis B, Lutscher F (2007) Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods. Ecol Lett 10:551–563
46. Lele S, Nadeem K, Schmuland B (2010) Estimability and likelihood inference for generalized linear mixed models using data cloning. J Am Stat Assoc 105:1617–1625
47. Liu J (2001) Monte Carlo strategies in scientific computing. Springer, Berlin
48. Luenberger D, Ye Y (2008) Linear and nonlinear programming, 3rd edn. Springer Science and Business Media, New York
49. McNeil J, Cox-Foster D, Gardner M, Slavicek J, Thiem S, Hoover K (2010) Pathogenesis of Lymantria dispar multiple nucleopolyhedrovirus (LdMNPV) in L. dispar and mechanisms of developmental resistance. J Gen Virol 91:1590–1600
50. Meynell G (1957) The applicability of the hypothesis of independent action to fatal infections in mice given Salmonella typhimurium by mouth. J Gen Microbiol 16:396–404
51. Miller G (2010) Markov chain Monte Carlo calculations allowing parallel processing using a variant of the Metropolis algorithm. Open Numer Methods J 2:12–17
52. Morgan B (1992) Analysis of quantal response data. Chapman & Hall, London
53. Mudholkar G, Srivastava D, Kollia G (1996) A generalization of the Weibull distribution with application to the analysis of survival data. J Am Stat Assoc 91:1575–1583
54. Plummer M, Best N, Cowles K, Vines K. (2009) coda: Output analysis and diagnostics for MCMC. R package version 0.13-4
55. Ponciano J, Burleigh J, Braun E, Taper M (2012) Assessing parameter identifiability in phylogenetic models using data cloning. Syst Biol 61:955–972
56. Press W, Teukolsky S, Vetterling W, Flannery B (1992) Numerical recipes in C. Cambridge University Press, Cambridge
57. Development Core Team R (2009) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0
58. Robert C, Cornuet J, Marin J, Pillai N (2011) Lack of confidence in approximate Bayesian computation model choice. Proc Natl Acad Sci USA 108:15112–15117
59. Roberts G, Gelman A, Gilks W (1997) Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann Appl Probab 7:110–120
60. Rosenthal J (2000) Parallel computing and Monte Carlo algorithms. Far East J Theor Stat 4:207–236
61. Saaty T (1961) Some stochastic-processes with absorbing barriers. J R Stat Soc Ser B Stat Methodol 23:319–334
62. Schmid-Hempel P (2005) Evolutionary ecology of insect immune defenses. Annu Rev Entomol 50:529–551
63. Schölkopf B, Smola A, Müller K (1998) Nonlinear component analysis as a kernel eigenvalue problem. Neural Comput 10:1299–1319
64. Shapiro M, Farrar R Jr, Domek J, Javaid I (2002) Effects of virus concentration and ultraviolet irradiation on the activity of corn earworm and beet armyworm (Lepidoptera:Noctuidae) nucleopolyhedroviruses. J Econ Entomol 95:243–249
65. Shapiro M, Robertson J, Bell R (1986) Quantitative and qualitative differences in gypsy moth (Lepidoptera: Lymantriidae) nucleopolyhedrosis virus produced in different-aged larvae. J Econ Entomol 79:1174–1177
66. Shortley G (1965) A stochastic model for distributions of biological response times. Biometrics 21:562–582
67. Solonen A, Ollinaho P, Laine M, Haario H, Tamminen J, Jarvinen H (2012) Efficient MCMC for climate model parameter estimation: parallel adaptive chains and early rejection. Bayesian Anal 7:715–736
68. Strid I (2010) Efficient parallelisation of Metropolis-Hastings algorithms using a prefetching approach. Comput Stat Data Anal 54:2814–2835
69. Trudeau D, Washburn J, Volkman L (2001) Central role of hemocytes in Autographa californica M nucleopolyhedrovirus pathogenesis in Heliothis virescens and Helicoverpa zea. J Virol 75:996–1003
70. Turchin P (2003) Complex population dynamics: a theoretical/empirical synthesis. Princeton University Press, Princeton
71. van Beek N, Flore P, Wood H, Hughes P (1990) Rate of increase of Autographa californica nuclear polyhedrosis virus in Trichoplusia ni larvae determined by DNA-DNA hybridization. J Invertebr Pathol 55:85–92
72. van Beek N, Hughes P, Wood H (2000) Effects of incubation temperature on the dose-survival time relationship of Trichoplusia ni larvae infected with Autographa californica nucleopolyhedrovirus. J Invertebr Pathol 76:185–190
73. van Beek N, Wood H, Hughes P (1988) Quantitative aspects of nuclear polyhedrosis virus infections in Lepidopterous larvae: the dose-survival time relationship. J Invertebr Pathol 51:58–63
74. van den Berg S, Beem L, Boomsma D (2006) Fitting genetic models using Markov chain Monte Carlo algorithms with BUGS. Twin Res Hum Genet 9:334–342
75. Vaughan T, Drummond P, Drummond A (2012) Within-host demographic fluctuations and correlations in early retroviral infection. J Theor Biol 295:86–99
76. Wilkinson D (2005) Handbook of Parallel computing and statistics, chapter parallel Bayesian computation. Dekker/CRC Press, New York
77. Yan J, Cowles M, Wang S, Armstrong M (2007) Parallelizing MCMC for Bayesian spatiotemporal geostatistical models. Stat Comput 17:323–335
78. Zwart M, Hemerik L, Cory J, de Visser J, Bianchi F, Van Oers M, Vlak J, Hoekstra R, Van der Werf W (2009) An experimental test of the independent action hypothesis in virus-insect pathosystems. Proc R Soc Lond Ser B-Biol Sci 276:2233–2242