第八课征服一般线性模型general linear model菜单详解上

1、SPSS 第八课：征服一般线性模型General Linear M菜单详解（上）（医学统计：）上次更新日期：两方差分析univarate框界面说明8.1.2结果解释协方差分析分析步骤8.2.2结果解释8.3 其他较简单的方差分析问题8.48.4.18.4.28.58.5.18.5.2多元方差分析分析步骤 结果解释重复测量的方差分析 Repeated measures结果解释框界面说明请注意，本章的标题用了一些修辞手法，一般线性模型可不是用一章就可以说清楚的，因为它包括的内容实在太多了。那么，究竟用到的哪些分析会包含在其中呢？简而言之：凡是和方差分析粘边的都可以用他来做。比如成组设计的方差分析（

2、即单方差分析）、配伍设计的方差分析（即两方差分析）、交叉设计的方差分析、析因设计的方差分析、重复测量的方差分析、协方差分析等等。因此，能真正掌握 GLM 菜单的用使大家的统计分析能力有极大地提高。实际上一般线性模型包括的统计模型还不止这些，我这里举出来的只是从用SPSS 作统计分析的角度而言的一些。好了，既然一般线性模型的能力如此强大，那么下属的四个子菜单各自的功能是什么呢？请看：Univariate 子菜单：四个菜单中的大哥大，绝大部分的方法分析都在这里面进行。Multivariate 子菜单：当结果变量（应变量）不止一个时，当然要用他来分析啦！Repeted Measures 子菜单：顾名

3、思义，重复测量的数据就要用他来分析，这一点我可能要强调一下，用前两个菜单似乎都可以分析出来结果，但在许多情况下该结果是不正确的，应该用重复测量的分析方法才对（不能再讲了，再讲下去就会扯到多水平模型去了）。Variance Components 子菜单：用于作方差成份模型的，这个模型实在太深，不是一时半会说的请的，所以我在这里就干脆不讲了。出于模型复杂性、篇幅、应用范围及乱七八糟一系列的理由，当然主要是 我懒得一一解释，我决定本章采用举例讲解的方式，及讲解一些常见的分析实例，通过这种方法来熟悉那些最为常用的分析方法。对统计分析的数据格式不太熟悉的朋友，请一定先去看看统计第一课： 论统计中的数据录

4、入格式，会大有帮助的。8.1两方差分析下面的这个例子来自卫生统计学第四版，书还没有出来，大家先尝尝鲜。例 8.1 对小白鼠喂以A、B、C 三种不同的营养素，目的是了解不同营养素增重的效果。采用随机区组设计方法，以窝别作为划分区组的特征，以消除遗传对体重增长的影响。现将同品系同体重的 24 只小白鼠分为 8 个区组，每个区组 3 只小白鼠。三周后体重增量结果(克)列于下表，问小白鼠经三种不同营养素喂养后所增体重有无差别？区组号A 营养素B 营养素C 营养素1234567850.1047.8053.1063.5071.2041.4061.9042.2058.2048.5053.8064.2068.

5、4045.7053.0039.8064.5062.4058.6072.5079.3038.4051.2046.20根据统计分析的要求，建立了三个变量来包括上述信息，即 group表示区组，food 代表使用的营养素，weight 表示最终的重量，即：group food weight111250.0158.20依此类推。8.1.1univarate框界面说明这里只有一个结果变量 weight，要采用 univarate框，如下所示：在上面的这些框框中，最常用的有：DependentVariable 框、FixedFactors 框、M钮、t Hoc 钮，下面来一一解释。【Dependent V

6、ariable 框】选入需要分析的变量（应变量），只能选入一个。这里weight，将他选入即可。的应变量为【Fixed Factors 框】即固定，说的通俗一些，就是哎呀，我都不知道怎么解释好了，这样，如果你搞不明白，那么绝大多数要分析的都应该往里面选。这里要分析的是 group 和 food 两个变量，把他们全都给我抓进去！固定指的是在样本中它所有可能的取值都出现了，比如例中的 food，只可能有 1、2、3 这三个值，并且都出现了，就被称作固定效应；而相对应的随机效应的指的是所有可能的取值在样本中没有都出现，或不可能都出现，如本例中的 group，实际上总体中当然不可能只有这 8 窝，因此

