if (FALSE) {
set.seed(1)
m <- MTDmodel(Lambda = c(1, 3), A = c(0, 1), lam0 = 0.05)
print(m) # compact display: lags (Z^-) and state space
s <- summary(m)
str(s)
coef(m) # list(lambdas = ..., pj = ..., p0 = ...)
transitP(m) # global transition matrix P
pj(m); p0(m); lambdas(m); lags(m); Lambda(m); states(m)
X <- perfectSample(m, N = 400)
logLik(m, X)
}
Run the code above in your browser using DataLab