統(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

2、. 2160. 1170.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

3、.9176.4164.2177.2168.3177.8173.0167.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. 11

4、69. 5173.4170. 3數(shù)據(jù)格式如下:X1164.42175.53171. 74171.8|5172.26176.47164. 38169. 99175.810168. 311168.812170. 3|13169.214170. 315168.816174.917172.818168.4|19166.920172. 121171. 122177. 123177.624167. 525174. 726172. 327171.82818029178.6|3017131176.432170.433174.834172. 535175.2|36171. 23717538168.239174.

5、 540173. 541172|42169. 143164.244174.845172. 746166. 347171. 548170.44916550166. 751171. 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.8

6、80170. 581174. 182162. 383172.684164. 185160. 186171.68717488181. 389166. 390165.2911739216593177.294171. 31 95 1184.896 |169.597170.998165.299169100173.8101176. 1102168103167.4104172.8105170106177107168. 3108 |173.4109170. 1110172111174.5112164.7113169.2|114 |172.4115173.2116166.8117172.7118163.511

7、9182.8j 120 1170. 3為制作頻數(shù)表，鍵入Stata命令：.gen f=in t (x-160)/2) *2+160產(chǎn)生用以作頻數(shù)表的新變量“ f ”,tab f對變量“ f ”作頻數(shù)表“ge n”命令產(chǎn)生新變量“ f ”，將各觀察值轉(zhuǎn)換成相應該組的下限值。int 為取整函數(shù)，結(jié)果為括號內(nèi)函數(shù)值的整數(shù)部分，如int(3.24)=3 o “ 160”為第一組的下限，“ 2 ”為組距。以第一例觀察值160.1cm為例，f=int(160. 1- 160)/2)*2+160=160,則它應歸入 “1606” 組。結(jié)果如下：f Freq. Perce nt Cum160 |10. 830

8、. 83162 |32. 503. 33164 |108. 3311.67166119. 1720. 83168 |1613. 3334. 17170 |2218. 3352.50172 |2218. 3370. 83174 |1411.6782.50176108. 3390. 83178 |54. 1795. 0018032.5097. 50182 |21.6799. 17184 |10. 83100. 00+Total |120100. 00作頻數(shù)圖hist f, start(160) width(2)結(jié)果如下:I 含I2 午.0 -IIIIiI>6i 5>7» 8例

9、2-2某醫(yī)生測定230名正常成年男子的空腹血清胰島素樣生長因子-1（IGF-1F）水平，整理后編制為頻數(shù)分布表（表2-2 ），請根據(jù)該頻數(shù)分布表作頻數(shù)圖。表2-2 230名正常人空腹血清胰島素樣生長因子水平中位數(shù)的計算IGF-1F(1)人數(shù)f頻率(%) (3)累計頻數(shù)f（A累計頻率（）（5） = （ 4） / n503013.043013.0150"7130. 8710143.92504921. 3015065.2350"2812. 1717877.4450、146. 0919283.5550"125. 2220488. 7650"104. 352149

10、3.0750"83.4822296.5850"52. 1722798.7950105031. 30230100.0合計230100. 00Stata數(shù)據(jù)格式如下:Xf150302150713250494350285450146550127650108750898505109503tab x weight二f制作匯總表的頻數(shù)表tab x weight=f制作匯總表的頻數(shù)表在本例中，x的取值為各組的下限，而f是各組的頻數(shù)，所以需要用Stata制作頻數(shù)圖的頻數(shù)選項weight=。結(jié)果:x IFreq. +Perce ntCum.3013.0413. 04150 |7130. 87

11、43.912504921. 3065. 22350 |2812. 1777. 39450 |146. 0983. 48550 |125. 2288. 70650 |104. 3593. 04750 |83. 4896. 52850 |52. 1798.7095031. 30100. 00hist x weight=f, start(50) width(100)Total 230100.00制作頻數(shù)圖hist x weight=f, start(50) width(100)amol /L例2-3隨機測量某地10名2030歲健康男性居民血清鐵含量（）,測量值分別為 6. 58 , 7. 42 ,

