特殊设计的方差分析.ppt

,SPSS软件应用 潘发明 2010年10月 安徽医科大学流行病与卫生统计学系,第二讲 方差分析及卡方检验,主要内容 1.单因素的方差分析（one-way anova）； 2.双因素的方差分析（two-way anova）； 3.析因设计的方差分析（factorial Design）； 4.重复设计的方差分析（repeated measures anova）；,One-way Anova (单因素方差分析),作用： 1、检验单一因素影响的多个样本是否来自均数相同的总体。要求各样本相互独立、均来自正态总体、总体方差齐。 2、均数之间两两比较,同时可以精确设定均数比较方式。 3、为什么多个均数之间的比较不能用两两t检验比较？ 可选入多个应变量；只能选入一个分组变量。 以lung.sav为例，分析三组石棉矿工的用力肺活量之间有无统计学意义。 结果包括： 各组变量值的描述统计；方差齐性检验； 方差分析表；各组均数的两两比较表。,例题,某职业病防治院对31名石棉矿工中的石棉肺患者、可疑患者及非患者进行了用力肺活量测定,问三组石棉矿工的用力肺活量有无差别? 具体数据见LUNG.SAV数据库。 问题1：如何建立数据库？（两个变量） 问题2：如何进行两两比较（ SNK ）？,Analyze compare means one-way anova Dependent list: lung 要分析的变量 Factor: group 分组变量 Options: Homogeneity-of variance 方差齐性检验 Continue Post Hoc: S-N-K:continue 两两比较的方法 Ok,两两比较方法的选择,两两方法的选择是个令人头疼的问题？ 一般可以参考以下标准： (1)如果存在明确的对照组，要进行的是验证性的研究，即计划好的某两个或几个组间（和对照组）的比较，宜用（Dunnett-t,LSD）法. (2)若需要进行的是多个均数的两两比较(探索性的研究),且各组人数相等,适宜用Tukey,SNK法;其他情况宜用 Scheffe法. (3)看看几种方法检验的结果是否一致然后选择比较的方法. (4)国内外都以SNK法最为常用,但不适宜于组数较多的情况.(因为假阳性 率较高) (5)常用比较方法很多，有多重比较，线形对比，正交对比等,配伍设计的方差分析,又称为双因素设计的方差分析（Randomized Block Design） 可以分为两种情况: (1)对同一受试对象在同一处理不同水平间的比较; (2)将几个受试对象按一定条件划分成配伍组,再将每一配伍组的各受试者 随机分配到各个处理组中去; *:需要特别指出的是：如果配伍因素是希望研究的因素，或者配伍因素与研究因素间的交互作用不能忽略，则不应当配伍设计，而应当采用下面的各种设计方法（如析因设计等）。,例题,某厂医务室测定了10名氟作业工人工前、工中及工后4小时的尿氟浓度，问氟作业工人在这三个不同时间的尿氟浓度有无差别？ 具体见数据库 twoway.sav 建立数据库 三个变量:其中主要的因素是time; 测定时间（time） 区组(worker) 尿氟浓度(fu),SPSS分析过程,Analyze General Lineal Model Univariate Dependent List: fu Fixed Factor: group work Model: Custom Build Terms:Main effects 模型中准备纳入主效应 Model: time、 worker 模型中只纳入主效应:时间和区组因素 Contiune: Post Hoc: Post Hoc Test For: time SNK Continue OK,析因设计,当几个因素间存在交互作用时，析因设计是一种非常理想的设计。 （1）检验两因素各水平之间是否有统计学意义； （2）可以检验两因素间的交互作用； 其优点： 可以节约样本含量，（1/21/3）； #：注意：当各因素间存在交互作用，此时各因子的主效应检验结果已经无实际意义，应当按各因素各水平的组合来研究。 具体的例题见数据库析因设计.SAV,Univariate（单变量过程）,介绍析因设计的双因素方差分析： 原始资料为： 整理后的表格为： 整理后的SPSS数据库格式见析因设计.sav（变量值为红细胞增加数）。,Univariate,问：A和B两药单独使用的治疗效果如何？两药联合使用的治疗效果如何？ Analyze GLM Univariate ，打开对话框，把变量放入应变量栏中，两药因素放入自变量栏中；单击“Plots”按钮，打开对话框，设置两因素的交叉作用。 结果中包括两药的单独作用以及联合作用的情况。,因为A、B两药存在交互作用，因此他们各自的主效应无实际意义 需要按以下组合来分析： （1）不使用乙药时甲药有无作用 （2）使用乙药时甲药有无作用 （3）不使用甲药时乙药有无作用 （4）使用甲药时乙药有无作用 这样才能真正准确地回答两种药物疗效如何的问题,重复测量数据的方差分析 （repeated measures anova）,重复测量数据的定义：是指某一应变量在同一时间点上被重复多次测量所得到的数据。 为什么不能采用简单模型因素方差分析（即使设置时间变量） 主要的原因： （1）同一个体的不同次测量间往往存在相关性（检验方法是球型检验-Mauchlys test of Sphericity）； （2）由于实验设计的要求或数据本身的特点:即使符合球型检验,采用单因素方差分析模型也非常困难;特别是在重复测量数据中研究者关心的是随着时间的变化,测量指标的变动趋势如何,而不关系各次测量间两两比较的结果，因此对不同次测量间它只提供趋势检验。,例题,见数据库anxity2.sav; Anxiety和 tension对实验结果(trial1trial4)有无影响;四次实验间有无差异;实验次数和两个变量间有无交互作用.,具体的操作过程 Analyze General lineal model Repeated measures Within-subject factor name: trial 定义重复测量的变量名 Number of levels 4 add 重复测量的次数为4次 Define Within-subject variables:trial1trial4 代表4次测量结果 Between subjects factor: anxiety tension Model Custom Within-subject Model :trial 分析4次测量间有无趋势 Between subjects Model : anxiety tension 只分析主效应 Continue OK,多元方差分析结果:主要是针对组内变量TRIAL以及和引入的另外两个自变量的交互作用是否有统计学意义,系统给了四种统计学方法,但一般结果是同向的,以上是球型检验的结果:卡方=9.383 P=0.097 符合球形分布,以一元方差结果为准 如果不符合球形分布,后面有三种EPSILON三种校正方法,应该以多元方差分析为准,上表结果是采用一元方差分析对组内因素进行调整,本例题是以该分析结果为准。,分类资料统计分析,Analyze Descriptive Statistics Crosstabs,Crosstables 过程 （列联表分析）,1 确定行变量、列变量、频数变量； 2 必须指定频数变量，Data Weight Cases； 3 根据不同资料、不同的研究目的，在 “Statistics”对话框中选入不同的分析指标。 4 在结果中，根据不同的研究目的，或根 据四格表资料的样本例数和各格子中的理 论频数大小，选择所需要的指标。,四格表资料2检验的适用条件,T为理论频数， n为样本含量。 1、T5，n40 ，选择不校正的2值； 2、140 ，选择连续性校正后的2值； 3、T1 或 n40，选择确切概率法的2值。,一般四格表,问：两种药的愈合率是否相同？在“Statistics”对话框中选中“Chi-square”即可数据库PPPP.SAV,配对分类资料,问：甲、乙两法检验结果有无差别？ 在“Sta