巢式嵌套设计的方差分析

《巢式嵌套设计的方差分析》由会员分享，可在线阅读，更多相关《巢式嵌套设计的方差分析（4页珍藏版）》请在装配图网上搜索。

1、巢式嵌套设计的方差分析巢式嵌套设计的方差分析嵌套设计被称为巢式设计（nested design）或系统分组(hierarchal classification)的资料，有些教科书上称这类资料为组内又分亚组的分类资料。根据因素数的不同，套设计可分为二因素（二级）、三因素（三级）等套设计。将全部 k 个因素按主次排列，依次称为 1 级，2 级 k 级因素，再将总离差平和及自由度进行分解，其基本思想与一般方差分析相同。所不同的是分解法有明显的区别，它侧重于主要因素，并且，第 i 级因素的显著与否，是分别用第 i 级与第 i+1 级因素的均方为分子和分母来构造 F 统计量，并以 F测验为其理论根据的。

2、嵌套试验设计应用：(1)情形一受试对象本身具有按其隶属关系进行分组再分组的各种因素。比如：选取某种植物 3 个品种（A），在每一株内选取两片叶子（B）（嵌套在植株因素下的第二个因素），用取样器从每一片叶子上选取同样面积的两个样本（两次重复），称取湿重。不是把 B 因素的 2 个水平简单地看作是与 A 因素 3 个水平的全面组合，而是分别嵌套在 A1、A2、A3三个水平之下，相当于 B 因素有 6 个水平，但它们所产生的离差平和中又包含了 A 因素的作用，一般来说，它大于模型所提供的总误差，用它作为度量 A 因素作用大小的误差项，是严格考核 A 因素的一种措施。方差分析表：变异来源AB（A）E总

3、和自由度df1=a-1df2=a(b-1)dfe=(n-1)abnab-1平方和SS1=Q1-CSS2=Q2-Q1SSe=W-Q2SSt=W-C均方MS1=SS1/df1MS2=SS2/df2MSe=SSe/dfeF 值MS1/MS2MS2/MSeProbF注：A 因素有 a 个水平，B 因素有 b 个水平，重复 n 次，总样本含量为 nab，全部数据之总和为 T，全部数据之平方和为 W，校正数为 C，Qi 为系统分组分到第 i 级因素时，各小组数据之和的平方均值的和。W=(2)情形二受试对象本身并非具有分组再分组的各种分组因素，而是各之间在专业上有主次之分。区分嵌套设计与析因设计的关键是看因

4、素之间的地位是否平等，因素的地位平等则属于析因设计，不平等则属于嵌套设计。比如：为了研究某种抗菌药的效果，对小白鼠进行试验。考虑 3 个试验因素，因素 A(用此抗菌药与否)可分为 A1(对照组不用抗菌药)、A2(试验组用抗菌药)；因素 B(小白鼠代次)可分为 B1(第 1 代)、B2(第2 代)、B3(第 3 代)；因素 C(性别)可分为 C1(雄性)、C2(雌性)。让第 1 代小白鼠被这种细菌感染，按雌雄分别统计对照组和试验组小白鼠的存活率，让存活的小白鼠分别在各自的组内(指对照组和试验组)生产第二代，对于第 2 代重复上面的试验，3 因素各水平搭配下都重复次试验，每次都有足够数量的小白鼠(

5、因为观测指标是存活率)，有人进行了如下的设计并收集到如下的试验数据。由专业知识得知：3 因素的主次顺序为 ABC。试分析 3 因素对小白鼠存活率的影响。方差分析表：21TC=nabQ1=nbxijk2T2i.Q2=1nT2ij.变异来源AB（A）C（A B）E总和自由度df1=a-1df2=a(b-1)df3=ab(c-1)dfe=(n-1)abcnabc-1平方和SS1=Q1-CSS2=Q2-Q1SS3=Q3-Q2SSe=W-Q3SSt=W-C均方MS1=SS1/df1MS2=SS2/df2MS3=SS3/df3MSe=SSe/dfeF 值MS1/MS2MS2/MS3MS3/MSeProbF

