blowflies: Model for Nicholson's blowflies.

In pomp: Statistical Inference for Partially Observed Markov Processes

Description

blowflies1 and blowflies2 are pomp objects encoding stochastic delay-difference models.

Details

The data are from "population I", a control culture in one of A. J. Nicholson's experiments with the Australian sheep-blowfly Lucilia cuprina. The experiment is described on pp. 163–4 of Nicholson (1957). Unlimited quantities of larval food were provided; the adult food supply (ground liver) was constant at 0.4g per day. The data were taken from the table provided by Brillinger et al. (1980).

The models are discrete delay equations:

R[t+1] ~ Poisson(P N[t-tau] exp(-N[t-tau]/N0) e[t+1] dt)

S[t+1] ~ binomial(N[t],exp(-delta eps[t+1] dt))

N[t]=R[t]+S[t]

where e[t] and eps[t] are Gamma-distributed i.i.d. random variables with mean 1 and variances sigma.p^2/dt, sigma.d^2/dt, respectively. blowflies1 has a timestep (dt) of 1 day, and blowflies2 has a timestep of 2 days. The process model in blowflies1 thus corresponds exactly to that studied by Wood (2010). The measurement model in both cases is taken to be

y[t] ~ negbin(N[t],1/sigma.y^2)

, i.e., the observations are assumed to be negative-binomially distributed with mean N[t] and variance N[t]+(sigma.y N[t])^2.

Do

file.show(system.file("examples","blowflies.R",package="pomp"))

to view the code that constructs these pomp objects.

References

A. J. Nicholson (1957) The self-adjustment of populations to change. Cold Spring Harbor Symposia on Quantitative Biology, 22, 153–173.

Y. Xia and H. Tong (2011) Feature Matching in Time Series Modeling. Statistical Science 26, 21–46.

E. L. Ionides (2011) Discussion of “Feature Matching in Time Series Modeling” by Y. Xia and H. Tong. Statistical Science 26, 49–52.

S. N. Wood (2010) Statistical inference for noisy nonlinear ecological dynamic systems. Nature 466, 1102–1104.

W. S. C. Gurney, S. P. Blythe, and R. M. Nisbet (1980) Nicholson's blowflies revisited. Nature 287, 17–21.

D. R. Brillinger, J. Guckenheimer, P. Guttorp and G. Oster (1980) Empirical modelling of population time series: The case of age and density dependent rates. in G. Oster (ed.), Some Questions in Mathematical Biology, vol. 13, pp. 65–90. American Mathematical Society, Providence.

See Also

pomp

Examples

pompExample(blowflies)
plot(blowflies1)
plot(blowflies2)

Example output

In 'pomp': the following unrecognized argument(s) will be stored for use by user-defined functions: 'y.init'
In 'pomp': the following unrecognized argument(s) will be stored for use by user-defined functions: 'y.init'
newly created object(s):
 blowflies1 blowflies2 

pomp documentation built on May 2, 2019, 4:09 p.m.