《医学统计》第八章 征服一般线性模型――General Linear Model菜单详解（1/2）（张文彤）

SPSS第八课:征服一般线性模型 - General linear model菜单详解(上) (医学统计之星:张文彤) 请注意,本章的标题用了一些修辞手法,一般线性模型可不是用一章就可 以说清楚的,因为它包括的内容实在太多了。 那么,究竟我们用到的哪些分析会包含在其中呢?简而言之:凡是和方差 分析粘边的都可以用他来做。比如成组设计的方差分析(即单因素方差分析) 配伍设计的方差分析(即两因素方差分析)、交叉设计的方差分析、析因设计的 方差分析、重复测量的方差分析、协方差分析等等。因此,能真正掌握GLM菜单 的用法,会使大家的统计分析能力有极大地提高 实际上一般线性模型包括的统计模型还不止这些,我这里举出来的只是从用 SPSS作统计分析的角度而言的一些。 好了,既然一般线性模型的能力如此强大,那么下属的四个子菜单各自的 功能是什么呢?请看: Univariate子菜单:四个菜单中的大哥大,绝大部分的方法分析都在这 里面进行 Multivariate子菜单:当结果变量(应变量)不止一个时,当然要用他 来分析啦! Repeted Measures子菜单:顾名思义,重复测量的数据就要用他来分析, 这一点我可能要强调一下,用前两个菜单似乎都可以分析出来结果,但在 许多情况下该结果是不正确的,应该用重复测量的分析方法才对(不能再 讲了,再讲下去就会扯到多水平模型去了) ariance Components子菜单:用于作方差成份模型的,这个模型实在 太深,不是一时半会说的请的,所以我在这里就干脆不讲了 出于模型复杂性、篇幅、应用范围及乱七八糟一系列的理由,当然主要是 我懒得一一解释,我决定本章采用举例讲解的方式,及讲解一些常见的分析实例, 通过这种方法来熟悉那些最为常用的分析方法 D多对统计分析的数据格式不太熟悉的朋友,请一定先去看看统计软件第一课 论统计软件中的数据录入格式,会大有帮助的。 §8.1两因素方差分析 下面的这个例子来自《卫生统计学》第四版,书还没有出来,大家先尝尝 鲜

SPSS 第八课：征服一般线性模型 ――General Linear Model 菜单详解（上） （医学统计之星：张文彤） 请注意，本章的标题用了一些修辞手法，一般线性模型可不是用一章就可 以说清楚的，因为它包括的内容实在太多了。 那么，究竟我们用到的哪些分析会包含在其中呢？简而言之：凡是和方差 分析粘边的都可以用他来做。比如成组设计的方差分析（即单因素方差分析）、 配伍设计的方差分析（即两因素方差分析）、交叉设计的方差分析、析因设计的 方差分析、重复测量的方差分析、协方差分析等等。因此，能真正掌握 GLM 菜单 的用法，会使大家的统计分析能力有极大地提高。 实际上一般线性模型包括的统计模型还不止这些，我这里举出来的只是从用 SPSS 作统计分析的角度而言的一些。 好了，既然一般线性模型的能力如此强大，那么下属的四个子菜单各自的 功能是什么呢？请看： • Univariate 子菜单：四个菜单中的大哥大，绝大部分的方法分析都在这 里面进行。 • Multivariate 子菜单：当结果变量（应变量）不止一个时，当然要用他 来分析啦！ • Repeted Measures 子菜单：顾名思义，重复测量的数据就要用他来分析， 这一点我可能要强调一下，用前两个菜单似乎都可以分析出来结果，但在 许多情况下该结果是不正确的，应该用重复测量的分析方法才对（不能再 讲了，再讲下去就会扯到多水平模型去了）。 • Variance Components 子菜单：用于作方差成份模型的，这个模型实在 太深，不是一时半会说的请的，所以我在这里就干脆不讲了。 出于模型复杂性、篇幅、应用范围及乱七八糟一系列的理由，当然主要是 我懒得一一解释，我决定本章采用举例讲解的方式，及讲解一些常见的分析实例， 通过这种方法来熟悉那些最为常用的分析方法。 对统计分析的数据格式不太熟悉的朋友，请一定先去看看统计软件第一课： 论统计软件中的数据录入格式，会大有帮助的。 §8.1 两因素方差分析 下面的这个例子来自《卫生统计学》第四版，书还没有出来，大家先尝尝 鲜

