《spss中一般线性模型解析》由会员分享,可在线阅读,更多相关《spss中一般线性模型解析(58页珍藏版)》请在金锄头文库上搜索。
1、spssspss中一般线性模型解析中一般线性模型解析一般一般线性模型性模型n一般一般线性模型性模型单变量分析的基本量分析的基本过程程n完全随机完全随机设计资料的方差分析料的方差分析n随机区随机区组(单位位组)设计资料的方差分析料的方差分析2024/8/112一、一般一、一般线性模型性模型单变量分析的基本量分析的基本过程程nGeneral Linear Model(GLM,一般,一般线性模型)性模型)n包括:包括:lUnivariate(单因因变量多因素方差分析),量多因素方差分析),lMultivariate(多因(多因变量方差分析),量方差分析),lRepeated Measures(重复(
2、重复测量方差分析),量方差分析),lVariance(方差分量分析)(方差分量分析)l GLM可完成多自可完成多自变量、多水平、多因量、多水平、多因变量、重复量、重复测量方差分析以及量方差分析以及协方差分析等。方差分析等。2024/8/113Univariate(单因因变量方差分析)基本量方差分析)基本过程程2024/8/1141 主主对话框框nDependent Variable:因因变量量nFixed Facter: 固定因子,所有可能的水固定因子,所有可能的水平都出平都出现在在样本中,如分本中,如分组等等nRandom Facter: 随机因子,所有可能随机因子,所有可能的取的取值并不都
3、在并不都在样本中出本中出现,如,如观察个察个体体nCovariates:协变量,量,协方差分析方差分析时用用nWLS Weight: WLS权重。用于加重。用于加权最最小二乘分析。小二乘分析。2024/8/1152 功能按功能按钮nModel:分析模型:分析模型nContrast:对照方法照方法nPlots:分布:分布图形形nPost Hoc:多重比:多重比较nSave:保存运算:保存运算值nOption:选择输出出项。2024/8/1162.1 Model按按钮、在、在Specify Model栏中指定模型中指定模型类型型nFull Factorial,全模型,系,全模型,系统默默认。包括所
4、有因素的主效。包括所有因素的主效应和所有的交互效和所有的交互效应。例如有三。例如有三个因素个因素变量,全模型包括三个量,全模型包括三个因素的主效因素的主效应、两两的交互效、两两的交互效应和三个因素的高和三个因素的高级交互效交互效应。nCustom,自定,自定义模型。模型。选择此此项激活下面各操作框激活下面各操作框 2024/8/117、建立自定、建立自定义模型模型nFactorsCovariates 框中框中自自动列出可以作列出可以作为因素的因素的变量名,其后面的括号中量名,其后面的括号中标有字母有字母“F”(固定因(固定因子)、子)、“R”(随机因子)(随机因子)或者或者“C”(协变量)。量
5、)。2024/8/118A、选择效效应类型型nInteractin:交互效:交互效应nMain effects:主效:主效应nAll2-way: 所有所有2维交互效交互效应nAll3-way:所有:所有3维交互效交互效应nAll4-Way:所有:所有4维交互效交互效应nAll5-Way:所有:所有5维交互效交互效应2024/8/119B、选择模型中的主效模型中的主效应 (Model)n首先定首先定义效效应类型型为Main effectsn鼠鼠标键单击某一个因素,某一个因素,该变量名量名背景将改背景将改变颜色(一般色(一般变为蓝色),色),单击Build Term(s)栏中下面的中下面的箭箭头,
6、该变量出量出现在在Mode1中。一中。一个个变量名占一行称量名占一行称为主效主效应项。欲。欲在模型中包括几个主效在模型中包括几个主效应项,就,就进行几次如上的操作。行几次如上的操作。2024/8/1110C、建立模型中的交互、建立模型中的交互项 n例如,因素有例如,因素有Light(F)、Device(F)、)、Target(F),若要),若要求模型中包括求模型中包括变量量Light与与Device交互效交互效应。n首先定首先定义效效应类型型为Interactin,n然后在然后在FactorsCovariates框内的框内的变量表中,用鼠量表中,用鼠标单击Device变量使其背景改量使其背景改
7、变颜色,再用鼠色,再用鼠标单击变量量Light变量使其背景改量使其背景改变颜色;色;单击Build Term(s)栏内残数框的内残数框的箭箭头按按钮,一个交互效,一个交互效应出出现在在Model框中。模型增加了一框中。模型增加了一个交互效个交互效应项:Device*Light。2024/8/1111 C. 建立模型中的交互建立模型中的交互项n模型中包括三个模型中包括三个变量的所有量的所有2维交互效交互效应项,n定定义效效应类型型为All2-way,n单击light、Device、Target三个三个变量名,量名,n单击箭箭头按按钮。nModel中出中出现三个三个 2维交互效交互效应项: Lig
8、ht*Device、 Light*Target、Device*Target。n模型中包括所有模型中包括所有3维效效应,n定定义效效应类型型为All3way,n单击变量量Llight、 Device、 Target。n单击箭箭头按按钮,nModel框中出框中出现3维交互效交互效应项:Ligh*Device*Target。 2024/8/1112、选择平方和分解的方法平方和分解的方法nSum of squares:nTYPE(嵌套(嵌套设计)、)、 TYPE(平衡(平衡设计、仅主效主效应)、)、n TYPE (系(系统默默认、最常、最常用)用)nTYPEIV(不完整数据)。(不完整数据)。 202
9、4/8/11132.2 Contrasts按按钮nFactors框中框中显示出所有在示出所有在主主对话框中框中选中的因素,中的因素,其后的括号中是当前的其后的括号中是当前的对比方法了;比方法了;nChange Contrast栏中改中改变对照方法。照方法。2024/8/1114可供可供选择的的对照方法照方法nNone:不:不进行均数比行均数比较;nDeviation:比:比较预测变量或因素的每个水平的效量或因素的每个水平的效应。选择Last或或First作作为参照的水平;参照的水平;nSimple:对预测变量或因素量或因素变量的每一水平都与参照水平量的每一水平都与参照水平进行比行比较。选择La
10、st或或First作作为参照水平;参照水平;nDifference:对预测变量或因素每一水平的效量或因素每一水平的效应,除第一水平以,除第一水平以外,都与其前面各水平的平均效外,都与其前面各水平的平均效应进行比行比较。与。与Helmert对照方照方法相反;法相反; 注:注: 只有只有Deviation和和Simple 需要需要选择参考水平,参考水平,Last(系(系统默默认)和)和First。 2024/8/11152.3 Plots按按钮nFactor:主:主对话框中所框中所选因素因素变量名;量名;nHorizontal:横坐:横坐标框框nSeparate Lines:确定分:确定分线变量量
11、nSeparate Plots:确定分:确定分图变量量2024/8/11162.4 Post Hoc按按钮n均数多重比均数多重比较(事后(事后检验)2024/8/11172.5 Save按按钮(选择保存运算保存运算值)n通通过在在对话框中的框中的选择,可以将所可以将所计算的算的预测值、残差残差和和诊断断值(回(回归分析分析时)作作为新的新的变量保存在量保存在编辑数据文件中。以便在数据文件中。以便在其他其他统计分析中使用分析中使用这些些值。 2024/8/11182.5 Save按按钮(选择保存运算保存运算值) Predicted Values(预测值)nUnstandardized:非:非标准
12、化准化预测值nWeighted:如果在主:如果在主对话框框选择了了WLS变量,量,选中中该复复选项将保存加将保存加权非非标准化准化预测值nStandard Error:预测值的的标准准误 Residuals(残差(残差栏)nUnstandarized:非非标准化残差准化残差nWeighted: 加加权非非标准化残差准化残差nStandardized:标准化残差准化残差nStudentized:学生化残差学生化残差nDeleted:剔除残差剔除残差2024/8/11192.5 Save按按钮(选择保存运算保存运算值) Diagnostics(诊断断值栏)nCooks distance:Cook距
13、离;距离;nLeverage values:非中心化:非中心化Leverage值;Save to new filen将参数将参数协方差矩方差矩阵保存到一个保存到一个新文件中新文件中2024/8/11202.6 Option按按钮(选择输出出项)nDisplay Means for:显示分示分组因素因素nDisplay:指定:指定输出的出的统计量量nDescriptive statistics:描述:描述统计量,均量,均值、标准差,准差,样本量本量nEstimates Of effect size:效:效应量估量估计。nObserved power:检验假假设的功效。的功效。nParameter
14、 estimates:各因素:各因素变量的模型参量的模型参数估数估计、标准准误、t 检验的的t值、P值和和 95的置信区的置信区间。nSigificance level :指定:指定Confidence intervals的的显著性水平著性水平2024/8/1121nDescriptive statistics:描述:描述统计量,均量,均值、标准差,准差,样本量本量nEstimates Of effect size:效:效应量估量估计。nObserved power:检验假假设的功效。的功效。nParameter estimates:各因素:各因素变量的模型参数估量的模型参数估计、标准准误、t
15、 检验的的t值、P值和和95的置信区的置信区间。Contrast coefficient matrix:变换系数矩系数矩阵或或L矩矩阵。nHomogeneity tests:方差:方差齐性性检验。nSpread Vs level plot:绘制制观测量均量均值-标准差准差图、观测量均量均值-方方差差图。nResiduals plot:绘制残差制残差图。nLack of fit:检查因素和因因素和因变量量间的关系是否被充分描述。的关系是否被充分描述。nGeneral estimable function:可以根据一般估:可以根据一般估计函数自定函数自定义假假设检验。对比系数矩比系数矩阵的行与一般
16、估的行与一般估计函数是函数是线性性组合的。合的。2024/8/1122n例例1 为研究多酚保健研究多酚保健饮料料对急性缺氧的影响,将急性缺氧的影响,将60只只Wistar小白鼠随机分小白鼠随机分为低、中、高低、中、高三个三个剂量量组和和一个一个对照照组,每,每组15只小白鼠。只小白鼠。对照照组给予蒸予蒸馏水水0.25ml灌胃,低、中、高灌胃,低、中、高剂量量组分分别给予予2.0、4.0、8.0g/kg的的饮料溶于料溶于0.20.3ml蒸蒸馏水后灌胃,每天一水后灌胃,每天一次。次。40天后,天后,对小白鼠小白鼠进行耐缺氧存活行耐缺氧存活时间实验,结果如表果如表1。试比比较不同不同剂量的茶多酚保健
17、量的茶多酚保健饮料料对延延长小白鼠的平均耐缺氧存活小白鼠的平均耐缺氧存活时间有无差有无差别。二、完全随机二、完全随机设计资料的方差分析料的方差分析2024/8/1123对照组对照组低剂量组低剂量组中剂量组中剂量组高剂量组高剂量组21.31 23.4820.16 26.1335.07 24.3330.23 38.4723.14 20.3424.49 25.2428.11 33.9736.84 35.1027.48 26.9821.32 20.2324.74 21.8638.61 28.0119.54 19.5619.46 22.4729.79 28.6527.13 23.3718.03 17.3
18、925.63 29.3822.68 25.1328.79 28.4424.03 24.3728.81 20.1623.01 34.4433.24 34.2222.82 16.0118.74 22.5128.32 31.6931.68 35.0818.7218.4229.0428.29表表1 各组小白鼠耐缺氧时间各组小白鼠耐缺氧时间/min2024/8/1124n建立数据文件:耐缺氧建立数据文件:耐缺氧时间.sav.n定定义变量量2024/8/1125n建立数据文件:耐缺氧建立数据文件:耐缺氧时间.sav.n定定义变量量n输入数据入数据2024/8/1126n建立数据文件:耐缺氧建立数据文件:耐
19、缺氧时间.sav.n定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroupFixed Factors2024/8/1127n建立数据文件:耐缺氧建立数据文件:耐缺氧时间.sav.n定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroupFixed FactorsnPost Hoc : Group Post Hoc Tests for LSD,SNK
20、, Bonferroni2024/8/1128n建立数据文件:耐缺氧建立数据文件:耐缺氧时间.sav.n定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroupFixed FactorsnPost Hoc :nOptions:Group Display Means for Descriptive Statistics, Homogeneity tests2024/8/1129n主要主要结果果-描述性描述性统计量量2024/8/1130n主要主要结果果方差方差齐性性检验的的
21、P值方差方差齐性性检验的的统计量量F值结论:尚不能:尚不能认为各各总体方差不体方差不齐2024/8/1131方差分析方差分析的的P值方差分方差分析析F值方差分方差分析析MS值(均方)均方)方差分方差分析析自由自由度度方差分方差分析析SS值(平方和)平方和)总变异异组间变异异组内内变异异(误差差变异)异)结论:各:各组总体均数不等。体均数不等。2024/8/1132各各组总体均体均数的数的95%置置信区信区间按方差分析中按方差分析中的的误差均方差均方计算的算的标准准误2024/8/1133Bonferroni方法的方法的P值,是校正了是校正了a 后的后的2024/8/1134均数均数标在同一列的
22、在同一列的组间差异无差异无统计学意学意义,在不同,在不同列的差异有列的差异有统计学意学意义,即,即对照照组与低与低剂量差异无量差异无统计学意学意义,其他均有,其他均有统计学意学意义SNK检验2024/8/1135n用用Compare Means One-Way ANOVAn建立数据同前建立数据同前n开始分析:开始分析:analyze Compare Means One-Way ANOVAnYDependent ListnGroupFactor2024/8/1136n用用Compare Means One-Way ANOVAn建立数据同前建立数据同前n开始分析:开始分析:analyze Comp
23、are Means One-Way ANOVAnYDependent ListnGroupFactornPost Hoc : LSD,Bonferroni,SNK2024/8/1137n用用Compare Means One-Way ANOVAn建立数据同前建立数据同前n开始分析:开始分析:analyze Compare Means One-Way ANOVAnYDependent ListnGroupFactornPost Hoc : LSD,Bonferroni,SNKnOptions Statistics Descriptive, Homogeneity tests2024/8/1138
24、各各组总体均体均数的数的95%置置信区信区间按各自的方按各自的方差差计算的算的标准准误2024/8/1139n主要主要结果果方差方差齐性性检验的的P值方差方差齐性性检验的的统计量量F值2024/8/1140方差分析方差分析的的P值方差分方差分析析F值方差分方差分析析MS值(均方)均方)方差分方差分析析自由自由度度方差分方差分析析SS值(平方和)平方和)组内变异组内变异(误差变异)(误差变异)组间变异异总变总变异异2024/8/1141Bonferroni方法的方法的P值,是校正了是校正了a 后的后的2024/8/1142均数均数标在同一列的在同一列的组间差异无差异无统计学意学意义,在不同,在不
25、同列的差异有列的差异有统计学意学意义,即,即对照照组与低与低剂量差异无量差异无统计学意学意义,其他均有,其他均有统计学意学意义SNK检验2024/8/1143n练习1 某人研究北京机关工作人某人研究北京机关工作人员血脂水平,随机抽取不同年血脂水平,随机抽取不同年龄男性各男性各10名受名受试者,者,检测他他们的的总胆固醇胆固醇(TC)的含量的含量(mmol/L),其,其结果如下表:果如下表:青年组青年组中年组中年组老年组老年组5.004.785.125.145.245.234.855.185.135.165.265.214.934.894.894.985.234.985.185.075.205.
26、165.105.154.955.214.995.255.315.192024/8/1144三、随机区三、随机区组(单位位组)设计资料的方差分析料的方差分析n例例2 研究者欲比研究者欲比较生物蛋白粉生物蛋白粉饲料料、血血浆蛋白粉蛋白粉饲料料和和普通普通饲料料喂养喂养断奶仔猪断奶仔猪的的增重效果增重效果。为了消除了消除和控制其他因素的影响,研究者将断奶仔猪配成若和控制其他因素的影响,研究者将断奶仔猪配成若干干区区组(block),每个区,每个区组3只仔猪,并且只仔猪,并且满足同一区足同一区组的仔猪是的仔猪是同同窝别、同性同性别、同日同日龄、体重接近体重接近,共配成共配成10个区个区组。然后在每个区
27、。然后在每个区组内随机将内随机将3只仔猪只仔猪分配到各分配到各实验组。比。比较喂养喂养10天后各天后各实验仔猪的仔猪的平平均体重增加量均体重增加量(kg),结果果见表表2。试比比较各种各种饲料的料的增重效果有无差异。增重效果有无差异。2024/8/1145表表2 生物蛋白粉、血生物蛋白粉、血浆蛋白粉和普通蛋白粉和普通饲料料饲养仔猪增重量养仔猪增重量/Kg2024/8/1146n建立数据文件:仔猪增重量建立数据文件:仔猪增重量.savn定定义变量量2024/8/1147n建立数据文件:仔猪增重量建立数据文件:仔猪增重量.savn定定义变量量n输入数据入数据2024/8/1148n建立数据文件:仔
28、猪增重量建立数据文件:仔猪增重量.savn定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroup,BlockFixed Factors2024/8/1149n建立数据文件:仔猪增重量建立数据文件:仔猪增重量.savn定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroup,BlockFixed FactorsnModel Custom Main e
29、ffects(Group,Block)2024/8/1150n建立数据文件:仔猪增重量建立数据文件:仔猪增重量.savn定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroup,BlockFixed FactorsnModel Custom Main effects(Group,Block)nPost Hoc : Group Post Hoc Tests for SNK,Dunnett(first)2024/8/1151n建立数据文件:仔猪增重量建立数据文件:仔猪增重量.
30、savn定定义变量量n输入数据入数据n开始分析:开始分析:analyze General Linear Model UnivariatenYDependent VariablenGroup,BlockFixed FactorsnModel Custom Main effects(Group,Block)nPost Hoc :nOptions:Group,Block Display Means for 2024/8/1152总变异异处理理组间变异异误差差变异异区区组间变异异n主要主要结果果2024/8/1153各各组总体均体均数的数的95%置置信区信区间2024/8/11542024/8/115
31、5均数均数标在同一列的在同一列的组间差异无差异无统计学意学意义,在不同,在不同列的差异有列的差异有统计学意学意义,即血,即血浆蛋白与生物蛋白的蛋白与生物蛋白的增重效果差异无增重效果差异无统计学意学意义,其他均有,其他均有统计学意学意义2024/8/1156n练习2. 四个种系未成年雌性大白鼠各三只,每只按四个种系未成年雌性大白鼠各三只,每只按一种一种剂量注射雌激素,一段量注射雌激素,一段时间后,解剖称子后,解剖称子宫重量。重量。结果果见下表,下表,试分析注射不同雌激素分析注射不同雌激素剂量,子量,子宫重量重量是否不同。是否不同。种系种系剂剂 量量0.2(1)0.4(2)0.8(3)A(1)106116145B(2)4268115C(3)70111133D(4)4263872024/8/1157结束!结束!