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
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
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).
3 matches
Mail list logo