例8.1对小白鼠喂以A、B、C三种不同的营养素,目的是了解不同营养素 增重的效果。采用随机区组设计方法,以窝别作为划分区组的特征,以消除遗传 因素对体重增长的影响。现将同品系同体重的24只小白鼠分为8个区组,每个 区组3只小白鼠。三周后体重增量结果(克)列于下表,问小白鼠经三种不同营养 素喂养后所增体重有无差别? 区组号A营养素B营养素C营养素 150.1058.2064.50 8048.5062.40 353.1053.8058.60 463.5064.2072.50 571.2068.4079.30 641 45.7038.40 761.9053.0051.20 842.2039.8046.20 根据统计分析的要求,我们建立了三个变量来包括上述信息,即 group 表示区组,food代表使用的营养素, weight表示最终的重量,即 group food weight 1150.01 依此类推 8.1.1 uni karate对话框界面说明 这里只有一个结果变量 weight,要采用 univariate对话框,如下所示

例 8.1 对小白鼠喂以 A、B、C 三种不同的营养素，目的是了解不同营养素 增重的效果。采用随机区组设计方法，以窝别作为划分区组的特征，以消除遗传 因素对体重增长的影响。现将同品系同体重的 24 只小白鼠分为 8 个区组，每个 区组 3 只小白鼠。三周后体重增量结果(克)列于下表，问小白鼠经三种不同营养 素喂养后所增体重有无差别？ 区组号 A 营养素 B 营养素 C 营养素 1 50.10 58.20 64.50 2 47.80 48.50 62.40 3 53.10 53.80 58.60 4 63.50 64.20 72.50 5 71.20 68.40 79.30 6 41.40 45.70 38.40 7 61.90 53.00 51.20 8 42.20 39.80 46.20 根据统计分析的要求，我们建立了三个变量来包括上述信息，即 group 表示区组，food 代表使用的营养素，weight 表示最终的重量，即： group food weight 1 1 50.01 1 2 58.20 依此类推。 8.1.1 univarate 对话框界面说明 这里只有一个结果变量 weight，要采用 univarate 对话框，如下所示：

4 Univariate 参 group Dependent variable: Model Food ght Contrasts Fi× ed Factor(s Post hoc Random Factor[s]: Save Options Covariate[s]: YLS Weig OK Paste Reset Cancel Help 在上面的这些框框钮钮中,最常用的有: Dependent variable框、 Fixed Factors框、 Model钮、 Post hoc钮,下面我们来一一解释 【 Dependent variable框】 选入需要分析的变量(应变量),只能选入一个。这里我们的应变量为 weight,将他选入即可。 【 Fixed Factors框】 即固定因素,说的通俗一些,就是一一哎呀,我都不知道怎么解释好了, 这样,如果你搞不明白,那么绝大多数要分析的因素都应该往里面选。这里我们 要分析的是 group和food两个变量,把他们全都给我抓进去! 固定因素指的是在样本中它所有可能的取值都出现了,比如例中的food 只可能有1、2、3这三个值,并且都出现了,就被称作固定效应:而相对应的随 机效应的因素指的是所有可能的取值在样本中没有都出现,或不可能都出现,如 本例中的 group,实际上总体中当然不可能只有这8窝,因此要用样本中 group 的情况来推论总体中 group未出现的那些取值的情况时就会存在误差,因此被称 为随机因素。我这里让 group也选入固定框是基于下面的事实:这样做统计分析 的结论是完全相同的。不同的只是推论的那部分

在上面的这些框框钮钮中，最常用的有：Dependent Variable 框、Fixed Factors 框、Model 钮、Post Hoc 钮，下面我们来一一解释。 【Dependent Variable 框】 选入需要分析的变量（应变量），只能选入一个。这里我们的应变量为 weight，将他选入即可。 【Fixed Factors 框】 即固定因素，说的通俗一些，就是－－哎呀，我都不知道怎么解释好了， 这样，如果你搞不明白，那么绝大多数要分析的因素都应该往里面选。这里我们 要分析的是 group 和 food 两个变量，把他们全都给我抓进去！ 固定因素指的是在样本中它所有可能的取值都出现了，比如例中的 food， 只可能有 1、2、3 这三个值，并且都出现了，就被称作固定效应；而相对应的随 机效应的因素指的是所有可能的取值在样本中没有都出现，或不可能都出现，如 本例中的 group，实际上总体中当然不可能只有这 8 窝，因此要用样本中 group 的情况来推论总体中 group 未出现的那些取值的情况时就会存在误差，因此被称 为随机因素。我这里让 group 也选入固定框是基于下面的事实：这样做统计分析 的结论是完全相同的。不同的只是推论的那部分

【 Random factors框】 用于选入随机因素,如果你弄不明白,假装没看见他就是了。 【 Covariate框】 用于选入协方差分析时的协变量,现在还用不到,不过下一个例子我们就 要给他送礼了。 【 WLS Weight框】 即用于选入最小二乘法权重系数。别理他,根据我的理解,只有统计分析 的变态狂才会想起来用他(如有雷同,纯属巧合)! 【 Model钮】 单击后出现一个对话框,用于设置在模型中包含哪些主效应和交互因子 默认情况为 Full factorial,即分析所有的主效应和交互作用。我们这里没有 交互作用可分析,所以要改一下,否则将作不出结果来。将按钮切换到右侧的 custum,这时中部的 Build term下拉列表框就变黑可用,该框用于选择进入模 型的因素交互作用级别,即是分析主效应、两阶交互、三阶交互、还是全部分析 这里我们只能分析主效应:选择main,再用黑色箭头将 group和food选入右侧 的 model框中,如果对这段叙述不太清楚,请参考下面的动画。 Specify Model F Full factorial C Custom Continue Factors covariates Model Cancel group) food[F] Help Build Term(sI Interaction Sum of squares vre-]lmeudetetmde

