tests/simulate_data_frame.R

In kidusasfaw/pomp: Statistical Inference for Partially Observed Markov Processes

options(digits=3)

library(pomp2)
library(dplyr)

set.seed(1041414791L)

data.frame(t=0:10,y=runif(11)) %>%
  simulate(times="t",t0=0,
    rprocess=discrete_time(
      step.fun=function(r,x,...) {
        c(x=r*x*exp(-x))
      }),
    rmeasure=function(x,...) {
      c(y=rpois(n=1,100*x))
    },
    params=list(r=15,x.0=0.1)
  )

try(data.frame(t=0:10,y=runif(11)) %>%
    simulate(times="t",t0=0,
      rprocess=discrete_time(
        step.fun=function(r,x,...) {
          c(x=r*x*exp(-x))
        }),
      rmeasure=function(x,...) {
        c(y=rpois(n=1,100*x))
      },
      params=list(r=15,x.0=0.1),
      format="bob"
    )
)

try(data.frame(t=0:10,y=runif(11)) %>%
  simulate(times=1:50))
try(data.frame(t=0:10,y=runif(11)) %>%
    simulate(times="t"))
try(data.frame(t=0:10,y=runif(11)) %>%
    simulate(times="t",t0=0))
try(data.frame(t=0:10,y=runif(11)) %>%
    simulate(times="t",t0=0))
try(data.frame(t=0:10,y=runif(11)) %>%
    simulate(times="t",t0=0,params=c(x.0=1)))