1、關于醫學統計學列聯表資料分析1第一張，PPT共五十八頁，創作于2022年6月2第七章 列聯表資料分析列聯表：由兩個或兩個以上的變量進行交叉分類的頻數分布表。第二張，PPT共五十八頁，創作于2022年6月3列聯表（行列表）資料模式橫向變量R類, 縱向變量C類, 簡稱RC表第三張，PPT共五十八頁，創作于2022年6月4教學目的與要求 掌握：2檢驗的基本思想和用途四格表獨立性檢驗配對四格表獨立性檢驗了解：Fisher精確概率法的計算單變量擬合度檢驗頻數分布的擬合度檢驗 第四張，PPT共五十八頁，創作于2022年6月5教學內容提要 重點講解：2檢驗的基本思想四格表獨立性檢驗配對四格表獨立性檢驗介紹：

2、Fisher精確概率法的計算擬合度檢驗 第五張，PPT共五十八頁，創作于2022年6月6第一節 列聯表 2檢驗2檢驗及原理四格表資料的2檢驗雙向無序行列表的2檢驗配對計數資料的2檢驗四格表確切概率法第六張，PPT共五十八頁，創作于2022年6月7雙向無序列聯表資料單向有序列聯表資料雙向有序且屬性相同列聯表資料雙向有序且屬性不同列聯表資料一、列聯表（行列表）資料類型第七張，PPT共五十八頁，創作于2022年6月8三組消化性潰瘍住院病人四種證型的構成 1. 雙向無序列聯表資料第八張，PPT共五十八頁，創作于2022年6月92. 單向有序列聯表資料第九張，PPT共五十八頁，創作于2022年6月103

3、. 雙向有序且屬性不同列聯表資料第十張，PPT共五十八頁，創作于2022年6月114. 雙向有序且屬性相同列聯表資料第十一張，PPT共五十八頁，創作于2022年6月12四 格 表 資 料第十二張，PPT共五十八頁，創作于2022年6月13二、四格表資料的2 檢驗 2檢驗是一種檢驗實際頻數（Actual frequecy）與理論頻數（Theoretical frequecy ）之間偏離度的方法。 偏離度用2值（統計量）表示。2 =（A-T）2/T1. 2 檢驗的基本思想第十三張，PPT共五十八頁，創作于2022年6月142. 四格表資料的2檢驗 比較兩個樣本率、兩種治療方法之間的差異有無顯著性。

4、這類資料由4個基本數據構成：兩種治療方法的有效例數和無效例數，統計學稱這類資料為四格表資料。第十四張，PPT共五十八頁，創作于2022年6月152. 四格表資料的2檢驗 四格表資料的模式 組 別陽 性陰 性合計甲aba+b乙cdc+d合計a+cb+dnA11 (T11)第十五張，PPT共五十八頁，創作于2022年6月16【例7-2】為了解鉛中毒病人是否有尿棕色素增加現象，分別對病人組和對照組的尿液作尿棕色素定性檢查，結果如下，問鉛中毒病人和對照人群的尿棕色素陽性率有無差別？兩組人群尿棕色素陽性率比較第十六張，PPT共五十八頁，創作于2022年6月 建立檢驗假設和確定檢驗水準 H0：總體陽性率1

5、=2 H1：12 =0.05 計算檢驗統計量其中為Ai實際頻數，Ti為理論頻數, nR是行合計，nC是列合計，n是四格之和。【例7-2】假設檢驗步驟：第十七張，PPT共五十八頁，創作于2022年6月T11T12 理論頻數是在H0成立的條件下每個格子理論上的數據。 （1=2= 兩樣本合計率） 若H0成立，則合計的陽性率為 38/73=52.05% 則鉛中毒病人呈陽性的理論上應有 T11=3638/73=18.74人第十八張，PPT共五十八頁，創作于2022年6月2值的大小決定于： （1）A-T的差值； （2）格子數（嚴格地說是df）的多少。 T12=17.26 T21=19.26 T22=17.

6、74 =23.12 確定p值 =(R-1)(C-1)=1 , 查2界值表（附表6）得 P40，T值 2=4.125確定P值作推論1486100309012044176220第二十三張，PPT共五十八頁，創作于2022年6月【例7-4】 某礦石粉廠當生產一種礦石粉石時，在數天內即有部分工人患職業性皮膚炎，在生產季節開始，隨機抽取15名車間工人穿上新防護服，其余仍穿原用的防護服，生產進行一個月后，檢查兩組工人的皮膚炎患病率，結果如下：穿新舊兩種防護服工人的皮膚炎患病率比較第二十四張，PPT共五十八頁，創作于2022年6月問兩組工人的皮膚炎患病率有無差別？ H0：1=2；H1：12。=0.05 計算

