Re: [R] can I get same results using lme and gls?

From: <toby909_at_gmail.com>
Date: Mon, 21 May 2007 11:30:01 -0700

Douglas Bates wrote:
> On 5/20/07, toby909_at_gmail.com <toby909_at_gmail.com> wrote:
>
>
>>I was wondering how to get the same results with gls and lme. In my lme, the
>>design matrix for the random effects is (should be) a identity matrix and
>>therefore G should add up with R to produce the R matrix that gls would report
>>(V=ZGZ'+R). Added complexity is that I have 3 levels, so I have R, G and say H
>>(V=WHW'+ZGZ'+R). The lme is giving me the correct results, I am having trouble
>>finding the right corresponding specification for the gls.
>
>
> Thanks for including a reproducible example. However, I'm a bit at a
> loss as to why you would want to try to create a gls model that fits a
> mixed-effects model that has random effects for intercept and slope at
> two nested levels. I don't think that corCompSymm will do what you
> want but, to tell the truth, I have difficulty in thinking of the
> model in that form. I much prefer the mixed-effects form.

Thanks for your response. I was too vague, sorry. I think I want to fit a gls model at the lowest level and a mixed effects model with random effects for intercept and slope the the highest level. I dont know if I express that right:
would know how to do using GLS: V=R, R is a "double"block diagonal matrix would also know how to do w/LME: V=WHW'+ZGZ'+R, R is diagonal wanted: V=WHW'+R, R is block diagonal.

This is just for my learning purposes to replicate a textbookexample. This example has 3 levels. The lowest are 2, possibly more, multivatiate measurements (not longitudinal, just to measures, math score and language score in the example, ie unstructured covariance matrix). These are nested within individuals and the individuals are nested within schools. For the lowest level I would choose a gls model for the higher level I would want a lme model both at the same time.

Thanks for your thought.

Toby

>
>
>
>>Thanks for your help.
>>
>>Toby
>>
>>
>>dtaa =
>>read.table("http://www.ats.ucla.edu/stat/mplus/examples/ma_snijders/mlbook1.dat",
>>sep=",")
>>dta1 = reshape(dtaa, list(c("V10","V12")), "score", direction="long",
>>drop=c("V2","V3","V4","V5","V6","V7","V8","V9","V11","V13","V14","V15","V16","V17","V18","V19","V20","V21","V22","V23","V24","V25"))
>>colnames(dta1)[1] = "schoolNR"
>>dta2 = dta1[order(dta1$id),]
>>head(dta2)
>>timef = factor(dta2$time)
>>
>>summary(mdl1l <- lme(score~timef-1, dta2, ~timef-1|schoolNR/id,,,,"ML"))
>>summary(mdl1g <- gls(score~timef-1, dta2, corCompSymm(, ~timef|schoolNR/id),
>>varIdent(, ~1|id*timef),,"ML"))
>>
>>______________________________________________
>>R-help_at_stat.math.ethz.ch 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.
>>
>
>
> ______________________________________________
> R-help_at_stat.math.ethz.ch 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.
>



R-help_at_stat.math.ethz.ch 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. Received on Mon 21 May 2007 - 18:36:12 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 21 May 2007 - 19:30:58 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.