1、第五章 参数估计和假设检验的 Stata 实现本章用到的 Stata 命令有可信区间计算 立即命令为正态分布： cii 样本量 样本均数 样本标准差率： cii 样本量 阳性数Poisson 分布均数： cii 观察单位数 观察值,p原始数据的命令正态分布：ci 变量二项分布：ci 变量,b假设检验 均数的检验 ：立即命令为：ttesti 样本量 样本均数 样本标准差 （待检验的总体均数）0原始数据的命令：ttest 变量名 = 0例 51 随机抽取某地 25 名正常成年男子，测得其血红蛋白含量如下：146 139 153 138 137 125 142 134 133 122 137 128
2、 140137 139 128 131 158 138 151 147 144 151 117 118该样本的均数为 137.32g/L，标准差为 10.63g/L，求该地正常成年男子血红蛋白含量总体均数的 95%可信区间。数据格式为x1 1462 1393 1534 1385 1376 1257 1428 1349 13310 12211 13712 12813 14014 13715 13916 12817 13118 15819 13820 15121 14722 14423 15124 11725 118计算 95可信区间的 Stata 命令为：ci x 计算 95可信区间结果为Var
3、iable | 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)，标准差Xs=4.75(cm)，计算该市 7 岁男童总体均数 90%的可信区间。在 Stata 中有即时命令可以直接计算仅给出均数和标准差时的可信区间。cii 120 123.62 4.75,level(90)结果为：Variable | Obs Mean
4、Std. Err. 90% Conf. Interval-+-| 120 123.62 .4336137 122.9012 124.3388 该市 7 岁男童总体均数 90%的可信区间（122.90124.34） 。例 53 为研究铅暴露对儿童智商(IQ) 的影响，某研究调查了 78 名铅暴露(其血铅水平 40g/100ml)的 6 岁儿童，测得其平均 IQ 为 88.02，标准差为12.21；同时选择了 78 名铅非暴露的 6 岁儿童作为对照，测得其平均 IQ 为92.89，标准差为 13.34。试估计铅暴露的儿童智商 IQ 的平均水平与铅非暴露儿童相差多少，并估计两个人群 IQ 的总体均数
5、之差的 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.88771y | 78 88.02 1.38251 12.21 85.26707 90.77293-+-combined | 156 90.45 1.039008 12.97721 88.3975
6、6 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 .52
7、91687肺癌人群中 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,
8、 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.58
9、3477 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。为研究纺机噪声对纺织女工的听力是否有影响，随机调查
10、了 20 名工龄在 2 年以上的纺织女工，测得其听阈值(dB)如下：10 11 12 13 14 14 16 17 18 1818 18 19 20 20 23 22 23 24 26研究者的问题是：纺织女工的听阈值是否与正常成年女子不同？Stata 数据格式为：x1 102 113 124 135 146 147 168 179 1810 1811 1812 1813 1914 2015 2016 2317 2218 2319 2420 26H0： =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 = 19Ha: mean 15Pr(T |t|) = 0.0119 Pr(T t) = 0.0060统计量 t=2.7810，Pr(|T| |t|) = 0.0119，所以可以拒绝 H0，可以认为纺织女工与正常成年女子的平均听阈值的差异有统计学意义。
