[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