Nothing
## extra tests of correctness requiring addition Monte Carlo intensity
set.seed(42)
library(spatPomp)
print("Test he10 with towns_selected argument")
h4 <- he10(U=4,towns_selected=c(1,2,11,12),
basic_params = c(
alpha =0.99, iota=0, R0=30,
cohort=0.5, amplitude=0.3, gamma=52,
sigma=52, mu=0.02, sigmaSE=0.05,
rho=0.5, psi=0.1, g=800,
S_0=0.036, E_0=0.00007, I_0=0.00006
))
s4 <- simulate(h4,seed=27)
obs(s4)[,1:2]
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.