spss中一般线性模型课件

1、一般线性模型（一）,1,学习交流PPT,一般线性模型,一般线性模型单变量分析的基本过程 完全随机设计资料的方差分析 随机区组（单位组）设计资料的方差分析,2,学习交流PPT,一、一般线性模型单变量分析的基本过程,General Linear Model（GLM，一般线性模型） 包括： Univariate（单因变量多因素方差分析）， Multivariate（多因变量方差分析）， Repeated Measures（重复测量方差分析）， Variance（方差分量分析） GLM可完成多自变量、多水平、多因变量、重复测量方差分析以及协方差分析等。,3,学习交流PPT,Univariate（单因变

2、量方差分析）基本过程,4,学习交流PPT,1 主对话框,Dependent Variable:因变量 Fixed Facter: 固定因子，所有可能的水平都出现在样本中，如分组等 Random Facter: 随机因子，所有可能的取值并不都在样本中出现，如观察个体 Covariates:协变量，协方差分析时用 WLS Weight: WLS权重。用于加权最小二乘分析。,5,学习交流PPT,2 功能按钮,Model：分析模型 Contrast：对照方法 Plots：分布图形 Post Hoc：多重比较 Save：保存运算值 Option：选择输出项。,6,学习交流PPT,2.1 Model按钮、

3、在Specify Model栏中指定模型类型,Full Factorial，全模型，系统默认。包括所有因素的主效应和所有的交互效应。例如有三个因素变量，全模型包括三个因素的主效应、两两的交互效应和三个因素的高级交互效应。 Custom，自定义模型。选择此项激活下面各操作框,7,学习交流PPT,、建立自定义模型,FactorsCovariates 框中自动列出可以作为因素的变量名，其后面的括号中标有字母“F”（固定因子）、“R”（随机因子）或者“C”（协变量）。,8,学习交流PPT,A、选择效应类型,Interactin：交互效应 Main effects：主效应 All2-way： 所有2维交

4、互效应 All3-way：所有3维交互效应 All4-Way：所有4维交互效应 All5-Way：所有5维交互效应,9,学习交流PPT,B、选择模型中的主效应 （Model）,首先定义效应类型为Main effects 鼠标键单击某一个因素，该变量名背景将改变颜色（一般变为蓝色），单击Build Term（s）栏中下面的箭头，该变量出现在Mode1中。一个变量名占一行称为主效应项。欲在模型中包括几个主效应项，就进行几次如上的操作。,10,学习交流PPT,C、建立模型中的交互项,例如，因素有Light(F)、Device（F）、Target（F），若要求模型中包括变量Light与Device交互

5、效应。 首先定义效应类型为Interactin， 然后在FactorsCovariates框内的变量表中，用鼠标单击Device变量使其背景改变颜色，再用鼠标单击变量Light变量使其背景改变颜色；单击Build Term（s）栏内残数框的箭头按钮，一个交互效应出现在Model框中。模型增加了一个交互效应项：Device*Light。,11,学习交流PPT,C. 建立模型中的交互项,模型中包括三个变量的所有2维交互效应项， 定义效应类型为All2-way， 单击light、Device、Target三个变量名， 单击箭头按钮。 Model中出现三个 2维交互效应项： Light*Device、

6、 Light*Target、Device*Target。 模型中包括所有3维效应， 定义效应类型为All3way， 单击变量Llight、 Device、 Target。 单击箭头按钮， Model框中出现3维交互效应项：Ligh*Device*Target。,12,学习交流PPT,、选择平方和分解的方法,Sum of squares: TYPE（嵌套设计）、 TYPE（平衡设计、仅主效应）、 TYPE （系统默认、最常用） TYPEIV（不完整数据）。,13,学习交流PPT,2.2 Contrasts按钮,Factors框中显示出所有在主对话框中选中的因素，其后的括号中是当前的对比方法了；

