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

第十六章析因设计和交叉设计资料分析的Stata实现
本章应用的Stata命令为:
方差分析
anova varname [term [/] [term [/] ...]]
重组交叉设计数据
pkshape id sequence period1 period2
分析交叉设计数据
pkcross e
例16-3 治疗白细胞降低病人16例,随机分为4组,给予相应治疗措施:基础药物(对照组)、基础药物+甲药、基础药物+乙药、基础药物+甲药+乙药,3个月后观察白细胞增加数(109/L),以评价甲药和乙药以及它们联合用药的作用,得到结果如表16-3。
表16-3 治疗白细胞降低四种不同疗法3个月后白细胞增加数
第一组(g=1)
(基础治疗)
第二组(g=2)
(基础治疗+甲药)
第三组(g3)
(基础治疗+乙药)
第四组(g=4)
(基础治疗+甲药+乙药)
















Stata数据如下:
a
b
g
x
0
0
1
1
0
0
1

0
0
1

0
0
1

1
0
2

1
0
2

1
0
2

1
0
2

0
1
3

0
1
3

0
1
3

0
1
3

1
1
4

1
1
4

1
1
4

1
1
4

Stata命令为:
anova x a b a*b
结果为:
Number of obs = 16 R-squared =
Root MSE = .181968 Adj R-squared =
Source | Partial SS df MS F Prob > F
-----------+----------------------------------------------------
Model | 3
|
a | 1
b | 1
a*b | .336400032 1 .336400032
|
Residual | .397350008 12 .033112501
-----------+----------------------------------------------------
Total | 15 .338544999
计算残差
predict e,residual
对残差进行正态性检验()
swilk e
Shapiro-Wilk W test for normal data
Variable |