Re: [R] lme to determine if there is a group effect

2016-08-25 Thread Thierry Onkelinx
Dear John, lme() not longer requires a GroupedData object. You can directly use a data.frame which is easier to specify different models. You want something like lme(value ~ time * group, random = ~ time|SS, data = data1) PS Note that the R-Sig-mixedmodels is more suited for this kind of

Re: [R] lme to determine if there is a group effect

2016-08-24 Thread Bert Gunter
I never used the groupedData structure, precisely because I found it confusing, but I think: 1. group is *not* a (random) grouping variable; it's a fixed effect covariate. 2. so I believe your groupedData call should be: GD<- groupedData(value~time|SS,data=data1,outer = group) Of course, as you

[R] lme to determine if there is a group effect

2016-08-24 Thread John Sorkin
I apologize for sending this message again. The last time I sent it, the subject line was not correct. I have corrected the subject line. I am trying to run a repeated measures analysis of data in which each subject (identified by SS) has 3 observations at three different times (0, 3, and 6).