析因设计与方差分析.docxVIP

析因设计与方差分析1. 基本概念完全随机设计(单因素)随机区组设计(两因素, 无重复)拉丁方设计(三因素, 无重复)析因设计(两因素以上, 至少重复2次以上)析因设计的意义在评价药物疗效时，除需知道A药和B药各剂量的疗效外（主效应），还需知道两种药同时使用的协同疗效。析因设计及相应的方差分析能分析药物的单独效应、主效应和交互效应。例：A因素食物中蛋白含量;B因素食物中脂肪含量A因素缺乏(a1)正常(a2)B因素缺乏(b1)a1b1a2b1正常(b2)a1b2a2b2BA平均a2-a1a1a223644408平均333835.55b2-b16129(1)单独效应: 在每个B水平, A的效应。或在每个A水平，B的效应。(2)主效应：某因素各水平的平均差别。(3)交互效应：某因素各水平的单独效应随另一因素水平变化而变化，则称两因素间存在交互效应。如果，存在交互效应。如果，协同作用。如果，拮抗作用。如果不存在交互效应，则只需考虑各因素的主效应。在方差分析中，如果存在交互效应，解释结果时，要逐一分析各因素的单独效应，找出最优搭配。在两因素析因设计时，只需考虑一阶交互效应。三因素以上时，除一阶交互效应外，还需考虑二阶、三阶等高阶交互效应，解释将更复杂。析因设计的优点：用相对较小样本，获取更多的信息，特别是交互效应分析。析因设计的缺点：当因素增加时，实验组数呈几何倍数增加。实际工作中部分交互效应，特别是高阶交互效应可以根据临床知识排除，这时可选用正交设计。2. 析因设计与结果的方差分析(1)实验设计设有k个因素，每个因素有L1, L2, …, Lk个水平，那么共有G= L1×L2×…×Lk个处理组。例如有三个因素，分别是A，B，C。A因素有两水平，B因素有3水平和C因素有2水平，则共有G=2×3×2=12个处理组。大家可以自己回去将这12种组合排列出来。确定了处理组数后，将实验对象分配到各组的方法可以采用完全随机设计、随机区组设计或拉丁方设计。注意：析因设计的基本要求，各组例数相等，每组例数必须2例以上。(2)析因设计资料的方差分析第一步：与一般的方差分析一样，将总变异分离成组间变异和组内变异。如果是随机区组设计还需从组内变异分离出单位组间变异和误差变异。方差来源DFSSMS总变异(T)N-1组间变异(B)G-1SSB/(G-1)组内变异(E)N-GSST-SSBSSE/(N-G)Tk(k=1,2,…,G)为各处理组观察值小计，r为各处理组例数，C=(ΣX)2/N第二步：将组间变异分解出主效应项和交互效应项，以两因素析因设计为例，i和j分别是因素A和因素B的水平数，Ai和Bj分别是各水平观察值的小计。方差来源DFSSMSF主效应Ai-1SS(A)/dfaMS(A)/MSEBj-1SS(B)/dfbMS(B)/MSE交互效应AB(i-1)(j-1)SSB-SS(A)-SS(B)SS(AB)/dfabMS(AB)/MSE两个因素以上的析因设计，计算原理类似，但手工计算较繁琐。当有计算机后和统计软件的帮助后，已完全没必要手工计算。但是掌握变异来源分解的基本思想很重要，应该将哪项变异作为误差项，如何解释结果都与此有关。下面用实例介绍计算过程：A(缝合方法)外膜缝合(a1)束膜缝合(a2)总计B(缝合时间)1月(b1)2月(b2)1月(b1)2月(b2)110301050210302050340703070450605060510303030Tk120220140260740Σx244001120048001440034800C=(740)2/20=27380方差分析表SSDFMSFPT74201932.9110.05E480016300A1=120+220=340, A2=400, B1=260, B2=480600.05B2420124208.070.05AB201200.070.05结论：缝合时间（B）的主效应有统计学意义，即B的主效应=48-26=22。第五节裂区设计与结果的方差分析(1)基本概念裂区设计与一般析因设计的区别在于每种处理因素分别作用于不同级别的实验单位。如眼科实验中，兔子为一级实验单位，每只兔子的两只眼睛为二级实验单位。当处理因素分别作用于一级实验单位和二级实验单位时，称裂区设计。如果将作用于二级实验单位的处理因素称为二级处理，作用于一级实验单位的处理因素称为一级处理。显然前者为区组设计，后者为完全随机设计，两种处理的设计精度不同。因此又称这类设计为不完整析因设计。一般在设计时，常选最感兴趣的主要研究因素为二级处理因素或出于区组的限制，选水平数少的研究因素为二级处理因素(2)实验设计方法设一级处理因素有i个水平，二级处理因素有j个水平。一级实验单位的分配方式可采用完全