data(polio)
## ARMA(2,1)
m <- gcmr(y~., data=polio, marginal=nb.marg(), cormat=arma.cormat(2,1),
options=list(seed=71271, nrep=100))
m
## various type of standard errors
res <- cbind(se(m,"he"),se(m,"v"),se(m,"s"),se(m,"hac"))
colnames(res) <- c("hessian", "vscore", "sandwich", "HAC")
resRun the code above in your browser using DataLab