Kappa一致性相关分析中经常应用的Kappa系数有三种,即简单Kappa系数,加权Kappa系数和总Kappa系数及标准误和检验统计量的计算公式,并针对Kappa系数仅适用于行数和列数相等的方表的问题,给出了用SPSS软件实现对行列数不等资料的Kappa检验方法。

1简单Kappa系数的计算公式[1]

K=P0-Pe1-Pe

其中P0=∑ipii,称为观测一致率,Pe=∑ipi.pi,称为期望一致率,即两次检验结果由于偶然机会所造成的一致率,其中pi.=RiN,pi=CiN,Ri,Ci分别为第i个格点所对的行合计和列合计,N为总例数。当两个诊断完全一致时,P0=1,此时Kappa值为1。当观测一致率大于期望一致率时,Kappa值为正数,且Kappa值越大,说明一致性越好。当观察一致率小于期望一致率时,Kappa值为负数,这种情况一般来说比较少见。根据边缘概率的计算,Kappa值的范围值应在-1~1之间。Kappa≥075两者一致性较好;0.75Kappa≥0.4两者一致性一般;Kapp.4两者一致性较差。

Kappa系数标准误的计算公式为:

S=Pe+P2e-∑ipi.piN

其95%的置信区间为:

由于Kappa值是一个样本统计量,作是否有统计学意义的假设检验时,应选用统计量:

U=KappaS

2加权的Kappa系数[2]

加权的Kappa系数是简单Kappa系数的推广,是用加权的方法对两个评价结果进行量化。对于四格表来说,简单Kappa系数与加权的Kappa系数是相等的,对于一般的行列表,加权的Kappa系数的计算公式为:

Kw=P0-Pe1-Pe

P0=∑i∑jwijpij

Pe=∑i∑jwijpi.pj

其中0≤wij=wji1,i≠j,wij=1。

加权Kappa系数的标准误计算公式为:

Skw=∑i∑jpi.pj[wij-]2-P2e)2N

95%的置信区间为:

假设检验的统计量为U=kWSkw

Kappa的权系数一般使用CicchettiAllison和FleissCohen两种权值类型[2],CicchettiAllison的计算公式为:

wij=1- Ci-Cj Ck-Ci

FleissCohen的计算公式为:

wij=1-22

其中,Ci表示第i列的评价分值,k表示列数。假如是数值型变量,评价分值Cij就是第i行第j列对应的具体数值;假如是分类变量,可按照相应级别进行赋值。由于wij=1,而当i≠j时,0≤wij1,所以由Kappa系数的计算公式可知加权的Kappa系数大于简单Kappa系数。

3总Kappa系数

假设列联表为多向列联表,且每个变量有两个水平,即为2×2×…×2列联表,令ki表示第i个变量的Kappa系数,Ski表示第i个Kappa系数的标准误,则总Kappa系数的计算公式为:

K总=∑qi=1kiSki/∑qi=11Ski

若要检验各变量Kappa系数是否都相等,可采用自由度为q-1的χ2检验,计算公式为:

χ2=∑qi=12Ski

、两个公式均适用于加权的Kappa系数。

4行列数不等时Kappa系数的计算

Kappa系数的计算适用于两个评价人分级水平数相同的情况,即数据格式为行数和列数相等的方表。而在实际操作中,经常会出现分级水平数不一致,即行列数不等的情况。我们来看一个实例:两名医生按照某项指标的1~4个等级来评价8个病人。一个医生用全部4个等级进行评价,而另一医生只有3个等级进行评价。此时,对于两个医生来说,他们评价的级别范围不同。数据见表1。

表1甲乙两医生对病人的评价

下面我们使用SPPS软件获得Kappa值及检验结果。SPSS110或更低的版本在这种情况下均无获得Kappa统计量。SPSS115以上的版本可以计算出Kappa值。首先进入数据编辑器并给甲医生添加额外的观测值0001。数据录入见图1。选择Analyze→DescriptiveStatistics→Crosstabs,将变量甲,乙分别放入对应的行列框内,选择Statistics按钮,得到如下对话框,选择Kappa复选框,按continue即可输出Kappa值、标准误和P值。

相应的结果为:

Kappa=0.810,P=0.001,S=0.176。

5讨论

在诊断试验的研究中,数据资料多为双向有序的列联表资料,即两个变量都是有序变量,而且属性相同。属性相同分为三种情况,一种情况是属性、分级水平数和分级水平都完全相同。如甲医生和乙医生都把病人的检查结果分为1、2、3、4四个等级。此时可直接作Kappa检验。当这两个变量都只有2个水平时,就成为配对设计的四格表资料,可使用配对χ2检验,即McNemar检验。第二种是属性相同的分级水平数相同,但分级水平不全相同。如甲医生和乙医生都把病人的检查结果分为四个等级,但甲医生的分级为1、2、3、4,而乙医生的分级为2、3、4、5。在这种情况下,由于列联表的行数和列数仍然是一致的,即列联表仍为方表,所以也可计算出相应的Kappa统计量。第三种是属性相同,但分级水平数和分级水平不全相同。这种情况就是我们所说的列联表的行列数不一致。由于收集上来的数据不能轻易删除掉,所以我们考虑添加行或列使联表成为方表。如行数为n,例数为n-1,则我们只需要添加第n列,在第n行第n列的格点中添加权值0001,而第n行的其它格点均设为0,就可以命名其成为方表,并计算Kappa统计量了。由于权值系数很小,所以不会影响Kappa值的计算结果。

另一方面,假如两个变量中有一个变量是金标准,那么我们不但能分析出检验结果的一致性,还可以计算出敏感度、特异度、误诊率和漏诊率等指标。假如有不同的诊断分界点,还可以绘制出ROC曲线。

