选择的钙磷含量水平来看，有一种随着饲料中钙含量的减少，要求磷含量增加的趋势。当磷含量0.8%时，钙以0.4%为好，但除显著高于钙为1.0%的水平外，与钙为0.6%、0.8%的差异不显著；当磷的水平为0.6%时，钙的水平也以0.6%为好，但除显著高于钙为1.0%的水平外，与钙为0.4%、0.8%的差异不显著；磷含量0.4%时，钙含量以0.8%为好；磷含量为0.2%时，钙水平达到1.0%效果较好，但与钙为0.8%的差异不显著。同样也呈现一种随着磷含量降低，钙水平应提高的趋势。
综观全试验，以A2B3(钙0.8%，磷0.4%)效果最好，钙磷含量均高或均低效果都差。
二、系统分组资料的方差分析
在生物科学的研究中，实际问题是多种多样的，有些涉及多因素问题的研究或试验用交叉分组是困难的。例如，要比较a头公畜的种用价值，就必须考虑到与配的母畜。这是因为公畜的种用价值是通过后代的表现来评定的，而后代的表现除受公畜的影响外还要受到母畜的影响。但是在同期，公畜和母畜这两个因素的不同水平(不同公畜和不同母畜)是不能交叉的，即同一头母畜不能同时与不同的公畜交配产生后代。合理的方法是，选择一些生产性能大体一致的同胎次母畜随机分配与a头公畜交配，即公畜A1与一组母畜交配，公畜A2与另一组母畜交配??。然后通过后代的性能表现来判断这些公畜的种用价值有无显著差异。又如，为了比较利用同一设备生产同一种饲料的不同班组产品质量有无差异，我们可从每班组所生产的饲料中随机抽取若干样品，每个样品作若干次测定，根据测定结果判断不同班组的产品质量有无差异。
在安排多因素试验方案时，将A因素分为a个水平，在A因素每个水平Ai下又将B因素分成b个水平，再在B因素每个水平Bij下将C因素分c个水平??，这样得到各因素水平组合的方式称为系统分组(hierarchical classification)或称多层分组、套设计、窝设计。
在系统分组中，首先划分水平的因素(上述的不同公畜、不同班组)叫一级因素(或一级样本)，其次划分水平的因素(如上述的母畜、抽取的样品)叫二级因素(二级样本，次级样本)，类此有三级因素??。在系统分组中，次级因素的各水平会套在一级因素的每个水平下，它们之间是从属关系而不是平等关系，分析侧重于一级因素。
由系统分组方式安排的多因素试验而得到的资料称为系统分组资料。根据次级样本含量是否相等，系统分组资料分为次级样本含量相等与不等两种。最简单的系统分组资料是二因素系统分组资料。
如果A因素有a个水平；A因素每个水平Ai下，B因素分b个水平；B因素每个水平Bij
下有n个观测值，则共有abn个观测值，其数据模式如表6-35所示。
表6-35 二因素系统分组资料数据模式
(i?1,2,?,a;j?1,2,?,b;j?1,,2,?,n)
一级因素
A
A1
二级因素
B B11 B12 ?
x111 x121 ?
观测值C xijl x112 x122 ?
二级因素 总和xij.
平均x
ij.一级因素 总和xi.. x1..
平均x
i..? ? ? x11n x12n ? x11. x12. ?
x11. x12. ?x1..
107
??B1b
A2
B21 B22 ? B2b ? Ba1 Ba2 ? Bab
x1b1 x211 x221 ? x2b1
?x1b2 x212 x222 ? x2b2
??
Aa 合计
xa11 xa21 ? xab1 aa12 xa22 ? xab2
? ? ? ? ? ? ? ? ? ?
x1bn x21n x22n ? x2bn
?x1b. x21. x22. ? x2b.
?x1b. x21.
x22.
?
x2.. xa.. x…
x2..
x2b. ?
xa.. x...xa1n xa2n ? xabn xa1. xa2. ? xab.
xa1. xa2. ?
xab.
表6-35中，
xij.?xi..??x xijll?1bnnij.?xij./n
??j?1al?1bxijln xi..?xi../bn?x.../abn x数学模型为
...????x xijli?1j?1l?1...xijl????i??ij??ijl(i?1,2,?,a;j?1,2,?,b;j?1,,2,?,n)iij (6-36)
??i??,?ij??ij??i;?i、
2
??为Ai内Bij的效应，式中μ为总体平均数，?为Ai的效应，?ijijli分别为Ai、Bij观测值总体平均数。?为随机误差，相互独立，且都服从N(0,σ)。 系统分组资料的数学模型与交叉分组的不同，其中不包含交互作用项；并且因素B的效
应?是随着A的水平的变化而变化的，这就是说次级因素的同一水平在一级因素不同水平
ij中有不同的效应。因此，须把一级因素不同水平中的次级因素同一水平看作是不同水平。至于??，??是否一定为零，应视?,?是固定还是随机而定。
iijiij表6-35数据的总变异可分解为A因素各水平(Ai)间的变异(一级样本间的变异)，A因素各水平(Ai)内B因素各水平(Bij)间的变异(一级样本内二级样本间的变异)和试验误差(B因素各水平内观测值间的变异)。对两因素系统分组资料进行方差分析，平方和与自由度的剖分式为：
SST=SSA+SSB(A)+SSe (6-37) dfT=dfA+dfB(A)+dfe
各项平方和与自由度计算公式如下：
108
总平方和及其自由度C?x.../abnabn2(xijlabn2ijl2SST????i?1j?1l?1n?x...)2????xi?1j?1l?1?CdfT?abn?1一级因素间平方和及其SSdfA自由度a2?bn?(xi?1ai..?x...)?1bn?xi?12i..?CA?a?1方和及其自由度aij.b2xijlab (6-38)
1n一级因素内二级因素平SSdfB(A)?n??(xi?1j?1a?xi..)2???i?1j?1?1bn?xi?1n2i..B(A)?a(b?1))平方和及其自由度nijlaba2xijlb2ij.b误差(二级因素内三级因素SSe?SSC(B)????(xi?1j?1l?1?xij.)2????i?1j?1l?1?1n??xi?1j?1dfe?dfC(B)?ab(n?1)各项均方如下：
一级因素的均方
MSA?SS?SS?SSA/dfA B(A)一级因素内二级因素的均方 MS误差(二级因素内三级因素)均方 MSF检验时F值的计算： 当检验一级因素时，用MSB(A)B(A)B(A)/df C(B)e/dfe 作分母，即: F?MSeA/MSB(A) /MSe当检验一级因素内二级因素时，用MS作分母，即: F由数学模型决定的)，有关这方面的内容将在第四节介绍。
?MSB(A) 实际上，计算F值时分母项的选择是由有关因素的效应是固定还是随机所决定的(即是
(一)次级样本含量相等的系统分组资料的方差分析
【例6.7】为测定3种不同来源的鱼粉的蛋白质消化率，在不含蛋白质的饲料里按一定比例分别加入不同的鱼粉A1,A2,A3，配制成饲料，各喂给3头试验动物(B)。收集排泄物、风干、粉碎、混和均匀。分别从每头动物的排泄物中各取两份样品作化学分析。测定结果(xijl)列于表6-36，试分析不同来源鱼粉的蛋白质消化率是否有显著差异。
表6-36 蛋白质的消化率 鱼粉A A1
个体B
B11 B12 B13 B21 B22 B23 B31 B32 B33
测定结果C(xijl) 82.5 82.4 87.1 86.5 84.0 83.9 86.6 86.2 87.0 82.0 80.0 79.5
85.8 85.7 87.6 81.5 80.5 80.3
xij.
xij.
xi..
xi..
164.9 173.6 167.9 172.4 171.9 174.6 163.5 160.5 159.8
82.5 86.8 84.0 86.2 86.0 87.3 81.8 80.3 79.9
506.4 84.4
A2 518.9 86.5
A3 483.8 80.6
109
Σ
x…=1509.1
这是一个二因素系统分组资料，A因素的水平数a=3，Ai内B因素的水平数b=3，Bij
内重复测定次数n=2，共有abn=3×3×2=18个观测值，方差分析如下。
1、计算各项平方和与自由度
矫正数 C总平方和及其自由度
SST??x.../abn?1509.1/18?126521.267222
???xijl?C?(82.522?82.42???79.52?80.3)?126521.26722?126653.6100?126521.2672?132.3428dfT?abn?1?3?3?2?1?17
鱼粉间平方和及其自由度
SSA?1bn?xi..?C?213?2(506.4?518.9?483.8)?126521.2672222?126626.7683?126521.2672?105.5011dfA?a?1?3?1?2
鱼粉内个体间的平方和及其自由度
SS?dfB(A)?1n2??xij.?1bn22?xi..?122(164.9?173.6???160.5?159.8)222213?2B(A)(506.4?518.9?483.8)?126652.2250?126626.7683?25.45672
?a(b?1)?3?(3?1)?6误差（个体内分析样品间)平方和及其自由度
SSe?SSC(B)????xijl?21??nxij.?126653.6100?126652.2250?1.38502
dfe?dfC(B)?ab(n?1)?3?3?(2?1)?92、列出方差分析表，进行F检验
表6-37 不同来源鱼粉蛋白质消化率方差分析表
变异来源
鱼粉间A
鱼粉内个体间B(A)
误差C(B) 总变异
平方和 105.5011 25.4567 1.3850 132.3428 自由度 2 6 9 17 均方 52.7506 4.2428 0.1539 F值 12.43** 27.57**
查临界F值：F0.01(2,6)=10.92，F0.01(6,9)=5.80，因为鱼粉间的F＞F0.01(2,6)，鱼粉内个体间的F＞F0.01(6,9)，表明不同来源的鱼粉蛋白质消化率差异极显著，即3种鱼粉的质量差异极显著；喂同一鱼粉的不同个体对鱼粉的消化利用能力差异也极显著。
3.三种鱼粉平均消化率的多重比较（SSR法） 因为对一级因素（鱼粉）进行F检验时是以鱼粉内个体间均方作为分母，鱼粉的重复数为bn，所以鱼粉的标准误为：
Sx?MSB(A)/bn?4.2428?0.8409
x以dfB(A)=6，查附表6得k=2,3时SSR0.05和SSR0.01的值与S相乘求出相应的LSR0.05和LSR0.01的值，得:
110