統(tǒng)計描述的Stata實現(xiàn)

1、第二章 統(tǒng)計描述的Stata實現(xiàn)本章主要用到的Stata主要命令*描述性統(tǒng)計： summary 變量名 weight,detail頻數(shù)表: tabulate 變量名 weight計算均數(shù)： means 變量名 weight頻數(shù)分布圖： histogram 變量名,bin() start()*Stata的命令可以用前幾個字母代替，如summary可以由su代替，tabulate可以由tab代替。下面就直接用縮寫，不再贅述。例2-1 某市1982年調(diào)查120名20歲男子的身高（cm）資料如下，試編制頻數(shù)表和頻數(shù)圖。164.4169.2174.7175.0165.0162.8170.2160.117

2、0.9170.1175.5170.3172.3168.2166.7171.7166.8171.6165.2172.0171.7168.8171.8174.5171.7172.7166.3174.0169.0174.5171.8174.9180.0173.5178.1169.7176.1181.3173.8164.7172.2172.8178.6172.0182.5172.0173.5166.3176.1169.2176.4168.4171.0169.1166.9181.1170.4165.2168.0172.4164.3166.9176.4164.2177.2168.3177.8173.016

3、7.4173.2169.9172.1170.4174.8172.1164.4170.5165.0172.8166.8175.8171.1174.8172.7169.4178.2174.1177.2170.0172.7168.3177.1172.5166.3175.1174.4162.3171.3177.0163.5168.8177.6175.2171.5172.5179.1172.6184.8168.3182.8170.3167.5171.2170.4166.9178.5164.1169.5173.4170.3數(shù)據(jù)格式如下：x1164.42175.53171.74171.85172.26176

4、.47164.38169.99175.810168.311168.812170.313169.214170.315168.816174.917172.818168.419166.920172.121171.122177.123177.624167.525174.726172.327171.82818029178.63017131176.432170.433174.834172.535175.236171.23717538168.239174.540173.54117242169.143164.244174.845172.746166.347171.548170.44916550166.7511

5、71.752178.153182.554166.955177.256172.157169.458175.159172.560166.961162.862171.763172.764169.76517266181.167168.368164.469178.270174.471179.172178.573170.274166.875166.376176.177173.578170.479177.880170.581174.182162.383172.684164.185160.186171.68717488181.389166.390165.2911739216593177.294171.3951

6、84.896169.597170.998165.299169100173.8101176.1102168103167.4104172.8105170106177107168.3108173.4109170.1110172111174.5112164.7113169.2114172.4115173.2116166.8117172.7118163.5119182.8120170.3 為制作頻數(shù)表，鍵入Stata命令：.gen f=int(x-160)/2)*2+160產(chǎn)生用以作頻數(shù)表的新變量“f”.tab f對變量“f”作頻數(shù)表“gen”命令產(chǎn)生新變量“f”，將各觀察值轉(zhuǎn)換成相應(yīng)該組的下限值。 i

7、nt為取整函數(shù)，結(jié)果為括號內(nèi)函數(shù)值的整數(shù)部分，如int(3.24)=3。“160”為第一組的下限，“2”為組距。以第一例觀察值160.1cm為例，f=int(160.1-160)/2)*2+160=160，則它應(yīng)歸入“160”組。結(jié)果如下： f | Freq. Percent Cum.-+- 160 | 1 0.83 0.83 162 | 3 2.50 3.33 164 | 10 8.33 11.67 166 | 11 9.17 20.83 168 | 16 13.33 34.17 170 | 22 18.33 52.50 172 | 22 18.33 70.83 174 | 14 11.67

8、 82.50 176 | 10 8.33 90.83 178 | 5 4.17 95.00 180 | 3 2.50 97.50 182 | 2 1.67 99.17 184 | 1 0.83 100.00-+- Total | 120 100.00hist f,start(160) width(2)作頻數(shù)圖結(jié)果如下：例2-2某醫(yī)生測定230名正常成年男子的空腹血清胰島素樣生長因子-1（IGF-1F）水平，整理后編制為頻數(shù)分布表（表2-2），請根據(jù)該頻數(shù)分布表作頻數(shù)圖。表2-2 230名正常人空腹血清胰島素樣生長因子-1水平中位數(shù)的計算 IGF-1F （1） 人數(shù) （2）頻率(%)(3)累計頻

9、數(shù)（4） 累計頻率(%) （5）=（4）/ 50 30 13.04 30 13.0 150 71 30.87 101 43.9 250 49 21.30 150 65.2 350 28 12.17 178 77.4 450 14 6.09 192 83.5 550 12 5.22 204 88.7 650 10 4.35 214 93.0 750 8 3.48 222 96.5 850 5 2.17 227 98.7 9501050 3 1.30 230100.0 合計230100.00Stata數(shù)據(jù)格式如下： xf150302150713250494350285450146550127650

10、108750898505109503tab x weight=f制作匯總表的頻數(shù)表在本例中，x的取值為各組的下限，而f是各組的頻數(shù)，所以需要用Stata的頻數(shù)選項weight=。結(jié)果： x | Freq. Percent Cum.-+- 50 | 30 13.04 13.04 150 | 71 30.87 43.91 250 | 49 21.30 65.22 350 | 28 12.17 77.39 450 | 14 6.09 83.48 550 | 12 5.22 88.70 650 | 10 4.35 93.04 750 | 8 3.48 96.52 850 | 5 2.17 98.70

