Modeling Residuals vs. Modeling Random-effects Covariances

Dear all,

I had a quick question. The variance-covariance structure from the
"glmmTMB" package and the "nlme" package are not the same thing?

Is there any other package capable of modeling the Level-1 residual
variance-covariance structure like in "nlme"? *-- Many Thanks, Simon *




*glmmTMB::glmmTMB(y ~ times + ar1(times | subjects), data = data) ## DON'T
RUN   nlme::lme    (y ~ times, random = ~ times | subjects, correlation =
corAR1(), data = data) ## DON'T RUN *