tests/bake.R

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

### test of reproducibility utilities

library(pomp2)
options(digits=2)

set.seed(5499)
w1 <- runif(2)
w4 <- freeze({runif(5)},seed=c(499586,588686,39995866))
w2 <- runif(2)
w5 <- freeze(runif(5),seed=499586)
set.seed(5499)
w3 <- runif(4)
stopifnot(all.equal(c(w1,w2),w3))
stopifnot(all(w4==w5))

set.seed(32765883)
x1 <- bake({runif(4)},file=file.path(tempdir(),"bake1.rds"))
x2 <- bake({runif(4)},file=file.path(tempdir(),"bake2.rds"),seed=32765883)
x3 <- bake({runif(4)},file=file.path(tempdir(),"bake1.rds"))
stopifnot(all.equal(as.numeric(x1),as.numeric(x2)))
stopifnot(all.equal(as.numeric(x1),as.numeric(x3)))
rm(.Random.seed)
x4 <- bake({runif(4)},file=file.path(tempdir(),"bake3.rds"),seed=59566)

set.seed(113848)
stew({y1 <- runif(4)},file=file.path(tempdir(),"stew1.rds"))
stew({y2 <- runif(4)},file=file.path(tempdir(),"stew2.rds"),seed=113848)
print(stew({y3 <- runif(4)},file=file.path(tempdir(),"stew1.rds")))
stopifnot(all.equal(y1,y2))
stopifnot(!exists("y3"))
rm(.Random.seed)
stew({y4 <- runif(4)},file=file.path(tempdir(),"stew3.rds"),seed=59566)

window(sir2(),end=0.5) -> po
simulate(po,seed=1347484107L) -> x
freeze(simulate(po),seed=1347484107L) -> y
attr(y,"seed") <- NULL
attr(y,"kind") <- NULL
attr(y,"normal.kind") <- NULL
stopifnot(identical(x,y))

stopifnot(
  is.null(freeze({rnorm(5); NULL},seed=3494995)),
  is.character(bake({rnorm(5); NULL},seed=3494995,
    file=file.path(tempdir(),"bake4.rds")))
)

rm(.Random.seed)
invisible(bake({runif(4)},file=file.path(tempdir(),"b99.rds"),seed=32765883))
rm(.Random.seed)
invisible(stew({runif(4)},file=file.path(tempdir(),"s99.rda"),seed=32765883))
rm(.Random.seed)
invisible(freeze({runif(4)},seed=32765883))
invisible(freeze({runif(4)}))

detach("package:pomp2", unload=TRUE)