Hello, I was under the impression that in a call to gls, if i specify
gls(.., cor=AR1(form= ~ 1 | subject)) and suppose there are N subjects, then i would have obtained estimates for N Phi (s) (the AR1 coefficient) - but I only get 1 estimate. I then checked "Linear Mixed-Effects Models Using R" and it i must have confused it for the weights(e.g. varIdent) where there are indeed group(strata) specific estimates. Is there a way to include correlations and have group specific estimates? Regards Saptarshi ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
