第十六章析因设计和交叉设计资料分析的SAS实现

第十六章 析因设计和交叉设计资料分析的 SAS 实现例16-3 治疗白细胞降低病人16例，随机分为4组，给予相应治疗措施：基础药物（对比组）、基础药物+甲药、基础药物+乙药、基础药物+甲药+乙药，3个月后观察白细胞增加数（109/L），以评价甲药和乙药以及它们联合用药的作用，得到结果如表16-3。表16-3 治疗白细胞降低四种不同疗法3个月后白细胞增加数第一组（基础治疗）第二组(基础治疗+甲药)第三组(基础治疗+乙药)第四组(基础治疗+甲药+乙药)1.000.600.880.701.481.421.111.301.421.280.900.982.322.382.122.21样本均数()0.801.331.152.26总体均数(μ)μ00μ10μ01μ11SAS 分析程序data ch16_3; 定义数据集名do a=1 to 2; 定义变量a，a=1表示不用甲药，a=2表示使用甲药do b=1 to 2; 定义变量b，b=1表示不用乙药，b=2表示用乙药do i=1 to 4;input x@@; output;end;end;end;cards;1.00 0.60 0.88 0.70 1.42 1.28 0.90 0.981.48 1.42 1.11 1.302.32 2.38 2.12 2.21;proc anova; 利用anova过程进行方差分析class a b; 指定分类变量model x=a b a*b; 指定模型，考虑a,b主效应以及交互作用run;SAS 软件输出结果The ANOVA ProcedureDependent Variable: x Sum of Source DF Squares Mean Square F Value Pr > F Model 3 4.68082500 1.56027500 47.12 <.0001 Error 12 0.39735000 0.03311250 Corrected Total 15 5.07817500 R-Square Coeff Var Root MSE x Mean 0.921753 13.17418 0.181968 1.381250 Source DF Anova SS Mean Square F Value Pr > F a 1 2.70602500 2.70602500 81.72 <.0001 b 1 1.63840000 1.63840000 49.48 <.0001 a*b 1 0.33640000 0.33640000 10.16 0.0078The ANOVA ProcedureSourceDFSum of squaresMean squareF valuePr > F变异来源自由度平方和均方和F 值P 值a(甲药)12.706025002.7060250081.72<.0001b(乙药)11.638400001.6384000049.48<.0001a*b(交互作用)10.336400000.3364000010.160.0078Error(误差)120.397350000.03311250SAS 软件输出结果解释甲药与乙药的交互作用的统计检验 P＝0.0078<0.05，因此有统计学意义，可以认为甲药与乙药对提高白细胞数有交互作用。教材中的说明甲药与乙药的交互作用的统计检验P＝0.008<0.05，因此有统计学意义，可以认为甲药与乙药对提高白细胞数有交互作用（推断...