(完整word版)STATA第四章t检验和单因素方差分析命令输出结果说明
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第四章 t检验和单因素方差分析命令与输出结果说明
·单因素方差分析
单因素方差分析又称为Oneway ANOVA,用于比较多组样本的均数是否相同,并假定:每组的数据服从正态分布,具有相同的方差,且相互独立,则无效假设。
:各组总体均数相同。
原假设:H
在STATA中可用命令:
oneway 观察变量分组变量[, means bonferroni]
其中子命令bonferroni是用于多组样本均数的两两比较检验。
例:测定健康男子各年龄组的淋巴细胞转化率(%),结果见表,问:各组的淋巴细胞转化率的均数之间的差别有无显著性?
健康男子各年龄组淋巴细胞转化率(%)的测定结果:
11-20 岁组:58 61 61 62 63 68 70 70 74 78
41-50 岁组:54 57 57 58 60 60 63 64 66
61-75 岁组:43 52 55 56 60
用变量x 表示这些淋巴细胞转化率以及用分组变量group=1,2,3分别表示
则用 STATA 命令:
oneway x group, mean bonferroni
| Summary of x
group | Mean ①
-------------+------------
1 | 66.5
2 | 59.888889
3 | 53.2
------+------------
Total | 61.25 ②
Analysis of Variance
Source SS df MS F Prob > F
------------------------------------------------------------------------------- Between groups 616.311111③ 2 ④ 308.155556⑤ 9.77⑥ 0.0010⑦Within groups 662.188889⑧ 21⑨ 31.5328042⑴
------------------------------------------------------------------------------- Total 1278.50 23 55.586956
(2)Bartlett's test for equal variances:chi2(2) = 2.1977 (3)Prob>chi2=0.333
Comparison of x by group
(Bonferroni)
Row Mean- |
Col Mean | 1 2
-------------- --|--------------------------------------
2 | -6.61111 (4)
| 0.054 (5)
|
3 | -13.3 (6) -6.68889(8)
| 0.001 (7) 0.134 (9)
①对应三个年龄组的淋巴细胞转化率的均数;②三组合并在一起的总的样本
均数;③组间离均差平方和;④组间离均差平方和的自由度;⑤组间均方和(即:
⑤=③/④);⑧组内离均差平方和;⑨组内离均差平方和的自由度;(1)组内均
方和(即:(1)=⑧/⑨);⑥为F 统计值(即为⑤/(1));⑦为相应的p值;(2)
为方差齐性的Bartlett检验;(3)方差齐性检验相应的p值;(4)第二组的淋
巴细胞转化率样本均数—第一组的淋巴细胞转化率的样本均数的差;(5)第二和
第一组均数差的显著性检验所对应p 值;(6)第三组的淋巴细胞转化率样本均数—第一组的淋巴细胞转化率的样本均数的差;(7)第三和第一组均数差的显著
性检验所对应的 p 值;(8)第三组的淋巴细胞转化率样本均数—第二组的淋巴
细胞转化率的样本均数的差;(9)第三和第二组均数差的显著性检验所对应的p 值。
由上述结果可知:三组方差无显著地齐性,因此若三组数据近似服从正态
分布,无效假设Ho检验所对应的p值<0.01,可以认为这三组均数有显著差异。
由 Bonferroni统计检验结果表明:第一组淋巴细胞转化率显著地高于第三组淋
巴细胞转化率(p<0.005),其它各组之间均数无显著性差异。
附1: STATA运行结果
. oneway x group, mean bonferroni
Summary of
x
group Mean
1 66.5
2 59.888889
3 53.2
Total 61.25
Analysis of Variance
Source SS df MS F Prob > F
Between groups 616.311111 2 308.155556 9.77 0.0010
Within groups 662.188889 21 31.5328042
Total 1278.5 23 55.5869565
Bartlett's test for equal variances: chi2(2) = 2.1977 Prob>chi2 = 0.333 Comparison of x by group
(Bonferroni)
Row Mean-
Col Mean 1 2
2 -6.61111
0.054
3 -13.3 -6.68889
0.001 0.134
附录2:如何输入分组数据
将数据导入EXCEL 如图
再将文件转化为CSV文件,导入STA TA即可,如果不懂请下载第一章。