STATA 第四章 t检验和单因素方差分析命令输出结果说明

1第四章t检验和单因素方差分析命令与输出结果说明·单因素方差分析单因素方差分析又称为OnewayANOVA,用于比较多组样本的均数是否相同，并假定：每组的数据服从正态分布,具有相同的方差，且相互独立，则无效假设。原假设：H0:各组总体均数相同。在STATA中可用命令：oneway观察变量分组变量[,meansbonferroni]其中子命令bonferroni是用于多组样本均数的两两比较检验。例：测定健康男子各年龄组的淋巴细胞转化率(%)，结果见表，问:各组的淋巴细胞转化率的均数之间的差别有无显著性？健康男子各年龄组淋巴细胞转化率(%)的测定结果:11-20岁组：5861616263687070747841-50岁组：54575758606063646661-75岁组：4352555660用变量x表示这些淋巴细胞转化率以及用分组变量group=1，2，3分别表示11-20岁组，41-50岁组和61-75岁组,即：数据表示为：x586161626368707074785457group111111111122x575860606364664352555660group222222233333则用STATA命令：onewayxgroup,meanbonferroni|Summaryofxgroup|Mean①-------------+------------1|66.52|59.8888893|53.2------+------------Total|61.25②2AnalysisofVarianceSourceSSdfMSFProbF-------------------------------------------------------------------------------Betweengroups616.311111③2④308.155556⑤9.77⑥0.0010⑦Withingroups662.188889⑧21⑨31.5328042⑴-------------------------------------------------------------------------------Total1278.502355.586956（2）Bartlett'stestforequalvariances:chi2(2)=2.1977（3）Probchi2=0.333Comparisonofxbygroup(Bonferroni)RowMean-|ColMean|12----------------|--------------------------------------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统计检验结果表明：第一组淋巴细胞转化率显著地高于第三组淋巴细胞转化率(p0.005)，其它各组之间均数无显著性差异。3附1：STATA运行结果附录2：如何输入分组数据将数据导入EXCEL如图再将文件转化为CSV文件，导入STATA即可，如果不懂请下载第一章。0.0010.1343-13.3-6.688890.0542-6.61111ColMean12RowMean-(Bonferroni)ComparisonofxbygroupBartlett'stestforequalvariances:chi2(2)=2.1977Probchi2=0.333Total1278.52355.5869565Withingroups662.1888892131.5328042Betweengroups616.3111112308.1555569.770.0010SourceSSdfMSFProbFAnalysisofVarianceTotal61.25353.2259.888889166.5groupMeanxSummaryof.onewayxgroup,meanbonferroni