【Random Factors 框】 用于选入随机因素，如果你弄不明白，假装没看见他就是了。 【Covariate 框】 用于选入协方差分析时的协变量，现在还用不到，不过下一个例子我们就 要给他送礼了。 【WLS Weight 框】 即用于选入最小二乘法权重系数。别理他，根据我的理解，只有统计分析 的变态狂才会想起来用他（如有雷同，纯属巧合）！ 【Model 钮】 单击后出现一个对话框，用于设置在模型中包含哪些主效应和交互因子， 默认情况为 Full factorial，即分析所有的主效应和交互作用。我们这里没有 交互作用可分析，所以要改一下，否则将作不出结果来。将按钮切换到右侧的 custum，这时中部的 Build Term 下拉列表框就变黑可用，该框用于选择进入模 型的因素交互作用级别，即是分析主效应、两阶交互、三阶交互、还是全部分析。 这里我们只能分析主效应：选择 main，再用黑色箭头将 group 和 food 选入右侧 的 model 框中，如果对这段叙述不太清楚，请参考下面的动画

该对话框中还有两个元素:左下方的 Sum of squares框用于选择方差分析 模型类别,有1型到4型四种,如果你搞不清他们之间的区别,使用默认的3 型即可:中下部有个 Include intercept in model复选框,用于选择是否在模 型中包括截距,不用改动,默认即可。 【 Contrast钮】 弹出 Contrast对话框,用于对精细趋势检验和精确两两比较的选项进行定 义,在这里,该对话框比单因素方差分析的时候还要专业,使用频率也更少,反 正我都没用过,就干脆就不介绍了。 【 Plots钮】 用于指定用模型的某些参数作图,比如用food和 group来作图,用的也比 较少(指国内,因为它主要是用来做模型诊断用的)。 【 Post hoc钮】 该按钮弹出的两两比较对话框和第7章单因素方差分析中的一模一样,不 再重复。本题对food作两两比较,方法为SNK法 【Save钮】 将模型拟合时产生的中间结果或参数保存为新变量供继续分析时用,可保 存的东东有预测值、残差、诊断用指标等。 【 Options钮】 当然是定义选项啦!可以定义输出哪些指标的估计均数、并做所选择的两 两比较,还有其他一些输出,如常用描述指标、方差齐性检验等。 好了,都解释完了,再重复以下,我们所作的操作为: 1. Analyze==>General Lineal mode l==>Univariate 2. Dependent variable框:选入 weight 3. Fixed Factors框:选入 group和food 4. Model钮:单击 5. Custom单选钮:选中 6. Model框:选入 group和food 7.单击OK 8. Post hoc钮:单击 9. Post hoc test for框:选入food 10.SNK复选框:选中 11.单击OK 12.单击OK