12、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é)果：VariableOb s+Me a n Std. Dev.MinMaxX 11015.6175.0752546.5822.64所以平均平均血清鐵含量為15. 617例2-4計算例2-1的頻數(shù)表（表2-3 ）中120名男子的平均身高顯然，如果用Stata,對于例2-1的

13、資料不需要事先轉(zhuǎn)化為頻數(shù)表，然后再計算均數(shù)。不妨用例2-2作為例子展示Stata計算頻數(shù)表資料的均數(shù)。su x weight=f計算頻數(shù)表資料的均數(shù)結(jié)果:Variable+ObsWeightMea n Std. Dev.MinMaxX 110230290223.412350950例2-5 8名麻疹易感兒接種麻疹疫苗3周后，其血凝抑制抗體滴度分別為1:4, 1:8, 1:16, 1:32, 1:64, 1:128, 1:256, 1:512。試求其平均抗體滴度。數(shù)據(jù)格式：316432564612872568gen logx=log(x)計算X的自然對數(shù)su logx計算自然對數(shù)的均數(shù)，得到3.8

14、1231disp exp(3. 81231)計算均數(shù)的反對數(shù)，即指數(shù)得到幾何均數(shù)為45. 254857另外，Stata還有一個直接計算幾何均數(shù)的命令：mea ns x結(jié)果為：Variable Type Obs Mea n 95% Con f. I ntervalXArithmetic8127.5-20.53203275.532Geometric845. 2548310. 94481 187. 1206Harm on ic816. 06275 Mi ssi ng values in con fide nee in terval (s) for harm onic me a n in dicate

15、 that con fide nee in terval isun defi ned for corresp onding variable(s). Con suit Refere nee Man ual for details.Means命令計算算術(shù)均數(shù)、幾何均數(shù)以及調(diào)和均數(shù)，其中幾何均數(shù)為 45.25483,和前面的結(jié)果相同。例2-6某地區(qū)50名麻疹易感兒童接種麻疹疫苗3周后，測其血凝抑制抗體滴度，如表2-4中第(1)欄和第(2)欄，求平均抗體滴度。表2-4 50名麻疹易感兒童平均抗體滴度計算表抗體滴度(1)人數(shù)f(2)滴度倒數(shù)X (3)1g X (4)fig X(5) = ( 2)X(

16、4)1:4140. 60210. 60211:8280. 90311. 80621:166161.20417. 22461:3210321. 505115. 05101:6416641. 806228. 89921:12881282. 107216. 85761:25652562. 408212. 04101:51225122. 70935.4186合計5087. 9003數(shù)據(jù)如下：x1231643256461281016256512gen logx=log(x)計算X的自然對數(shù)su logx weight二f)計算自然對數(shù)的均數(shù)，得到 4. 04798disp exp (4. 04798)計

17、算均數(shù)的反對數(shù)，即指數(shù)得到幾何均數(shù)為57.2816312也可以用mea ns x weight=f結(jié)果:(an alytic weights assumed)Variable Type ObsMea n95% Co nf. In tervalx I Arithmetic 895.76-2. 30755 193. 8275Geometric 857.281622.59723 145.2029Harmo nic 8 32.82051.Missi ng values in con fide nee in terval (s) for harm onic mea n in dicate that co

18、n fide nee in terval is un defi ned for corresp ending variable(s). Con suit Refere nee Man ual for details.例2-7 5名成年男子的體重(kg)為60, 70, 75, 80, 90，求中位數(shù)。例2-8某醫(yī)生測定了 6名正常成年男子的空腹血清胰島素樣生長因子-1(IGF-1F )水平為 150, 170, 185, 245, 265, 280,求中位數(shù)。注：這兩個例子用軟件算實在是浪費，可以用Isu x, detail來處理。對于 summary加上detail后可以輸出百分位數(shù)，其中的P50就是中位數(shù)。PercentilesSmallest1%60605