《Content： 析因设计、正交设计、重复测量》由会员分享，可在线阅读，更多相关《Content： 析因设计、正交设计、重复测量（15页珍藏版）》请在金锄头文库上搜索。
1、*=* Title: 多样本均数比较的方差分析 2多因素试验资料的 SAS 方差分析程序 Content: 析因设计、正交设计、重复测量 Department: Health Statistics Teacher: 吴骋 *=* *=* 多因素试验：处理因素不只一个，研究目的不仅要比较处理因素之间的差别，还要分析不同处理因素 水平组合之间的差别。 *概念回顾： 1、单独效应：其它因素的水平固定时，同一因素不同水平间的差别。 2、主效应：某一因素各水平间的平均差别。 3、交互作用：当某因素的各个单独效应随另一因素变化而变化时，称这两个因素间存在交互作用。 *=* *=* 第五节 析因设计资料的方
2、差分析 析因设计(factorial design) 既考虑各因素的效应，又考虑到因素间的交互作用。实验中将各因素的 所有水平相互交叉进行组合，每种组合看作一种处理，然后在每种处理中进行实验。 析因设计资料的方差分析不仅要分析单个因素的作用，又要分析交互作用。 *=* *=* 一、两因素两水平的析因设计 例 6.5 将20只家兔随机等分 4 组，每组 5 只，进行神经损伤后的缝合试验。处理由两个因素组合而成，A因素为 缝合方法，有两水平，一为外膜缝合，另一为束膜缝合，B因素为缝合后的时间，有两水平，一为缝合后 1 月， 另一为缝合后 2 月? 试验结果为家兔神经缝合后的轴突通过率(%)（注：测
3、量指标，视为计量资料） ，见表 6.5，试做析因分析。 表 6.5 家兔神经缝合后的轴突通过率(%) - 缝合后时间(B) 缝合方法(A) - 外膜缝合 束膜缝合 1 月 10 10 10 20 40 30 50 50 10 30 2 月 30 50 30 50 70 70 60 60 30 30 - *=*; data prg6_5; do a=1 to 2; /*变量a表示 A 因素：缝合方法，1 代表外膜缝合，2 代表束膜缝合*/ do b=1 to 2; /*变量b 表示B 因素：缝合后时间，1 代表1个月，2 代表 2 个月*/ do i=1 to 5; /*变量 i 表示每次实验的
4、观测值次序*/ input x ; /*变量 x 表示家兔神经缝合后的轴突通过率(%)*/ output; end; end; end; cards; 10 10 40 50 10 30 30 70 60 30 10 20 30 50 30 50 50 70 60 30 ; proc anova;/*采用anova 过程*/ class a b;/*分类变量为两个因素：a缝合方法、b缝合后时间*/ model x=a b a*b; /*在anova过程中定义模型时， 用*表示两因素的交互作用， 本例考虑 a因素和 b 因素的交互作用*/ run; *=* 二、两因素三水平的析因分析 例 6.6
5、 观察 A,B 两种镇痛药物联合运用在产妇分娩时的镇痛效果。A 药取 3 个剂量：1mg，2.5mg，3mg； B 药也取 3 个剂量：5ug，15ug，30ug；共 9 个处理组。将 27 名产妇随机等分 9 组，记录分娩时的镇痛时间见表 6.6。 试分析 A、B两药联合运用的镇痛效果。 表 6.6 A、B 两药联合运用的镇痛时间（min） - A 药物剂量 B 药物剂量 - 5 15 30 - 105 115 75 1.0 mg 80 105 95 65 80 85 75 125 135 2.5 mg 115 130 120 80 90 150 85 65 180 5.0 mg 120 1
6、20 190 125 100 160 - *=*; data prg6_6; do a=1 to 3; /*变量a 表示A 药物剂量*/ do b=1 to 3; /*变量 B 表示B药物剂量*/ do i=1 to 3; /*变量 i 表示每次实验的观测值次序*/ input x ; /*变量x 表示镇痛时间*/ output; end; end; end; cards; 105 80 65 115 105 80 75 95 85 75 115 80 125 130 90 135 120 150 85 120 125 65 120 100 180 190 160 ; proc anova;
7、/*采用anova 过程*/ class a b; /*定义分组变量*/ model x=a b a*b; /*模 型考虑 a、b主效应，以及 a*b 交互作用*/ means a*b; /*计算a*b不同水平组合下镇痛时间 x 的描述性统计量（均数、标准差）*/ run; *=* 三、三因素不同水平的析因分析 例 6.7 用 522 析因设计研究 5种类型的军装在两种环境、两种活动状态下的散热效果，观察指标是受试者 的主观热感觉(从冷到热按等级评分)，结果见表 6.7。试进行方差分析。 表 6.7 不同条件下战士的主观热感觉 - 活动环境 活动状态 军 装 类 型 - 1 2 3 4 5 -
8、 干燥 静坐 0.25 0.30 0.75 0.20 -0.10 -0.25 0.10 0.50 -1.00 0.00 1.25 0.50 -0.60 0.85 2.50 -0.75 -0.35 0.40 -0.50 0.10 0.40 0.05 -0.20 0.90 -0.10 活动 4.75 4.60 4.55 4.25 4.725 3.45 4.80 3.50 3.10 4.30 4.00 4.00 4.25 4.00 4.10 4.85 5.20 4.10 5.00 4.80 4.55 4.30 4.40 4.20 3.60 潮湿 静坐 0.50 1.50 0.75 -0.75 1.75 2.10 1.50 2.65 0.90 2.40 2.75 1.25 3.00 0.95 1.75 1.00 1.37 0.05 0.625 3.05 2.35 2.55 1.17 1.05 2.75 活动 3.75 4.00 4.10 3.275 4.80 4.00 4.05 5.00 4.25 4.02 4.00 4.15 4.20 4.00 4.15 4.25 4.10 4.15 4.25 4.75 4.60 4.25 4.175 4.25 4.80 - *=*; data prg6_7; do a=1 to 5;/*变量 a: 军装类型*/ do b=1
