裂区设计(析因设计特例)(医学统计学卫生统计学).pptVIP

裂区(分割)设计与统计分析 随机区组的裂区设计 例1 为研究3种小麦种子的产量，以及2种化肥对不同农作物产量的影响，在田地里进行了如下实验： 小麦种子（因素A）：A1,A2,A3 化肥（因素B）：B1,B2 表1 农作物种子和化肥对产量的影响 完全随机的裂区设计 例2 试验一种全身注射抗毒素对皮肤损伤的保护作用，将10只家兔随机等分两组，一组注射抗毒素，一组注射生理盐水作对照。之后，每只家兔取甲、乙两部位，随机分配分别注射低浓度毒素和高浓度毒素，观察指标为皮肤受损直径。 表2 家兔皮肤损伤直径（mm） 什么是裂区设计？ 例1：因素在施加上有先后顺序 因素A先施加（如小麦种子） 因素B后施加（如化肥） 例2：实验单位具有隶属关系（分级） 因素A作用于一级实验单位（如家兔） 因素B作用于二级实验单位（某部位） 按一级实验单位是否形成区组，可采用完全随机的裂区设计和随机区组的裂区设计。 什么是裂区设计？ 和析因设计不同 析因设计的处理因素是同时施加的； 析因设计的处理全部作用于同一级别的实验单位，各因素在专业上的地位是等价的。 方差分析的策略 检验一级单位和二级单位因素作用时检验的精度不同，检验二级单位因素作用时的检验精度高于一级单位。 分析一级单位的因素作用时，用和该因素有关的某种变异来源作为误差项，即用个体因素（区组因素）×一级单位因素的交互作用作为误差项。 完全随机裂区设计的方差分析 SAS程序 Proc glm; Class no a b; Model y=a no*a b a*b; Test h=a e=no*a; *检验a因素作用时，用 no*a作为误差项； Run; 随机区组裂区设计的方差分析 SAS Proc glm; Class s a b; Model y=a s s*a b a*b; Test h=a e=s*a; *检验a因素作用时，用 no*a作为误差项； Run; * * 1 3 2 4 A1 B1 B2 A2 A3 A3 A2 A1 4 A3 A2 A1 3 A3 A2 A1 2 87 98 A3 75 93 A2 95 99 A1 1 B2 B1 化肥B 种子 A 区组 10 9 8 7 6 5 4 3 2 1 家兔编号 no 23.75 20.75 24.75 21.50 23.50 19.75 21.50 18.50 22.25 18.25 生理盐水A2 20.00 16.50 20.50 17.00 18.50 15.50 20.75 15.50 19.00 15.75 抗毒素A1 高浓度B2 低浓度B1 毒素浓度（因素B） 注射药物 （因素A） 8 二级单位误差 1 A×B 1 因素B 10(19－9) 二级单位合计（个体内变异） 8 一级单位误差 1 因素A 9 一级单位合计（个体间变异） 19 二级单位总计（总变异） F MS 自由度 SS 变异来源 二级单位误差 A×B 因素B 二级单位合计 一级单位误差(S×A) 因素A 一级单位合计 区组S 二级单位总计 F MS 自由度 SS 变异来源