使用R中的irr包计算Fleiss的kappa

时间:2015-11-19 14:59:36

标签: r categorical-data

我正在研究一个受试者在早年和成年期的标准体重(轻度年龄或年龄偏重)之间是否达成一致。我已经将“年龄轻”重新编码为“1”,重度年龄为“2”,并使用R中的irr包计算Fleiss'kappa。我的数据如下:

phase1=c(2,1,2,1,2,2,2,1,2,2,1,2,2,2,2,2,2,2,2,2,2,2,2,1,1,1,2,1,1,1,2,2,2,2,2,2,2,2,2,2,1,2,2,2,2,2,2,2,2,1,2,2,2,2)
phase3=c(2,1,1,NA,2,1,2,1,NA, NA, 1, NA, 1, NA, 1, 2, NA, NA, 1, 2, NA, NA, 2, 1, 2, 2, 1, 1, 1, 1, NA, NA, 2, 2, 1, 2, NA, NA, 1, NA, NA, 1, 2, 2, 2, 1, 1, 2, 1, 1, 2, NA, NA, NA)
ydata=data.frame(phase1, phase3)

'第1阶段'是早期生命,'第3阶段'是成年期。

这是我的程序:

library(irr)
phases=kappam.fleiss(ydata, exact=TRUE, detail=FALSE)
phases$value

接下来,我想看看早年生活中轻微的年龄('1')是否会在成年期保持轻度年龄:

phases2=kappam.fleiss(ydata[ydata$phase1==1, ], exact=TRUE, detail=FALSE)
phases2$value

根据第二次分析中的数据,我预计会达成非常高的一致性,但是我得到了极低的kappa。所以我想知道我是否可能错误地实现了这一点。

0 个答案:

没有答案