I have a (lme) mixed model of this type:
model=lme(trait1~tl_s+year+month_s+I(month_s^2)+I(month_s^3),random=~1|ID,correlation=corAR1(form=~tim),data=depm3,method="REML")
I'm interested in estimating the repeatability of the trait I'm modelling (trait1 in this case). I get a point estimate with this code:
R=as.numeric(VarCorr(mod))[1]/(as.numeric(VarCorr(mod))[1]+as.numeric(VarCorr(mod))[2])
But what would be the code to extract confidence interval for R estimate? I read that parametric bootstrapping would be an option here, but I don't really know how to implement it. I have check some recent posts on the topic like this one...
How to calculate estimated proportions and their confidence intervals from a mixed model?
...but I haven't really found the answer to my question.
Any advise would be appreciated
Thanks