1、第八章 多组独立定量资料的统计分析本章应用的 Stata 命令为：单因素方差分析 oneway 应变量 分组变量, bonferroni 非参数多组比较 kwallis 应变量 , by(分组变量)例 8-1 为研究茶多酚保健饮料对急性缺氧的影响，某研究者将 60 只小白鼠随机分为低、中、高三个剂量组和一个对照组，每组 15 只小白鼠。对照组给予蒸馏水 0.25ml 灌胃，低、中、高剂量组分别给予2.0gkg、4.0gkg、8.0gkg 的饮料溶于 0.20.3ml 蒸馏水后灌胃。每天一次，40 天后，对小白鼠进行耐缺氧存活时间试验，结果如表 8-1。试比较不同剂量的茶多酚保健饮料对延长小白鼠
2、的平均耐缺氧存活时间有无差别。表 8-1 各组小白鼠耐缺氧时间组别 耐缺氧时间 （min）ijY20.79 22.91 27.21 19.34 17.85 23.79 22.60 18.53对照组（ ）1i23.23 20.14 26.71 19.36 17.22 24.13 15.8522.22 24.74 21.53 19.66 25.89 29.10 18.93 18.64低剂量组（ ）226.39 25.49 20.43 22.69 29.67 20.36 22.7428.56 28.67 25.28 30.38 23.13 23.47 28.88 29.62中剂量组（ ）3i24.8
3、2 34.64 22.29 29.22 25.63 35.12 32.3231.93 37.94 39.76 27.94 29.65 34.23 32.63 29.13高剂量组（ ）439.62 36.15 28.85 24.07 29.29 35.24 36.13Stata 数据格式为：g x1 20.791 22.911 27.211 19.341 17.851 23.791 22.61 18.531 23.231 20.141 26.711 19.361 17.221 24.131 15.852 22.222 24.742 21.532 19.662 25.892 29.12 18.93
4、2 18.642 26.392 25.492 20.432 22.692 29.672 20.362 22.743 28.563 28.673 25.283 30.383 23.133 23.473 28.883 29.623 24.823 34.643 22.293 29.223 25.633 35.123 32.324 31.934 37.944 39.764 27.944 29.654 34.234 32.634 29.134 39.624 36.154 28.854 24.074 29.294 35.244 36.13：四个总体均数相等0H：四个总体均数不等或不全相等105.Stata
5、 命令为：oneway x g结果为：Analysis of VarianceSource SS df MS F Prob F-Between groups 1205.79582 3 401.931941 26.09 0.0000Within groups 862.796912 56 15.4070877-Total 2068.59273 59 35.0608938Bartletts test for equal variances: chi2(3) = 1.7256 Probchi2 = 0.631检验统计量 F26.09，p F-Between groups .173146778 2 .0
6、86573389 5.90 0.0075Within groups .396399936 27 .014681479-Total .569546715 29 .019639542Bartletts test for equal variances: chi2(2) = 1.9926 Probchi2 = 0.369Comparison of x by g(Bonferroni)Row Mean-|Col Mean | 1 2-+-2 | .098| 0.245|3 | .186 .088| 0.006 0.3483 个组比较 F=5.9，p0.0075 ，差别有统计意义，然后用 Bonferr
7、oni 法进行两两比较，其中 1 比 2 的 p0.245，1 比 3 为 0.006，2 比 3 为 0.348。例 8-4 对按完全随机设计分组的四组大白鼠，给予不同剂量的某种激素后，测量耻骨间隙宽度的增加量（） ，结果如表 8-9 所示，试分析给予不同剂量的某种激素后大白鼠耻骨间隙宽度的平均增加量有无差异？表 8-9 耻骨间隙宽度的增加量（）1 组 0.15 0.30 0.40 0.502 组 1.20 1.35 1.40 1.50 1.90 2.303 组 0.50 1.20 1.40 2.00 2.20 2.204 组 1.50 1.50 2.50 2.50对于表8-9数据，如果对四
8、组数据进行Levene方差齐性检验，P0.025，按 水准，表明该资料不满足方差分析要求方差齐性的基本条件，所以0.1这里采用Kruskal-Wallis检验。1） 建立检验假设：接受不同剂量（4 种）激素的大白鼠耻骨间隙宽度的增加量总体分布0H相同：接受不同剂量（4 种）激素的大白鼠耻骨间隙宽度的增加量总体1H分布不全相同=0.05Stata 数据格式为g x1 1 0.152 1 0.33 1 0.44 1 0.55 2 1.26 2 1.357 2 1.48 2 1.59 2 1.910 2 2.311 3 0.512 3 1.213 3 1.414 3 215 3 2.216 3 2.
9、217 4 1.518 4 1.519 4 2.520 4 2.5Stata 命令为：kwallis x,by( g)结果为：Test: Equality of populations (Kruskal-Wallis test)+-+| g | Obs | Rank Sum |-+-+-| 1 | 4 | 10.50 | 2 | 6 | 68.00 | 3 | 6 | 68.50 | 4 | 4 | 63.00 |+-+chi-squared = 10.501 with 3 d.f.probability = 0.0148chi-squared with ties = 10.572 with 3 d.f.probability = 0.0143按 =0.05 水准，拒绝 ，接受 ，故可认为接受不同剂量（4 种）激素0H1的大白鼠平均耻骨间隙宽度的增加量总体分布不全相同。
