第七章两独立样本定虽资料的统计分析的Stata实现例7-1为研究金届镉中蠹对大鼠肝脏中锌含量的影响，随机抽取20只小鼠，随机分为实验组和对照组。实验组小鼠每日经饮水染蠹，对照组正常饮水。1个月后，测量小鼠肝脏中锌含量（ug/ml）（表7-1），试分析两种饮用水小鼠的 肝脏中锌平均含量有无差异。表7-1实验组和对照组小鼠肝脏中锌含量（ug/ml）实验组7.145.957.108.2610.087.919.079.308.648.51对照组6.617.317.206.597.655.596.396.577.917.86例7-2为研究接触某重金届对人体血胰岛素水平有无影响，研究者从接触某重金届的职业工人中随机抽取 14人，从非接触工人中随机抽取14人，测量每个工人的血胰岛素水平（ulu/ml）（表7-2），试分析上述两个人群的血胰岛素 平均水平有无差异。表7-2某重金届接触工人和非接触工人血胰岛素水平（ulu/ml）7.84桂触页金属工人8.118.238.458.929.383.0210.224.5412.287.1111.576.0411.3817.947.429.129.457.8616.4824.12非接触工人8.578.926.749.2316.295.738.51例7-3随机抽取30名在校大学生，男性16人，女性14人，检测血红蛋白含量（g/ul）,问男性和女性大学生的血红蛋白平均含量有无差别？（数据 略）例7-4 例7-1中要检验两种饮水小鼠的肝脏锌平均含量有无差别，需要先检验实验组和对照组所来自的总体方差是否相等，即对例7-1中资料进行方差齐 性检验。例7-1的Stata数据格式如下:groupx1实验组7.142实验组5.953实验组7.14实验组8.265实验组10.086实验组7.917实验组9.078实验组9.39实验组8.6410实验组8.5111对照组6.6112对照组7.3113对照组7.214对照组6.5915对照组7.6516对照组5.5917对照组6.3918对照组6.5719对照组7.9120对照组7.86正态性检验的程序和结果如下：a =0.05swilk x if group=实验组对实验组资料进行正态性检验Shapiro-Wilk W test for normal dataVariable | ObsWVz Probz+x |100.979840.311-1.8030.96432P=0.96432 %按正态性假定进行统计分析swilk x if group=对照组对对照组资料进行正态性检验Shapiro-Wilk W test for normal dataVariable | +ObsWVzProbzx |100.935021.0010.0020.49903P=0.49903 a,按正态性假定进行统计分析方差齐性检验程序和结果如下：(: = 0.10)sdtest x,by(group)方差齐性检验Variance ratio testGroup |+ObsMeanStd. Err. Std. Dev.95% Conf. Interval对照组|106.968.2352625.74396536.4357997.500201实验组|+108.196.38405211.2144797.3272149.064786combined |207.582.26054561.1651957.036672Ho: 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.6658.127328Ha: sd(1) sd(2)P F_obs=0.0802P F_U = 0.1604P F_obs = 0.9198方差齐性检验的P= 0.1604,可以按方差齐性的假定作进一步统计分析例7-5对例7-1资料，试比较两种饮水小鼠肝脏平均锌含量有无差别？相应的Stata程序为ttest x,by(group)两独立样本 t检验输出结果为Two-sample t test with equal variancesGroup |ObsMeanStd. Err.-Std. Dev.95% Conf. Interval+-对照组|106.968.2352625.74396536.4357997.500201实验组|108.196.38405211.2144797.3272149.064786+-combined |207.582.26054561.1651957.0366728.127328+-diff |-1.228.4503826-2.174219-.2817813Degrees of freedom: 18-Ho:mean（对照组）-mean（实验组）=diff =0Ha: diff 0t =-2.7266t = -2.7266t = -2.7266P |t| =0.0138P t =0.9931t检验的P=0.01380.05 ,可以认为饮用污染的水使的小鼠肝脏中锌平均含 量高丁饮用未污染水的小鼠。例7-6对例7-2中的资料进行方差齐性检验，得 F= 4.13, P0.02 （具体计算略），可认为两总体方差不等。试比较接触重金届工人和非接触工人的血胰岛素平均水平有无差别？例7-2的Stata数据为:groupx1接触重金属工人7.842接触重金属工人8.113接触重金属工人8.454接触重金属工人9.385接触重金属工人10.226接触重金属工人12.287接触重金属工人11.578接触重金属工人11.389接触重金属工人8.2310接触重金属工人8.9211接触重金属工人3.0212接触重金属工人4.5413接触重金属工人7.1114接触重金属工人6.0415非接触工人17.9416非接触工人7.4217非接触工人9.1218非接触工人9.4519非接触工人7.8620非接触工人16.4821非接触工人24.1222非接触工人8.5723非接触工人8.9224非接触工人6.7425非接触工人9.2326非接触工人16.2927非接触工人5.7328非接触工人8.51由丁方差不齐，所以采用方差不齐的t检验。程序和结果为:ttest x,by(group) unequal两独立样本 t检验，方差不齐Two-sample t test with unequal variancesGroup |ObsMeanStd. Err.Std. Dev.95% Conf. Interval+非接触工|1411.171.4296815.3493758.08136314.25864接触重金|148.363571.70326992.6313956.8442499.882894+combined |289.766786.82708214.3765078.06975311.46382+diff |2.8064291.593291Satterthwaites degrees of freedom:18.9433 -.5290436.1419Ho: mean（非接触工）-mean（接触重金）=diff = 0Ha: diff 0t =1.7614t =1.7614t =1.7614P |t| =0.0943P t =0.0472t检验的P=0.09430.05 ,可以认为接触重金届工人和非接触工人的血胰岛 素平均水平没有差别。例7-7为研究孕妇患有妊娠合并症对葡萄糖耐受水平有无影响，对17名孕妇作葡萄糖耐受水平试验，其中有 9名孕妇曾患有妊娠合并症，而另外 8名没 有患妊娠合并症（表7-4中第（1）、（3）列）。试比较患与未患妊娠合并症的孕 妇葡萄糖平均耐受水平是否相同？表7-4患与未患妊娠合并症的孕妇葡萄糖耐受水平试验比较患妊娠合并症未患妊娠合并症葡萄糖耐受水平秩次葡萄糖耐受水平秩次(1)(2)(3)(4)11011204.51193140111339162141276184171411213281172128713510177151204.514313