比值比(Odds Ratio) ——jamovi软件实现

发布于 2023年1月1日星期日 11:44:23 浏览：1037

原创不易，转载请注明来源，感谢！

附件下载：

比值比.zip 请勿重复点击，如无响应请耐心等待或稍后再试。

比值比(Odds Ratio, OR)是医学研究中常用的效应指标，本文将实例演示在jamovi软件中实现比值比计算的操作步骤。

关键词：jamovi; 比值比; 数据加权; OR值

一、案例介绍

某研究员欲探究脑卒中的危险因素，在某地区抽样调查了345名35岁以上居民脑卒中的患病情况，现将所有调查对象按照年龄是否≥60岁分为两组[中年组(35~59岁)和老年组(≥60岁)]，数据见图1。问老年组的脑卒中患病风险是中年组的多少倍？

二、问题分析

本案例的分析目的是计算老年组的脑卒中患病风险是中年组的多少倍，也就是计算老年组与中年组脑卒中患病的比值比(Odds Ratio，OR)。

本案例以独立样本χ²检验为例计算OR值，但需要满足2个条件：

条件1：自变量和因变量均为二分类变量。本案例中的年龄分组和脑卒中患病情况均为二分类变量，满足该条件。

条件2：独立性，即观测值之间互不干扰。本案例的观测值间相互独立，满足该条件。

三、软件操作及结果解读

(一) 适用条件判断

本案例数据满足上述条件。

(二) 统计描述及推断

1. 软件操作

选择“分析”—“频数”—“独立样本χ²关联检验”(图2)。

将自变量“年龄”、因变量“脑卒中”和“频数”选入右侧对应的变量框内(图3)。

在“单元格”下的“计数”中勾选“实际频数”和“理论频数”，在“百分比”中勾选“列”(图4)。

在“统计”的“检验”中勾选“χ²”，在“比较方法(仅限于2×2四格表)”中勾选“比值比”和“置信区间”(图5)。

在“绘图”中勾选“条形图”“堆积”“百分比”“列内”“列”(图6)。

2. 结果解读

(1) 统计描述

“列联表”(图7)中给出了患脑卒中和未患脑卒中两组人群中老年组(≥60岁)和中年组(35~59岁)的例数及百分比。由结果可知，患脑卒中的人群中老年组占70.7%，未患脑卒中的人群中老年组占41.5%。将两组人群中不同年龄构成可视化为堆积的百分比柱状图，见图8。

(2) 统计学推断

解读χ²检验结果时，遵循以下原则：

①当总例数≥40，且所有期望频数(理论频数)≥5时，可使用Pearson χ²检验。

②当总例数≥40，有1个期望频数≥1且<5，可使用连续校正χ²检验。

③当总例数＜40或有1个期望频数＜1，需要使用Fisher确切概率法。

④无论总例数的大小及期望频数分布情况如何，Fisher确切概率法均可使用。在计算机能够满足运行负荷的情况下，推荐使用Fisher确切概率法。

⑤如果使用Pearson χ²检验或连续校正χ²检验得出的P值比较接近检验水准，建议采用Fisher确切概率法。

由图7可知没有单元格的期望频数少于5，最小期望计数为70.991。同时本案例样本数>40，所以检验结果可采用图9第一行的Pearson χ²检验。由“χ²检验”结果(图9)可知，患脑卒中和未患脑卒中两组人群中年龄构成的差异有统计学意义(χ²=29.471，P<0.001)。由“比较方法”结果(图10)可知，中年组脑卒中的患病风险是老年组的0.294倍[OR=(78×46)/(111×110)]，OR值的95% CI为0.187~0.461。反之，可计算老年组脑卒中的患病风险是中年组的3.403倍[OR=(111×110)/(78×46)]，OR值的95% CI为2.170~5.336。