第五章参数估计和假设检验Stata实现.docx

1、第五章参数估计和假设检验Stata实现第五章参数估计和假设检验的Stata实现本章用到的Stata命令有可信区间计算立即命令为正态分布： cii 样本量 样本均数 样本标准差率： cii 样本量 阳性数Poisson分布均数： cii 观察单位数 观察值,p 原始数据的命令正态分布：ci 变量二项分布：ci 变量,b假设检验均数的检验 ：立即命令为： ttesti 样本量 样本均数 样本标准差 （待检验的总体均数）原始数据的命令：ttest 变量名 =例51 随机抽取某地25名正常成年男子，测得其血红蛋白含量如下：1461391531381371251421341331221371281401

2、37139128131158138151147144151117118该样本的均数为137.32g/L，标准差为10.63g/L，求该地正常成年男子血红蛋白含量总体均数的95%可信区间。数据格式为x11462139315341385137612571428134913310122111371212813140141371513916128171311815819138201512114722144231512411725118计算95可信区间的Stata命令为：ci x计算95可信区间结果为 Variable | Obs Mean Std. Err. 95% Conf. Interval-+-

3、 x | 25 137.32 2.126594 132.9309 141.7091该地正常成年男子血红蛋白含量总体均数的95%可信区间为（132.93141.71）例52 某市2005年120名7岁男童的身高=123.62(cm)，标准差s=4.75(cm)，计算该市7岁男童总体均数90%的可信区间。在Stata中有即时命令可以直接计算仅给出均数和标准差时的可信区间。cii 120 123.62 4.75,level(90)结果为： Variable | Obs Mean Std. Err. 90% Conf. Interval-+- | 120 123.62 .4336137 122.901

4、2 124.3388 该市7岁男童总体均数90%的可信区间（122.90124.34）。例53 为研究铅暴露对儿童智商(IQ)的影响，某研究调查了78名铅暴露(其血铅水平40 g/100ml)的6岁儿童，测得其平均IQ为88.02，标准差为12.21；同时选择了78名铅非暴露的6岁儿童作为对照，测得其平均IQ为92.89，标准差为13.34。试估计铅暴露的儿童智商IQ的平均水平与铅非暴露儿童相差多少，并估计两个人群IQ的总体均数之差的95%可信区间。本题也可以应用Stata的即时命令：ttesti 78 92.88 13.34 78 88.02 12.21结果：Two-sample t tes

5、t with equal variances- | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- x | 78 92.88 1.510458 13.34 89.87229 95.88771 y | 78 88.02 1.38251 12.21 85.26707 90.77293-+-combined | 156 90.45 1.039008 12.97721 88.39756 92.50244-+- diff | 4.86 2.047637 .8149179 8.905082-差值为4.86，差值的可信区间为0.818.90。例5-4 为

6、研究肿瘤标志物癌胚抗原(CEA)对肺癌的灵敏度，随机抽取140例确诊为肺癌患者，用CEA进行检测，结果呈阳性反应者共62人，试估计肺癌人群中CEA的阳性率。Stata即时命令为cii 140 62结果为 - Binomial Exact - Variable | Obs Mean Std. Err. 95% Conf. Interval-+- | 140 .4428571 .0419808 .3590149 .5291687肺癌人群中CEA的阳性率为44.28，可信区间为35.9052.82。例5-5 某医生用A药物治疗幽门螺旋杆菌感染者10人，其中9人转阴，试估计该药物治疗幽门螺旋杆菌感染者

7、人群的转阴率。Stata即时命令为cii 10 9结果为 - Binomial Exact - Variable | Obs Mean Std. Err. 95% Conf. Interval-+- | 10 .9 .0948683 .5549839 .9974714例56 某市区某年12个月发生恶性交通事故的次数分别为：5, 4, 6, 12, 7, 8, 10, 7, 6, 11, 3, 5假设每个月恶性交通事故的次数服从Poisson分布，试估计该市平均每个月恶性交通事故的次数的95%可信区间。将1个月视为一个单位时间，因Poisson分布具有可加性，我们先计算12个单位时间内平均脉冲数

8、估计值及其95%可信区间。X=5+4+6+12+7+8+10+7+6+11+3+5=84由于在12个单位时间内总的发生次数为84，所以可以用cii 12 84,poisson得到结果： - Poisson Exact - Variable | Exposure Mean Std. Err. 95% Conf. Interval-+- | 12 7 .7637626 5.583477 8.666438例57 续例34。从某纯净水生产厂家生产的矿泉水中随机取1升水样，检出3个大肠菌群。试估计该家生产的矿泉水中每升水中大肠杆菌数的95可信区间。Stata的命令为：cii 1 3,poisson结果为

9、： - Poisson Exact - Variable | Exposure Mean Std. Err. 95% Conf. Interval-+- | 1 3 1.732051 .6186721 8.767273例58 大规模调查表明正常成年女子的双耳在4kHz频率时的纯音气传导听阈值平均为15dB。为研究纺机噪声对纺织女工的听力是否有影响，随机调查了20名工龄在2年以上的纺织女工，测得其听阈值(dB)如下：10 11 12 13 14 14 16 17 18 1818 18 19 20 20 23 22 23 24 26研究者的问题是：纺织女工的听阈值是否与正常成年女子不同？Stata

10、数据格式为：x11021131241351461471681791810181118121813191420152016231722182319242026H0： = 0， H1： 0。即 H0： =15(dB)， H1： 15(dB)。Stata命令为：ttest x=15结果为：One-sample t test-Variable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval-+- x | 20 17.8 1.006819 4.502631 15.6927 19.9073- mean = mean(x) t = 2.7810Ho: mean = 15 degrees of freedom = 19 Ha: mean 15 Pr(T |t|) = 0.0119 Pr(T t) = 0.0060统计量t=2.7810，Pr(|T| |t|) = 0.0119，所以可以拒绝H0，可以认为纺织女工与正常成年女子的平均听阈值的差异有统计学意义。