[R-sig-ME] Calculating repeatability with ordinal data
Federico Calboli
f.calboli at imperial.ac.uk
Wed May 2 12:38:24 CEST 2012
On 2 May 2012, at 11:08, Samantha Patrick wrote:
>
> Hi
>
> Firstly I sent this message last week but can find no evidence that it actually sent.
> However, if this is a double posting I am very sorry.
>
> I am currently working with repeated measures for individuals and I am
> trying to quantify individual repeatability. Normally, for continuous
> distributions, I would use a mixed model and calculate the variance
> explained by individual divided by the total variance. However my individual scores
> are ordinal and I have been using the clmm function in the Ordinal package:
>
> example of data:
> ID Bird Sex scoremax year
> 622 BS8831 M 2 2008
> 623 BS8831 M 1 2010
> 624 BS8831 M 1 2011
> 625 BS9065 M 1 2010
> 626 BS9065 M 3 2011
> 627 BS19724 F 4 2010
> 628 BS19724 F 5 2010
> 629 BS21302 F 1 2010
> 630 BS25376 F 1 2011
> 631 BS9184 F 2 2009
> 632 BS19989 M 3 2011
> 633 BS21617 M 4 2008
> 634 BS21617 M 2 2009
> 635 BS21617 M 1 2010
>
> where scoremax ranges from 1-5, and there are 1188 birds and 1638
> observations.
>
> scoremax<-as.factor(scoremax)
does
scoremax = as.ordered(scoremax)
make any difference in your results?
I ask this because as.factor() does not, strictly speaking, create an ordered factor.
BW
F
> bird<-as.factor(bird)
> year<-as.factor(year)
> fmm1<- clmm(scoremax~year+ (1|bird), link = c("probit"), Hess =TRUE)
>
> summary(fmm1)
>
> but this only gives the variance estimate for bird, with no residual
> estimate. Some investigations reveal that using an ordinal regression
> in MCMCglmm will also not estimate the residual variance, and it seems
> you need to constrain this value. I have been unable to find any posts
> about repeatability in ordinal data.
>
> My questions I guess are:
> Is using a mixed model appropriate for calculating the repeatability of
> ordinal data (and if not does anyone know any other methods)?
>
> If it is, does anyone have any hints on how to calculate the residual variance,
> to enable repeatability estimates to be calculated.
>
> Many Thanks
>
> Sam
>
> --
>
> Dr Samantha Patrick
> Post Doctoral Fellow
> Centre d'Etudes Biologiques de Chizé - CNRS
> 79360 Villiers-en-Bois
> France
> T:+33 549 097 846
> M:+33 675 603 451
> Skype: sammy_patrick
> http://www.cebc.cnrs.fr/Fidentite/patrick/patrick.htm
> http://www.researchgate.net/profile/Samantha_Patrick/
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
--
Federico C. F. Calboli
Neuroepidemiology and Ageing Research
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG
Tel +44 (0)20 75941602 Fax +44 (0)20 75943193
f.calboli [.a.t] imperial.ac.uk
f.calboli [.a.t] gmail.com
More information about the R-sig-mixed-models
mailing list