7、要用样本中 group的情况来推论总体中 group 未出现的那些取值的情况时就会存在误差，因此被称为随机。我这里让 group 也选入固定框是基于下面的事实：这样做统计分析的结论是完全相同的。不同的只是推论的那部分。【Random Factors 框】用于选入随机，如果你弄不明白，假装没看见他就是了。【Covariate 框】用于选入协方差分析时的协变量，现在还用不到，不过下一个例子要给他送礼了。就【WLS Weight 框】即用于选入最小二乘法权重系数。别理他，根据理解，只有统计分析的狂才会想起来用他（雷同，纯属巧合）！【M钮】单击后出现一个框，用于设置在模型中包含哪些主效应和交互因子，默

8、认情况为 Full factorial，即分析所有的主效应和交互作用。这里没有交互作用可分析，所以要改一下，否则将作不出结果来。将按钮切换到右侧的custum，这时中部的 Build Term 下拉列表框就变黑可用，该框用于选择进入模型的这里的m交互作用级别，即是分析主效应、两阶交互、三阶交互、还是全部分析。只能分析主效应：选择 main，再用黑色箭头将 group 和 food 选入右侧框中，如果对这段叙述不太清楚，请参考下面的动画。该框中还有两个元素：左下方的 Sum of squares 框用于选择方差分析模型类别，有 1 型到 4 型四种，如果你搞不清他们之间的区别，使用默认的 3型即

9、可；中下部有个 Includeercept复选框，用于选择是否在模型中包括截距，不用改动，默认即可。【Contrast 钮】弹出 Contrast义，在这里，该框，用于对精细趋势检验和精确两两比较的选项进行定框比单方差分析的时候还要专业，使用频率也更少，反正我都没用过，就干脆就不介绍了。【Plots 钮】用于指定用模型的某些参数作图，比如用 food 和 group 来作图，用的也比较少（指国内，因为它主要是用来型用的）。【t Hoc 钮】该按钮弹出的两两比较框和第 7方差分析中的一模一样，不再重复。本题对 food 作两两比较，方法为 SNK 法。【Save 钮】将模型拟合时产生的中间结果或

10、参数保存为新变量供继续分析时用，可保存的东东有值、残差、用指标等。【Options 钮】当然是定义选项啦！可以定义输出哪些指标的估计均数、并做所选择的两两比较，还有其他一些输出，如常用描述指标、方差齐性检验等。好了，都解释完了，再重复以下，所作的操作为：1.2.3.4.5.6.7.8.9.yze=General Lineal m=UnivariateDependent Variable 框：选入 weight Fixed Factors 框：选入 group 和 foodM钮：单击Custom 单选钮：选中M框：选入 group 和 food单击 OKt Hoc 钮：单击t Hoc test

11、for 框：选入 food SNK 复选框：选中单击 OK10.11.12.单击 OK8.1.2结果解释按照上题的操作，结果输出如下：UnivariateysisofVariance这是一个所分析的取值情况列表，没不好懂的。现在大家看到的是一个典型的方差分析表，只不过是两的而已，我来解释一下：首先是所用方差分析模型的检验，F 值为 00.517，P 小于 0.05，因此所用的模型有统计学意义，可以用它来判断模型中系数有无统计学意义；第二行是截距，它在的分析中没有实际意义，忽略即可；第三行是变量 GROUP，可见它也有统计学意义，不过关心的也不是他；第四行是真正要分析的FOOD，非常遗憾，它的

12、P 值为 0.084，还没有统计学意义。尽管不太愿意，的结论也只能是：尚不能认为三种营养素喂养的小白鼠体重增量有差别。上表的标题内容翻译如下：t Hoc TestsFOODHomogeneoubsets变异来源III 型方差 SS度均方 MS统计量FP 值校正的模型2521.2949280.14411.517.000截距74359.534174359.5343056.985.000GROUP2376.3767339.48213.956.000FOOD144.917272.4592.979.084误差340.5431424.324合计77221.37024校正的合计2861.83623现在是两两

