第十六章 析因设计和交叉设计资料分析.docVIP

第十六章 析因设计和交叉设计资料分析的Stata实现 本章应用的Stata命令为： 方差分析 anova varname [term [/] [term [/] ...]] 重组交叉设计数据 pkshape id sequence period1 period2 分析交叉设计数据 pkcross outcome 例16-3 治疗白细胞降低病人16例，随机分为4组，给予相应治疗措施：基础药物（对照组）、基础药物+甲药、基础药物+乙药、基础药物+甲药+乙药，3个月后观察白细胞增加数（109/L） 表16-3 治疗白细胞降低四种不同疗法3个月后白细胞增加数 第一组(g=1) （基础治疗） 第二组(g=2) (基础治疗+甲药) 第三组(g3) (基础治疗+乙药) 第四组(g=4) (基础治疗+甲药+乙药) 1.00 0.60 0.88 0.70 1.48 1.42 1.11 1.30 1.42 1.28 0.90 0.98 2.32 2.38 2.12 2.21 Stata数据如下： a b g x 0 0 1 1 0 0 1 0.6 0 0 1 0.88 0 0 1 0.7 1 0 2 1.48 1 0 2 1.42 1 0 2 1.11 1 0 2 1.3 0 1 3 1.42 0 1 3 1.28 0 1 3 0.9 0 1 3 0.98 1 1 4 2.32 1 1 4 2.38 1 1 4 2.12 1 1 4 2.21 Stata命令为： anova x a b a*b 结果为： Number of obs = 16 R-squared = 0.9218 Root MSE = .181968 Adj R-squared = 0.9022 Source | Partial SS df MS F Prob F -----------+---------------------------------------------------- Model | 4 3 1 47.12 0.0000 | a | 2 1 2 81.72 0.0000 b | 1 1 1 49.48 0.0000 a*b | .336400032 1 .336400032 10.16 0.0078 | Residual | .397350008 12 .033112501 -----------+---------------------------------------------------- Total | 5 15 .338544999 计算残差 predict e,residual 对残差进行正态性检验() swilk e Shapiro-Wilk W test for normal data Variable | Obs W V z Probz -------------+------------------------------------------------- e | 16 0.94820 1.049 0.096 0.46179 P=0.46179(，按正态分布假定进行统计分析。 由于交付作用项a*b的P=0.0078，所以可以认为存在交互作用。 对例16-3作简单效应检验（Bonferroni校正）和方差齐性检验（）。 oneway x g , b t 结果为： |