Nothing
library(spatPomp)
set.seed(3)
l <- lorenz(U=5,N=2)
try(lorenz(U=2,N=2))
l_pf <- pfilter(l,Np=10)
paste("lorenz pfilter loglik: ", round(logLik(l_pf),10))
## A call to igirf using the moment-based guide function can test compiled code for eunit_measure, munit_measure, vunit_measure, dunit_measure, runit_measure, rprocess, skeleton, rinit and partrans.
l_igirf <- igirf(l,
Ngirf = 2,
rw.sd = rw_sd(F=0.02, tau=0.02,X1_0=ivp(0),X2_0=ivp(0)),
cooling.type = "hyperbolic",
cooling.fraction.50 = 0.5,
Np=5,
Ninter = 2,
lookahead = 1,
Nguide = 5,
kind = 'bootstrap',
verbose = FALSE
)
paste("lorenz igirf bootstrap hyperbolic loglik:", round(logLik(l_igirf),10))
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.