Skip to content
Prev 7115 / 20628 Next

MIXED MODEL WITH REPEATED MEASURES

Good suggestions; however, there is inherent value in the temporal
progression of the repeated measures, so I need to capture that in some way.
For similar reasons, averaging the values of the independent variables is
problematic, as they progress over time to a final, actual value, which
presumably should be weighted more heavily. In other words, truth is known
on the final repeated measure, but I wish to make accurate predictions much
earlier than the final repeated measure.

As for your first concern, I actually have other dependent variables - I
have just omitted them for simplicity.

I ran your version of the model and unfortunately I obtained the same
singularity error.

Thanks,
Erin

-----Original Message-----
From: r-sig-mixed-models-bounces at r-project.org
[mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf Of Ben Bolker
Sent: Thursday, December 08, 2011 5:01 PM
To: r-sig-mixed-models at r-project.org
Subject: Re: [R-sig-ME] MIXED MODEL WITH REPEATED MEASURES

Erin Ryan <erin at ...> writes:
[snip]
subjects.
I think you're right that DepVar is fixed per individual.
Technical details aside, I'm having trouble seeing how you're going to
estimate the effects of predictor variables that vary within subject when
you've only got one response per subject.
Furthermore, I think what you're terming "RandomVar1" and "RandomVar2"
are probably *not* random variables, but rather are variables
that vary within subject.   For this response variable, I would
suggest averaging the values of RandomVar1 and RandomVar2 per subject and
collapsing the data set to a simple linear model on subjects -- and get rid
of the correlation model at the same time.  For response variables that do
vary within subject, I would suggest

ModelFit <- lme(fixed = DepVar ~FixedVar1+FixedVar2+
   RandomVar1 + RandomVar2, random = 1 | Subject,
  na.action = na.omit, data = dataset, corr = corAR())

_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models