学习R:vcd包Kappa()函数实现一致性检验
我需要引用一个案例,但是我的英语不太行,所以我就不翻译了,直接给案例背景的英文介绍。如下:
Data from Hout et al. (1987) given by Agresti (1990) summarizing the responses of married couples to the questionnaire item: Sex is fun for me and my partner: (a) never or occasionally, (b) fairly often, (c) very often, (d) almost always.
通俗讲就是91对已婚夫妇,丈夫和妻子对sex的打分数据,我想了解一下丈夫和妻子的评价一致性。
先录入数据:
mydata <- matrix(c(7,2,1,2,7,8,5,8,2,3,4,9,3,7,9,14),nrow = 4,ncol = 4,dimnames = list(A=c(1:4),B=c(1:4)))
(a) never or occasionally, (b) fairly often, (c) very often, (d) almost always。咱们一次定义1234,显然是有有序分类变量,有序分类资料研究一致性,主要考虑加权的kappa一致性检验。
现在有米好做饭,直接计算kappa一致性系数,本案例认为1234重要性是等间距的,考虑线性加权。
Kappa(mydata,weights = "Equal-Spacing")
来看结果吧:
第一行是未加权的kappa系数,
第二行是线性加权的kappa系数,也是本例的结果。
结果显示,丈夫和妻子对sex看法的weighted kappa系数为0.2374,P=0.0024<0.01,他们的看法一致性很一般啊。
。
结果同SPSS。
本文完
文/图=数据小兵
更多R统计文章
练习R语言:blandr包实现Bland-Altman一致性评价