1、 第七章两独立样本定量资料的统计分析 的 Stata 实现 例 7-1 为研究金属镉中毒对大鼠肝脏中锌含量的影响，随机抽取 20 只小鼠，随机分为实验组和对照组。实验组小鼠每日经饮水染毒，对照组正常饮水。1 个月后，测量小鼠肝脏中锌含量（ ug/ml）（表 7-1），试分析两种饮用水小鼠的肝脏中锌平均含量有无差异。 表 7-1 实验组和对照组小鼠肝脏中锌含量 （ ug/ml） 例 7-2 为研究接触某重金属对人体血胰岛素水平有无影响，研究者从接触某重金属的职业工人中随机抽取 14 人，从非接触工人中随机抽取 14 人，测量每个工人的血胰岛素水平（ ulu/ml）（表 7-2），试分析上述两个人
2、群的血胰岛素平均水平有无差异。 表 7-2 某重金属接触工人和非接触工人血胰岛素水平（ ulu/ml） 接触重金属工人 7.84 8.11 8.45 9.38 10.22 12.28 11.57 11.38 8.23 8.92 3.02 4.54 7.11 6.04 非接触工人 17.94 7.42 9.12 9.45 7.86 16.48 24.12 8.57 8.92 6.74 9.23 16.29 5.73 8.51 例 7-3 随机抽取 30 名在校大学生，男性 16 人，女性 14 人，检测血红蛋白含量（ g/ul），问男性和女性大学生的血红蛋白平均含量有无差别？（数据略） 例 7-
3、4 例 7-1 中要检验两种饮水小鼠的肝脏锌平均含量有无差别，需要先检验实验组和对照组所来自的总体方差是否相等，即对例 7-1 中资料进行方差齐性检验 。 例 7-1 的 Stata 数据格式如下： group x 1 实验组 7.14 2 实验组 5.95 3 实验组 7.1 4 实验组 8.26 实验组 7.14 5.95 7.10 8.26 10.08 7.91 9.07 9.30 8.64 8.51 对照组 6.61 7.31 7.20 6.59 7.65 5.59 6.39 6.57 7.91 7.86 5 实验组 10.08 6 实验组 7.91 7 实验组 9.07 8 实验组
4、9.3 9 实验组 8.64 10 实验组 8.51 11 对照组 6.61 12 对照组 7.31 13 对照组 7.2 14 对照组 6.59 15 对照组 7.65 16 对照组 5.59 17 对照组 6.39 18 对照组 6.57 19 对照组 7.91 20 对照组 7.86 正态性检验的程序 和结果 如下： 0.05 swilk x if group=“实验组 “ 对实验组资料进行正态性检验 Shapiro-Wilk W test for normal data Variable | Obs W V z Probz -+- x | 10 0.97984 0.311 -1.803
5、 0.96432 P=0.96432，按正态性假定进行统计分析。 swilk x if group=“对照组 “ 对对照组资料进行正态性检验 Shapiro-Wilk W test for normal data Variable | Obs W V z Probz -+- x | 10 0.93502 1.001 0.002 0.49903 P=0.49903，按正态性假定进行统计分析。 方差齐性检验程序和结果如下： ( 0.10) sdtest x,by(group) 方差齐性检验 Variance ratio test - Group | Obs Mean Std. Err. Std.
6、Dev. 95% Conf. Interval -+- 对照组 | 10 6.968 .2352625 .7439653 6.435799 7.500201 实验组 | 10 8.196 .3840521 1.214479 7.327214 9.064786 -+- combined | 20 7.582 .2605456 1.165195 7.036672 8.127328 - Ho: sd(对照组 ) = sd(实验组 ) F(9,9) observed = F_obs = 0.375 F(9,9) lower tail = F_L = F_obs = 0.375 F(9,9) upper
7、 tail = F_U = 1/F_obs = 2.665 Ha: sd(1) sd(2) P F_U = 0.1604 P F_obs = 0.9198 方差齐性检验的 P= 0.1604，可以按方差齐性的假定作进一步统计分析。 例 7-5 对例 7-1 资料，试比较两种饮水小鼠肝脏平均锌含量有无差别？ 相应的 Stata 程序为 ttest x,by(group) 两独立样本 t 检验 输出结果为 : Two-sample t test with equal variances - Group | Obs Mean Std. Err. Std. Dev. 95% Conf. Interva
8、l -+- 对照组 | 10 6.968 .2352625 .7439653 6.435799 7.500201 实验组 | 10 8.196 .3840521 1.214479 7.327214 9.064786 -+- combined | 20 7.582 .2605456 1.165195 7.036672 8.127328 -+- diff | -1.228 .4503826 -2.174219 -.2817813 - Degrees of freedom: 18 Ho: mean(对照组 ) - mean(实验组 ) = diff = 0 Ha: diff 0 t = -2.726
9、6 t = -2.7266 t = -2.7266 P |t| = 0.0138 P t = 0.9931 t 检验的 P=0.0138 0 t = 1.7614 t = 1.7614 t = 1.7614 P |t| = 0.0943 P t = 0.0472 t 检验的 P=0.09430.05，可以认为 接触重金属工人和非接触工人的血胰岛素平均水平没有差别 。 例 7-7 为研究孕妇患有妊娠合并症对葡萄糖耐受水平有无影响，对 17 名孕妇作葡萄糖耐受水平试验，其中有 9 名孕妇曾患有妊娠合并症，而另外 8 名没有患妊娠合并症（表 7-4 中第（ 1）、（ 3）列）。试比较患与未患妊娠合并
10、症的孕妇葡萄 糖平均耐受水平是否相同？ 表 7-4 患与未患妊娠合并症的孕妇葡萄糖耐受水平试验比较 患妊娠合并症 未患妊娠合并症 葡萄糖耐受水平 （ 1） 秩次 （ 2） 葡萄糖耐受水平 （ 3） 秩次 （ 4） 110 1 120 4.5 119 3 140 11 133 9 162 14 127 6 184 17 141 12 132 8 117 2 128 7 135 10 177 15 120 4.5 143 13 181 16 n1=8 T1=47.5 n2=9 T2=105.5 例 7-7 的 Stata 数据如下： group x 1 患妊娠合并症 110 2 患妊娠合并症 11
11、9 3 患妊娠合并症 133 4 患妊娠合并症 127 5 患妊娠合并症 141 6 患妊娠合并症 117 7 患妊娠合并症 135 8 患妊娠合并症 120 9 未患妊娠合并症 120 10 未患妊娠合并症 140 11 未患妊娠合并症 162 12 未患妊娠合并症 184 13 未患妊娠合并症 132 14 未患妊娠合并症 128 15 未患妊娠合并症 177 16 未患妊娠合并症 143 17 未患妊娠合并症 181 Stata 程序和结果如下： ranksum x, by(g) Two-sample Wilcoxon rank-sum (Mann-Whitney) test g | o
12、bs rank sum expected -+- 1 | 8 47.5 72 2 | 9 105.5 81 -+- combined | 17 153 153 unadjusted variance 108.00 adjustment for ties -0.13 - adjusted variance 107.87 Ho: x(g=1) = x(g=2) z = -2.359 Prob |z| = 0.0183 t 检验的 P=0.0183 chi2 = 0.0000 Log likelihood = -7.612088 Pseudo R2 = 0.8264 - x | Coef. Std.
13、 Err. z P|z| 95% Conf. Interval -+- group | -.6711683 .0810558 -8.28 0.000 -.8300348 -.5123018 _cons | 6.780416 .1150509 58.93 0.000 6.55492 7.005912 - Poisson 回归分组的 P chi2 = 0.0000 Log likelihood = -6.1989818 Pseudo R2 = 0.6759 - x | Coef. Std. Err. z P|z| 95% Conf. Interval -+- group | -.8102747 .1618094 -5.01 0.000 -1.127415 -.493134 _cons | 5.174376 .2389773 21.65 0.000 4.70599 5.642763 W | (exposure) - Poisson 回归分组的 P0.001，可以认为 甲镇原发性肝癌的死亡率比乙镇高。
