应用STATA软件包计算Kappa统计量

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

　Kappa统计量常用于评价两次调查应答结果的重现性，研究不同检验(或测定，诊断)方法的结果间，或不同评定者评定结果间的一致性。加权Kappa值和多个评定者评定结果一致性的Kappa值用手工计算时，计算繁杂，多数统计软件无计算该统计量的程序。STATA软件包提供了计算该统计量的功能。但国内编译的说明书中未介绍，本文简要介绍应用STATA软件包计算Kappa统计量。
　　STATA软件包关于Kappa统计量的命令有三个：kap、kapwgt和kappa。kap命令用于计算两种检验(或测定)方法结果一致性的Kappa统计量；若计算加权Kappa统计量，需先用kapwgt命令为列联表中各格子设定权重值；kappa命令用于两个或两个以上评定者之间评定结果一致性的Kappa统计量的计算，各命令格式如下：
kap　varname1 varname2［weight］［if exp］［in rang］［,tab wgt(wgtid)］
　　kapwgt wgtid［1＼#1［＼##1…］］
　　kappa varlist［if exp］［in range］
　　kap命令中varname1和varname2分别表示两种检验(或测定)方法的变量名称，选择项tab要求输出列联表，选择项wgt要求按照一定的权重计算Kappa统计量，权重值需在kap命令之前由kapwgt命令指定。
　　例1　分别用血清法和滤纸片法测定某抗体血凝试验〔1〕，结果如表1。
表1　某抗体间接血凝试验结果

纸滤片法 　　(m1)	血清法(m2)
	-	+	++	+++
-	8	0	1	0
+	0	10	1	0
++	0	1	26	2
+++	0	1	0	6

　　假定数据已输入计算机内存，m1、m2分别表示滤纸片法和血清法的测定结果，f表示各格子中的例数，输入以下命令：kap m1 m2［freq=f］,tab
　　即可得实际一致率=89.29%、期望一致率=34.18%、Kappa统计量=0.8372、z值=10.12和P值=0.0000。
　　若计算加权Kappa统计量，则需为每个格子设定权数，命令如下：
　　kapwgt weight 1＼0.6667 1＼0.3333 0.6667 1＼0 0.3333 0.6667 1
　　可再用一个kapwgt语句要求STATA输出权重阵，以检查权重是否正确，然后计算加权Kappa统计量，如：
　　kapwgt weight
　　kap m1 m2［freq=f］,tab wgt(weight)
　　输出结果为实际一致率=95.24%，期望一致率=68.37%，Kappa统计量=0.8495，z=9.22，P=0.0000。
　　如果评定者为两人以上，可用Kappa命令计算Kappa值，以下例说明。
　　例2　6名评定者评定13名残疾人上下楼梯的功能〔2〕，残疾人功能分为：不能上下楼梯、能上下但需别人帮助及能独立上下楼梯三个的等级。分别以cat0、cat1和cat2表示对每个残疾人功能评定时，在三个等级上评定结果一致的评定者人数，见表2。
表2　6名评定者评定13名残疾人上下楼梯功能结果

残疾人编号	评定结果一致者人数
	不能做(cat0)	需帮助(cat1)	独立做(cat2)
1	0	0	6
2	6	0	0
12	4	1	1
13	6	0	0

　　假定数据已输入计算机内存，可用以下命令：
　　kappa cat0 cat1 cat2
　　STATA将分别给出各等级及总的一致性Kappa统计量及相应的z值和P值：

结果	Kappa	z	Pr＞Z
cat0	0.8663	12.10	0.0000
cat1	-0.0263	-0.37	0.6434
cat2	0.8445	11.79	0.0000
合计	0.8132	12.20	0.0000

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享1 收藏1 回帖

关键词：kappa统计量 stata软件 kappa Stata tata 说明书 软件包 Kappa kappa 统计

楼主万岁，请教楼主关于kappa值CI的计算方法

正好要用到，太好了～～～