Feng Jingyu wrote: > > I used gls and it still does not provide me different estimates of > variance for each treatment group. Did I do anything wrong? > > lm3<-gls(GSI~treatment,data=z,weights=varIdent(form=~treatment),method="ML") >
try weights = varIdent(form~1|treatment) See the example in library/nlme/scripts/ch05.r ,fm1Orth.gls Dieter -- View this message in context: http://www.nabble.com/Fit-unequal-variance-model-in-R-tp22829549p22873236.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list 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.
