嵌套裂区方差分析

嵌套裂区方差分析第一页，共二十五页，2022年，8月28日嵌套设计又被称为巢式设计（nesteddesign）、系统分组设计(hierarchalclassification)或组内又分亚组的设计。根据因素数的不同，嵌套设计可分为二级嵌套（二因素）、三级嵌套（三因素）等套设计。一、嵌套设计的方差分析第二页，共二十五页，2022年，8月28日

将全部k个因素按主次排列，依次称为1级，2级…

k级因素，再将总离差平和及自由度进行分解，其基本思想与一般方差分析相同。所不同的是分解法有明显的区别，它侧重于主要因素，并且，第i级因素的显著与否，是用第i级与第i+1级因素的均方分别做为分子和分母来构造F统计量，并以F检验为其理论根据的。嵌套设计的方差分析第三页，共二十五页，2022年，8月28日(1)情形一受试对象本身具有按其隶属关系进行分组再分组的各种因素。(2)情形二受试对象本身并非具有分组再分组的各种分组因素，而是各之间在专业上有主次之分。区分嵌套设计与析因设计的关键是看因素之间的地位是否平等，因素的地位平等则属于析因设计，不平等则属于嵌套设计。嵌套设计应用第四页，共二十五页，2022年，8月28日例1（嵌套关系）：选取某种植物3个品种(A)的植株，在每一株内选取2片叶子(B)，用取样器从每一片叶子上选取同样大小的两块（重复）进行检测。不能把B因素的2个水平简单地看作是与A因素3个水平的全面组合，而是分别嵌套在A1、A2、A3三个水平之下，相当于B因素有6个水平，但它们所产生的离差平方和中又包含了A因素的作用，一般来说，用它作为度量A因素作用大小的误差项，是严格考核A因素的一种措施。嵌套设计举例1第五页，共二十五页，2022年，8月28日例2（因素分主次）：为了研究某种抗菌药的效果，考虑3个因素对小白鼠进行试验。因素A可分为A1(对照组不用抗菌药)、A2(试验组用抗菌药)；因素B(小白鼠代次)可分为B1(第1代)、B2(第2代)、B3(第3代)；因素C(性别)可分为C1(雄性)、C2(雌性)。让第1代小白鼠被这种细菌感染，按雌雄分别统计对照组和试验组小白鼠的存活率，对于第2代、第2代重复上面的试验，……，观测小白鼠存活率。由专业知识得知：3因素的主次顺序为A→B→C。嵌套设计举例2第六页，共二十五页，2022年，8月28日二级嵌套的方差分析表嵌套设计分析表变异来源自由度平方和均方F值Prob>FAdf1=a-1SS1=Q1-CMS1=SS1/df1MS1/MS2

B(A)df2=a(b-1)SS2=Q2-Q1MS2=SS2/df2MS2/MSe

Edfe=(n-1)abSSe=W-Q2MSe=SSe/dfe

总和dfT=nab-1SSt=W-C

A因素有a个水平，B因素有b个水平，重复n次，总样本含量为nab，全部数据之总和为T，全部数据之平方和为W，校正数为C=T2/nab，Qi为系统分组分到第i级因素时，各小组数据之和的平方均值的和。

W=C=

Q1=Q2=第七页，共二十五页，2022年，8月28日注意：由于在嵌套设计中，两个因素之间不能自由交错地组成各种处理组合，因而，不能考察因素之间地交互作用。所以，凡是交互作用比较重要的试验，都不应采用这种设计。

嵌套设计注意事项第八页，共二十五页，2022年，8月28日进行嵌套设计的方差分析的过程：nested、anova或GLM，1、过程nested过程更简单一些，但它受数据中读取分组变量先后顺序的影响很大，应先用Sort过程按因素重要程度进行排序整理；2、过程anova或GLM不需要排序整理，且容易进行均值比较。平衡资料用过程nested、anova或GLM三种过程做结果是一样的；若不是平衡资料，须用过程GLM。

