data(polio)
## ARMA(2,1)
m <- gcmr( y ~ . , data = polio, marginal = negbin.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