主处理按拉丁方、副处理按随机区组安排的裂区设计.docVIP

主处理按拉丁方，副处理按随机区组安排的裂区设计试验资料的分析 N P K N B K D N C P A K P K N P N K P N A P B N D K C P K P N K P K P N C N A P B N D K K N K P N P K N D P C K A N B P K N P K 图2.15给出了一个主处理按拉丁方，副区按随机区组安排的裂区设计的例子。如果在一个裂区设计的试验中，主(整)区因素(A)有a（i＝1, 2, …, a）个水平，按拉丁方安排了a（j＝1, 2, … , a）个横行区组、a（k＝1, 2, … , a）个直行区组。共有a2个主区；又，副(裂)区因素(B)有b（l＝1, 2, …, b）个水平，每个整区再被划分为b个副(裂)区，于是整个试验共有a2b个副(裂)区，即整个试验共有a 2b个观察值。此类试验资料中各观察值的数学模型为： （i＝1，2，…，a；j＝1，2，…，a；k＝1，2，…，a；l＝1，2，…，b） 其中为整区ij(k)的效应值，为第i横行区组的效应值，为第j直行区组的效应值，为第k个A水平的效应值，其下标加括号是为了表示这个下标与横行、直行的下标有重叠作用，为主区误差效应；为第l个B水平的效应值，为Ak与Bl之间的交互作用，为裂区误差。仿照主副处理都按随机区组排列的裂区设计资料的情况，将分析分为主区部分和裂区部分两步进行，方差分析表分别如表8.60、表8.61所示。 表8.60 第一步：将整区变异分解为横行变异、直行变异、主处理(A)间变异和主区误差 变异来源 自由度 平方和 均方 横行间 dfR＝a－1 SSR＝((各横行之和) 2/每横行中观察值数目－C － 直行间 dfC＝a－1 SSC＝((各直行之和) 2/每直行中观察值数目－C － A间 dfA＝a－1 SSA＝((各A水平之和) 2/每A水平中观察值数目－C MSA 误差A dfea＝(a－1) (a－2) SSea＝SSTA―SSR―SSC―SSA MSea 整区变异 dfTA＝a2－1 SSTA＝((各整区之和) 2/每整区中观察值数目－C － 表8.61 第二步：按模型将总变异进行分解 变异来源 自由度 平方和 均方 整区变异 dfTA＝a2－1 SSTA＝((各整区之和) 2/每整区中观察值数目－C － B间 dfB＝b－1 SSB＝((各B水平之和) 2/每B水平中观察值数目－C MSB AB互作 dfAB＝(a－1) (b－1) SSAB＝处理组合平方和－SSA－SSB MSAB 误差B dfeb＝dfT－dfTA－dfB－dfAB SSeb＝SST－SSTA－SSB－SSAB MSeb 总变异 dfT ＝a2b－1 SST ＝所有观察值的平方之和－C MST 其中处理平方和要从AB二向表来计算，并利用它计算出AB间交互作用的平方和。将这两步合并便可以得到总的方差分析表。各项的期望均方与表8.46相同。 整个分析过程与 相同，不同之处只在于将第一步中的区组变异一项变成了横行变异和直行变异两项。这里就不拟再举例了。