第一节 Linear过程

8.1.1 主要功能

调用此过程可完成二元或多元的线性回归分析。在多元线性回归分析中，用户还可根据需要，选用不同筛选自变量的方法（如：逐步法、向前法、向后法，等）。

8.1.2 实例操作

［例8.1］某医师测得10名3岁儿童的身高（cm）、体重（kg）和体表面积（cm²）资料如下。试用多元回归方法确定以身高、体重为自变量，体表面积为应变量的回归方程。

8.1.2.1 数据准备

激活数据管理窗口，定义变量名：体表面积为Y，保留3位小数；身高、体重分别为X1、X2，1位小数。输入原始数据，结果如图8.1所示。

8.1.2.2 统计分析

激活Statistics菜单选Regression中的Linear...项，弹出Linear Regression对话框（如图8.2示）。从对话框左侧的变量列表中选y，点击Ø钮使之进入Dependent框，选x1、x2，点击Ø钮使之进入Indepentdent(s)框；在Method处下拉菜单，共有5个选项：Enter（全部入选法）、Stepwise（逐步法）、Remove（强制剔除法）、Backward（向后法）、Forward（向前法）。本例选用Enter法。点击OK钮即完成分析。

用户还可点击Statistics...钮选择是否作变量的描述性统计、回归方程应变量的可信区间估计等分析；点击Plots...钮选择是否作变量分布图（本例要求对标准化Y预测值作变量分布图）；点击Save...钮选择对回归分析的有关结果是否作保存（本例要求对根据所确定的回归方程求得的未校正Y预测值和标准化Y预测值作保存）；点击Options...钮选择变量入选与剔除的α、β值和缺失值的处理方法。

8.1.2.3 结果解释

在结果输出窗口中将看到如下统计数据：

结果显示，本例以X1、X2为自变量，Y为应变量，采用全部入选法建立回归方程。回归方程的复相关系数为0.94964，决定系数（即r²）为0.90181，经方差分析，F=34.14499，P=0.0003，回归方程有效。回归方程为Y=0.0687101X1+0.183756X2-2.856476。

本例要求按所建立的回归方程计算Y预测值和标准化Y预测值（所谓标准化Y预测值是指将根据回归方程求得的Y预测值转化成按均数为0、标准差为1的标准正态分布的Y值）并将计算结果保存入原数据库。系统将原始的X1、X2值代入方程求Y值预测值（即库中pre_1栏）和标准化Y预测值（即库中zpr_1栏），详见图8.3。

图8.3 计算结果的保存

本例还要求对标准化Y预测值作变量分布图，系统将绘制的统计图送向Chart Carousel窗口，双击该窗口可见下图显示结果。

图8.4 对标准化Y预测值所作的正态分布图

8.2.1 主要功能

调用此过程可完成下列有关曲线拟合的功能：

1、Linear：拟合直线方程（实际上与Linear过程的二元直线回归相同，即Y = b₀+ b₁X）；

2、Quadratic：拟合二次方程（Y = b₀+ b₁X+b₂X²）；

3、Compound：拟合复合曲线模型（Y = b_0×b₁^X）；

4、Growth：拟合等比级数曲线模型（Y = e^(b0+b1X)）；

5、Logarithmic：拟合对数方程（Y = b₀+b₁lnX）

6、Cubic：拟合三次方程（Y = b₀+ b₁X+b₂X²+b₃X³）；

7、S：拟合S形曲线（Y = e^(b0+b1/X)）；

8、Exponential：拟合指数方程（Y = b₀ e^b1X）；

9、Inverse：数据按Y = b₀+b₁/X进行变换；

10、Power：拟合乘幂曲线模型（Y = b₀X ^b1）；

11、Logistic：拟合Logistic曲线模型（Y = 1/（1/u + b_0×b₁^X）。

8.2.2 实例操作

[例8.2]某地1963年调查得儿童年龄（岁）X与锡克试验阴性率（%）Y的资料如下，试拟合对数曲线。

8.2.2.1 数据准备

激活数据管理窗口，定义变量名：锡克试验阴性率为Y，年龄为X，输入原始数据。

8.2.2.2 统计分析

激活Statistics菜单选Regression中的Curve Estimation...项，弹出Curve Estimation对话框（如图8.5示）。从对话框左侧的变量列表中选y，点击Ø钮使之进入Dependent框，选x，点击Ø钮使之进入Indepentdent(s)框；在Model框内选择所需的曲线模型，本例选择Logarithmic模型（即对数曲线）；选Plot models项要求绘制曲线拟合图；点击Save...钮，弹出Curve Estimation:Save对话框，选择Predicted value项，要求在原始数据库中保存根据对数方程求出的Y预测值，点击Continue钮返回Curve Estimation对话框，再点击OK钮即可。

