Blowfly model

Description

synlik object containing the blowfly model proposed by Wood (2010). The main components are the simulator blowSimul and the statistics blowStats, described in the same reference.

Author(s)

Simon Wood and Matteo Fasiolo <matteo.fasiolo@gmail.com>.

References

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

See Also

blowSimul

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
data(blow_sl)

plot(blow_sl)
simulate(blow_sl, stats = TRUE)

slik(blow_sl,
     param  = log( c( "delta" = 0.16, "P" = 6.5, "N0" = 400,
                      "var.p" = 0.1, "tau" = 14, "var.d" = 0.1)  ),
     nsim    = 1e3)

# Using Nicholson's data
data(bf1)

plot(blow_sl)

blow_sl@data <- bf1$pop
blow_sl@extraArgs$obsData <- bf1$pop #Important: blow_sl@blowStats uses the observed data

slik(blow_sl,
     param  = log( c( "delta" = 0.16, "P" = 6.5, "N0" = 400,
                      "var.p" = 0.1, "tau" = 14, "var.d" = 0.1)  ),
     nsim    = 1e3)

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.