13、比较的结果，方法为 SNK 法，由于前面总的比较无差异，所以这里三种食物均在一个亚组内，检验无差异，P 值为 0.121。前面方差分析 FOOD 的 P 值不是 0.084 吗？这里又是 0.121，究竟哪个为准？两两比较只是近似的比较结果，应以前面方差分析的 P 为准，不过这两个 P 值不会在检验结果上发生质的，一般只是大小不同而已。好了，上面是正确的结果，如果 m会得出方差分析表如下：选择是采用 Full factor 又如何呢？看到了吗？由于所谓的交互作用将用于统计分析，什么结果也做不出来。度给全部“了，没有误差可8.2协方差分析例 8.2某医生欲了解成年重正常者与超重者的胆固醇是否不同

14、。而胆固醇含量与有关，资料见下表。正常组(X1) 胆固醇(Y1)超重组(X2) 胆固醇(Y2)483351434463494240474141563.54.65.85.84.98.73.65.54.95.14.14.65.1584171764933546539524558677.34.78.48.85.14.96.76.46.07.56.46.89.2该题选自医学统计学第二版第七章。考虑到统计分析对数据格式的要求，醇。这里建立三个变量：GROUP 表示组别，AGE 代表，CHOL 则表示胆固8.2.1分析步骤由于协方差分析涉及到许多较深的统计理论，这里我只好采用照本宣科的方法，告诉大家如何作，

15、而不作过多解释，欲进一步了解原理的朋友请参考医学统计学原书。首先应进行预分析，了解资料是否符合协方差分析的要求，最重要的一点就是看 age 的影响在两组中是否相同，这可以用 age 与 group 是否存在交互作用来表示。对该问题，粗糙的方法可以是作分组散点图，差不多就可以，也可以进行预分析，看交互作用有无统计学意义，这里用后讲解。法中最为精确的步骤来预分析步骤：1.yze=General Lineal m=Univariate2. Dependent Variable 框：选入 chol3.4.5.6.Fixed Factors 框：选入 groupM钮：单击Custom 单选钮：选中M可选

16、入)框：选入 group、age 和 group*age(后者用eraction 方法就7.8.9.Sum of squares 列表框：改为M单击 OK单击 OKI该步骤用于判断 group 和 age 间是否存在交互作用，如存在，则协方差分析的条件不满足，分析不能继续。注意这里选择了 MI，从而拟合结果和模型中变量的引入顺序有关，即侧重点在 group 对 chol 的影响大小和交互作用上。8.2.2结果解释预分析步骤的结果如下：Univariateysis of Variance上表显示交互作用无统计学意义，而且 P 值非常大，因此交换 group 和 age多半交互作用也无统计学意义，

17、因此可以不继续作预分析了，当然，严格的步骤应当交换两者的顺序继续进行预分析。正式分析步骤：1.2.3.4.5.6.7.8.9.yze=General Lineal m=UnivariateDependent Variable 框：选入 chol Fixed Factors 框：选入 groupM钮：单击Custom 单选钮：选中M框：选入 group、ageSum of squares 列表框：改为M单击 OK Options 钮：单击IIIDisplsy means for 框：选入 groupCompare mean effects 复选框：选中（下面的区间调整方法就用LSD(none)即

18、可）单击 OK 13.单击 OKUnivariateysis of Variance这是正式的统计分析结果，显示 group 和 age 都对胆固醇含量有影响，P 值分别为 0.038 和小于 0.001。Estimated Marginal Means这是两组的修正均数及相应的下方的提示表明该修正均数是按区间，显然超重组的胆固醇均值较高。为 50.2308 岁的情形计算的。8.3其他较简单的方差分析问题其他各种不太复杂的方差分析，如交叉设计的方差分析、析因设计的方差分析等的菜单选择和统计结果的解释我就不一一详细讲解了，大家举一反三，类似上面的方法就可以作出来。这里只是列举对于初学者来说可能有用的几个问题：需要分析的影响可以都选入