## beta regression with ARMA(1,3) errors
data(HUR)
trend <- scale(time(HUR))
m <- gcmr(HUR ~ trend | trend, marginal = beta.marg,
cormat = arma.cormat(1, 3))
## 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