嵌套设计分析的sas实现第九页，共二十五页，2022年，8月28日例1的SAS程序DataNested;Inputplant$leafwt@@;Cards;a112.1a112.1b114.4b114.4c123.1c123.4a212.8a212.8b214.7b214.5c228.1c228.8;例1Sas程序第十页，共二十五页，2022年，8月28日过程步NestedProcSort; Byplantleaf;Run;ProcNested; Classplantleaf; Varwt;Run;例1Sas程序第十一页，共二十五页，2022年，8月28日过程步Anova/GLMProcAnova; Classplantleaf; Modelwt=plantleaf(plant); TestH=plantE=leaf(plant); Meansplant/DuncanE=leaf(plant);Run;例1Sas程序第十二页，共二十五页，2022年，8月28日裂区设计与两因素随机区组设计近似，但是两者不一样的。不同点之一是后者在每一区组内A、B两因素的ab次处理是完全随机化的。而裂区设计的每一区组内A因素先分为a个处理，在每一处理内B因素再分为b个处理。随机化过程只能在A因素的a个处理和B因素的b个处理之间进行。由A因素所划分的A个部分称为主区或整区，每一主区再划分的b个部分称为裂区或副区。二、裂区设计的方差分析ˉ第十三页，共二十五页，2022年，8月28日不同点之二是方差分析计算时F值时误差项的选择，裂区设计方差分析时有两个误差项，区组和整区是同一个误差，而裂区和交互作用则用另一个误差项，而二因素随机区组设计方差分析时用一个误差项。裂区设计的方差分析第十四页，共二十五页，2022年，8月28日进行裂区试验设计时首先要分清主要因子和次要因子，主要因子是想要获得较高精确度的因子，次要因子是精确度可以低些的因子。主要因子的各个水平随机安排在裂区，次要因子的各个水平随机安排在整区，只有这样，主要因子的各水平的重复数才会大大的多于次要因子的各个水平的重复数，才能获得较高的精确度。裂区设计的原则第十五页，共二十五页，2022年，8月28日适用范围：１．复因子试验中，两个因子要求的精确度不一时，可用裂区设计。２．各个因子的各个水平需要的面积大小不一时，亦可用裂区设计。３．在原有的试验的基础上，临时加入一个研究因子时，可用裂区设计。适用范围

第十六页，共二十五页，2022年，8月28日优点：１．某个因子可获得较高的精确度。２．能利用原有的试验地及试验材料，进行深一步的研究。缺点：１．资料的统计分析比较复杂，不易掌握。２．次要因子的精确度较低。优缺点

第十七页，共二十五页，2022年，8月28日以两个因素的裂区设计进行方差分析：设有A、B两个试验因素，A因素有a个水平，安排在整区，B因素有b个水平，安排在裂区，整个试验有n个重复区组。总变异分解为整区部分和裂区部分，整区部分总平方和（SS1）可分解为区组平方和（SSr）、A因素水平间平方和（SSA）和整区误差平方和（SSeA）；裂区部分分解为B因素水平间平方和（SSB）、交互作用AxB平方和（SSAB）和裂区误差平方和（SSeB）。平方和分解

第十八页，共二十五页，2022年，8月28日

平方和分解

变异来源自由度平方和均方F值Prob整区部分

Blockdfr=n-1SSrMSr=SSr/dfrMSr/MSeA

Adfa=a-1SSAMSA=SSA/dfaMSA/MSeA

EAdfea=(n-1)(a-1)SSeA=SSI-SSrSSAMSeA=SSeA/dfea

裂区部分

BdfB=b-1SSBMSB=SSB/dfBMSB/MSeB

AxBdfAB=(a-1)(b-1)SSAB=SSK-SSA-SSBMSAB=SSAB/dfABMSAB/MSeB

EBdfeB=a(b-1)(n-1)SSeB=SST-SSB-SSABMSeB=SSeB/dfeB

总和dfT=nab-1SST=W-C

注：全部数据之总和为T，全部数据之平和为W，校正数为C=T2/nab，Tj、Tm、Tl、Tml、Tjm分别为各区组、A各水平、B各水平、A和B各水平组合、区组n和A各水平组合的总和数。

SSI=SSr=SSA=SSB=SSk=第十九页，共二十五页，2022年，8月28日例题1

江苏某地在不整地和整地（A）的基础上，比较对照、猪牛粪、绿肥、堆肥和草塘泥等五种施肥（B1-B5），对早稻产量的影响，采用裂区设计，整区处理为不整地（A1）和整地（A2），裂区处理为施用不同农家肥，设4个重复区组，试作分析。

裂区设计举例B2B4B5B1B4B5B1B2B3B3B3B5B5B4B1B2B4B3B2B1B1B3B3B2B5B4B2B1B4B5B4B2B1B3B3B1B5B4B2B5区组IIIIIIIVA1A2A2A1A2A1A1A2第二十页，共二十五页，2022年，8月28日试验数据见下表：kg/667.7㎡

裂区设计举例裂区处理IIIIIIIVA1A2A1A2A1A2A1A2不施肥B1176445192445192448304524猪牛粪B2352592256504246520388500绿

肥B3416604325604406640486650堆

肥B4280548240485320584320524草塘泥B5405640444565366660456616第二十一页，共二十五页，2022年，8月28日Datasplit;Doblock=1to4;Doa=1to2;Dob=1to5;Inputyield@@;Output;End;End;End;裂区设计SAS分析－编程第二十二页，共二十五页，2022年，8月28日Cards;176352416280405445592604548640192256325240444445504604485565192246406320366448520640584660304388486320456524500650524616;

裂区设计SAS分析第二十三页，共二十五页，2022年，8月28日ProcAnova;Classblockab;Modelyield=blockabblock*aa*b;Testh=blockae=block*a;Meansa/Duncane=block*a;Meansb/Duncan;Run;

裂区设计SAS分析第二十四页，共二十五页，2022年，8月28日方差分析结果：SourceDFAnovaSS