中一般线性模型.ppt

*,1,一般线性模型（一）,*,2,一般线性模型,一般线性模型单变量分析的基本过程完全随机设计资料的方差分析随机区组（单位组）设计资料的方差分析,*,3,一、一般线性模型单变量分析的基本过程,General Linear Model（GLM，一般线性模型）包括：Univariate（单因变量多因素方差分析），Multivariate（多因变量方差分析），Repeated Measures（重复测量方差分析），Variance（方差分量分析）GLM可完成多自变量、多水平、多因变量、重复测量方差分析以及协方差分析等。,*,4,Univariate（单因变量方差分析）基本过程,*,5,1 主对话框,Dependent Variable:因变量Fixed Facter:固定因子，所有可能的水平都出现在样本中，如分组等Random Facter:随机因子，所有可能的取值并不都在样本中出现，如观察个体Covariates:协变量，协方差分析时用WLS Weight:WLS权重。用于加权最小二乘分析。,*,6,2 功能按钮,Model：分析模型Contrast：对照方法Plots：分布图形Post Hoc：多重比较Save：保存运算值Option：选择输出项。,*,7,2.1 Model按钮、在Specify Model栏中指定模型类型,Full Factorial，全模型，系统默认。包括所有因素的主效应和所有的交互效应。例如有三个因素变量，全模型包括三个因素的主效应、两两的交互效应和三个因素的高级交互效应。Custom，自定义模型。选择此项激活下面各操作框,*,8,、建立自定义模型,FactorsCovariates 框中自动列出可以作为因素的变量名，其后面的括号中标有字母“F”（固定因子）、“R”（随机因子）或者“C”（协变量）。,*,9,A、选择效应类型,Interactin：交互效应Main effects：主效应All2-way：所有2维交互效应All3-way：所有3维交互效应All4-Way：所有4维交互效应All5-Way：所有5维交互效应,*,10,B、选择模型中的主效应（Model）,首先定义效应类型为Main effects鼠标键单击某一个因素，该变量名背景将改变颜色（一般变为蓝色），单击Build Term（s）栏中下面的箭头，该变量出现在Mode1中。一个变量名占一行称为主效应项。欲在模型中包括几个主效应项，就进行几次如上的操作。,*,11,C、建立模型中的交互项,例如，因素有Light(F)、Device（F）、Target（F），若要求模型中包括变量Light与Device交互效应。首先定义效应类型为Interactin，然后在FactorsCovariates框内的变量表中，用鼠标单击Device变量使其背景改变颜色，再用鼠标单击变量Light变量使其背景改变颜色；单击Build Term（s）栏内残数框的箭头按钮，一个交互效应出现在Model框中。模型增加了一个交互效应项：Device*Light。,*,12,C.建立模型中的交互项,模型中包括三个变量的所有2维交互效应项，定义效应类型为All2-way，单击light、Device、Target三个变量名，单击箭头按钮。Model中出现三个 2维交互效应项：Light*Device、Light*Target、Device*Target。模型中包括所有3维效应，定义效应类型为All3way，单击变量Llight、Device、Target。单击箭头按钮，Model框中出现3维交互效应项：Ligh*Device*Target。,*,13,、选择平方和分解的方法,Sum of squares:TYPE（嵌套设计）、TYPE（平衡设计、仅主效应）、TYPE（系统默认、最常用）TYPEIV（不完整数据）。,*,14,2.2 Contrasts按钮,Factors框中显示出所有在主对话框中选中的因素，其后的括号中是当前的对比方法了；Change Contrast栏中改变对照方法。,*,15,可供选择的对照方法,None：不进行均数比较；Deviation：比较预测变量或因素的每个水平的效应。选择Last或First作为参照的水平；Simple：对预测变量或因素变量的每一水平都与参照水平进行比较。选择Last或First作为参照水平；Difference：对预测变量或因素每一水平的效应，除第一水平以外，都与其前面各水平的平均效应进行比较。与Helmert对照方法相反；,注：只有Deviation和Simple 需要选择参考水平，Last（系统默认）和First。,*,16,2.3 Plots按钮,Factor：主对话框中所选因素变量名；Horizontal：横坐标框Separate Lines：确定分线变量Separate Plots：确定分图变量,*,17,2.4 Post Hoc按钮,均数多重比较（事后检验）,*,18,2.5 Save按钮(选择保存运算值),通过在对话框中的选择，可以将所计算的预测值、残差和诊断值（回归分析时）作为新的变量保存在编辑数据文件中。以便在其他统计分析中使用这些值。,*,19,2.5 Save按钮(选择保存运算值),Predicted Values（预测值）Unstandardized：非标准化预测值Weighted：如果在主对话框选择了WLS变量，选中该复选项将保存加权非标准化预测值Standard Error：预测值的标准误,Residuals（残差栏）Unstandarized:非标准化残差Weighted:加权非标准化残差Standardized:标准化残差Studentized:学生化残差Deleted:剔除残差,*,20,2.5 Save按钮(选择保存运算值),Diagnostics（诊断值栏）Cooks distance：Cook距离；Leverage values：非中心化Leverage值；Save to new file将参数协方差矩阵保存到一个新文件中,*,21,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的显著性水平,*,22,Descriptive 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：检查因素和因变量间的关系是否被充分描述。General estimable function：可以根据一般估计函数自定义假设检验。对比系数矩阵的行与一般估计函数是线性组合的。,*,23,例1 为研究多酚保健饮料对急性缺氧的影响，将60只Wistar小白鼠随机分为低、中、高三个剂量组和一个对照组，每组15只小白鼠。对照组给予蒸馏水0.25ml灌胃，低、中、高剂量组分别给予2.0、4.0、8.0g/kg的饮料溶于0.20.3ml蒸馏水后灌胃，每天一次。40天后，对小白鼠进行耐缺氧存活时间实验，结果如表1。试比较不同剂量的茶多酚保健饮料对延长小白鼠的平均耐缺氧存活时间有无差别。,二、完全随机设计资料的方差分析,*,24,表1 各组小白鼠耐缺氧时间/min,*,25,建立数据文件：耐缺氧时间.sav.定义变量,*,26,建立数据文件：耐缺氧时间.sav.定义变量输入数据,*,27,建立数据文件：耐缺氧时间.sav.定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroupFixed Factors,*,28,建立数据文件：耐缺氧时间.sav.定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroupFixed FactorsPost Hoc：Group Post Hoc Tests for LSD，SNK,Bonferroni,*,29,建立数据文件：耐缺氧时间.sav.定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroupFixed FactorsPost Hoc：Options：Group Display Means for Descriptive Statistics，Homogeneity tests,*,30,主要结果-描述性统计量,*,31,主要结果,方差齐性检验的P值,方差齐性检验的统计量F值,结论：尚不能认为各总体方差不齐,*,32,方差分析的P值,方差分析F值,方差分析MS值(均方）,方差分析自由度,方差分析SS值(平方和）,总变异,组间变异,组内变异（误差变异）,结论：各组总体均数不等。,*,33,各组总体均数的95%置信区间,按方差分析中的误差均方计算的标准误,*,34,Bonferroni方法的P值，是校正了a 后的,*,35,均数标在同一列的组间差异无统计学意义，在不同列的差异有统计学意义，即对照组与低剂量差异无统计学意义，其他均有统计学意义,SNK检验,*,36,用Compare Means One-Way ANOVA建立数据同前开始分析：analyze Compare Means One-Way ANOVAYDependent ListGroupFactor,*,37,用Compare Means One-Way ANOVA建立数据同前开始分析：analyze Compare Means One-Way ANOVAYDependent ListGroupFactorPost Hoc：LSD，Bonferroni，SNK,*,38,用Compare Means One-Way ANOVA建立数据同前开始分析：analyze Compare Means One-Way ANOVAYDependent ListGroupFactorPost Hoc：LSD，Bonferroni，SNKOptions Statistics Descriptive，Homogeneity tests,*,39,各组总体均数的95%置信区间,按各自的方差计算的标准误,*,40,主要结果,方差齐性检验的P值,方差齐性检验的统计量F值,*,41,方差分析的P值,方差分析F值,方差分析MS值(均方）,方差分析自由度,方差分析SS值(平方和）,组内变异（误差变异）,组间变异,总变异,*,42,Bonferroni方法的P值，是校正了a 后的,*,43,均数标在同一列的组间差异无统计学意义，在不同列的差异有统计学意义，即对照组与低剂量差异无统计学意义，其他均有统计学意义,SNK检验,*,44,练习1 某人研究北京机关工作人员血脂水平，随机抽取不同年龄男性各10名受试者，检测他们的总胆固醇(TC)的含量(mmol/L)，其结果如下表：,*,45,三、随机区组（单位组）设计资料的方差分析,例2 研究者欲比较生物蛋白粉饲料、血浆蛋白粉饲料和普通饲料喂养断奶仔猪的增重效果。为了消除和控制其他因素的影响，研究者将断奶仔猪配成若干区组(block)，每个区组3只仔猪，并且满足同一区组的仔猪是同窝别、同性别、同日龄、体重接近，共配成10个区组。然后在每个区组内随机将3只仔猪分配到各实验组。比较喂养10天后各实验仔猪的平均体重增加量(kg)，结果见表2。试比较各种饲料的增重效果有无差异。,*,46,表2 生物蛋白粉、血浆蛋白粉和普通饲料饲养仔猪增重量/Kg,*,47,建立数据文件：仔猪增重量.sav定义变量,*,48,建立数据文件：仔猪增重量.sav定义变量输入数据,*,49,建立数据文件：仔猪增重量.sav定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroup，BlockFixed Factors,*,50,建立数据文件：仔猪增重量.sav定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroup，BlockFixed FactorsModel Custom Main effects（Group，Block）,*,51,建立数据文件：仔猪增重量.sav定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroup，BlockFixed FactorsModel Custom Main effects（Group，Block）Post Hoc：Group Post Hoc Tests for SNK，Dunnett（first）,*,52,建立数据文件：仔猪增重量.sav定义变量输入数据开始分析：analyze General Linear Model UnivariateYDependent VariableGroup，BlockFixed FactorsModel Custom Main effects（Group，Block）Post Hoc：Options：Group，Block Display Means for,*,53,总变异,处理组间变异,误差变异,区组间变异,主要结果,*,54,各组总体均数的95%置信区间,*,55,*,56,均数标在同一列的组间差异无统计学意义，在不同列的差异有统计学意义，即血浆蛋白与生物蛋白的增重效果差异无统计学意义，其他均有统计学意义,*,57,练习2.四个种系未成年雌性大白鼠各三只，每只按一种剂量注射雌激素，一段时间后，解剖称子宫重量。结果见下表，试分析注射不同雌激素剂量，子宫重量是否不同。,