第八章spss之征服线性模型2(张文彤)

SPSS第八课：征服一般线性模型――GeneralLinearModel菜单详解（下）（医学统计之星：张文彤）上次更新日期：8.1两因素方差分析8.1.1univarate对话框界面说明8.1.2结果解释8.2协方差分析8.2.1分析步骤8.2.2结果解释8.3其他较简单的方差分析问题8.4多元方差分析8.4.1分析步骤8.4.2结果解释8.5重复测量的方差分析8.5.1Repeatedmeasures对话框界面说明8.5.2结果解释§8.4多元方差分析所谓的多元方差分析，就是说存在着不止一个应变量，而是两个以上的应变量共同反映了自变量的影响程度。比如要研究某些因素对儿童生长的影响程度，则身高、体重等都可以作为生长程度的测量因子，即都应作为应变量。8.4.1分析步骤为了方便起见，我们这里直接利用SPSS自带的数据集plastic.sav，假设tear_res、gloss和opacity都使反应橡胶质量的指标（不要笑，是假设），现在要研究extrusn和additive对橡胶的质量影响如何，则应采用多元方差分析。选择Analyze==GeneralLinearModel==Multivariate，则弹出Multivariate对话框，请注意，除了没有randomeffect外，它的所有元素都是和univariate对话框相同的，里面的内容也相同，因此我们这里就不再重复了。按照我们的分析要求，对话框操作步骤如下：1.Analyze==GeneralLinealmodel==Multivariate2.DependentVariable框：选入tear_res、gloss和opacity3.FixedFactors框：选入extrusn和additive4.单击OK此处两个自变量均是二分类变量，故无需选择两两比较方法。8.4.2结果解释按上面的选择，分析结果如下：GeneralLinearModel这是引入模型的自变量的取值情况列表。上表是针对模型中的自变量间及其交互作用所做的检验，采用的是四种多元检验方法。一般他们的结果都是相同的，如果不同，一般以Hotelling'sTrace方法的结果为准。可见在所用的模型中，extrusn和additive对结果变量是有统计学意义的，但交互作用无统计学意义。上表实际上是四个一元方差分析表的合并，即分别考虑四个应变量时的方差分析结果。上面的多元方差分析已经得知两自变量对应变量有影响，从现在的分析表就可以更清楚的知道是对那些自变量影响较大。对照可知，extrusn和additive对tearresistance和gloss都有较大影响，而他们的交互作用对gloss有影响，他们（及交互作用）对Opacity都没有影响。§8.5重复测量的方差分析重复测量的方差分析指的是一个应变量被重复测量好几次，从而同一个个体的几次观察结果间存在相关，这样就不满足普通分析的要求，需要用重复测量的方差分析模型来解决。8.5.1Repeatedmeasures对话框界面说明实际上，如果对普通方差分析模型作出正确的设置，两者的分析结果是完全相同的，即都正确，那么，重复测量的方差分析过程有何优势呢？我们通过下面的例子来看看：例8.3在数据集anxity2.sav中判断：anxiety和tension对实验结果（即trial1~trial4）有无影响；四次试验间有无差异；试验次数和两个变量有无交互作用。anxity2.sav和anxity.sav实际上是同一个数据，但根据不同的分析目的采用了不同的数据排列方式。如果采用anxity.sav进行分析，我们可以分析四次试验间有无差异的问题，但对另两个问题就无能为力了，因为用普通的方差分析模型，anxity和tension的影响被合并到了subject中，根本就无法分解出来进行分析，这时，我们就只能求助于重复测量的方差分析模型。在菜单中选择Analyze==GeneralLinealmodel==Repeatedmeasures，系统首先会弹出一个重复测量因子定义对话框如下：因为是重复测量的模型，应变量被重复测量了几次，分别存放在几个变量中，所以我们这里要自行定义应变量。默认的名称为factor1，我们将其改为trail，下面的因素等级数填入4(因一共测量了四次)。单击Add钮，则该变量被加入，我们就完成了模型设置的第一步：应变量名称和测量次数定义。单击define，我们开始进行下一个步骤：具体重复测量变量定义及模型设置，对话框如下：这个对话框和我们以前看到的方差分析对话框不太一样：它没有应变量框，而是改为了组内效应框，实际上是一回事，上面我们定义了trial有四次测量，此处就给出了四个空让你填入相应代表四次测量的变量，选中trial1~trial4，将其选入；然后要选择自变量了（这里又将其称为了betweensubjectsfactor），将剩下的三个都选入即可。最后，根据题意，不需要检验anxity与tension的交互作用对试验次数有无交互作用，所以要在model中作相应设置，把那个东东拉出来。详细的操作步骤如下：1.Analyze==GeneralLinealmodel==Repeatedmeasures2.Within-subjectfactorname框：键选入trial3.numberoflevels框：键入44.单击ADD钮5.单击DEFINE钮6.Within-subjectvariables(trial)框：选入trial1~trial47.betweensubjectsfactor框：选入subject、anxity和tension8.单击MODEL钮9.Custom单选钮：选中10.Within-subjectModel框：选入trial11.betweensubjectsModel框：选入anxity和tension12.单击CONTINUE13.单击OK请注意，这里没有选入变量subject，因为它实际上在这里成为了一个记录ID，要是将它选入，则什么都检验不了了。8.5.2结果解释本题的分析结果如下：GeneralLinearModel上表给出了所定义的4次测量的变量名，在模型中它们都代表一个应变量trial，只是测量的次数不同而已。这是引入模型的其它自变量的情况列表。上表是针对所检验的结果变量trial，以及他和另两个引入模型的自变量间的交互作用是否存在统计学意义，采用的是四种多元检验方法。一般他们的结果都是相同的，如果不同，我一般以Hotelling'sTrace方法的结果为准。可见在所用的模型中，trial的四次测量间的确是存在着统计学差异的，但它和另两个变量间的交互作用无统计学意义。上表是球形检验，因为重复测量的方差分析模型要求所检验的应变量服从一种叫做球形分布的东东。上面可能有些内容不好懂，不过没关系，只要看到近似卡方为9.383，自由度为5，P值为0.097就可以了。因此trial是勉强服从球形分布的，可以进行重复测量的方差分析。来源医学统计之星2004-7-1222:06:00SPSS10.0高级教程十：征服一般线性模型（2）张文彤分享|收藏热点推荐生物谷盘点20092009年，国内生物医药的突破之年。不仅有干细胞发现的新突破，还有转基因作物政策的新举措。更多专题相关阅读推荐学术会议学术研讨会、峰会论坛、展会更多会议§8.4多元方差分析所谓的多元方差分析，就是说存在着不止一个应变量，而是两个以上的应变量共同反映了自变量的影响程度。比如要研究某些因素对儿童生长的影响程度，则身高、体重等都可以作为生长程度的测量因子，即都应作为应变量。8.4.1分析步骤为了方便起见，我们这里直接利用SPSS自带的数据集plastic.sav，假设tear_res、gloss和opacity都使反应橡胶质量的指标（不要笑，是假设），现在要研究extrusn和additive对橡胶的质量影响如何，则应采用多元方差分析。选择Analyze==GeneralLinearModel==Multivariate，则弹出Multivariate对话框，请注意，除了没有randomeffect外，它的所有元素都是和univariate对话框相同的，里面的内容也相同，因此我们这里就不再重复了。按照我们的分析要求，对话框操作步骤如下：1.Analyze==GeneralLinealmodel==Multivariate2.DependentVariable框：选入tear_res、gloss和opacity3.FixedFactors框：选入extrusn和additive4.单击OK此处两个自变量均是二分类变量，故无需选择两两比较方法。8.4.2结果解释按上面的选择，分析结果如下：GeneralLinearModel这是引入模型的自变量的取值情况列表。上表是针对模型中的自变量间及其交互作用所做的检验，采用的是四种多元检验方法。一般他们的结果都是相同的，如果不同，一般以Hotelling'sTrace方法的结果为准。可见在所用的模型中，extrusn和additive对结果变量是有统计学意义的，但交互作用无统计学意义。上表实际上是四个一元方差分析表的合并，即分别考虑四个应变量时的方差分析结果。上面的多元方差分析已经得知两自变量对应变量有影响，从现在的分析表就可以更清楚的知道是对那些自变量影响较大。对照可知，extrusn和additive对tearresistance和gloss都有较大影响，而他们的交互作用对gloss有影响，他们（及交互作用）对Opacity都没有影响。§8.5重复测量的方差分析重复测量的方差分析指的是一个应变量被重复测量好几次，从而同一个个体的几次观察结果间存在相关，这样就不满足普通分析的要求，需要用重复测量的方差分析模型来解决。8.5.1Repeatedmeasures对话框界面说明实际上，如果对普通方差分析模型作出正确的设置，两者的分析结果是完全相同的，即都正确，那么，重复测量的方差分析过程有何优势呢？我们通过下面的例子来看看：例8.3在数据集anxity2.sav中判断：anxiety和tension对实验结果（即trial1~trial4）有无影响；四次试验间有无差异；试验次数和两个变量有无交互作用。anxity2.sav和anxity.sav实际上是同一个数据，但根据不同的分析目的采用了不同的数据排列方式。如果采用anxity.sav进行分析，我们可以分析四次试验间有无差异的问题，但对另两个问题就无能为力了，因为用普通的方差分析模型，anxity和tension的影响被合并到了subject中，根本就无法分解出来进行分析，这时，我们就只能求助于重复测量的方差分析模型。在菜单中选择Analyze==GeneralLinealmodel==Repeatedmeasures，系统首先会弹出一个重复测量因子定义对话框如下：因为是重复测量的模型，应变量被重复测量了几次，分别存放在几个变量中，所以我们这里要自行定义应变量。默认的名称为factor1，我们将其改为trail，下面的因素等级数填入4(因一共测量了四次)。单击Add钮，则该变量被加入，我们就完成了模型设置的第一步：应变量名称和测量次数定义。单击define，我们开始进行下一个步骤：具体重复测量变量定义及模型设置，对话框如下：这个对话框和我们以前看到的方差分析对话框不太一样：它没有应变量框，而是改为了组内效应框，实际上是一回事，上面我们定义了trial有四次测量，此处就给出了四个空让你填入相应代表四次测量的变量，选中trial1~trial4，将其选入；然后要选择自变量了（这里又将其称为了betweensubjectsfactor），将剩下的三个都选入即可。最后，根据题意，不需要检验anxity与tension的交互作用对试验次数有无交互作用，所以要在model中作相应设置，把那个东东拉出来。详细的操作步骤如下：1.Analyze==GeneralL