[R] lme command

From: Enrico Colosimo <enricoc57_at_gmail.com>
Date: Mon, 14 Jun 2010 10:37:43 -0300


Hi,

I am doing a longitudinal data set fit using lme. I used two forms of the lme command and I am getting two different outputs.

FIRST
out<-lme(Altura~Idade+Idade2+sexo+status+Idade:sexo+Idade:status+Idade2:sexo+Idade2:status, random=(list(ident=~Idade+Idade2)))

SECOND
out<-lme(Altura~Idade+Idade2+sexo+status+Idade:sexo+Idade:status+Idade2:sexo+Idade2:status, random= ~Idade+Idade2|ident,data=dados)

I got weird results from the first one and could not understand the reason of it. All the results are
exactly the same but the intercetp, and the two main terms sexo (gender) and status (treatment).
That differences made a lot of difference in the final results.

Anybody can tell me what is the differences between them? Thanks.

Enrico.



R-help_at_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. Received on Mon 14 Jun 2010 - 15:34:43 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 14 Jun 2010 - 16:10:32 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.

list of date sections of archive