8.2.2.3 结果解释

在结果输出窗口中将看到如下统计数据：

在以X为自变量、Y为应变量，采用对数曲线拟合方法建立的方程，决定系数R²=0.913（接近于1），作拟合优度检验，方差分析表明：F=52.32，P=0.001，拟合度很好，对数方程为：Y=61.3259+20.6704lnX。

本例要求绘制曲线拟合图，结果如图8.6所示。

图8.6 对数曲线拟合情形

图8.7 计算结果的保存

第三节 Logistic过程

8.3.1 主要功能

调用此过程可完成Logistic回归的运算。所谓Logistic回归，是指应变量为二级计分或二类评定的回归分析，这在医学研究中经常遇到，如：死亡与否（即生、死二类评定）的概率跟病人自身生理状况和所患疾病的严重程度有关；对某种疾病的易感性的概率（患病、不患病二类评定）与个体性别、年龄、免疫水平等有关。此类问题的解决均可借助逻辑回归来完成。

特别指出，本节介绍的Logistic过程，应与日常所说的Logistic曲线模型（即S或倒S形曲线）相区别。用户如果要拟合Logistic曲线模型，可调用本章第二节Curve Estimation过程，系统提供11种曲线模型，其中含有Logistic曲线模型（参见上节）。

在一般的多元回归中，若以P（概率）为应变量，则方程为P=b0+b1X1+b2X2+…+bkXk,

但用该方程计算时，常会出现P>1或P<0的不合理情形。为此，对P作对数单位转换，即logitP=ln(P/1-P)，于是，可得到Logistic回归方程为：

e^{b0+b1X1+b2X2+…+bkXk}

P = ———————————

1+ e^{b0+b1X1+b2X2+…+bkXk}

8.3.2 实例操作

[例8.3]某医师研究男性胃癌患者发生术后院内感染的影响因素，资料如下表，请通过Logistic回归统计方法对主要影响因素进行分析。

8.3.2.1 数据准备

激活数据管理窗口，定义变量名：术后感染为Y（字符变量，有输入Y、无输入N），年龄为X1，手术创伤程度为X2，营养状态为X3，术前预防性抗菌为X4（字符变量，有输入Y、无输入N），白细胞数为X5，癌肿病理分度为X6。按要求输入原始数据。

8.3.2.2 统计分析

激活Statistics菜单选Regression中的Logistic...项，弹出Logistic Regression对话框（如图8.8示）。从对话框左侧的变量列表中选y，点击Ø钮使之进入Dependent框，选x1、x2、x3、x4、x5和x6，点击Ø钮使之进入Covariates框；点击Method处的下拉按钮，系统提供7种方法：

1、Enter：所有自变量强制进入回归方程；

2、Forward: Conditional：以假定参数为基础作似然比概率检验，向前逐步选择自变量；

3、Forward: LR：以最大局部似然为基础作似然比概率检验，向前逐步选择自变量；

4、Forward: Wald：作Wald概率统计法，向前逐步选择自变量；

5、Backward: Conditional：以假定参数为基础作似然比概率检验，向后逐步选择自变量；

6、Backward: LR：以最大局部似然为基础作似然比概率检验，向后逐步选择自变量；

7、Backward: Wald：作Wald概率统计法，向后逐步选择自变量。

本例选用Forward: Conditional法，以便选择有主要作用的影响因素；点击Options...钮，弹出Logistic Regression:Options对话框，在Display框中选取At last step项，要求只显示最终计算结果，点击Continue钮返回Logistic Regression对话框，再点击OK钮即可。

8.3.2.3 结果解释

在结果输出窗口中将看到如下统计数据：

系统先对字符变量进行重新赋值，对于应变量Y，回答是（Y）的赋值为0，回答否（X）的赋值为1；对于应变量X4，回答是（Y）的赋值为-1，回答否（X）的赋值为1。

结果表明，第一步自变量X3入选，方程分类能力达80.00%；第二步自变量X6入选，方程分类能力达93.33%（参见结果中的分类分析表）；方程有效性经χ²检验，χ²=15.276，P=0.0005。

Logistic回归的分类概率方程为：

e^{123.4053-30.5171X3-10.2797X6}

P = ——————————————

1+ e^{123.4053-30.5171X3-10.2797X6}

根据该方程，若一胃癌患者营养状态评分（X3）为3，癌肿病理分度（X6）为9，则其P=4.5×10^-27≈0，这意味着术后将发生院内感染；另一胃癌患者营养状态评分（X3）为1，癌肿病理分度（X6）为4，则其P=0.98105≈1，这意味着术后将不会发生院内感染。

