Message-ID: <22873995.post@talk.nabble.com>
Date: 2009-04-03T18:33:15Z
From: Feng Jingyu
Subject: Fit unequal variance model in R
In-Reply-To: <22873236.post@talk.nabble.com>
Thanks a lot. The problem is solved. It took me a while to understand the
output from the R. With little calculation, I am able to match results from
R to SAS.
Dieter Menne wrote:
>
>
> 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-tp22829549p22873995.html
Sent from the R help mailing list archive at Nabble.com.