学习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()函数实现趋势性卡方检验

(0)

相关推荐