computeICCrater: Compute intraclass correlation coefficients ICC(2, 1) and...

In felix-hof/biostatUZH: Misc Tools of the Department of Biostatistics, EBPI, University of Zurich

Compute intraclass correlation coefficients ICC(2, 1) and ICC(3, 1)

Description

Interrater reliability for continuous measurements presented to two or more raters.

Usage

computeICCrater(dat)

Arguments

dat	Data frame that contains the columns score, pat, rater.

Details

Interrater reliability, or more specifically an intraclass correlation coefficient (ICC), is computed when the same continuous measurements are presented to two or more raters. Depending on the underlying model, several different ICCs can be computed. Here, we provide ICC(2, 1) which is derived from a regression model with random effects for rater and patients, and ICC(3, 1) where raters are considered to be fixed.

Value

A list containing the following elements:

sig.pat	Variance of random subject effect.
sig.rater	Variance of the random rater effect.
sig.res	Residual variance.
icc(2, 1)	ICC(2, 1): Intraclass correlation from a two-random effects model.
icc(3, 1)	ICC(3, 1): Intraclass correlation from a model with fixed rater effect.

Note

Since this function relies on a mixed-model formulation using package lme4 to compute ICCs, it can also handle unbalanced data.

Author(s)

Kaspar Rufibach
[e-mail]

References

Shrout, P.E. and Fleiss, J.L. (1979). Intraclass correlations: uses in assessing rater reliability. Psychological Bulletin, 36, 420–428.

Rousson, V., Gasser, T., and Seifert, B. (2002). Assessing intrarater, interrater and test-retest reliability of continuous measurements. Statist. Med., 21, 3431–3446.

See Also

confIntICC

Examples

## generate dataset (balanced data)
set.seed(1977)
n <- 40
r1 <- round(runif(n, 1, 20))
dat <- data.frame(
    "score" = c(r1, r1 + abs(round(rnorm(n, 1, 3))),
        r1 + abs(round(rnorm(n, 1, 3)))), 
    "pat" = rep(c(1:n), 3),
    "rater" = rep(1:3, each = n))
computeICCrater(dat)

## also works for unbalanced data
dat2 <- dat[sort(sample(1:(3 * n))[1:100]), ]
computeICCrater(dat2)

felix-hof/biostatUZH documentation built on Sept. 27, 2024, 1:48 p.m.