7、 最小的理論頻數T11=1511/43=3.84, 1T1140，所以宜用2值的校正公式 =2.94 查2界值表，得0.05p 0.10 按=0.05水平不拒絕H0，尚不能認為穿不同防護服的皮膚炎患病率有差別。 若不校正，2=4.33 p0.05第二十五張，PPT共五十八頁，創作于2022年6月 【例7-5】 研究復方哌唑嗪對高血壓病治療效果的臨床試驗并與復方降壓片和安慰劑作對照，結果如下表，問三種藥物效果有無差別？組別 有效 無效 合計 有效率% 復方哌唑嗪30104075.00復方降壓片20103066.67安慰劑7253221.88合計574510255.88三種藥物治療高血壓病的有效率

8、比較第二十六張，PPT共五十八頁，創作于2022年6月27 【例7-6】某省為了解花生黃曲霉素污染，隨機觀察了三個地區的花生受黃曲霉素B1污染的情況，結果見下表。請問這三個地區花生的黃曲霉素B1污染率是否不同？某省三個地區花生的黃曲霉素B1污染率調查地區受檢樣品合計污染率（）未污染污染甲 6232979.3乙30144431.8丙 8 31127.3合計44408447.6第二十七張，PPT共五十八頁，創作于2022年6月28三、雙向無序列聯表的2檢驗 1. 列聯表的2值計算 2 = = =(R-1)(C-1) 第二十八張，PPT共五十八頁，創作于2022年6月29【例7-5】三種藥物治療高血

9、壓病有效率比較 2. 多個率的比較組別 有效 無效 合計 有效率% 復方哌唑嗪30104075.00復方降壓片20103066.67安慰劑7253221.88合計574510255.88第二十九張，PPT共五十八頁，創作于2022年6月H0：三種藥物的總體有效率相等，即1=2=3 H1：1、2、3不相等或不全相等 =0.05 2=(2-1)(3-1)=2 查2界值表得 p0.005 按=0.05水平拒絕H0，認為三種藥物的治療效果不相同或不全相同。 第三十張，PPT共五十八頁，創作于2022年6月31【例7-7】 1979年某地爆發松毛蟲病，333例患者按年齡以14歲為界分成2組，資料如表，試

10、考察兩組病人病變類型的構成比有無差別？ 某地兩組松毛蟲病患者型別構成比較 3. 構成比的比較第三十一張，PPT共五十八頁，創作于2022年6月32a. H0：兩組患者中患各型松毛蟲病的構成比相同 H1：兩組患者中患各型松毛蟲病的構成比不同 =0.05b. 2=70.14c. 查表得p0.005d.按=0.05水平拒絕H0，認為兩組患者中患各型松毛蟲病的構成比不同。3. 構成比的比較第三十二張，PPT共五十八頁，創作于2022年6月行列表的2檢驗要求理論頻數T不宜太小，要求不宜有1/5以上的格子的T5，或有一個T1。如果以上條件不能滿足，可采用： 增加樣本含量 刪去某行或某列 合理地合并部分行或

11、列 用精確概率法多個率或構成比比較的2檢驗，結論為拒絕H0時，僅表示幾組有差別，并非任2組之間都有差別。若要了解之，可進行多個率之間的兩兩比較(如2表的分割法)。4. 行列表的2檢驗注意事項 第三十三張，PPT共五十八頁，創作于2022年6月345. 多個率(構成比)兩兩比較行列表的分割（略） 【例7-5】 三種藥物治療高血壓病的有效率比較 2=組別 有效 無效 合計 有效率% 復方哌唑嗪30104075.00復方降壓片20103066.67安慰劑7253221.88合計574510255.88第三十四張，PPT共五十八頁，創作于2022年6月35分割表1 分割表2組別有效無效合計組別有效無效

12、合計哌唑嗪301040用藥組502070降壓片201030安慰劑72532合計502070合計574510221=0.583 22=21.873【例7-8】（續） 哌唑嗪組與降壓片組的有效率接近，將此2組合并，再與另一個樣本比較，得到2個分割表：20.05,1=3.84第三十五張，PPT共五十八頁，創作于2022年6月36 可以得出21+22=22.456約等于原來的2=22.36，兩者不相等是由于分割表中理論頻數的數值與總表中理論頻數不同所致。 結論：哌唑嗪與降壓片總體有效率相同,但與安慰劑的有效率不相同。第三十六張，PPT共五十八頁，創作于2022年6月37例題： 考察一種新型透析療法的效

13、果病人序號透析前透析后 131.618.2 220.7 7.3 336.426.5 433.123.7 529.522.6 620.710.7 750.325.1 831.220.9 936.623.71028.116.5四、22列聯表資料的2檢驗 第三十七張，PPT共五十八頁，創作于2022年6月38【例7-8】用兩種檢驗方法對某食品作沙門氏菌檢驗，結果如下，試比較兩種方法的陽性結果是否有差別？ 配對設計兩種方法結果比較熒光抗體法（甲法）常規培養法（乙法）合計+-+160(a)26(b)186-5(c)48(d)53合計16574239第三十八張，PPT共五十八頁，創作于2022年6月39配

