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

计算机应用 应用STATA软件包计算Kappa统计量 山东医科大学(210012) 刘言训 王洁贞 王 萍 韩 兢 Kappa统计量常用于评价两次调查应答结果的重 现性,研究不同检验(或测定,诊断)方法的结果间,或 不同评定者评定结果间的一致性加权Kappa值和多 个评定者评定结果一致性的Kappa值用手工计算时, 计算繁杂,多数统计软件无计算该统计量的程序 STA TA软件包提供了计算该统计量的功能但国内 编译的说明书中未介绍,本文简要介绍应用STA TA 软件包计算Kappa统计量 STA TA软件包关于Kappa统计量的命令有三 个: kap、kapw gt和kappakap命令用于计算两种检验 (或测定)方法结果一致性的Kappa统计量;若计算加 权Kappa统计量,需先用kapw gt命令为列联表中各 格子设定权重值; kappa命令用于两个或两个以上评 定者之间评定结果一致性的Kappa统计量的计算,各 命令格式如下: kap varname1 varname2 [w eight ] [ if exp ] [ in rang][, tab w gt(w gtid ) ] kapw gt w gtid[1ł# 1[ł# # 1⋯]] kappa varlist[if exp ][in range] kap命令中varname1和varname2分别表示两种 检验(或测定)方法的变量名称,选择项tab要求输出 列联表,选择项w gt要求按照一定的权重计算Kappa 统计量,权重值需在kap命令之前由kapw gt命令指 定。

例1 分别用血清法和滤纸片法测定某抗体血凝 试验 〔1〕 ,结果如表1 表1 某抗体间接血凝试验结果 纸滤片法 (m1) 血清法(m2) -++ ++ + + -8010 +01010 + +01262 + + +0106 假定数据已输入计算机内存,m1、m2分别表示滤 纸片法和血清法的测定结果,f表示各格子中的例数, 输入以下命令: kapm1m2[freq= f], tab 即可得实际一致率= 89129%、 期望一致率= 34118%、 Kappa统 计 量= 018372、z值= 10112和P值= 010000 若计算加权Kappa统计量,则需为每个格子设定 权数,命令如下: kapw gt w eight 1ł016667 1ł013333 016667 1ł0 013333 016667 1 可再用一个kapw gt语句要求STA TA输出权重 阵,以检查权重是否正确,然后计算加权Kappa统计 量,如: kapw gt w eight kapm1m2[freq=f], tab w gt(w eight) 输出结果为实际一致率= 95124% ,期望一致率= 68137% , Kappa统 计 量= 018495,z= 9122,P= 010000。

如果评定者为两人以上,可用Kappa命令计算 Kappa值,以下例说明 例2 6名评定者评定13名残疾人上下楼梯的功 能 〔2〕 ,残疾人功能分为:不能上下楼梯、 能上下但需别 人帮助及能独立上下楼梯三个的等级分别以cat0、 cat1和cat2表示对每个残疾人功能评定时, 在三个等 级上评定结果一致的评定者人数,见表2 表2 6名评定者评定13名残疾人上下楼梯功能结果 残疾人编号 评定结果一致者人数 不能做( cat0) 需帮助( cat1) 独立做( cat2) 1006 2600 ???? 12411 13600 假定数据已输入计算机内存,可用以下命令: kappa cat0 cat1 cat2 STA TA将分别给出各等级及总的一致性 Kappa统计量及相应的z值和P值: 结果KappazP rZ cat0 01866312110010000 cat1- 010263 - 0137016434 cat2 01844511179010000 合计01813212120010000 参 考 文 献 11 王洁贞.赵跃进等. Kappa统计量及其应用.中国卫生统计, 1995, 12 (6): 49. 21 栾建安,王纪宪等.多类别多评估者的Kappa分析.中国卫生统计, 1995, 12(6): 20. 381 中国卫生统计1999年6月第16卷第3期 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved. 。