blow_sl: Blowfly model

Description Author(s) References See Also Examples

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 <[email protected]>.

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: [email protected] 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)

Example output

This is "synlik"  0.1.1 
         [,1]          [,2]          [,3]       [,4]         [,5]          [,6]
[1,] 1.091458 -0.0001588608 -3.511022e-08 -0.2045743 1.591921e-05 -2.011993e-09
         [,7]          [,8]     [,9]    [,10]   [,11]   [,12]   [,13]   [,14]
[1,] 1.125509 -2.207958e-05 2428.776 356.9915 3343398 3234499 3033722 2739543
       [,15]   [,16]   [,17]   [,18]  [,19]    [,20]     [,21]     [,22] [,23]
[1,] 2385304 1983366 1534151 1051698 555459 62508.81 -410824.3 -857623.4    26
[1] -63.19726
[1] -629.1512

synlik documentation built on May 29, 2017, 12:16 p.m.