【原】学习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语言：bruceR包Freq()函数频率百分比统计

练习R语言：CrossTable()函数实现交叉表卡方检验

练习R语言：blandr包实现Bland-Altman一致性评价

练习R语言：ls()函数列出变量名称

学习R：matrix()函数录入矩阵数据

学习R语言：CATT()函数实现趋势性卡方检验