该对话框中还有两个元素：左下方的 Sum of squares 框用于选择方差分析 模型类别，有 1 型到 4 型四种，如果你搞不清他们之间的区别，使用默认的 3 型即可；中下部有个 Include intercept in model 复选框，用于选择是否在模 型中包括截距，不用改动，默认即可。 【Contrast 钮】 弹出 Contrast 对话框，用于对精细趋势检验和精确两两比较的选项进行定 义，在这里，该对话框比单因素方差分析的时候还要专业，使用频率也更少，反 正我都没用过，就干脆就不介绍了。 【Plots 钮】 用于指定用模型的某些参数作图，比如用 food 和 group 来作图，用的也比 较少（指国内，因为它主要是用来做模型诊断用的）。 【Post Hoc 钮】 该按钮弹出的两两比较对话框和第 7 章单因素方差分析中的一模一样，不 再重复。本题对 food 作两两比较，方法为 SNK 法。 【Save 钮】 将模型拟合时产生的中间结果或参数保存为新变量供继续分析时用，可保 存的东东有预测值、残差、诊断用指标等。 【Options 钮】 当然是定义选项啦！可以定义输出哪些指标的估计均数、并做所选择的两 两比较，还有其他一些输出，如常用描述指标、方差齐性检验等。 好了，都解释完了，再重复以下，我们所作的操作为： 1. Analyze==>General Lineal model==>Univariate 2. Dependent Variable 框：选入 weight 3. Fixed Factors 框：选入 group 和 food 4. Model 钮：单击 5. Custom 单选钮：选中 6. Model 框：选入 group 和 food 7. 单击 OK 8. Post Hoc 钮：单击 9. Post Hoc test for 框：选入 food 10. SNK 复选框：选中 11. 单击 OK 12.单击 OK

8.1.2结果解释 按照上题的操作,结果输出如下 Univariate Analysis of Variance Between-Subjects Factors GROUP 1.0 2.00 4.00 3 5.00 6.00 3 00D1.00 这是一个所分析因素的取值情况列表,没有什么不好懂的。 Tests of Between-Subjects Effects Dependent Variable: WEIGHT Type I Sum Source of Squares df Mean Square Corrected Mode|2521.2944928014411517000 Intercept 74359534 1743595343056985000 GROUP 2376.376 33948213956000 144.917 2.979 340.54314 77221.37024 2861.836 a R Squared=. 881(Adjusted R Squared=. 805) 现在大家看到的是一个典型的方差分析表,只不过是两因素的而已,我来 解释一下:首先是所用方差分析模型的检验,F值为00.517,P小于0.05,因此 所用的模型有统计学意义,可以用它来判断模型中系数有无统计学意义;第二行 是截距,它在我们的分析中没有实际意义,忽略即可;第三行是变量 GROUP,可 见它也有统计学意义,不过我们关心的也不是他;第四行是我们真正要分析的

8.1.2 结果解释 按照上题的操作，结果输出如下： Univariate Analysis of Variance 这是一个所分析因素的取值情况列表，没有什么不好懂的。 现在大家看到的是一个典型的方差分析表，只不过是两因素的而已，我来 解释一下：首先是所用方差分析模型的检验，F 值为 00.517，P 小于 0.05，因此 所用的模型有统计学意义，可以用它来判断模型中系数有无统计学意义；第二行 是截距，它在我们的分析中没有实际意义，忽略即可；第三行是变量 GROUP，可 见它也有统计学意义，不过我们关心的也不是他；第四行是我们真正要分析的

