USL20.data <- simulate(USL20,sampleT=50,start=1950)
m <- l(USL20,USL20.data)
# tfplot(m)
AIC(m)
USL20.f <- forecast(m,horizon=50)
USL20.fx <- merge.forecast(USL20.f)
tfplot(USL20.f)
#
# USL20c Growth Controlled
# by stabilizing growth rates
#
f[1,1] <- 1.02
f[2,2] <- 0.97
f[3,3] <- 0.97
#
USL20c <- SS(F=f,H=h,K=k,z0=c( 0.166789178, 0.004293523, -0.006797188,
1.0000000),
output.names=c("US1","US2","US3"))
stability(USL20c)
shockDecomposition(toSSChol(USL20c))
USL20c.data <- simulate(USL20c,sampleT=50,start=1950)
m <- l(USL20c,USL20c.data)
# tfplot(m)
AIC(m)
USL20c.f <- forecast(m,horizon=50)
USL20c.fx <- merge.forecast(USL20c.f)
tfplot(USL20c.f)