tests/gompertz.R

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

options(digits=3)
png(filename="gompertz-%02d.png",res=100)

library(pomp2)

gompertz() -> po

set.seed(1438408329L)

rinit(po)
coef(po)

stopifnot(all.equal(coef(po),partrans(po,coef(po,transform=TRUE),dir="from")))
plot(simulate(po,seed=1438408329L))
pf <- freeze(pfilter(po,Np=1000),seed=1438408329L)
plot(pf)
tj <- trajectory(po,params=c(K=1,r=0.1,sigma=0.1,tau=0.1,X_0=3))
plot(time(po),tj[,,],type="l")

dev.off()