title: A multiscale maximum entropy moment closure for locally regulated space-time point process models of population dynamics.
abstract: The prevalence of structure in biological populations challenges fundamental assumptions at the heart of continuum models of population dynamics based only on mean densities (local or global). Individual-based models (IBMs) were introduced during the last decade in an attempt to overcome this limitation by following explicitly each individual in the population. Although the IBM approach has been quite useful, the capability to follow each individual usually comes at the expense of analytical tract ability, which limits the generality of the statements that can be made. For the specific case of spatial structure in populations of sessile (and identical) organisms, space-time point processes with local regulation seem to cover the middle ground between analytical tract ability and a higher degree of biological realism. This approach has shown that simplified representations of fecundity, local dispersal and density-dependent mortality weighted by the local competitive environment are sufficient to generate spatial patterns that mimic field observations. Continuum approximations of these stochastic processes try to distill their fundamental properties, and they keep track of not only mean densities, but also higher order spatial correlations. However, due to the non-linearities involved they result in infinite hierarchies of moment equations. This leads to the problem of finding a 'moment closure'; that is, an appropriate order of (lower order) truncation, together with a method of expressing the highest order density not explicitly modelled in the truncated hierarchy in terms of the lower order densities. We use the principle of constrained maximum entropy to derive a closure relationship for truncation at second order using normalisation and the product densities of first and second orders as constraints, and apply it to one such hierarchy. The resulting 'maxent' closure is similar to the Kirkwood superposition approximation, or 'power-3' closure, but it is complemented with previously unknown correction terms that depend mainly on the avoidance function of an associated Poisson point process over the region for which third order correlations are irreducible. This domain of irreducible triplet correlations is found from an integral equation associated with the normalisation constraint. This also serves the purpose of a validation check, since a single, non-trivial domain can only be found if the assumptions of the closure are consistent with the predictions of the hierarchy. Comparisons between simulations of the point process, alternative heuristic closures, and the maxent closure show significant improvements in the ability of the truncated hierarchy to predict equilibrium values for mildly aggregated spatial patterns. However, the maxent closure performs comparatively poorly in segregated ones. Although the closure is applied in the context of point processes, the method does not require fixed locations to be valid, and can in principle be applied to problems where the particles move, provided that their correlation functions are stationary in space and time.
authors.affiliation: Department of Mathematics, University of Glasgow, Glasgow, UK. mraghib@princeton.edu
authors.forename: Michael
authors.initials: M
authors.lastname: Raghib
authors.forename: Nicholas A
authors.initials: NA
authors.lastname: Hill
authors.forename: Ulf
authors.initials: U
authors.lastname: Dieckmann
mesh_terms: D000465:Algorithms; D003198:Computer Simulation; D005298:Fertility; D016015:Logistic Models; D008954:Models, Biological; D009026:Mortality; D018521:Plant Physiological Phenomena; D011156:Population Density; D011157:Population Dynamics; D011158:Population Growth; D011336:Probability; D016008:Statistical Distributions; D013269:Stochastic Processes
references.citation: J Theor Biol. 1999 Feb 21;196(4):397-454
references.pmid: 10036198
references.citation: Science. 2000 Dec 1;290(5497):1758-61
references.pmid: 11099413
references.citation: J Chem Phys. 2004 Aug 22;121(8):3657-66
references.pmid: 15303932
references.citation: Proc Biol Sci. 1999 Apr 22;266(1421):859-67
references.pmid: 10343409
references.citation: Theor Popul Biol. 2009 Mar-May;75(2-3):133-41
references.pmid: 19344628
references.citation: Proc Natl Acad Sci U S A. 1920 Jun;6(6):275-88
references.pmid: 16576496
references.citation: Nature. 2007 Sep 13;449(7159):209-12
references.pmid: 17851523
references.citation: Proc Biol Sci. 2005 Jun 22;272(1569):1279-85
references.pmid: 16024393
references.citation: Theor Popul Biol. 2006 Feb;69(1):13-33
references.pmid: 16246386
references.citation: Ecol Lett. 2008 Sep;11(9):929-36
references.pmid: 18485001
references.citation: J Math Biol. 2000 Nov;41(5):387-429
references.pmid: 11151706
references.citation: J Theor Biol. 2004 Aug 7;229(3):421-32
references.pmid: 15234208
references.citation: Proc Biol Sci. 1997 Aug 22;264(1385):1149-56
references.pmid: 9308191
references.citation: J Math Biol. 1988;26(3):263-98
references.pmid: 3411255
references.citation: Bull Math Biol. 2005 Jul;67(4):855-73
references.pmid: 15893556
references.citation: Bull Math Biol. 2006 Aug;68(6):1315-33
references.pmid: 17149819
references.citation: Bull Math Biol. 2001 Jul;63(4):603-24
references.pmid: 11497160
references.citation: Nature. 2001 Jul 19;412(6844):328-31
references.pmid: 11460162
references.citation: Proc Natl Acad Sci U S A. 2006 Aug 22;103(34):12781-6
references.pmid: 16912114
references.citation: Am Nat. 2003 Aug;162(2):135-48
references.pmid: 12858259
references.citation: J Chem Phys. 2005 Sep 8;123(10):104507
references.pmid: 16178610
references.citation: Bull Math Biol. 2004 Sep;66(5):1027-37
references.pmid: 15294417
references.citation: Math Biosci. 2003 May;183(1):15-35
references.pmid: 12604133
references.citation: Theor Popul Biol. 2006 Aug;70(1):26-42
references.pmid: 16442137
references.citation: Science. 2000 May 26;288(5470):1414-8
references.pmid: 10827950
references.citation: Theor Popul Biol. 1997 Dec;52(3):179-97
references.pmid: 9466960