7、Change Contrast栏中改变对照方法。,14,学习交流PPT,可供选择的对照方法,None：不进行均数比较； Deviation：比较预测变量或因素的每个水平的效应。选择Last或First作为参照的水平； Simple：对预测变量或因素变量的每一水平都与参照水平进行比较。选择Last或First作为参照水平； Difference：对预测变量或因素每一水平的效应，除第一水平以外，都与其前面各水平的平均效应进行比较。与Helmert对照方法相反；,注： 只有Deviation和Simple 需要选择参考水平，Last（系统默认）和First。,15,学习交流PPT,2.3 Plots

8、按钮,Factor：主对话框中所选因素变量名； Horizontal：横坐标框 Separate Lines：确定分线变量 Separate Plots：确定分图变量,16,学习交流PPT,2.4 Post Hoc按钮,均数多重比较（事后检验）,17,学习交流PPT,2.5 Save按钮(选择保存运算值),通过在对话框中的选择，可以将所计算的预测值、残差和诊断值（回归分析时）作为新的变量保存在编辑数据文件中。以便在其他统计分析中使用这些值。,18,学习交流PPT,2.5 Save按钮(选择保存运算值),Predicted Values（预测值） Unstandardized：非标准化预测值 W

9、eighted：如果在主对话框选择了WLS变量，选中该复选项将保存加权非标准化预测值 Standard Error：预测值的标准误,Residuals（残差栏） Unstandarized:非标准化残差 Weighted: 加权非标准化残差 Standardized:标准化残差 Studentized:学生化残差 Deleted:剔除残差,19,学习交流PPT,2.5 Save按钮(选择保存运算值),Diagnostics（诊断值栏） Cooks distance：Cook距离； Leverage values：非中心化Leverage值； Save to new file 将参数协方差矩阵保

10、存到一个新文件中,20,学习交流PPT,2.6 Option按钮(选择输出项),Display Means for：显示分组因素 Display：指定输出的统计量 Descriptive statistics：描述统计量，均值、标准差，样本量 Estimates Of effect size：效应量估计。 Observed power：检验假设的功效。 Parameter estimates：各因素变量的模型参数估计、标准误、t 检验的t值、P值和 95的置信区间。 Sigificance level ：指定Confidence intervals的显著性水平,21,学习交流PPT,Descr

11、iptive statistics：描述统计量，均值、标准差，样本量 Estimates Of effect size：效应量估计。 Observed power：检验假设的功效。 Parameter estimates：各因素变量的模型参数估计、标准误、t 检验的t值、P值和95的置信区间。Contrast coefficient matrix：变换系数矩阵或L矩阵。 Homogeneity tests：方差齐性检验。 Spread Vs level plot：绘制观测量均值-标准差图、观测量均值-方差图。 Residuals plot：绘制残差图。 Lack of fit：检查因素和因变量

12、间的关系是否被充分描述。 General estimable function：可以根据一般估计函数自定义假设检验。对比系数矩阵的行与一般估计函数是线性组合的。,22,学习交流PPT,例1 为研究多酚保健饮料对急性缺氧的影响，将60只Wistar小白鼠随机分为低、中、高三个剂量组和一个对照组，每组15只小白鼠。对照组给予蒸馏水0.25ml灌胃，低、中、高剂量组分别给予2.0、4.0、8.0g/kg的饮料溶于0.20.3ml蒸馏水后灌胃，每天一次。40天后，对小白鼠进行耐缺氧存活时间实验，结果如表1。试比较不同剂量的茶多酚保健饮料对延长小白鼠的平均耐缺氧存活时间有无差别。,二、完全随机设计资料的

13、方差分析,23,学习交流PPT,表1 各组小白鼠耐缺氧时间/min,24,学习交流PPT,建立数据文件：耐缺氧时间.sav. 定义变量,25,学习交流PPT,建立数据文件：耐缺氧时间.sav. 定义变量 输入数据,26,学习交流PPT,建立数据文件：耐缺氧时间.sav. 定义变量 输入数据 开始分析：analyze General Linear Model Univariate YDependent Variable GroupFixed Factors,27,学习交流PPT,建立数据文件：耐缺氧时间.sav. 定义变量 输入数据 开始分析：analyze General Linear Mod

