参数估计和假设检验的Stata实现

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

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

3、9309 141.7091该地正常成年男子血红蛋白含量总体均数的95%可信区间为(132.93141.71 )例5-2 某市2005年120名7岁男童的身高 X =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 +|120123.62.4336137122.9012124.3388该市7岁男童总体均数 90%的可信区间(

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

5、s| Obs Mean Std. Err. Std. Dev. 95% Conf. Interval+x |7892.881.510458y |7888.021.3825113.34 89.8722995.8877112.2185.2670790.77293+combined |15690.45 1.039008 12.97721 88.39756 92.50244+diff |4.862.047637.8149179 8.905082差值为4.86,差值的可信区间为0.818.90例5-4为研究肿瘤标志物癌胚抗原（CEA）对肺癌的灵敏度，随机抽取140例确 诊为肺癌患者， 用CEA进行检测，

6、结果呈阳性反应者共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.90 %52.82 %。例5-5某医生用A药物治疗幽门螺旋杆菌感染者10人，其中9人转阴，试估计该药物治疗幽门螺旋杆菌感染者人群的转阴率。Stata即时命令为cii 10 9结果为-Binomial Exact -Variable

7、 |Obs Mean Std. Err.95% Conf. Interval10.9.0948683.5549839 .9974714例5 - 6某市区某年12个月发生恶性交通事故的次数分别为：5, 4, 6, 12, 7, 8, 10, 7, 6, 11, 3, 5假设每个月恶性交通事故的次数服从Poisson分布，试估计该市平均每个月恶性交通事故的次数的95%可信区间。将1个月视为一个单位时间，因 Poisson分布具有可加性，我们先计算12个单位时间内平均脉冲数估计值及其95%可信区间。X=5+4+6+12+7+8+10+7+6+11+3+5=84由于在12个单位时间内总的发生次数为84

8、,所以可以用cii 12 84,poisson得到结果：-Poisson Exact -Variable | Exposure Mean Std. Err. 95% Conf. Interval +|127 .76376265.583477 8.666438例5 - 7续例3-4。从某纯净水生产厂家生产的矿泉水中随机取1升水样，检出3个大肠菌群。试估计该家生产的矿泉水中每升水中大肠杆菌数的95%可信区间。Stata的命令为：cii 1 3,poisson结果为:-Poisson Exact -Variable | ExposureMean Std. Err.95% Conf. Interval

9、31.732051.6186721 8.767273例5-8大规模调查表明正常成年女子的双耳在4kHz频率时的纯音气传导听阈值平均为 15dB。为研究纺机噪声对纺织女工的听力是否有影响，随机调查 了 20名工龄在2年以上的纺织女工，测得其听阈值（dB）如下：10 11 12 13 14 14 16 17 18 18 18 18 19 20 20 23 22 23 2426研究者的问题是：纺织女工的听阈值是否与正常成年女子不同?Statax1102II31241351461471681791810181118121813191420152016231722182319242026H :0即H0 :0=15(dB)H ： 丰。10H115(dB)Stata命令为:ttest x=15结果为:One-sample t testVariable | Obs Mean Std. Err. Std. Dev, 95% Conf. Interval +t = 2.7810x |2017.81,0068194.50263115.692719.9073mean = mean(x)Ho: mean = 15degrees of freedom =19Ha: mean < 15Pr(T < t) = 0,9940Ha: