tests/slice_design.R

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

png(filename="slice_design-%02d.png",res=100)

library(pomp2)
library(dplyr)

## A single 11-point slice through the point c(A=3,B=8,C=0) along the B direction.
x <- sliceDesign(center=c(A=3,B=8,C=0),B=seq(0,10,by=1))
x %>% count(slice) %>% as.data.frame()
plot(x)

## Two slices through the same point along the A and C directions.
x <- sliceDesign(c(A=3,B=8,C=0),A=seq(0,5,by=1),C=seq(0,5,length=11))
x %>% count(slice) %>% as.data.frame()
plot(x)

try(x <- sliceDesign(center=c(A=3,C=0),B=seq(0,10,by=1)))
try(x <- sliceDesign(center=c(A=3),B=seq(0,10,by=1),C=c(1,2,3)))
try(x <- sliceDesign(center=c(3),B=seq(0,10,by=1),C=c(1,2,3)))
try(x <- sliceDesign(center=c(A="3"),B=seq(0,10,by=1),C=c(1,2,3)))
try(x <- sliceDesign(center=c(A=3,B=2),seq(1,3),A=seq(0,5)))

dev.off()