SAS第四讲.pptVIP

例：医学统计学 P253 例11-6 试验三种催化剂在不同温度下对某化合物的转化作用。由于各催化剂所要求的温度范围不同，将催化剂作为一级试验单位（A1、A2、A3），3种温度作为二级试验因素（B1、B2、B3），采用嵌套设计，每个处理重复2次，试验结果见下表4。 表4 化合物的转化率（%） 85 89 B1 (100) 71 67 B1 (90) A3 75 78 56 60 62 59 65 61 85 83 91 88 82 84 试验结果 X B1 (95) B3 (75) B2 (65) B1 (55) B3 (90) B2 (80) B1 (70) 温度（℃） A2 A1 催化剂 data nested; do a=1 to 3; do b=1 to 3; do i=1 to 2; input x @@; output; end; end; end; cards; 82 84 91 88 85 83 65 61 62 59 56 60 71 67 75 78 85 89 ; proc glm; class a b; model x=a b(a)/ss1; test h=a e=b(a); means a b(a); run; proc nested; class a b; var x; run; 裂区试验设计方差分析 裂区试验设计（split-plot design）也称为分割实验设计。裂区设计是析因设计的一种特殊形式，其安排的处理也是各因素各水平的全面组合，所不同的是裂区设计一般将试验单位分成两级，一级试验单位包含着二级试验单位，且一级处理因素作用于一级试验单位，二级处理因素作用于二级试验单位，而析因设计的全部处理因素作用于一级试验单位。 裂区设计的特点： 1. 因为二级试验单位包含在一级试验单位中，故可节省试验材料。 表5 完全随机分组3×2裂区试验设计（n=2) a3b2 a3b2 a3b1 a3b1 3 6 3水平 a2b2 a2b2 a2b1 a2b1 2 5 2水平 a1b2 a1b2 a1b1 a1b1 1 4 1水平 2水平 1水平 B因素 受试对象 A因素 （裂区组） 每个组合至少2个受试对象，共6个受试对象，A因素作用于一级试验单位，B因素作用于二级试验单位。 2. 一级处理与一级试验单位混杂，而二级处理与一级单位不混杂，即二级单位同质性好，二级单位误差小于一级单位误差，因此设计时可将最感兴趣、主要的因素、差异较小和要求精度较高的因素作为二级因素。 表6 完全随机分组3×2析因试验设计（n=2) a3b2 a3b2 a3b1 a3b1 a2b2 a2b2 a2b1 a2b1 a1b2 a1b2 a1b1 a1b1 b2因素 b1因素 b2因素 b1因素 b2因素 b1因素 a3因素 a2因素 a1因素 每个组合至少2个受试对象，共12个受试对象，全部为一级试验单位。 完全随机 I×J 裂区设计资料方差分析 1. 设计方案 （1） 将一级试验单位随机等分为I组(裂区组），每组例数为n(≥2),分别接受a1 , a2, ┄，aI 各水平的处理。 （2）分别将各一级试验单位内的二级试验单位随机分配接受b1，b2，┄，bJ的处理。（见表1）。 2. 方差分析 完全随机 I×J 裂区设计资料方差分析由两部分组成，表示 A因素主效应的FA以MS一级单位误差作为误差均方，表示B因素主效应和AB交互作用的FB、FAB以MS二级单位误差作为误差均方，（见书中表7-6）。 例7-4 试验一种全身注射抗毒素对皮肤损伤的保护作用，将10只家兔随机等分两组，一组注射抗毒素，一组注射生理盐水作对照。分组后，每只家兔取甲、乙两部位，分别注射低浓度毒素和高浓度毒素，观察指标为皮肤受损范围直径（mm），结果如表7-7所示。试作方差分析。 表7-7 家兔皮肤损伤直径（mm） 22.25 21.50 23.50 24.75 23.75 18.25 18.50 19.75 21.50 20.75 2 3 5 8 9 生理盐水（a2) 19.00 20.75 18.50 20.50 20.00 15.75 15.50 15.50 17.00 16.50 1 4 6 7 10 抗毒素（a1） 高浓度（b2） 低浓度（b1） 毒素浓度（B因素） 家兔编号 注射药物 （A因素） data ex16; do