巢式嵌套设计的方差分析.docVIP

巢式嵌套设计的方差分析 嵌套设计被称为巢式设计（nested design）或系统分组(hierarchal classification)的资料，有些教科书上称这类资料为组内又分亚组的分类资料。根据因素数的不同，套设计可分为二因素（二级）、三因素（三级）等套设计。 将全部k个因素按主次排列，依次称为1级，2级 … k级因素，再将总离差平和及自由度进行分解，其基本思想与一般方差分析相同。所不同的是分解法有明显的区别，它侧重于主要因素，并且，第i级因素的显著与否，是分别用第i级与第i+1级因素的均方为分子和分母来构造F统计量，并以F为其理论根据的。试验设计(1)情形一 受试对象本身具有按其隶属关系分组再分组的各种因素。3个品种（A），在每一株内选取两片叶子（B）（嵌套在植株因素下的第二个因素），用取样器从每一片叶子上选取同样面积的两个样本（两次重复），称取湿重。 不是把B因素的2个水平简单地看作是与A因素个水平的全面组合，而是分别嵌套在A1、A2水平之下，相当于B因素有个水平，但它们所产生的离差平和中又包含了A因素的作用，一般来说， 它大于模型所提供的总误差，用它作为度量A因素作用大小的误差项，是严格考核A因素的一种措施。 变异来源 自由度 平方和 均方 F值 ProbF A df1=a-1 SS1=Q1-C MS1= SS1/df1 MS1/ MS2 B（A） df2=a(b-1) SS2=Q2-Q1 MS2= SS2/df2 MS2/ MSe E dfe= (n-1)ab SSe=W-Q2 MSe= SSe/dfe 总和 nab-1 SSt=W-C 注：A因素有a个水平，B因素有b个水平，重复n次，总样本含量为，全部数据之总和为T，全部数据之平和为W，校正数为CQi为系统分组分到第i级因素时各小组数据之和的平均值的和 C= Q1= Q2= (2)情形二 受试对象本身并非具有分组再分组的各种分组因素，而是各之间在专业上有主次之分。区分设计与析因设计的关键是看因素之间的地位是否平等因素平等则属于析因设计，不平等则属于设计为了研究某种抗菌药的效果，对小白鼠进行试验。考虑3个试验因素， 因素A(用此抗菌药与否)可分为A1(对照组不用抗菌药)、A2(试验组用抗菌药)因素B(小白鼠代次)可分为B1(第1代)、B2(第2代)、B3(第3代)因素C(性别)可分为C1(雄性)、C2(雌性)。让第1代小白鼠被这种细菌感染，按雌雄分别统计对照组和试验组小白鼠的存活率，让存活的小白鼠分别在各自的组内(指对照组和试验组)，对于第2代重复上面的试验，……， 3因素各水平搭配下都重复２次试验，每次都有足够数量的小白鼠(因为观测指标是存活率)，有人进行了如下的设计并收集到如下的试验数据。由专业知识得知3因素的主次顺序为A→B→C。试分析3因素对小白鼠存活率的影响 变异来源 自由度 平方和 均方 F值 ProbF A df1=a-1 SS1=Q1-C MS1= SS1/df1 MS1/ MS2 B（A） df2=a(b-1) SS2=Q2-Q1 MS2= SS2/df2 MS2/ MS3 C（A B） df3=ab (c-1) SS3=Q3-Q2 MS3= SS3/df3 MS3/ MSe E dfe=(n-1)abc SSe=W-Q3 MSe= SSe/dfe 总和 nabc-1 SSt=W-C 注：A因素有a个水平，B因素有b个水平，C因素有c个水平，重复n次，总样本含量为，全部数据之总和为T，全部数据之平和为W，校正数为C=Qi为系统分组分到第i级因素时各小组数据之和的平均值的和但它受数据中读取分组变量先后顺序的影响很大，先用SORT过程整理整理平衡资料若不是平衡资料，必须用过程。 选取某种植物3个品种（Plant），在每一株内选取两片叶子（Leaf）（嵌套在植株因素下的第二个因素），用取样器从每一片叶子上选取同样面积的两个样本（两次重复），称取湿重。对以上结果进行方差分析。（固定效应Plant+随机效应Leaf） 解： SAS程序如下： Data Nested; Input plant $ leaf wt @@; Cards; a 1 12.1 a 1 12.1 b 1 14.4 b 1 14.4 c 1 23.1 c 1 23.4 a 2 12.8 a 2 12.8 b 2 14.7 b 2 14.5 c 2 28.1 c 2 28.8 ; Proc Anova; Class plant leaf ; Model wt = plant leaf(plant); Test H = plant E =