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

第五章参数估计和假设检验的Stata实现本章用到的Stata命令有可信区间计算立即命令为正态分布： cii 样本量 样本均数 样本标准差率： cii 样本量 阳性数Poisson分布均数： cii 观察单位数 观察值,p 原始数据的命令正态分布：ci 变量二项分布：ci 变量,b假设检验均数的检验 ：立即命令为： ttesti 样本量 样本均数 样本标准差 （待检验的总体均数）原始数据的命令：ttest 变量名 =例51 随机抽取某地25名正常成年男子，测得其血红蛋白含量如下：146139153138137125142134133122137128140137139128131158138151147144151117118该样本的均数为137.32g/L，标准差为10.63g/L，求该地正常成年男子血红蛋白含量总体均数的95%可信区间。数据格式为x11462139315341385137612571428134913310122111371212813140141371513916128171311815819138201512114722144231512411725118计算95可信区间的Stata命令为：ci x计算95可信区间结果为 Variable | Obs Mean Std. Err. 95% Conf. Interval-+- 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.9012 124.3388 该市7岁男童总体均数90%的可信区间（122.90124.34）。例53 为研究铅暴露对儿童智商(IQ)的影响，某研究调查了78名铅暴露(其血铅水平40mg/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. 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 为研究肿瘤标志物癌胚抗原(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人转阴，试估计该药物治疗幽门螺旋杆菌感染者人群的转阴率。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个单位时间内平均脉冲数估计值及其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结果为： - Poisson Exact - Variable | Exposure Mean Std. Err. 95% Conf. Interval-+- | 1 3 1.732051 .6186721 8.767273例58 大规模调查表明正常成年女子的双耳在4kHz频率时的纯音气传导听阈值平均为15dB。为研究纺机噪声对纺织女工的听力是否有影响，随机调查了20名工龄在2年以上的纺织女工，测得其听阈值(dB)如下：1011121314141617181818181920202322232426研究者的问题是：纺织女工的听阈值是否与正常成年女子不同？Stata数据格式为：x11021131241351461471681791810181118121813191420152016231722182319242026H0：m =m0， H1：mm0。即 H0：m =15(dB)，H1：m15(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: