第七章两独立样本定量资料的统计分析的Stata实现Word文档格式.docx

1、7.86例7-2 为研究接触某重金属对人体血胰岛素水平有无影响，研究者从接触某重金属的职业工人中随机抽取14人，从非接触工人中随机抽取14人，测量每个工人的血胰岛素水平（ulu/ml）（表7-2），试分析上述两个人群的血胰岛素平均水平有无差异。表7-2 某重金属接触工人和非接触工人血胰岛素水平（ulu/ml）接触重金属工人7.848.118.459.3810.2212.2811.5711.388.238.923.024.547.116.04非接触工人17.947.429.129.4516.4824.128.576.749.2316.295.73 例7-3 随机抽取30名在校大学生，男性16人，

2、女性14人，检测血红蛋白含量（g/ul），问男性和女性大学生的血红蛋白平均含量有无差别？（数据略）例7-4 例7-1中要检验两种饮水小鼠的肝脏锌平均含量有无差别，需要先检验实验组和对照组所来自的总体方差是否相等，即对例7-1中资料进行方差齐性检验。例7-1的Stata数据格式如下：groupx1237.1456789.39101112137.214151617181920正态性检验的程序和结果如下：swilk x if group=实验组 对实验组资料进行正态性检验 Shapiro-Wilk W test for normal data Variable | Obs W V z Probz-+

3、- x | 10 0.97984 0.311 -1.803 0.96432P=0.96432 ，按正态性假定进行统计分析。对照组 对对照组资料进行正态性检验 Shapiro-Wilk W test for normal data x | 10 0.93502 1.001 0.002 0.49903P=0.49903方差齐性检验程序和结果如下：sdtest x,by（group） 方差齐性检验Variance ratio test- Group | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- 对照组 | 10 6.968 .2352625

4、 .7439653 6.435799 7.500201 实验组 | 10 8.196 .3840521 1.214479 7.327214 9.064786combined | 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 tail = F_U = 1/F_obs = 2.665 Ha: sd（1） sd（2） P F_obs = 0.08

5、02 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 diff | -1.228 .4503826 -2.174219 -.2817813Degrees of freedom: 18 mean（对照组） - mean（实验组） = diff = 0 diff 0 t = -2.72

6、66 t = -2.7266 t = -2.7266 t = 0.0069 P |t| = 0.0138 P t = 0.9931t检验的P=0.01380.05，可以认为饮用污染的水使的小鼠肝脏中锌平均含量高于饮用未污染水的小鼠。例7-6 对例7-2中的资料进行方差齐性检验，得F4.13，P |t| = 0.0943 P t = 0.0472t检验的P=0.09430.05，可以认为接触重金属工人和非接触工人的血胰岛素平均水平没有差别。例7-7 为研究孕妇患有妊娠合并症对葡萄糖耐受水平有无影响，对17名孕妇作葡萄糖耐受水平试验，其中有9名孕妇曾患有妊娠合并症，而另外8名没有患妊娠合并症（表7

7、-4中第（1）、（3）列）。试比较患与未患妊娠合并症的孕妇葡萄糖平均耐受水平是否相同？表7-4 患与未患妊娠合并症的孕妇葡萄糖耐受水平试验比较患妊娠合并症未患妊娠合并症葡萄糖耐受水平（1）秩次（2）（3）（4）1101204.5119140133162127184141132117128135177143181n1=8T1=47.5n2=9T2=105.5例7-7的Stata数据如下：Stata程序和结果如下：ranksum x, by（g）Two-sample Wilcoxon rank-sum （Mann-Whitney） test g | obs rank sum expected-+-

8、 1 | 8 47.5 72 2 | 9 105.5 81 combined | 17 153 153unadjusted variance 108.00adjustment for ties -0.13 -adjusted variance 107.87Ho: x（g=1） = x（g=2） z = -2.359 Prob |z| = 0.0183t检验的P=0.0183|z| 95% Conf. Interval-+- group | -.6711683 .0810558 -8.28 0.000 -.8300348 -.5123018 _cons | 6.780416 .1150509 5

9、8.93 0.000 6.55492 7.005912Poisson回归分组的P0.001，可以认为甲饮料中的大肠杆菌比乙饮料多。例7-10 为了解2005年某地区甲乙两镇原发性肝癌的死亡情况，该地区慢性病控制中心在甲乙两镇进行了抽样调查。甲镇调查了123443人，死于原发性肝癌的病例97人，乙镇调查了180274人，死于原发性肝癌的63人。试问甲乙两镇原发性肝癌的死亡率有无差别？Stata的数据为9712344363180274分析程序和结果为：poisson x group ,exposure（ w） LR chi2（1） = 25.86Log likelihood = -6.1989818 Pseudo R2 = 0.6759 group | -.8102747 .1618094 -5.01 0.000 -1.127415 -.493134 _cons | 5.174376 .2389773 21.65 0.000 4.70599 5.642763 W | （exposure）0.001，可以认为甲镇原发性肝癌的死亡率比乙镇高。