醫(yī)學統(tǒng)計學實習課：STATA軟件的應用3)

1、醫(yī)學統(tǒng)計學實習課STATA軟件及其應用III課件下載賬號：stata_密碼：stata123456課件：stata3（4次）分類資料的統(tǒng)計分析分類變量的統(tǒng)計描述分類資料的統(tǒng)計推斷總體率和事件數的可信區(qū)間估計兩個率的假設檢驗卡方檢驗四格表卡方檢驗RC表卡方檢驗配對四格表卡方檢驗確切概率法分類變量的統(tǒng)計描述tabulate 變量1 變量2 , cell column row nofreqtab 分組變量 ， sum(數值變量)tab 分組變量1 分組變量2 ， sum(數值變量)cell /* 打印每個格子的頻數占總頻數的百分比column /* 打印每個格子的頻數占相應列合計的百分比row /*

2、 打印每個格子的頻數占相應行合計的百分比nofreq /* 不打印頻數例題（ex1） 有三組(group)患者，男女(sex)若干人，sex=1表示男性，sex=0表示女性。測得其血紅蛋白濃度(x,%) x group sex 1. 3.9 1 0 2. 4.2 1 0 3. 3.7 1 0 4. 4 1 0 5. 4.4 1 0 例題tab grouptab sextab sex grouptab group,sum(x)tab sex,sum(x)tab group sex,sum(x)tab groupgroup | Freq. Percent Cum.-+- 1 | 12 40.00

3、40.00 2 | 11 36.67 76.67 3 | 7 23.33 100.00-+- Total | 30 100.00tab sex sex | Freq. Percent Cum.-+- 0 | 15 50.00 50.00 1 | 15 50.00 100.00-+- Total | 30 100.00tab group sex | sex group | 0 1 | Total-+-+- 1 | 8 4 | 12 2 | 4 7 | 11 3 | 3 4 | 7 -+-+- Total | 15 15 | 30tab group,sum(x) | Summary of x gr

4、oup | Mean Std. Dev. Freq.-+- 1 | 3.8166667 .93889033 12 2 | 4.8727273 .52932203 11 3 | 3.7714286 .58513326 7-+- Total | 4.1933333 .88236879 30tab sex,sum(x) | Summary of x sex | Mean Std. Dev. Freq.-+- 0 | 4.1266667 .8224238 15 1 | 4.26 .9627342 15-+- Total | 4.1933333 .88236879 30tab group sex,sum

5、(x)Means, Standard Deviations and Frequencies of x | sex group | 0 1 | Total-+-+- 1 | 3.8125 3.825 | 3.8166667 | .90307009 1.1528949 | .93889033 | 8 4 | 12-+-+- 2 | 4.85 4.8857143 | 4.8727273 | .3872984 .6256425 | .52932203 | 4 7 | 11-+-+- 3 | 4 3.6 | 3.7714286 | .43588989 .68313003 | .58513326 | 3

6、4 | 7-+-+- Total | 4.1266667 4.26 | 4.1933333 | .8224238 .9627342 | .88236879 | 15 15 | 30率和事件計數的可信區(qū)間僅有原始資料，統(tǒng)計量未知時 ci 變量名統(tǒng)計量已知時cii 樣本含量 樣本率（發(fā)生數），選項level(#) /* 指定可信度，缺失時為 95(%)例題矯正25名學生的近視眼，3人有效，求總體率的95的可信區(qū)間；cii 25 3 樣本量 發(fā)生數cii 25 3/25 樣本量 發(fā)生率某醫(yī)師治療10名患者，治愈8人，求治愈率的90可信區(qū)間cii 10 8, level(90)四格表的卡方檢驗tabu

7、late var1 var2 fw=頻數變量，選擇項tabi a bc dchi2 /* (Pearson) 2檢驗lrchi2 /* 似然比2檢驗gamma /*Goodman-Kruskal的系數 taub /*Kendall的相關系數bV /* Cramer的列聯(lián)系數Vall /* 同時給出以上結果exact /* Fisher的確切概率cell /* 打印每個格子的頻數占總頻數的百分比 column /* 打印每個格子的頻數占相應列合計的百分比 row /* 打印每個格子的頻數占相應行合計的百分比 expected /* 計算期望頻數卡方檢驗方法的選擇卡方檢驗：N40，T 5(stat

8、a有)校正的卡方檢驗： N40，1T 5確切概率法：N40或者T 1，無須關心理論頻數(stata有)兩藥治療貧血有效率的比較組別有效人數無效人數合計試驗組431053對照組401656合計8326109Tmin=T12=53*26/109=12.64數據錄入格式：rcf1143121021402216例題(ex2)tab r c fw=f,row all | c r | 1 2 | Total-+-+- 1 | 43 10 | 53 | 81.13 18.87 | 100.00 -+-+- 2 | 40 16 | 56 | 71.43 28.57 | 100.00 -+-+- Total |

9、 83 26 | 109 | 76.15 23.85 | 100.00 Pearson chi2(1) = 1.4115 Pr = 0.235 likelihood-ratio chi2(1) = 1.4231 Pr = 0.233 Cramers V = 0.1138 gamma = 0.2647 ASE = 0.213 Kendalls tau-b = 0.1138 ASE = 0.094tabi 43 1040 16,row chi2 exact | col row | 1 2 | Total-+-+- 1 | 43 10 | 53 | 81.13 18.87 | 100.00 -+-+

10、- 2 | 40 16 | 56 | 71.43 28.57 | 100.00 -+-+- Total | 83 26 | 109 | 76.15 23.85 | 100.00 Pearson chi2(1) = 1.4115 Pr = 0.235 Fishers exact = 0.267 1-sided Fishers exact = 0.168mcci #a #b #c #d配對四個表的應用條件：1，b+c40時，使用該公式：2，20b+c40時，使用該公式：3，b+c20，使用確切概率法：配對四格表甲種培養(yǎng)法+-乙種培養(yǎng)法 +3624-10135配對四格表 mcci 36 24 10

11、135有205份咽喉涂抹樣本，把每份標本依同樣的條件分別接種于甲、乙兩種白喉桿菌培養(yǎng)基上，觀察白喉桿菌生長的情況，結果如下，問兩種培養(yǎng)基的陽性率有無差別？R*C表的卡方檢驗tabulate var1 var2 fw=頻數變量 ，選擇項tabi多個率的比較 多組構成比的比較 多個率的比較分 組 陽性例數 陰性例數 合 計陽性率(%)a=11881620492.3a=21023 3330.3a=34933338212.8合 計24737261939.9 tabi 188 1610 2349 333 , row chi2 lrchi2多組構成比的比較 A B O AB 合計 (xx=1) (xx=2

12、) (xx=3) (xx=4) 傣 族 (mz=1) 112 150 205 40 507佤 族 (mz=2) 200 112 135 73 520土家族(mz=3) 362 219 310 69 960合 計 674 481 650 182 1987tabi 112 150 205 40200 112 135 73362 219 310 69, row chi2 lrchi2 無效有效合計甲法617乙法112031合計172138tabi 6 111 20,row all exact tab r c fw=freq ,row all exact確切概率法練習1某醫(yī)生用甲乙兩法檢測同一組患者,具體結果如下,試比較這兩種方法的陽性率是否相同.甲法 乙法+-合計+351853-533285合計8850138練習2比較男性和女性暴露率是否相同。暴露未暴露合計男性13518女性21728合計341246知識點考察形式例：某醫(yī)院用A，B兩種儀器測了16名健康男青年的血紅蛋白含量（g/L），檢測結果見下表。問兩種血紅蛋白測定儀器的檢測結果是否有差別。1234567891011121314151611312512613015014513510512813510013011011512015514015013