14、對設計的計量資料（觀察值）比較，比較兩總體均數的差異，用配對 t檢驗配對設計的計數資料，我們用22列聯表資料的 2檢驗 第三十九張，PPT共五十八頁，創作于2022年6月401. 22列聯表資料的2檢驗 配對四格表資料分析(略)（）獨立性檢驗（是否有關聯）用四格表2檢驗法 （）優勢性檢驗（是否相同）用McNemar檢驗法 （）一致性檢驗用Kappa檢驗法第四十張，PPT共五十八頁，創作于2022年6月41【例7-8】用兩種檢驗方法對某食品作沙門氏菌檢驗，結果如下，試比較兩種方法的陽性結果是否有差別？ 配對設計兩種方法結果比較熒光抗體法（甲法）常規培養法（乙法）合計+-+160(a)26(b)1

15、86-5(c)48(d)53合計16574239著眼于兩法結果不一致的部分。第四十一張，PPT共五十八頁，創作于2022年6月42兩法結果不一致的部分： H0：總體B=C； H：總體BC檢驗統計量當b+c40時 2=當b+c40時 2= =1 樣本例數較少時還可以用精確概率法, 原理同配對設計的符號檢驗。第四十二張，PPT共五十八頁，創作于2022年6月43H0：甲法陽性乙法陰性例數與甲法陰性乙法陽性例數相同，即總體B=C；H1：總體BC =0.05 2= =12.90查2界值表 p0.005，拒絕H0，認為兩法檢驗結果不一樣，熒光抗體法陽性結果高于常規培養法。 熒光抗體法常規培養法合計+-+

16、160(a)26(b)186-5(c)48(d)53合計16574239第四十三張，PPT共五十八頁，創作于2022年6月44配對（列聯表）資料的2檢驗 【例7-9】 有93份咽喉涂抹標本，每份標本分別接種在甲乙兩種白喉桿菌培養基上，觀察其生長情況，結果如下表，試分析兩種培養基的效果是否不同？表5 甲乙兩種培養基的生長情況乙培養法甲培養法合計陽性陰性陽性37 1047陰性19 2746合計563793第四十四張，PPT共五十八頁，創作于2022年6月451. 建立假設 H0 ： B=C ；H1 ： BC 確定顯著性水平 =0.052. 計算統計量: b+c=10+19=2940， 校正公式：

17、2 =（|b-c|- 1）2/ （b+c）=2.207 =（2-1）（2-1）=13. 確定P值4. 作推斷第四十五張，PPT共五十八頁，創作于2022年6月46五、四格表的確切概率法【例7-10】 在某牧區觀察慢性布魯氏病患者植物血凝素皮膚試驗反應，資料如下，問活動型與穩定型布氏患者反應陽性率有無差別？第四十六張，PPT共五十八頁，創作于2022年6月47五、四格表的確切概率法適用條件 四格表資料 n40 時；四格表中有A= 0時； 四格表中有T1時；用其它檢驗方法所得概率接近檢驗水準時。第四十七張，PPT共五十八頁，創作于2022年6月48五、四格表的確切概率法【例7-10】 在某牧區觀察

18、慢性布魯氏病患者植物血凝素皮膚試驗反應，資料如下，問活動型與穩定型布氏患者反應陽性率有無差別？第四十八張，PPT共五十八頁，創作于2022年6月49五、四格表的確切概率法基本思想 第四十九張，PPT共五十八頁，創作于2022年6月50H0：兩型陽性率相等，即1=2 H1：12在將原四格表周邊數字（行和與列和）不變的情況下，變動a值，以得到不同情況：第五十張，PPT共五十八頁，創作于2022年6月51p為| A- T|大于等于原來表中的| A- T|值（0.14）的情況下的p(i)的和，即 p=p(0)+p(1)+p(4)=0.2668按=0.05水準接受H0,認為兩型陽性率相等。 計算每種情況

19、下的|A- T|和|A- T|大于等于原來表中的| A- T|值（0.14）的p(i)值，p(i)值按下列公式來計算：第五十一張，PPT共五十八頁，創作于2022年6月52【例7-11】 某醫師為研究乙肝免疫球蛋白預防胎兒宮內感染HBV的效果，將33例HBsAg陽性孕婦隨機分為預防注射組和非預防組，結果見下表。問兩組新生兒的HBV總體感染率有無差別？第五十二張，PPT共五十八頁，創作于2022年6月532211T11=6 A-T=2.0第五十三張，PPT共五十八頁，創作于2022年6月54第二節 分類資料的擬合優度檢驗 此例只有兩種分類。因此應有理論次數= 總數分類項數，即=700.5=35（人）= = 【例7-12】隨機抽取70名學生，調查他們對高中分文理科