R Under development (unstable) (2024-02-17 r85935 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > library(spatPomp) Loading required package: pomp > set.seed(3) > > l <- lorenz(U=5,N=3) > > l_pf <- pfilter(l,Np=10) > paste("lorenz pfilter loglik: ", round(logLik(l_pf),10)) [1] "lorenz pfilter loglik: -43.091453281" > > ## 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=10, + Ninter = 2, + lookahead = 1, + Nguide = 10, + kind = 'bootstrap', + verbose = FALSE + ) > paste("lorenz igirf bootstrap hyperbolic loglik:", round(logLik(l_igirf),10)) [1] "lorenz igirf bootstrap hyperbolic loglik: -26.2228972528" > > > > > proc.time() user system elapsed 1.32 0.17 5.14