[R-sig-ME] Is it possible to account for temporal autocorrelation in MCMCglmm?
Jarrod Hadfield
j.hadfield at ed.ac.uk
Fri Apr 3 09:19:38 CEST 2015
Hi David,
Sorry - this is not possible in MCMCglmm (yet).
Cheers,
Jarrod
Quoting David Villegas Ríos <chirleu at gmail.com> on Wed, 1 Apr 2015
14:34:19 +0200:
> 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]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>
--
The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.
More information about the R-sig-mixed-models
mailing list