两独立样本定量资料的统计分析的Stata实现.DOC

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.0417.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-4 例 7-1 中要检验两种

3、饮水小鼠的肝脏锌平均含量有无差别，需要先检验实验组和对照组所来自的总体方差是否相等，即对例 7-1 中资料进行方差齐性检验。例 7-1 的 Stata 数据格式如下：group x1 实验组 7.142 实验组 5.95实验组 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.863 实验组 7.14 实验组 8.265 实验组 10.086 实验组 7.917 实验组 9.078 实验组 9.39 实验组 8.6410 实验组 8.5111

4、对照组 6.6112 对照组 7.3113 对照组 7.214 对照组 6.5915 对照组 7.6516 对照组 5.5917 对照组 6.3918 对照组 6.5719 对照组 7.9120 对照组 7.86正态性检验的程序和结果如下： 0.5swilk x if group=“实验组“ 对实验组资料进行正态性检验Shapiro-Wilk W test for normal dataVariable | Obs W V z Probz-+-x | 10 0.97984 0.311 -1.803 0.96432P=0.96432，按正态性假定进行统计分析。swilk x if group=“

5、对照组“ 对对照组资料进行正态性检验Shapiro-Wilk W test for normal dataVariable | Obs W V z Probz-+-x | 10 0.93502 1.001 0.002 0.49903P=0.49903，按正态性假定进行统计分析。方差齐性检验程序和结果如下： (0.1)sdtest x,by(group) 方差齐性检验Variance ratio test-Group | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+-对照组 | 10 6.968 .2352625 .7439653 6.435

6、799 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.375F(9,9) lower tail = F_L = F_obs = 0.375F(9,9) upper tail = F_U = 1/F_obs = 2.665Ha: sd(1) sd(2)P F_U = 0.1604 P F_obs = 0.9198

7、方差齐性检验的 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. Interval-+-对照组 | 10 6.968 .2352625 .7439653 6.435799 7.500201实验组 | 10 8.196 .3840521 1.214479

8、 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: 18Ho: mean(对照组) - mean(实验组) = diff = 0Ha: diff 0t = -2.7266 t = -2.7266 t = -2.7266P |t| = 0.0138 P t = 0.9931t 检验的 P=0.0138 0t = 1.7614 t = 1.7614 t = 1.7614

9、P |t| = 0.0943 P t = 0.0472t 检验的 P=0.09430.05，可以认为接触重金属工人和非接触工人的血胰岛素平均水平没有差别。例 7-7 为研究孕妇患有妊娠合并症对葡萄糖耐受水平有无影响，对 17 名孕妇作葡萄糖耐受水平试验，其中有 9 名孕妇曾患有妊娠合并症，而另外 8 名没有患妊娠合并症（表 7-4 中第（1） 、 （3）列） 。试比较患与未患妊娠合并症的孕妇葡萄糖平均耐受水平是否相同？表 7-4 患与未患妊娠合并症的孕妇葡萄糖耐受水平试验比较患妊娠合并症 未患妊娠合并症葡萄糖耐受水平（1）秩次（2）葡萄糖耐受水平（3）秩次（4）110 1 120 4.5119

10、 3 140 11133 9 162 14127 6 184 17141 12 132 8117 2 128 7135 10 177 15120 4.5 143 13181 16n1=8 T1=47.5 n2=9 T2=105.5例 7-7 的 Stata 数据如下：group x1 患妊娠合并症 1102 患妊娠合并症 1193 患妊娠合并症 1334 患妊娠合并症 1275 患妊娠合并症 1416 患妊娠合并症 1177 患妊娠合并症 1358 患妊娠合并症 1209 未患妊娠合并症 12010 未患妊娠合并症 14011 未患妊娠合并症 16212 未患妊娠合并症 18413 未患妊娠合

11、并症 13214 未患妊娠合并症 12815 未患妊娠合并症 17716 未患妊娠合并症 14317 未患妊娠合并症 181Stata 程序和结果如下：ranksum x, by(g)Two-sample Wilcoxon rank-sum (Mann-Whitney) testg | obs rank sum expected-+-1 | 8 47.5 722 | 9 105.5 81-+-combined | 17 153 153unadjusted variance 108.00adjustment for ties -0.13-adjusted variance 107.87Ho: x

12、(g=1) = x(g=2)z = -2.359Prob |z| = 0.0183t 检验的 P=0.0183 chi2 = 0.0000Log likelihood = -7.612088 Pseudo R2 = 0.8264-x | Coef. Std. 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.0000Log 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.642763W | (exposure)-Poisson 回归分组的 P0.001，可以认为甲镇原发性肝癌的死亡率比乙镇高。