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

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

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

3、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 .122.9012 124.3388该市7岁男童总体均数90%的可信区间(122.90124.34)。例5 3

4、为研究铅暴露对儿童智商(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 test with equal variances| Obs Mean Std. Err.

5、 Std. Dev. 95% Conf. Interval +x |7892.881.51045813.3489.8722995.88771y |7888.021.3825112.2185.2670790.77293+combined |15690.45 1.039008 12.97721 88.39756 92.50244+diff |4.862.047637. 8.905082差值为4.86,差值的可信区间为0.818.90。例5-4为研究月中瘤标志物癌胚抗原（CEA）对肺癌的灵敏度，随机抽取 140 例确诊为肺癌患者，用CEA进行检测，结果呈阳性反应者共62人，试估计肺癌 人群中CEA的

6、阳性率。Stata即时命令为cii 140 62结果为-Binomial Exact -Variable | Obs Mean Std. Err. 95% Conf. Interval +|140.肺癌人群中CEA的阳性率为44.28%,可信区间为35.90%52.82%。例5-5某医生用A药物治疗幽门螺旋杆菌感染者10人，其中9人转阴， 试估计该药物治疗幽门螺旋杆菌感染者人群的转阴率。Stata即时命令为cii 10 9结果为-Binomial Exact -Variable | Obs Mean Std. Err. 95% Conf. Interval +|10.9.例5-6某市区某年12

7、个月发生恶性交通事故的次数分别为：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,所以可以用cii 12 84,poisson 得到结果：-Poisson Exact -Variable | Exposure Mean Std. Err. 95

8、% Conf. Interval +|127 .5.583477 8.666438例5- 7续例3-4。从某纯净水生产厂家生产的矿泉水中随机取 1升水样, 检出3个大肠菌群。试估计该家生产的矿泉水中每升水中大肠杆菌数的 95%可 信区间。Stata的命令为:cii 1 3,poisson结果为：-Poisson Exact -Variable | Exposure Mean Std. Err. 95% Conf. Interval +|13 1.732051. 8.767273例5- 8大规模调查表明正常成年女子的双耳在 4kHz频率时的纯音气传导听阈值平均为15dB。为研究纺机噪声对纺织女工

9、的听力是否有影响，随机调查 了 20名工龄在2年以上的纺织女工，测得其听阈值（dB）如下:10 11 12 13 14 14 16 17 18 1818 18 19 20 20 23 22 23 24 26研究者的问题是：纺织女工的听阈值是否与正常成年女子不同?Ho：= 0,Hi : w 0。H0：=15(dB),H1： w15(dB)。Stata数据格式为:x11021131241351461471681791810181118121813191420152016231722182319242026Stata命令为:ttest x=15结果为:One-sample t testVariable | Obs Mean Std. Err. Std. Dev. 95% Conf. Interval +x |2017.8 1.006819 4.50263115.692719.9073mean = mean(x)t = 2.7810Ho: mean = 15degrees of freedom =19Ha: mean < 15Ha: mean != 15Ha: mean &