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

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

2、1147144151117118該樣本的均數(shù)為137.32g/L，標準差為10.63g/L，求該地正常成年男子血紅蛋白含量總體均數(shù)的95%可信區(qū)間。數(shù)據(jù)格式為x11462139315341385137612571428134913310122111371212813140141371513916128171311815819138201512114722144231512411725118計算95可信區(qū)間的Stata命令為：ci x計算95可信區(qū)間結果為 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)，標準差s=4.75(cm)，計算該市7歲男童總體均數(shù)90%的可信區(qū)間。在Stata中有即時命令可以直接計算僅給出均數(shù)和標準差時的可信區(qū)間。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歲男童總體

4、均數(shù)90%的可信區(qū)間（122.90124.34）。例53 為研究鉛暴露對兒童智商(IQ)的影響，某研究調(diào)查了78名鉛暴露(其血鉛水平40mg/100ml)的6歲兒童，測得其平均IQ為88.02，標準差為12.21；同時選擇了78名鉛非暴露的6歲兒童作為對照，測得其平均IQ為92.89，標準差為13.34。試估計鉛暴露的兒童智商IQ的平均水平與鉛非暴露兒童相差多少，并估計兩個人群IQ的總體均數(shù)之差的95%可信區(qū)間。本題也可以應用Stata的即時命令：ttesti 78 92.88 13.34 78 88.02 12.21結果：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 為研究腫瘤標志物癌胚抗原(CEA)對肺癌

6、的靈敏度，隨機抽取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，可信區(qū)間為35.9052.82。例5-5 某醫(yī)生用A藥物治療幽門螺旋桿菌感染者10人，其中9人轉陰，試估計該藥物治療幽門螺旋桿菌感染者人群的轉陰率。Stata即時命令為ci

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

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

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

10、81118121813191420152016231722182319242026H0：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: mean = 15 degrees of freedom = 19 Ha: mean < 15 Ha: mean != 15 Ha: m