第五节 Nonlinear过程

8.5.1 主要功能

调用此过程可完成非线性回归的运算。所谓非线性回归，即为曲线型的回归分析，一些曲线模型我们已在本章第二节中述及。但在医学研究中经，还经常会遇到除本章第二节中述及的曲线模型，对此，SPSS提供Nonlinear过程让用户根据实际需要，建立各种曲线模型以用于研究变量间的相互关系。在医学中，如细菌繁殖与培养时间关系的研究即可借助Nonlinear过程完成。

下面一些曲线模型是在论文中较常见的，提供给用户应用时作参考：

模型名称 模型表达式

Asympt. Regression 1 Y = b1 + b2×exp( b3 ^X )

Asympt. Regression 2 Y = b1 -( b2 × ( b3 ^X ))

Density Y = ( b1 + b2 ×X ) ^{(-1/ b3 )}

Gauss Y = b1× (1- b3×exp( -b2 ×X ²))

Gompertz Y = b1×exp( -b2 ×exp( -b3 ×X ))

Johnson-Schumacher Y = b1×exp( -b2 / ( X + b3))

Log Modified Y = ( b1 + b3×X )^b2

Log-Logistic Y = b1 -ln(1+ b2 ×exp( -b3×X ))

Metcherlich Law of Dim. Ret. Y = b1 + b2×exp( -b3×X )

Michaelis Menten Y = b1×X /( X + b2 )

Morgan-Mercer-Florin Y = ( b1×b2 + b3×X ^b4 )/( b2 + X ^b4 )

Peal-Reed Y = b1 /(1+ b2 ×exp(-( b3×X + b4×X ²+ b5×X ³ )))

Ratio of Cubics Y = ( b1 + b2×X + b3×X ² + b4×X ³ )/( b5×X ³ )

Ratio of Quadratics Y = ( b1 + b2×X + b3×X² )/( b4×X ² )

Richards Y = b1 /((1+ b3×exp(- b2×X )) ^{(1/ b4 )} )

Verhulst Y = b1 /(1 + b3×exp(- b2×X ))

Von Bertalanffy Y = ( b1 ^{(1 - b4 )} - b2×exp( -b3×X )) ^{(1/(1 - b4 )})

Weibull Y = b1 - b2×exp(- b3×X ^b4 )

Yield Density Y = (b1 + b2×X + b3×X ² )^(-1)

8.5.2 实例操作

[例8.5]选取某地某年寿命表中40-80岁各年龄组的尚存人数资料如下表，请就该资料试拟合Gompertz曲线（Y = b1×b2^(b3X)）。

8.5.2.1 数据准备

激活数据管理窗口，定义变量名：年龄简化值为X，尚存人数为Y。输入原始数据。

8.5.2.2 统计分析

激活Statistics菜单选Regression中的Nonlinear...项，弹出Nonlinear Regression对话框（如图8.11示）。从对话框左侧的变量列表中选y，点击Ø钮使之进入Dependent框。由于SPSS系统尚无法智能地自动拟合用户所需的曲线，故一方面要求用户估计方程中常数项和各系数项进行叠代运算的起始值，另一方面要求用户列出方程模型。对此，可首先点击Nonlinear Regression对话框的Parameters...钮，弹出Nonlinear Regression: Parameters对话框（图8.12），在Name处定义系数名，在Start Value处输入起始值（这项工作是十分重要的，否则系统可能无法运算，甚至会因叠代次数过大导致SPSS系统的崩溃），本例定义b1=8500、b2=1、b3=1.5，每定义一个系数，即点击Add钮加以确定；若在后面的运算中出错，则还可修改系数项的起始值，修改后点击Change钮加以确定；然后点击Continue钮返回Nonlinear Regression对话框。在Model Expression处写出曲线方程表达式，用户可借助系统提供的数码盘和函数列表写出方程。本例要求计算根据回归方程求出的预测值，可点击Save钮，在Nonlinear Regression:Save New Variables对话框中选Predicted value项。最后点击OK钮即可。

8.12 系数项定义对话框

8.5.2.3 结果解释

在结果输出窗口中将看到如下统计数据：

经30次叠代运算后，相邻两次的方程剩余均方差值不大于规定的1×10^-8，满足要求；回归方程的决定系数R² = 0.99983，Gompertz曲线方程为：

Y = 85679.247671×0.94933859^{(1.469036403X)}

本例要求计算预测值，系统将结果存入原始数据库中（图8.13），系统以pred_作为预测值的变量名。由结果可见，预测值与实际值十分接近。

图8.13 原始数据及其预测值