11、950 | 3 1.30 100.00-+- Total | 230 100.00hist x weight=f,start(50) width(100)制作頻數(shù)圖例2-3 隨機測量某地10名2030歲健康男性居民血清鐵含量（），測量值分別為6.58，7.42，15.32，15.78，17.60，17.98，15.21，17.53，20.11，22.64，試求其平均血清鐵含量。Stata數(shù)據(jù)為x16.5827.42315.32415.78517.6617.98715.21817.53920.111022.64su x計算均數(shù)標準差結(jié)果： Variable | Obs Mean Std. Dev

12、. Min Max-+- x | 10 15.617 5.075254 6.58 22.64所以平均平均血清鐵含量為15.617。例2-4 計算例2-1的頻數(shù)表（表2-3）中120名男子的平均身高顯然，如果用Stata，對于例2-1的資料不需要事先轉(zhuǎn)化為頻數(shù)表，然后再計算均數(shù)。不妨用例2-2作為例子展示Stata計算頻數(shù)表資料的均數(shù)。su x weight=f計算頻數(shù)表資料的均數(shù)結(jié)果： Variable | Obs Weight Mean Std. Dev. Min Max-+- x | 10 230 290 223.4123 50 950例2-5 8名麻疹易感兒接種麻疹疫苗3周后，其血凝抑制

13、抗體滴度分別為1:4，1:8，1:16，1:32，1:64，1:128，1:256，1:512。試求其平均抗體滴度。數(shù)據(jù)格式：x1428316432564612872568512gen logx=log(x)計算x的自然對數(shù)su logx計算自然對數(shù)的均數(shù)，得到3.81231disp exp(3.81231)計算均數(shù)的反對數(shù)，即指數(shù)得到幾何均數(shù)為45.254857另外，Stata還有一個直接計算幾何均數(shù)的命令：means x結(jié)果為： Variable | Type Obs Mean 95% Conf. Interval-+- x | Arithmetic 8 127.5 -20.53203 2

14、75.532 | Geometric 8 45.25483 10.94481 187.1206 | Harmonic 8 16.06275 . . -Missing values in confidence interval(s) for harmonic mean indicate that confidence interval is undefined for corresponding variable(s).Consult Reference Manual for details.Means命令計算算術(shù)均數(shù)、幾何均數(shù)以及調(diào)和均數(shù)，其中幾何均數(shù)為45.25483，和前面的結(jié)果相同。例2

15、-6 某地區(qū)50名麻疹易感兒童接種麻疹疫苗3周后，測其血凝抑制抗體滴度，如表2-4中第（1）欄和第（2）欄，求平均抗體滴度。表2-4 50名麻疹易感兒童平均抗體滴度計算表 抗體滴度 （1） 人數(shù) （2） 滴度倒數(shù) （3）（4）（5）=（2）（4）1:4 1 40.6021 0.60211:8 2 80.9031 1.80621:16 6 161.2041 7.22461:3210 321.505115.05101:6416 641.806228.89921:128 81282.107216.85761:256 52562.408212.04101:512 25122.7093 5.4186 合

16、 計5087.9003數(shù)據(jù)如下：xf14128231664321056416612887256585122gen logx=log(x)計算x的自然對數(shù)su logx weight=f計算自然對數(shù)的均數(shù)，得到4.04798disp exp(4.04798)計算均數(shù)的反對數(shù)，即指數(shù)得到幾何均數(shù)為57.281631也可以用means x weight=f結(jié)果：(analytic weights assumed) Variable | Type Obs Mean 95% Conf. Interval-+- x | Arithmetic 8 95.76 -2.30755 193.8275 | Geom

17、etric 8 57.2816 22.59723 145.2029 | Harmonic 8 32.82051 . . -Missing values in confidence interval(s) for harmonic mean indicate that confidence interval is undefined for corresponding variable(s).Consult Reference Manual for details.例2-7 5名成年男子的體重（kg）為60，70，75，80， 90，求中位數(shù)。例2-8 某醫(yī)生測定了6名正常成年男子的空腹血清胰島

18、素樣生長因子-1（IGF-1F）水平為150，170，185，245，265，280，求中位數(shù)。注：這兩個例子用軟件算實在是浪費，可以用su x,detail來處理。對于summary，加上detail后可以輸出百分位數(shù)，其中的P50就是中位數(shù)。 x- Percentiles Smallest 1% 60 60 5% 60 7010% 60 75 Obs 525% 70 80 Sum of Wgt. 550% 75 Mean 75 Largest Std. Dev. 11.1803475% 80 7090% 90 75 Variance 12595% 90 80 Skewness 099% 90 90 Kurtosis 2.05例2-9 8名7歲男孩的身高（cm）分別為116，118，119，120，121，123，125，126。請計算這8名7歲男孩身高的第25百分位數(shù)和第90百分位數(shù)。Stata數(shù)據(jù)為：x116118119120121123125126Stata命令為：su x,d結(jié)果為： x- Percentiles Smallest 1% 116 116 5% 116 11810% 116 119 Obs 825% 118.5 120 Sum of Wgt. 850% 120.5 M