tests/pfilter_methods.R

In panelPomp: Inference for Panel Partially Observed Markov Processes

if (file.exists("_options.R")) source("_options.R")
library(panelPomp,quietly=TRUE)

TESTS_PASS <- NULL
test <- function(expr1,expr2,all="TESTS_PASS",env=parent.frame(),...)
  panelPomp:::test(expr1,expr2,all=all,env=env,...)

ppo <- panelRandomWalk(U=2,N=5)

##  check whether pfiltering a pfilterd.ppomp returns the same object
set.seed(21125715L)
ppf <- pfilter(ppo,Np=10)
logLik(ppf)
unitLogLik(ppf)
set.seed(21125715L)
ppf_ <- pfilter(ppf,Np=10)
logLik(ppf_)
unitLogLik(ppf_)

test(ppf,ppf_)

test(
  wQuotes("Error : in ''coef<-'': cannot change parameters of a filtered object.\n"),
  {coef(ppf_) <- c("test" = 1)}
)

test(
  wQuotes("Error : in ''shared<-'': cannot change parameters of a filtered object.\n"),
  {shared(ppf_) <- c("test" = 1)}
)

test(
  wQuotes("Error : in ''specific<-'': cannot change parameters of a filtered object.\n"),
  {specific(ppf_) <- c("test[unit1]" = 1)}
)

## check whether all tests passed
all(get(eval(formals(test))$all))
if (!all(get(eval(formals(test))$all))) stop("Not all tests passed!")