[R-sig-ME] unequal spacing in repeated measures
Ben Bolker
bolker at ufl.edu
Tue Apr 14 20:12:17 CEST 2009
Matthew P. Dekar wrote:
> I would appreciate any advice in regards to the handling of unequal
> spacing in repeated measures regression with mixed models. I sampled
> crayfish with funnel traps and measured environmental predictor
> variables monthly for two years in 12 stream pools. Data are
> presented in terms of the average number of crayfish per trap per
> pool (response variable = catch-per-unit-effort = cpue). However,
> the sampling interval was not fixed so I created a continuous time
> variable (day) indicating the number of days elapsed from the first
> sampling occasion (1,35,55,...,643). As an example, I modeled cpue
> in a repeated measures framework with a single predictor variable
> (stream temperature = temp), pool as the subject/random variable, and
> day was included using a spatial covariance structure (Exponential =
> corExp) in lme following:
>
> cpu1<-lme(cpue~temp, random=~1 | pool, data=CPU, method = "ML")
> cpu2<- update(cpu1, correlation = corExp(form=~day))
>
> Is this an appropriate usage of spatial covariance structures? Can
> the above analysis be replicated in lmer? Thanks very much for your
> time.
>
You should use correlation=corCAR1(form=~day) for a continuous
time covariate.
You can't do this in lmer (yet, or for a while) -- Doug Bates has
stated that implementing correlation structures a la nlme is lower
on his list than working out other issues.
good luck,
Ben Bolker
More information about the R-sig-mixed-models
mailing list