14、el Univariate YDependent Variable GroupFixed Factors Post Hoc ： Group Post Hoc Tests for LSD，SNK, Bonferroni,28,学习交流PPT,建立数据文件：耐缺氧时间.sav. 定义变量 输入数据 开始分析：analyze General Linear Model Univariate YDependent Variable GroupFixed Factors Post Hoc ： Options：Group Display Means for Descriptive Statistics， H

15、omogeneity tests,29,学习交流PPT,主要结果-描述性统计量,30,学习交流PPT,主要结果,方差齐性检验的P值,方差齐性检验的统计量F值,结论：尚不能认为各总体方差不齐,31,学习交流PPT,方差分析的P值,方差分析F值,方差分析MS值(均方）,方差分析自由度,方差分析SS值(平方和）,总变异,组间变异,组内变异（误差变异）,结论：各组总体均数不等。,32,学习交流PPT,各组总体均数的95%置信区间,按方差分析中的误差均方计算的标准误,33,学习交流PPT,Bonferroni方法的P值，是校正了a 后的,34,学习交流PPT,均数标在同一列的组间差异无统计学意义，在不同

16、列的差异有统计学意义，即对照组与低剂量差异无统计学意义，其他均有统计学意义,SNK检验,35,学习交流PPT,用Compare Means One-Way ANOVA 建立数据同前 开始分析：analyze Compare Means One-Way ANOVA YDependent List GroupFactor,36,学习交流PPT,用Compare Means One-Way ANOVA 建立数据同前 开始分析：analyze Compare Means One-Way ANOVA YDependent List GroupFactor Post Hoc ： LSD，Bonferron

17、i，SNK,37,学习交流PPT,用Compare Means One-Way ANOVA 建立数据同前 开始分析：analyze Compare Means One-Way ANOVA YDependent List GroupFactor Post Hoc ： LSD，Bonferroni，SNK Options Statistics Descriptive， Homogeneity tests,38,学习交流PPT,各组总体均数的95%置信区间,按各自的方差计算的标准误,39,学习交流PPT,主要结果,方差齐性检验的P值,方差齐性检验的统计量F值,40,学习交流PPT,方差分析的P值,方

18、差分析F值,方差分析MS值(均方）,方差分析自由度,方差分析SS值(平方和）,组内变异（误差变异）,组间变异,总变异,41,学习交流PPT,Bonferroni方法的P值，是校正了a 后的,42,学习交流PPT,均数标在同一列的组间差异无统计学意义，在不同列的差异有统计学意义，即对照组与低剂量差异无统计学意义，其他均有统计学意义,SNK检验,43,学习交流PPT,练习1 某人研究北京机关工作人员血脂水平，随机抽取不同年龄男性各10名受试者，检测他们的总胆固醇(TC)的含量(mmol/L)，其结果如下表：,44,学习交流PPT,三、随机区组（单位组）设计资料的方差分析,例2 研究者欲比较生物蛋白

19、粉饲料、血浆蛋白粉饲料和普通饲料喂养断奶仔猪的增重效果。为了消除和控制其他因素的影响，研究者将断奶仔猪配成若干区组(block)，每个区组3只仔猪，并且满足同一区组的仔猪是同窝别、同性别、同日龄、体重接近，共配成10个区组。然后在每个区组内随机将3只仔猪分配到各实验组。比较喂养10天后各实验仔猪的平均体重增加量(kg)，结果见表2。试比较各种饲料的增重效果有无差异。,45,学习交流PPT,表2 生物蛋白粉、血浆蛋白粉和普通饲料饲养仔猪增重量/Kg,46,学习交流PPT,建立数据文件：仔猪增重量.sav 定义变量,47,学习交流PPT,建立数据文件：仔猪增重量.sav 定义变量 输入数据,48,

20、学习交流PPT,建立数据文件：仔猪增重量.sav 定义变量 输入数据 开始分析：analyze General Linear Model Univariate YDependent Variable Group，BlockFixed Factors,49,学习交流PPT,建立数据文件：仔猪增重量.sav 定义变量 输入数据 开始分析：analyze General Linear Model Univariate YDependent Variable Group，BlockFixed Factors Model Custom Main effects（Group，Block）,50,学习交流PPT,建立数据文件：仔猪增重