诊断试验的评价[3]在医学研究中具有十分重要的意义,目前大多数文献都使用Kappa统计量来检验结果的一致性。所以本研究主要是对Kappa系数作一个探讨和分析。诊断试验评价的统计学方法还会随着更多问题的提出和解决而不断得到发展、修正和扩展。

matlab算kappa系数,[转载]Kappa系数的计算及应用相关推荐

  1. matlab布林线代码,[转载]布林指标的计算(附源码)

    鉴于通达信软件和其他国内证券分析软件布林带指标与国外Advanced GET 8.0 和 Gannalyst Professional 5.0 分析软件布林指标表示差距过大,特别写一篇博文帮助大家理解 ...

  2. FLeiss Kappa系数和Kappa系数的Python实现

    Kappa系数和Fleiss Kappa系数是检验实验标注结果数据一致性比较重要的两个参数,其中Kappa系数一般用于两份标注结果之间的比较,Fleiss Kappa则可以用于多份标注结果的一致性检测 ...

  3. matlab kappa计算,MSA系列9:在使用交叉表法计算Kappa时,期望值是怎么计算出来的?...

    原标题:MSA系列9:在使用交叉表法计算Kappa时,期望值是怎么计算出来的? 在使用如下的交叉表法计算计数型测量系统的Kappa时,需要计算期望值,例如表中的15.7,34.3等数值. 这些期望值是 ...

  4. 洛谷 P1313 计算系数 Label:杨辉三角形 多项式计算

    题目描述 给定一个多项式(by+ax)^k,请求出多项式展开后x^n*y^m 项的系数. 输入输出格式 输入格式: 输入文件名为factor.in. 共一行,包含5 个整数,分别为 a ,b ,k , ...

  5. matlab销量预测的数学模型,数学建模:酒店最优化问题.用matlab算出《酒店价格预测模型》...

    数学建模:酒店最优化问题.用matlab算出<酒店价格预测模型> 39人已浏览 时间 : 2020-10-22 08:23:25 导语:酒店客房管理的数据库模型1.表(客房号,楼,客房类型 ...

  6. 一元二次方程实根java_请依次输入一元二次方程的三个系数,并点击计算显示实根...

    [java]代码库import java.awt.*; import java.awt.event.*; import javax.swing.*; public class Chp13_1 exte ...

  7. 标准差 标准误 偏度系数和峰度系数的作用

    标准差和标准误是两个不同的概念,标准差就是一个变量的所有数据的离均差平方和再平均之后开平方,它是度量离散程度的指标.标准误用于估计参数的可信区间,进行假设检验等.它们与样本含量的关系不同:当样本含量 ...

  8. H.264的两个概念:DC系数和AC系数。 MV预测过程详解(附图)

    在做熵编码之前,先明确两个概念:DC系数和AC系数. 量化后得到的仍是64个系数,量化并没有改变系数的性质.大家知到DCT变换是将数据域从时(空)域变换到频域,在频域平面上变换系数是二维频域变量u和v ...

  9. matlab中右三角形方向,《有限元基础教程》_【MATLAB算例】4.7.1(2) 基于3节点三角形单元的矩形薄板分析(Triangle2D3Node)...

    [MATLAB 算例]4.7.1(2) 基于3节点三角形单元的矩形薄板分析(T riangle2D3Node) 如图4-20所示为一矩形薄平板,在右端部受集中力100 000F N =作用,材料常数为 ...

  10. 数据分析应用统计学之分散性与变异性的测量【极差、四分位差、偏态系数、峰态系数、统计指标】

    文章目录 1.极差与四分位差(R) 2.方差与标准差(西格玛.V) 3.偏态系数与峰态系数(SK.β) 4.统计指标类型 1.极差与四分位差(R) 1)极差:称为全距,极大值减去极小值所得:主要用于衡 ...

最新文章

  1. 怎么样才可批量导出word文档中的照片 不再一张一张的另存为!
  2. javaScript小问题,托拽元素,获得当前鼠标的坐标
  3. 程序员失业第一步?斯坦福研究员用AI从编译器反馈中学习改Bug
  4. 漫画:如何给女朋友解释什么是删库跑路?
  5. VTK与Qt整合的示例
  6. javascript的时间段选择
  7. (转)光照图的理论和实践
  8. watch监控,对比新值和旧值做出相应判断
  9. UI设计灵感|有声读物APP界面设计
  10. GridView 中如何给删除按钮添加提示
  11. Android开发布局 案例一
  12. 时空行为检测数据集 JHMDB UCF101_24 详解
  13. 宏基4752g linux驱动下载,宏基4752g显卡驱动
  14. 开根号的笔算算法图解_开根号手算方法
  15. 001-取整函数(ceil、floor、round)
  16. std::cout彩色输出
  17. 遥想大肠包小肠----python装饰器乱弹
  18. 【软考数据库】第五章 计算机网络
  19. Aria2-突破百度网盘下载限速
  20. The server time zone value xxxxxxx is unrecognized or represents more than one time zone.

热门文章

  1. CST启用GPU加速的调试笔记
  2. cve_2019_0708_bluekeep复现采坑
  3. 条件期望,重期望,相关知识点
  4. 优先队列优化迪杰斯特拉
  5. 极智AI | 目标检测 VOC 格式数据集制作
  6. Flink 开发环境部署和配置
  7. Asp.Net MVC 实用视频教程
  8. 学习笔记(1):Echarts+Asp.Net+Sql Server报表开发视频教程-折线图和曲线图
  9. xp提示计算机内存不足,在xp系统中,经常提示虚拟内存不足的原因是什么?
  10. 静坐常思己过,闲谈莫论人非