第五章參數(shù)估計(jì)和假設(shè)檢驗(yàn)的Stata實(shí)現(xiàn)

1、第五章參數(shù)估計(jì)和假設(shè)檢驗(yàn)的Stata實(shí)現(xiàn)本章用到的Stata命令有可信區(qū)間計(jì)算立即命令為正態(tài)分布： cii 樣本量 樣本均數(shù) 樣本標(biāo)準(zhǔn)差率： cii 樣本量 陽(yáng)性數(shù)Poisson分布均數(shù)： cii 觀察單位數(shù) 觀察值,p 原始數(shù)據(jù)的命令正態(tài)分布：ci 變量二項(xiàng)分布：ci 變量,b假設(shè)檢驗(yàn)均數(shù)的檢驗(yàn) ：立即命令為： ttesti 樣本量 樣本均數(shù) 樣本標(biāo)準(zhǔn)差 （待檢驗(yàn)的總體均數(shù)）原始數(shù)據(jù)的命令：ttest 變量名 =例51 隨機(jī)抽取某地25名正常成年男子，測(cè)得其血紅蛋白含量如下：14613915313813712514213413312213712814013713912813115813815

2、1147144151117118該樣本的均數(shù)為137.32g/L，標(biāo)準(zhǔn)差為10.63g/L，求該地正常成年男子血紅蛋白含量總體均數(shù)的95%可信區(qū)間。數(shù)據(jù)格式為x11462139315341385137612571428134913310122111371212813140141371513916128171311815819138201512114722144231512411725118計(jì)算95可信區(qū)間的Stata命令為：ci x計(jì)算95可信區(qū)間結(jié)果為 Variable | Obs Mean Std. Err. 95% Conf. Interval-+- x | 25 137.32 2.12

3、6594 132.9309 141.7091該地正常成年男子血紅蛋白含量總體均數(shù)的95%可信區(qū)間為（132.93141.71）例52 某市2005年120名7歲男童的身高=123.62(cm)，標(biāo)準(zhǔn)差s=4.75(cm)，計(jì)算該市7歲男童總體均數(shù)90%的可信區(qū)間。在Stata中有即時(shí)命令可以直接計(jì)算僅給出均數(shù)和標(biāo)準(zhǔn)差時(shí)的可信區(qū)間。cii 120 123.62 4.75,level(90)結(jié)果為： Variable | Obs Mean Std. Err. 90% Conf. Interval-+- | 120 123.62 .4336137 122.9012 124.3388 該市7歲男童總體

4、均數(shù)90%的可信區(qū)間（122.90124.34）。例53 為研究鉛暴露對(duì)兒童智商(IQ)的影響，某研究調(diào)查了78名鉛暴露(其血鉛水平40mg/100ml)的6歲兒童，測(cè)得其平均IQ為88.02，標(biāo)準(zhǔn)差為12.21；同時(shí)選擇了78名鉛非暴露的6歲兒童作為對(duì)照，測(cè)得其平均IQ為92.89，標(biāo)準(zhǔn)差為13.34。試估計(jì)鉛暴露的兒童智商IQ的平均水平與鉛非暴露兒童相差多少，并估計(jì)兩個(gè)人群IQ的總體均數(shù)之差的95%可信區(qū)間。本題也可以應(yīng)用Stata的即時(shí)命令：ttesti 78 92.88 13.34 78 88.02 12.21結(jié)果：Two-sample t test with equal varian

5、ces- | 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，差值的可信區(qū)間為0.818.90。例5-4 為研究腫瘤標(biāo)志物癌胚抗原(CEA)對(duì)肺癌

6、的靈敏度，隨機(jī)抽取140例確診為肺癌患者，用CEA進(jìn)行檢測(cè)，結(jié)果呈陽(yáng)性反應(yīng)者共62人，試估計(jì)肺癌人群中CEA的陽(yáng)性率。Stata即時(shí)命令為cii 140 62結(jié)果為 - Binomial Exact - Variable | Obs Mean Std. Err. 95% Conf. Interval-+- | 140 .4428571 .0419808 .3590149 .5291687肺癌人群中CEA的陽(yáng)性率為44.28，可信區(qū)間為35.9052.82。例5-5 某醫(yī)生用A藥物治療幽門螺旋桿菌感染者10人，其中9人轉(zhuǎn)陰，試估計(jì)該藥物治療幽門螺旋桿菌感染者人群的轉(zhuǎn)陰率。Stata即時(shí)命令為ci

7、i 10 9結(jié)果為 - Binomial Exact - Variable | Obs Mean Std. Err. 95% Conf. Interval-+- | 10 .9 .0948683 .5549839 .9974714例56 某市區(qū)某年12個(gè)月發(fā)生惡性交通事故的次數(shù)分別為：5, 4, 6, 12, 7, 8, 10, 7, 6, 11, 3, 5假設(shè)每個(gè)月惡性交通事故的次數(shù)服從Poisson分布，試估計(jì)該市平均每個(gè)月惡性交通事故的次數(shù)的95%可信區(qū)間。將1個(gè)月視為一個(gè)單位時(shí)間，因Poisson分布具有可加性，我們先計(jì)算12個(gè)單位時(shí)間內(nèi)平均脈沖數(shù)估計(jì)值及其95%可信區(qū)間。X=5+4+

8、6+12+7+8+10+7+6+11+3+5=84由于在12個(gè)單位時(shí)間內(nèi)總的發(fā)生次數(shù)為84，所以可以用cii 12 84,poisson得到結(jié)果： - Poisson Exact - Variable | Exposure Mean Std. Err. 95% Conf. Interval-+- | 12 7 .7637626 5.583477 8.666438例57 續(xù)例34。從某純凈水生產(chǎn)廠家生產(chǎn)的礦泉水中隨機(jī)取1升水樣，檢出3個(gè)大腸菌群。試估計(jì)該家生產(chǎn)的礦泉水中每升水中大腸桿菌數(shù)的95可信區(qū)間。Stata的命令為：cii 1 3,poisson結(jié)果為： - Poisson Exact -

9、 Variable | Exposure Mean Std. Err. 95% Conf. Interval-+- | 1 3 1.732051 .6186721 8.767273例58 大規(guī)模調(diào)查表明正常成年女子的雙耳在4kHz頻率時(shí)的純音氣傳導(dǎo)聽(tīng)閾值平均為15dB。為研究紡機(jī)噪聲對(duì)紡織女工的聽(tīng)力是否有影響，隨機(jī)調(diào)查了20名工齡在2年以上的紡織女工，測(cè)得其聽(tīng)閾值(dB)如下：1011121314141617181818181920202322232426研究者的問(wèn)題是：紡織女工的聽(tīng)閾值是否與正常成年女子不同？Stata數(shù)據(jù)格式為：x110211312413514614716817918101

10、81118121813191420152016231722182319242026H0：m =m0， H1：mm0。即 H0：m =15(dB)，H1：m15(dB)。Stata命令為：ttest x=15結(jié)果為：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 Ha: mean != 15 Ha: m