6、注：A 因素有 a 个水平，B 因素有 b 个水平，C 因素有 c 个水平，重复 n 次，总样本含量为 nab，全部数据之总和为 T，全部数据之平和为 W，校正数为 C=数据之和的平方均值的和。注意：由于在嵌套设计中，两个因素之间不能自由交错地组成各种处理组合，因而，不能考察因素之间地交互作用。所以，凡是交互作用比较重要的试验，都不应采用这种设计。进行嵌套设计的方差分析的过程：nested、anova 或 GLM，过程 nested 过程更简单一些，但它受数据中读取分组变量先后顺序的影响很大，应先用 SORT 过程按因素重要程度进行排序整理；过程 anova或 GLM 不需要排序整理，且容易进

7、行均值比较。平衡资料用过程 nested、anova 或 GLM 三种过程做结果是一样的；若不是平衡资料，必须用过程 GLM。例例 1 1选取某种植物 3 个品种（Plant），在每一株内选取两片叶子（Leaf）（嵌套在植株因素下的第二个因素），用取样器从每一片叶子上选取同样面积的两个样本（两次重复），称取湿重。对以上结果进行方差分析。（固定效应 Plant+随机效应 Leaf）解：解：SAS 程序如下：DataNested;Inputplant$leafwt;Cards;a1a1b1b1c1c1a2a2b2b2c2c2;ProcAnova;Proc Sort;By plant leaf;Ru

8、n;Proc Nested;Class plant leaf;Var wt;Classplant leaf;Modelwt=plant leaf(plant);TestH=plantE=leaf(plant);Means plant/Duncan E=leaf(plant);T2nabc，Qi 为系统分组分到第 i 级因素时，各小组Run;Run;输出结果见示例。例例 2 2在室内用 4 种培养液培养某植物，每种培养液培养 3 盆，每盆 4 株，一个月后测定其株高生长量（mm）。试对以上结果进行方差分析。（固定效应+随机效应）Data amm;Do a=1 to 4;Do b=1 to 3;D

9、o rep=1 to 4;Input y;Output;End;End;End;Cards;50 55 40 35 35 35 30 40 45 40 40 5050 45 50 45 55 60 50 50 55 45 65 5585 60 90 85 65 70 80 65 70 70 70 7060 55 35 70 60 85 45 75 65 65 85 75;Proc Sort;By a b rep;Run;Proc nested;Class a b;Var y;Run;Proc anova;/*均数的多重比较*/Class a b;Model y=a b(a);Test h=a

10、e=b(a);Means a/duncan e=b(a);Run;结果：例例 3 3 为了研究某种抗菌药的效果，对小白鼠进行试验。考虑 3 个试验因素，因素 A(用此抗菌药与否)可分为 A1(对照组不用抗菌药)、A2(试验组用抗菌药)；因素 B(小白鼠代次)可分为 B1(第 1 代)、B2(第 2代)、B3(第 3 代)；因素 C(性别)可分为 C1(雄性)、C2(雌性)。让第 1 代小白鼠被这种细菌感染，按雌雄分别统计对照组和试验组小白鼠的存活率，让存活的小白鼠分别在各自的组内(指对照组和试验组)生产第二代，对于第 2 代重复上面的试验，3 因素各水平搭配下都重复次试验，每次都有足够数量的小

11、白鼠(因为观测指标是存活率)，有人进行了如下的设计并收集到如下的试验数据。由专业知识得知：3 因素的主次顺序为 ABC。试分析 3 因素对小白鼠存活率的影响。DATA abc;DO r=1 to 2;DO a=1 to 2;DO b=1 to 3;DO c=1 to 2;INPUT p;y=ARSIN(SQRT(p/100);OUTPUT;END;END;END;END;CARDS;28 33 15 26 26 18 56 62 48 66 70 4919 27 11 22 20 11 51 60 44 61 65 55;PROC SORT;BY a b c r;RUN;PROC NESTED;CLASS a b c;VAR y;RUN;PROC GLM;CLASS a b c;MODEL y=a b(a)c(a b);TEST H=a E=b(a);TEST H=b(a)E=c(a b);Means a/duncan E=b(a);RUN;例例 4 4 鸡的试验 试验统计学P134区靖祥