全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 经管代码库
7587 6
2015-02-15
 Kappa统计量常用于评价两次调查应答结果的重现性,研究不同检验(或测定,诊断)方法的结果间,或不同评定者评定结果间的一致性。
      加权Kappa值和多个评定者评定结果一致性的Kappa值用手工计算时,计算繁杂,多数统计软件无计算该统计量的程序。STATA软件包提供了计算该统计量的功能。但国内编译的说明书中未介绍,本文简要介绍应用STATA软件包计算Kappa统计量。

  STATA软件包关于Kappa统计量的命令有三个:kap、kapwgt和kappa。
      kap命令用于计算两种检验(或测定)方法结果一致性的Kappa统计量;若计算加权Kappa统计量,需先用kapwgt命令为列联表中各格子设定权重值;kappa命令用于两个或两个以上评定者之间评定结果一致性的Kappa统计量的计算,各命令格式如下:

复制代码

  kap命令中varname1和varname2分别表示两种检验(或测定)方法的变量名称,选择项tab要求输出列联表,选择项wgt要求按照一定的权重计算Kappa统计量,权重值需在kap命令之前由kapwgt命令指定。
  


      例1 分别用血清法和滤纸片法测定某抗体血凝试验〔1〕,结果如表1。

      表1 某抗体间接血凝试验结果

纸滤片法
  (m1)
血清法(m2)
-++++++
-8010
+01010
++01262
+++0106
 

 假定数据已输入计算机内存,m1、m2分别表示滤纸片法和血清法的测定结果,f表示各格子中的例数,输入以下命令:
复制代码


  即可得实际一致率=89.29%、期望一致率=34.18%、Kappa统计量=0.8372、z值=10.12和P值=0.0000。


  若计算加权Kappa统计量,则需为每个格子设定权数,命令如下:
  
复制代码


  可再用一个kapwgt语句要求STATA输出权重阵,以检查权重是否正确,然后计算加权Kappa统计量,如:
  
复制代码


  输出结果为实际一致率=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)
1006
2600
12411
13600
  假定数据已输入计算机内存,可用以下命令:
  
复制代码

  STATA将分别给出各等级及总的一致性Kappa统计量及相应的z值和P值
结果KappazPr>Z
cat00.866312.100.0000
cat1-0.0263-0.370.6434
cat20.844511.790.0000
合计0.813212.200.0000
二维码

扫码加我 拉你入群

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

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

全部回复
2015-10-9 11:33:46
谢谢分享。
二维码

扫码加我 拉你入群

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

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

2019-3-29 17:14:21

谢谢分享。
二维码

扫码加我 拉你入群

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

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

2019-4-21 20:44:12
赞👍
二维码

扫码加我 拉你入群

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

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

2019-4-28 12:52:25


谢谢分享。
二维码

扫码加我 拉你入群

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

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

2019-12-12 23:16:38
请问为每个格子设定权数那步是怎么设定呢?
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群