第九章 配伍区组设计资料的统计分析的Stata实现.doc

第九章配伍区组设计资料的统计分析的Stata实现
配对t检验
ttest 变量1=变量2
配伍区组设计的方差分析
anova 应变量处理分组变量区组变量
方差分析的残差计算
在anova命令执行后,执行下列命令
predict 变量名,residual
配伍区组设计的Friedman检验
friedman 区组1 ……区组b
(见Stata7附加程序)
例9-1 某研究者用某药物治疗高血压患者10名,治疗前后舒张压的变化情况见表9-1。
表9-1 10名患者用某药物治疗后的舒张压测定值(mmHg)
患者编号
(1)
治疗前
(2)
治疗后
(3)
差值
(4)=(3)-(2)
1
2
3
4
5
6
7
8
9
10
115
110
129
109
110
116
116
116
120
104
116
90
108
89
92
90
110
120
88
96
-1
20
21
22
18
26
6
- 4
32
8
Stata数据为:

x1
x2
1
115
116
2
110
90
3
129
108
4
109
89
5
110
92
6
116
90
7
116
110
8
116
120
9
120
88
10
104
96
Stata命令为:
ttest x1= x2
结果为:
Paired t test
------------------------------------------------------------------------------
Variable | Obs Mean Std. Err. Std. Dev. [95% Conf. Interval]
---------+--------------------------------------------------------------------
x1 | 10
x2 | 10
---------+--------------------------------------------------------------------
diff | 10
------------------------------------------------------------------------------
Ho: mean(x1 - x2) = mean(diff) = 0
Ha: mean(diff) < 0 Ha: mean(diff) != 0 Ha: mean(diff) > 0
t = t = t =
P < t = P > |t| = P > t =
P=,治疗前后舒张压有差别,治疗后下降。
例9-2 某研究者对8名冻疮患者足部的两个冻疮部位(两个部位冻疮程度非常接近)用两种不同药物治疗,分别观测两个冻疮部位的痊愈时间,结果见表9-2。
表9-2 两种方法测定患者冻疮痊愈时间时间(天)结果
受试者编号
(1)
药物1
(2)
药物2
(3)
差值
(4)=(3)-(2)
1
8
12
4
2
10
9
-1
3
6
9
3
4
4
12
8
5
7
6
-1
6
10
10
0
7
8
11
3
8
10
11
1
Stata数据为:

x1
x2
1
8
12
2
10
9
3
6
9
4
4
12
5
7
6
6
10
10
7
8
11
8
10
11
建立检验假设,确定检验水准
,两种药物治疗的冻疮痊愈平均时间相同
,两种药物治疗的冻疮痊愈平均时间不同

Stata命令为:
ttest x1= x2
结果为:
Paired t test
------------------------------------------------------------------------------
Variable | Obs Mean Std. Err. Std. Dev.