SPSS实例教程：多重线性回归，你用对了么
在实际的医学研究中，一个生理指标或疾病指标往往受到多种因素的共同作用和影响，当研究的因变量为连续变量时，我们通常在统计分析过程中引入多重线性回归模型，来分析一个因变量与多个自变量之间的关联性。
一、多重线性回归的作用
多重线性回归模型在医学研究领域得到了广泛的应用，其作用主要体现在以下几个方面：
1、探索对于因变量具有影响作用的因素；
2、控制混杂因素，评价多个自变量对因变量的独立效应；
3、用已知的自变量来估计和预测因变量的值及其变化。
二、多重线性回归的形式
多重线性回归的一般表达形式为：
由表达式可以看出，每个因变量的实际测量值yi由两部分组成，即 和ei
之间的差值，表示不能由自变量决定的部分，而对于残差的分析是多重线性回归建模过程中需要重点关注的地方。
此外在多线性回归的表达式中，b0为常数项，表示当所有自变量取值为0时因变量
的估计值；bi为偏回归系数，表示当其他自变量不变时，xi每改变一个单位时所引起的的 变化量。
三、多重线性回归的适用条件
多重线性回归模型作为一种统计模型，它有严格的适用条件，在建模时也需要对这些适用条件进行判断。但是许多使用者往往忽视了这一点，在使用过程中只是单一的构建模型，最终很有可能得出错误的结论。因此在应用多重线性回归之前，我们应该了解它需要满足哪些前提条件呢？
总结起来可用4个词来概况：线性(Linear)，独立(Independence)，正态(Normality)，齐性(Equal variance)，缩写为LINE原则。
(1) 线性：各自变量xi与因变量yi之间存在线性关系，可以通过绘制散点图来进行判断；
(2) 独立：因变量yi的取值之间相互独立，反映到回归模型中，实际上就是要求残差ei之间相互独立；
(3) 正态性：构建多重线性回归模型后，残差ei服从正态分布；
(4) 方差齐性：残差ei的大小不随xi取值水平的变化而变化，即残差ei具有方差齐性。
只有准确把握了LINE核心原则，才能够保证构建符合统计学要求的多重线性回归模型。但是，由于多重线性回归模型具有一定的“抗偏倚性”，如果只是想通过构建方程来探讨自变量和因变量之间的关联性，而非对因变量进行预测，那么后面两个条件可以适当放宽。
此外，还应该注意以下几点：
(5) 因变量yi为连续性变量，而非分类变量；
(6) 自变量xi可以为连续变量，也可以为分类变量，当自变量为多分类无序变量时，则需要设置哑变量，当为有序变量时，则需要根据等级顺序进行赋值。
(7) 对于自变量xi的分布特征没有具体的限定，只要求自变量xi间相互独立，不存在多重共线性；
(8) 对于样本量的要求，根据经验一般要求样本量应当为纳入模型的自变量的20倍以上为宜，比如模型纳入5个自变量，则样本量应当在100以上为宜。
四、SPSS实例操作
1. 研究设计
某研究人员收集了100名研究者的最大摄氧量（VO2max），并记录了他们的年龄，体重，心率和性别，拟探讨年龄，体重，心率和性别对VO2
2. 录入数据
SPSS数据文件如图所示，分析数据结构：
因变量VO2max为连续变量，满足上述条件（5）；
自变量中年龄（age），体重（weight），心率（heart_rate）为连续变量，性别（gender）（女=0，男=1）为二分类变量，满足条件（6）；
样本量为100，纳入的自变量为4个，满足条件（8）中大于20倍的要求。在满足上述几个基本条件后，根据研究目的研究人员拟采用多重线性回归的方法来进行统计分析，而对于其他几点适用条件我们将在后面的介绍中进行一一验证。
3. 多重线性回归操作
(1)选择Analyze → Regression → Linear
在Linear Regression对话框中，将VO2max选入Dependent，将age，weight，heart_rate，gender选入Independent(s)中。
点击Method下拉列表，会出现Enter，Stepwise，Remove，Backward，Forward共5种方法可供选择，这里选择默认的Enter法，表示将所有的变量都纳入到回归模型中。（自变量筛选共有5种方法，每种方法的区别我们将会在以后的内容中进行详细介绍。）
(2)点击Statistic选项
在Regression Coefficients复选框中，勾选Estimates和Confidence Intervals Level（%）并设定为95，可输出自变量的回归系数及其95%可信区间。
选择Model fit，输出模型拟合优度的检验结果；选择Descriptive，输出对所有变量的基本统计描述；选择Part and partial correlations，输出自变量之间的相关系数；选择Collinearity diagnostics，输出对自变量进行共线性诊断的统计量。
在Residus复选框中，选择Durbin-Watson，输出值用于判断残差之间是否相互独立。选择Casewise Diagnositics，默认在3倍标准差内来判定离群点。一般来说，95%的值在 ± 2倍标准差内，99%的值在 ± 2.5倍标准差内，可根据具体情况来进行设定。
(3)点击Save选项
在Predicted Values复选框中选择Unstandardized，保存模型对因变量的原始预测值，在Residuals是复选框中选择Standardized，保存均数为0标准差为1的标准化残差值，在Prediction Intervals复选框中选择Individuals，设定Confidence Intervals为95%，保存个体预测值的95%可信区间。
（数据标准化：用观察值减去该变量的均数，然后除以标准差所得，标准化后数据的均数为0，标准差为1，经标准化的数据都是没有单位的纯数量。）
(4)点击Plot选项
在Plots对话框中将*ZRESID（标准化残差）放入Y轴，将*ZPRED（标准化预测值）放入X轴，绘制残差散点图；同时选择Histogram和Normal probability plot来绘制标准化残差图，考察残差是否符合正态分布；选择Produce all partial plots绘制每一个自变量与因变量残差的散点图。
(5)点击Continue回到Linear Regression主对话框，点击OK完成操作。
