前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >卡方检验简介

卡方检验简介

作者头像
一只羊
发布2022-11-30 14:18:57
1.7K0
发布2022-11-30 14:18:57
举报
文章被收录于专栏:生信了生信了

Chi square test(卡方检验)是用于评价两类变量之间是否存在相关性的统计检验方法。

医疗研究会产生大量不同类型的数据,最容易识别的是定量的数据。例如,直腿抬高 (SLR) 的受试者能够将腿抬高大于 0 度,这让我们可以计算两组的平均 SLR,并进行 t 检验。但并不是所有的数据都有这种定量特性。

例如,我们可能对两种治疗后患者的主观改善感兴趣(只使用“是”或“否”回答),而不是测量个体的 SLR。我们能够计算每组的平均改善程度,并做 t 检验吗?答案是否。处理这类数据最为常用的分析方法是 Chi Square 相关性检验。下面是最简单的一个例子。

坐骨神经痛的患者被分成两组,分别使用推拿(SMT)和电牵引(IMT)的方法进行了治疗,治疗的分组情况和病人反馈如下:

在这个例子中,我们的观测值是分类的而非定量的,所以我们应当关注比例而非均值。

注意:

p_1+p_2=q_1+q_2=1

.

我们感兴趣的统计假设总是无事发生(0 假设)。拓展到这个例子就是,

p_1=q_1

p_2=q_2

;即分组 2 中个体的分布不受分组 1 的影响。

为了测试这个假设,我们需要比较假设是真的情况下,期望值和我们实际观测值的差异。

在本例中,我们有 140 个患者认为自己改善了,相对于 390 个总患者来说,改善率为 36%。所以,如果治疗和改善之间没有联系(0 假设),那么对于每一个治疗分组,都应该有 36% 的改善率。

于是有:

注:括号中为 0 假设下的期望值。

获得了期望值之后,需要比较这些值和我们实际观测值之间的差距。

\chi^2=\sum_i \frac{(\mathrm{Observed}_i - \mathrm{Expected}_i)^2 }{\mathrm{Expected}_i}

计算表格如下:

此时,

\chi^2=32.53

根据

\chi^2

的计算公式我们知道,当零假设成立时,

\chi^2

的值会比较小,反之亦然。

接下来的问题是,当

\chi^2

多大时,我们会拒绝 0 假设?

\chi^2

值来自于 Chi Square distribution,这个分布由一个参数决定,即自由度。自由度取决于我们分析的表的大小,可用接下来的公式进行计算。

我们检测的 p-value(任何 2×2 table 的卡方检验),是计算出的卡方值到坐标最右侧曲线下的面积。

查表可知,当卡方值在 6.64 时,p-value 已经小于 0.01。由于我们的值是 32.53,其 p-value 自然小于 0.01。因此,我们拒绝了 0 假设并得出结论:患者接受两种治疗方式的受益是不一样的。

在很多实验中,改善会分多个 levels。例如,让我们对使用热包的脊椎按摩 (Trt 1) 和使用冷包的脊椎按摩 (Trt 2) 治疗急性腰痛进行比较试验。我们使用了 5 个分类来描述改善的状况:

零假设是,两种治疗方式没有差异。

下面计算零假设下的期望值以及最终的卡方值。

此时,自由度为:

(2-1)\times(5-1)=4

自由度为 4 的 Chi Square distribution 如下

卡方为 7.43 时,p-value 是 0.1148。如果我们的显著性水平定为 0.05,则我们无法拒绝零假设。此时,结论是两种治疗手段没有显著的区别。

要进一步解释这一点,请考虑表 8,其中的数据已转换为行百分比:

严格地讲,这些概率分布的比例并不相同。然而考虑到数据中的随机错误,我们没有足够的证据来说明观察到的差异表明了真正的潜在差异。

最后,在使用 \chi^2 检验时,需要遵循一些关键假设,包括了:

  1. 每个个体在表中只出现一次;
  2. 每个个体的结果独立于其他所有个体的结果;
  3. 期望值表中应该有 80% 的期望值大于 5。

参考文献

Ugoni A, Walker BF. The Chi square test: an introduction. COMSIG Rev. 1995 Nov 1;4(3):61-4. PMID: 17989754; PMCID: PMC2050386.

《生信了》2022年10月

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2022-10-21,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信了 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档