[R-sig-ME] Using lmer to determine inter-rater reliability
Kevin E. Thorpe
kevin.thorpe at utoronto.ca
Mon Aug 31 17:53:41 CEST 2015
Hello.
I have a data frame with following structure.
> str(vision)
'data.frame': 268 obs. of 9 variables:
$ Child : Factor w/ 67 levels "C01-05","C01-10",..: 43 43 43
43 44 44 44 44 42 42 ...
$ Test : Factor w/ 4 levels "1","2","3","4": 1 2 3 4 1 2 3 4
1 2 ...
$ Rater : Factor w/ 4 levels "F","L","P","S": 4 1 4 1 1 4 1 4
4 1 ...
$ Binoc : int 100 100 100 100 0 0 0 0 40 0 ...
$ Yield : int 100 100 100 80 100 20 50 30 100 100 ...
$ Tries : int 5 5 5 6 5 10 10 10 5 5 ...
$ Result : Factor w/ 5 levels "Pass","ReferBI",..: 1 1 1 1 3 2
3 2 3 3 ...
$ ResultCollapsed: Factor w/ 3 levels "Pass","Refer",..: 1 1 1 1 2 2 2
2 2 2 ...
$ Test1 : Factor w/ 16 levels "F:1","F:2","F:3",..: 13 2 15 4
1 14 3 16 13 2 ...
In these data, each subject is rated by 2 (of 4) raters twice. The Test1
variable was created from Test and Rater with
(Rater:Test)[drop=TRUE] to explicitly create the nesting.
I then fit the following model.
> binoc.lmer1 <- lmer(Binoc~1+(1|Child) + (1|Rater) +
(1|Test1),data=vision)
> binoc.lmer1
Linear mixed model fit by REML ['lmerMod']
Formula: Binoc ~ 1 + (1 | Child) + (1 | Rater) + (1 | Test1)
Data: vision
REML criterion at convergence: 2592.62
Random effects:
Groups Name Std.Dev.
Child (Intercept) 29.226
Test1 (Intercept) 2.292
Rater (Intercept) 5.823
Residual 26.330
Number of obs: 264, groups: Child, 66; Test1, 16; Rater, 4
Fixed Effects:
(Intercept)
51.68
Now my questions.
1. Have I fit the right model?
2. If so, would the right estimate for the rater ICC be
Rater/(Rater + Residual)
or
(Rater + Test1)/(Rater + Test1 + Residual)
3. Would Test1/(Test1 + Residual) give an estimate of intra-rater
reliability?
Thanks for you time.
Kevin
--
Kevin E. Thorpe
Head of Biostatistics, Applied Health Research Centre (AHRC)
Li Ka Shing Knowledge Institute of St. Michael's
Assistant Professor, Dalla Lana School of Public Health
University of Toronto
email: kevin.thorpe at utoronto.ca Tel: 416.864.5776 Fax: 416.864.3016
More information about the R-sig-mixed-models
mailing list