FOO,非常遗憾,它的P值为0.084,还没有统计学意义。尽管不太愿意,我们 的结论也只能是:尚不能认为三种营养素喂养的小白鼠体重增量有差别。 上表的标题内容翻译如下: 变异来源型方差SS自由度均方MS统计量FP 校正的模型2521.2949280.14411.517.000 74359.534 174359.5343056.985.000 GROUP 2376.3767332195900 FOOD 144.917 272.4592.979.084 误差 1424.324 合计 77221.37024 校正的合计2861.836 Post hoc test FOOD Homogeneous Subsets WEIGHT Subset FOOD N 853.9000 853.950 Means for groups in homogeneous subsets are displayed Based on Type Ill Sum of Squares The error term is Mean Square(Error=24. 324 a Uses Harmonic Mean Sample size =8.000 现在是两两比较的结果,方法为SN法,由于前面总的比较无差异,所以 这里三种食物均在一个亚组内,检验无差异,P值为0.121

FOOD，非常遗憾，它的 P 值为 0.084，还没有统计学意义。尽管不太愿意，我们 的结论也只能是：尚不能认为三种营养素喂养的小白鼠体重增量有差别。 上表的标题内容翻译如下： 变异来源 III 型方差 SS 自由度 均方 MS 统计量 F P 值 校正的模型 2521.294 9 280.144 11.517 .000 截距 74359.534 1 74359.534 3056.985 .000 GROUP 2376.376 7 339.482 13.956 .000 FOOD 144.917 2 72.459 2.979 .084 误差 340.543 14 24.324 合计 77221.370 24 校正的合计 2861.836 23 Post Hoc Tests FOOD Homogeneous Subsets 现在是两两比较的结果，方法为 SNK 法，由于前面总的比较无差异，所以 这里三种食物均在一个亚组内，检验无差异，P 值为 0.121

前面方差分析FOOD的P值不是0.084吗?这里又是0.121,究竟哪个为准? 两两比较只是近似的比较结果,应以前面方差分析的P为准,不过这两个P值不 会在检验结果上发生质的冲突,一般只是大小不同而己。 好了,上面是正确的结果,如果 model选择是采用 Full factor又如何呢? 会得出方差分析表如下: Tests of Between-Subjects Effects pendent Variable: WEIGHT Type Ill Sum df Mean Square F Sig Corrected Model 2861836323 124428 Intercept 74359534 74359534 GROUP 2376.376 339482 FOOD 144917 GROUP* FOOD 340.542 24.324 Error 000 Total 77221.370 a R Squared=1.000 (Adjusted R Squared =) 看到了吗?由于所谓的交互作用将自由度给全部“吃”掉了,没有误差可 用于统计分析,什么结果也做不出来 §8.2协方差分析 例8.2某医生欲了解成年人体重正常者与超重者的血清胆固醇是否不同 而胆固醇含量与年龄有关,资料见下表。 正常组 超重组 年龄(X1)胆固醇(Y1)年龄(X2)胆固醇(Y2) 3.5 4.6 51 5.8 8.4 5.8 4.9 8.7 33 3.6 .7 5.5

前面方差分析 FOOD 的 P 值不是 0.084 吗？这里又是 0.121，究竟哪个为准？ 两两比较只是近似的比较结果，应以前面方差分析的 P 为准，不过这两个 P 值不 会在检验结果上发生质的冲突，一般只是大小不同而已。 好了，上面是正确的结果，如果 model 选择是采用 Full factor 又如何呢？ 会得出方差分析表如下： 看到了吗？由于所谓的交互作用将自由度给全部“吃”掉了，没有误差可 用于统计分析，什么结果也做不出来。 §8.2 协方差分析 例 8.2 某医生欲了解成年人体重正常者与超重者的血清胆固醇是否不同。 而胆固醇含量与年龄有关，资料见下表。 正常组 超重组 年龄(X1) 胆固醇(Y1) 年龄(X2) 胆固醇(Y2) 48 3.5 58 7.3 33 4.6 41 4.7 51 5.8 71 8.4 43 5.8 76 8.8 44 4.9 49 5.1 63 8.7 33 4.9 49 3.6 54 6.7 42 5.5 65 6.4 40 4.9 39 6.0

5.1 52 4.1 4.6 6.8 5.1 该题选自《医学统计学》第二版第七章。考虑到统计分析对数据格式的要 求,我们这里建立三个变量: GROUP表示组别,AGE代表年龄,CHOL则表示胆固 8.2.1分析步骤 由于协方差分析涉及到许多较深的统计理论,这里我只好采用照本宣科的 方法,告诉大家如何作,而不作过多解释,欲进一步了解原理的朋友请参考《医 学统计学》原书。 首先应进行预分析,了解资料是否符合协方差分析的要求,最重要的一点 就是看age的影响在两组中是否相同,这可以用age与 group是否存在交互作用 来表示。对该问题,粗糙的方法可以是作分组散点图,差不多就可以,也可以进 行预分析,看交互作用有无统计学意义,这里用后一种方法中最为精确的步骤来 讲解 预分析步骤: 1. Analyze==>General Lineal mode l==>Univariate 2. Dependent variable框:选入chol 3. Fixed Factors框:选入 group 4. Model钮:单击 5. Custom单选钮:选中 6. Model框:选入 group、age和 group*age(后者用 interaction方法就 可选入) 7. Sum of squares列表框:改为 Model i 8.单击OK 9.单击OK 该步骤用于判断 group和age间是否存在交互作用,如存在,则协方差分 析的条件不满足,分析不能继续。注意这里选择了 Model i,从而拟合结果和模 型中变量的引入顺序有关,即侧重点在 group对chol的影响大小和交互作用上。 8.2.2结果解释 预分析步骤的结果如下

47 5.1 52 7.5 41 4.1 45 6.4 41 4.6 58 6.8 56 5.1 67 9.2 该题选自《医学统计学》第二版第七章。考虑到统计分析对数据格式的要 求，我们这里建立三个变量：GROUP 表示组别，AGE 代表年龄，CHOL 则表示胆固 醇。 8.2.1 分析步骤 由于协方差分析涉及到许多较深的统计理论，这里我只好采用照本宣科的 方法，告诉大家如何作，而不作过多解释，欲进一步了解原理的朋友请参考《医 学统计学》原书。 首先应进行预分析，了解资料是否符合协方差分析的要求，最重要的一点 就是看 age 的影响在两组中是否相同，这可以用 age 与 group 是否存在交互作用 来表示。对该问题，粗糙的方法可以是作分组散点图，差不多就可以，也可以进 行预分析，看交互作用有无统计学意义，这里用后一种方法中最为精确的步骤来 讲解。 预分析步骤： 1. Analyze==>General Lineal model==>Univariate 2. Dependent Variable 框：选入 chol 3. Fixed Factors 框：选入 group 4. Model 钮：单击 5. Custom 单选钮：选中 6. Model 框：选入 group、age 和 group*age(后者用 interaction 方法就 可选入) 7. Sum of squares 列表框：改为 Model I 8. 单击 OK 9. 单击 OK 该步骤用于判断 group 和 age 间是否存在交互作用，如存在，则协方差分 析的条件不满足，分析不能继续。注意这里选择了 Model I，从而拟合结果和模 型中变量的引入顺序有关，即侧重点在 group 对 chol 的影响大小和交互作用上。 8.2.2 结果解释 预分析步骤的结果如下：

Univariate Analysis of variance Between-Subjects Factors GROUP 100 13 Tests of Between-Subjects Effects Dependent variable: CHOL Type I Sum Source of Squares df Mean Square Sig Corrected Model 43.002 14.33414988000 Intercept 916.898 24.38 31111 916.898958.733000 18.61519465000 24.38025492000 GROUP*AGE 6.499E-03 6499E03 007935 Error 21.040 956 Total 980.940 a R Squared=671(Adjusted R Squared= 627) 上表显示交互作用无统计学意义,而且P值非常大,因此交换 group和age 多半交互作用也无统计学意义,因此可以不继续作预分析了,当然,严格的步骤 应当交换两者的顺序继续进行预分析。 正式分析步骤 1. Analyze==>General Lineal mode l==>Univariate 2. Dependent variable框:选入chol 3. Fixed factors框:选入 group 4. Model钮:单击 5. Custom单选钮:选中 Model框:选入 group、age 7. Sum of squares列表框:改为 Model lll 单击OK 9.0 ptions钮:单击 10. Display means for框:选入 group 11. Compare mean effects复选框:选中(下面的区间调整方法就用 LSD(none)即可) 2.单击OK 13.单击OK

Univariate Analysis of Variance 上表显示交互作用无统计学意义，而且 P 值非常大，因此交换 group 和 age 多半交互作用也无统计学意义，因此可以不继续作预分析了，当然，严格的步骤 应当交换两者的顺序继续进行预分析。 正式分析步骤： 1. Analyze==>General Lineal model==>Univariate 2. Dependent Variable 框：选入 chol 3. Fixed Factors 框：选入 group 4. Model 钮：单击 5. Custom 单选钮：选中 6. Model 框：选入 group、age 7. Sum of squares 列表框：改为 Model III 8. 单击 OK 9. Options 钮：单击 10. Displsy means for 框：选入 group 11. Compare mean effects 复选框：选中（下面的区间调整方法就用 LSD(none)即可） 12. 单击 OK 13.单击 OK