关键词:SPSS; 相关分析; Pearson相关分析; 连续变量相关分析; Pearson相关系数
一、案例介绍
某研究者研究健康成年人的体重和双肾脏总体积(ml)的关系,测得24名健康成年人的体重wt (kg)与双肾脏总体积volume (ml),拟探讨健康成年人的体重与双肾体积是否有关。
创建代表受试者编号的变量“ID”,测量尺度设为“名义”,创建代表体重的变量“wt”和代表双肾脏总体积的变量“volume ”,测量尺度均设为“标度”。部分数据见图1。本文案例可从“附件下载”处下载。
二、问题分析
本案例的分析目的是探索两个连续变量之间的相关性,可以使用Pearson相关分析。但需要满足5个条件:
条件1:两个变量均为连续变量。本研究中的成人体重和双肾总体积均为连续变量,该条件满足。
条件2:两个连续变量应当是配对的,即来源于同一个个体。本研究中同时检测了研究对象的体重和双肾总体积,变量之间是对应关系,该条件满足。
条件3:两个连续变量之间存在线性关系,该条件需要通过软件分析后判断。
条件4:两个变量都不存在明显的异常值,该条件需要通过软件分析后判断。
条件5:两个变量呈双变量正态分布或近似正态分布,该条件需要通过软件分析后判断。
三、软件操作及结果解读
(一) 适用条件判断
1. 条件3判断(线性关系分析)
(1) 软件操作
①选择“图形”—“旧对话框”—“散点图/点图…” (图2)。
② 在“散点图/点图…”对话框中选择第一种“简单散点”(图3),然后点击“定义”,出现主对话框。
③在“画简单散点”主对话框中将变量“wt”选入右侧X轴中,将变量“volume”选入右侧Y轴中(图4),点击“确定”则会输出结果。
(2) 结果解读
图5为体重和双肾总体积的散点图,散点大致呈一条直线,说明存在线性关系。该条件满足。
2. 条件4判断(异常值判断)
(1) 软件操作
①选择“分析”—“描述统计”—“探索”(图6)。
② 在“探索”对话框中将变量“wt”和“volume”选入右侧“因变量列表”框(图7)。
③ 点击“图”,在“图”对话框中取消勾选“茎叶图”,勾选“含检验的正态图”,点击“继续”(图8),回到主对话框后点击“确定”。
(2) 结果解读
图9“描述”表格中,列出了观察变量的最小值和最大值,依据专业可判断体重可能存在38 kg和85 kg,双肾总体积可能存在213 ml和348 ml的情况;此外,图10-1和图10-2两个箱线图也未提示任何异常值。综上,本案例未发现需要删除的异常值,满足条件4。
3. 条件5判断(正态性判断)
(1) 软件操作
Pearson相关分析要求变量呈双变量正态分布,而并非简单的要求两个自变量各自服从正态分布。有实验模拟显示,采用Shapiro-Wilk检验同时考察双变量正态分布和两变量分别的正态分布,在重复1万次服从两变量正态分布的实验中若以P=0.1为界值时,约有11%的实验不服从双变量正态分布;若以P=0.05为界值时,约有6%的实验不服从双变量正态分布。由于在SPSS软件中,尚无考察双变量正态分布的程序,且该条件比较宽松,因此可以考察两个自变量各自的正态分布情况代替分析双变量正态分布。正态性检验的操作在异常值判断过程中已完成。
(2) 结果解读
图11显示了两种正态性检验的结果,Kolmogorov-Smirnov (柯尔莫哥洛夫-斯米诺夫,K-S)检验和Shapiro-Wilk (夏皮罗-威尔克,S-W)检验。K-S检验适用于大样本资料,本案查看S-W检验结果,可见P值分别为0.337和0.312,均>0.1,提示两变量均服从正态分布。图12-1和图12-2分别列出了两个变量的Q-Q图,可见散点基本围绕对角线分布,也提示两个变量服从正态分布。综上,本案例满足条件5。
关于正态性检验的注意事项详见推文(医学统计学核心概念及重要假设检验的软件实现(2/4)——正态性假设检验的SPSS实现)。
(二) 统计描述及推断
1. 软件操作
①选择“分析”—“相关”—“双变量” (图13)。
②在“双变量相关性”对话框中将变量“wt”和“volume”选入右侧“变量”框,选择相关系数中的“皮尔逊(Pearson)”,其他设置保持默认不变,如图14所示,点击“确定”。
2. 结果解读
(1) 统计描述
从图9的“描述”表格可知,健康成年人的体重为59.83±14.39 kg,双肾总体积269.78±39.29 ml。
(2) 统计学推断
图15的“相关性”分析结果显示,两组Pearson相关系数r=0.947,P<0.001。可知本研究样本人群的双肾总体积与体重之间存在较强线性相关性。
四、结论
本研究采用Pearson相关分析判断健康成人的双肾总体积与体重是否有关。通过绘制散点图显示两变量之间存在正向线性关系;通过箱线图及专业知识判断,数据不存在需要处理的异常值;通过Q-Q图和Shapiro-Wilk检验,提示两个变量均服从正态分布。
描述性分析结果显示,健康成年人的体重为59.83±14.39 kg,双肾总体积为269.78±39.29 ml。Pearson相关分析显示,健康成人的双肾总体积与体重之间的Pearson相关系数r =0.947 ,P<0.001。综上,可知本研究样本人群的双肾总体积与体重之间存在较强正向线性相关。
五、分析小技巧
(一) 线性关系
- 要确定是否存在线性关系,研究者需要查看两个变量的散点图,如果散点图大致呈一条直线,说明有线性关系,如果不是一条直线则没有线性关系。常见散点图提示的相关类型如图16所示。
(二) 异常值处理
- 如果发现存在异常值,首先要考虑数据是否存在录入错误。如果存在,应修改为正确值。错误数值修订后,需要重新再做散点图,因为修订数值后,可能仍是异常值或者其他的观测值成为了异常值。如果不是录入错误,可以考虑是否为测量误差,如果存在,通常需剔除这些观测值,或将异常值替换为测量范围的上限值(研究者确认异常值超出仪器测量范围并且知道测量错误的方向)。如果异常值不是由于录入或者测量错误导致,这些异常值可能代表了真实的数据情况,尽管这些数据并不符合统计预期,但也没有理由删除,对这种情况的处理方法目前还没有统一标准。研究者可以根据实际情况选择合适的方法,也可以咨询相关专业人员。
(三) 相关系数
- Pearson相关系数r,是表示两个随机变量之间线性相关程度和方向的统计量,即r>0,为正相关;r=0为零相关。R的绝对值大小则表示两变量之间线性相关的密切程度,/r/越接近1,说明关联程度高,/r/越接近于0,说明关联程度越低。当0.9</r/<1,为高度相关;当0.7</r/<0.9,为强相关;0.4</r/<0.7,为中度相关;0.2</r/<0.4,为弱相关性;0</r/<0.2,为极弱相关或无相关性。
- 关于相关系数强弱的判断标准,对于不同专业存在稍微的差异,使用时需要结合专业背景综合决定。此外,Pearson相关系数多用于进行多重线性回归分析前的探索性分析,此时对于相关系数r的判断分级并不是很重要。