[R-sig-ME] Calculating repeatability with ordinal data
Samantha Patrick
spatrick at cebc.cnrs.fr
Wed May 2 13:11:09 CEST 2012
Hi
Thank you for the tip - I had editted my code! However, sadly it
doesn't solve the problem of how to estimate the residual variance.
Many Thanks
Sam
Le 02/05/2012 12:38, Federico Calboli a écrit :
> 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
>
>
>
--
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/
More information about the R-sig-mixed-models
mailing list