SPSS详细教程:相对危险度(RR)的计算
一、问题与数据
研究者想探索吸烟与肺癌间的关联,即吸烟者患肺癌的风险是否比不吸烟者高以及高多少。研究者从一般人群中随机抽样700名调查对象建立前瞻性队列研究。调查对象的吸烟状态为吸烟或不吸烟。在研究开始时,所有调查对象均未患肺癌。研究者随访10年记录调查对象是否患肺癌。
因此,研究者可以分别得到患肺癌的吸烟者、未患肺癌的吸烟者、患肺癌的不吸烟者和未患肺癌的不吸烟者的人数,据此可以计算吸烟组与不吸烟组的风险差异。
本例中代表吸烟状态的变量为smoking_status,吸烟赋值为1,不吸烟赋值为2;代表随访期肺癌发病与否的变量为lung_cancer,患肺癌赋值为1,未患肺癌赋值为2。
部分数据截图如下,左侧为原始数据,右侧为按不同吸烟状态和患肺癌与否统计的汇总数据。
二、对问题的分析
为计算相对危险度(RR),需要满足以下假设:
假设1:因变量和自变量均为二分类变量。
假设2:各观测间相互独立。
接下来将详细介绍如何计算相对危险度。
三、SPSS操作
1. 数据准备
如果研究者使用原始数据,则跳过数据准备步骤,直接计算相对危险度;如果使用按不同吸烟状态和患肺癌与否来统计得到汇总数据,则需要添加权重,步骤如下。
(1)点击主菜单Data > Weight Cases,如下图:
点击后出现Weight Cases对话框,如下图:
(2)勾选Weight cases by选项,激活 键和Frequency Variable: 框,如下图:
(3)将变量freq选入Frequency Variable框,如下图:
(4)点击OK键,为数据加权。
2. 相对危险度的SPSS操作
(1)点击主菜单Analyze > Descriptive Statistics > Crosstabs,如下图:
点击后出现Crosstabs对话框,如下图:
注意:如果使用频数统计的数据文件,Crosstabs对话框如下图:
(2)将自变量smoking_status选入Row(s)框,因变量lung_cancer选入Column(s)框,如下图:
注意:如果使用频数统计的数据文件,Crosstabs对话框如下图:
(3)点击Statistics键,出现Crosstabs:Statistics对话框,如下图:
(4)勾选Risk,如下图:
(5)点击Continue键,返回Crosstabs对话框。
(6)点击Cells键,出现Crosstabs:Cell Display对话框,如下图:
(7)勾选Percentanges区域的Row,如下图:
(8)点击Continue键,返回Crosstabs对话框。
(9)点击OK键,生成结果。
四、结果解释
1. 描述性统计分析
在报告相对危险度前,研究者应该先查看基本的一些统计量,了解数据特征。本例中在smoking_status*lung_cancer Crosstabulation表,如下图:
首先,表中可看到在350名吸烟者中患肺癌的人数,如下图高亮:
随访期间29名吸烟者患肺癌,即8.3%(29÷350*100%=8.3%)的吸烟者患肺癌。
注意:表中仅保留了一位小数,SPSS会自动计算相对危险度;如果研究者想得到更多小数位数手动计算的相对危险度,则双击smoking_status*lung_cancer Crosstabulation表,如下图。可观察到8.3%变为8.285714%。
其次,还能看到350名不吸烟者中患肺癌的人数,如下图:
随访期间9名不吸烟者患肺癌,即2.6%(9÷350*100%=2.6%)的不吸烟者患肺癌。研究者想得到更多小数位数手动计算的相对危险度,也可以如上操作,2.6%变为2.571429%,如下图:
从smoking_status*lung_cancer Crosstabulation表中可以得到初步结论:与不吸烟者相比,吸烟者患肺癌的风险更高。
研究者可以汇报:本研究随机抽样700名调查对象,吸烟者和不吸烟者各350名,吸烟者患肺癌的比例(8.3%)高于不吸烟者(2.6%)。
当然,研究者也可以选择分别汇报患肺癌的吸烟者、未患肺癌的吸烟者、患肺癌的不吸烟者和未患肺癌的不吸烟者的人数。
2. 相对危险度
相对危险度结果可以在Risk Table表中查看,如下图高亮显示。
注意:表中高亮显示的为“For cohort lung_cancer=Yes”行,不是“For cohort lung_cancer=No”,因为研究者感兴趣的是吸烟对患肺癌的影响。
吸烟者患肺癌的相对危险度是3.222,95%置信区间为1.548到6.707。95%置信区间说明研究者有95%的把握说明人群中真实的相对危险度在1.548和6.707间。
如果研究者手动计算相对危险度,则需要吸烟者患肺癌的比例即0.08285714和不吸烟者患肺癌的比例即0.02571429。将两个比例输入到相对危险度计算公式,如下图:
如果相对危险度大于1且95%置信区间不包括1,说明暴露组发生疾病的风险高于非暴露组,是结局的危险因素。本例中相对危险度3.222说明吸烟者患肺癌的风险高出不吸烟者222%。研究者还可以计算超额危险度,公式如下,其中RR为相对危险度。
本例中,超额危险度则为:
如果相对危险度小于1且95%置信区间不包括1,则意味着暴露组患病风险降低,是结局的保护因素。例如,相对危险度为0.34,代表与非暴露组相比,暴露组研究对象发生疾病的风险降低66%,计算如下。负号代表发生疾病的风险降低。
如果相对危险度的95%置信区间包括1,说明暴露组与非暴露组发生疾病的风险差异无统计学意义。
提示:研究者必须把自变量smoking_status选入Row(s),因变量lung_cancer选入Column(s)。如果选反,则会出现如下结果。
上表中相对危险度为1.574,95%置信区间上下限也分别变为1.297和1.910。因此,研究者需注意选对自变量和因变量。
综上,研究者可以汇报:本研究随机抽样700名调查对象,吸烟者和不吸烟者各350名,吸烟者患肺癌的比例(8.3%)高于不吸烟者(2.6%)。与不吸烟者相比,吸烟者患肺癌的相对危险度为3.222(95%置信区间为1.548-6.707)。
五、撰写结论
如果相对危险度95%置信区间不包括1,可以汇报:
本研究随机抽样700名调查对象,吸烟者和不吸烟者各350名,吸烟者患肺癌的比例(8.3%)高于不吸烟者(2.6%)。与不吸烟者相比,吸烟者患肺癌的相对危险度为3.222(95%置信区间为1.548-6.707),且有统计学意义。
如果相对危险度95%置信区间包括1,可以汇报(另外列举一个例子):
有900名永久性房颤病人随机分配至两组,每组450名,一组使用降低脑卒中风险的新药,另外一组使用常规药物。使用新药的房颤病人中11人(2.4%)患脑卒中,使用常规药物病人24人(4.7%)患脑卒中。与使用常规药物的房颤病人相比,使用新药的病人患脑卒中的相对危险度为0.524(95%置信区间为0.256-1.074),没有统计学意义。
数据分析咨询请扫描二维码
若不方便扫码,搜微信号:CDAshujufenxi
在当今这个数据驱动的时代,几乎每一个业务决策都离不开对数据的深入分析。而其中,指标波动归因分析更是至关重要的一环。无论是 ...
2025-02-18当数据开始说谎:那些年我们交过的学费 你有没有经历过这样的场景?熬了三个通宵做的数据分析报告,在会议上被老板一句"这数据靠 ...
2025-02-17数据分析作为一门跨学科领域,融合了统计学、编程、业务理解和可视化技术。无论是初学者还是有一定经验的从业者,系统化的学习路 ...
2025-02-17挖掘用户价值本质是让企业从‘赚今天的钱’升级为‘赚未来的钱’,同时让用户从‘被推销’变为‘被满足’。询问deepseek关于挖 ...
2025-02-17近来deepseek爆火,看看deepseek能否帮我们快速实现数据看板实时更新。 可以看出这对不知道怎么动手的小白来说是相当友好的, ...
2025-02-14一秒精通 Deepseek,不用找教程,不用买资料,更不用报一堆垃圾课程,所有这么去做的,都是舍近求远,因为你忽略了 deepseek 的 ...
2025-02-12自学 Python 的关键在于高效规划 + 实践驱动。以下是一份适合零基础快速入门的自学路径,结合资源推荐和实用技巧: 一、快速入 ...
2025-02-12“我们的利润率上升了,但销售额却没变,这是为什么?” “某个业务的市场份额在下滑,到底是什么原因?” “公司整体业绩 ...
2025-02-08活动介绍 为了助力大家在数据分析领域不断精进技能,我们特别举办本期打卡活动。在这里,你可以充分利用碎片化时间在线学习,让 ...
2025-02-071、闺女,醒醒,媒人把相亲的带来了。 我。。。。。。。 2、前年春节相亲相了40个, 去年春节相亲50个, 祖宗,今年你想相多少个 ...
2025-02-06在数据科学的广阔领域中,统计分析与数据挖掘占据了重要位置。尽管它们常常被视为有关联的领域,但两者在理论基础、目标、方法及 ...
2025-02-05在数据分析的世界里,“对比”是一种简单且有效的方法。这就像两个女孩子穿同一款式的衣服,效果不一样。 很多人都听过“货比三 ...
2025-02-05当我们只有非常少量的已标记数据,同时有大量未标记数据点时,可以使用半监督学习算法来处理。在sklearn中,基于图算法的半监督 ...
2025-02-05考虑一种棘手的情况:训练数据中大部分样本没有标签。此时,我们可以考虑使用半监督学习方法来处理。半监督学习能够利用这些额 ...
2025-02-04一、数学函数 1、取整 =INT(数字) 2、求余数 =MOD(除数,被除数) 3、四舍五入 =ROUND(数字,保留小数位数) 4、取绝对值 =AB ...
2025-02-03作者:CDA持证人 余治国 一般各平台出薪资报告,都会哀嚎遍野。举个例子,去年某招聘平台发布《中国女性职场现状调查报告》, ...
2025-02-02真正的数据分析大神是什么样的呢?有人认为他们能轻松驾驭各种分析工具,能够从海量数据中找到潜在关联,或者一眼识别报告中的数 ...
2025-02-01现今社会,“转行”似乎成无数职场人无法回避的话题。但行业就像座围城:外行人看光鲜,内行人看心酸。数据分析这个行业,近几年 ...
2025-01-31本人基本情况: 学校及专业:厦门大学经济学院应用统计 实习经历:快手数据分析、字节数据分析、百度数据分析 Offer情况:北京 ...
2025-01-3001专家简介 徐杨老师,CDA数据科学研究院教研副总监,主要负责CDA认证项目以及机器学习/人工智能类课程的研发与授课,负责过中 ...
2025-01-29