[R-sig-ME] Is it possible to account for temporal autocorrelation in MCMCglmm?
David Villegas Ríos
chirleu at gmail.com
Wed Apr 1 14:34:19 CEST 2015
Dear all,
For a number of individuals, I have measured several behavioral traits in the
wild. Those traits (e.g. home range) can be estimated on different temporal
scales, for example daily, weekly or monthly. I want to estimate repeatability
of those traits, assuming that the daily/weekly/monthly
measurements (averages) represent replicates. I have 3 months (90 days) of
data for each trait.
I have run a MCMCglmm model and extracted the residuals from it using:
res=data$trait-predict(model, marginal=NULL)
And then run:
acf(res)
And I found significant autocorrelation in lags 1 and 2.
Is there any way to account for this (temporal) autocorrelation in MCMCglmm?
Sorry for this pretty basic questions but I haven't found an answer so far.
Thanks!
David
[[alternative HTML version deleted]]
More information about the R-sig-mixed-models
mailing list