列聯(lián)表的SPSS分析概要課件

1、列聯(lián)表資料的SPSS分析 在實(shí)驗(yàn)研究與基礎(chǔ)實(shí)驗(yàn)研究中，所分析的指標(biāo)可以是定量的，也可以是定性的。其定量指標(biāo),有時(shí)也轉(zhuǎn)化成定性資料進(jìn)行分析。這些定性資料或由定量資料轉(zhuǎn)化而來的定性資料，一般都整理成列聯(lián)表形式，根據(jù)資料性質(zhì)和分析目的選擇恰當(dāng)?shù)姆治龇椒ㄟM(jìn)行統(tǒng)計(jì)分析，并將統(tǒng)計(jì)分析結(jié)果與專業(yè)知識(shí)相結(jié)合, 做出合理的解釋。 我們以實(shí)驗(yàn)研究列聯(lián)表資料實(shí)例，運(yùn)用SPSS軟件包進(jìn)行統(tǒng)計(jì)分析，對(duì)操作過程和結(jié)果解釋予以詳細(xì)說明，希望能給大家提供借鑒。列聯(lián)表資料的SPSS分析 一、一般四格表（2 2 列聯(lián)表）資料 實(shí)驗(yàn)研究一般四格表（22 列聯(lián)表）資料分析目的主要有 2 個(gè)，一是分析兩個(gè)比率總體的差別有無統(tǒng)計(jì)學(xué)意義或

2、兩樣本某指標(biāo)的分布（或構(gòu)成）總體是否相同，二是分析兩個(gè)分類特征是否有關(guān)聯(lián)。 例1. 某院欲比較異梨醇口服液（試驗(yàn)組）和氫氯噻嗪+地塞米松（對(duì)照組）降低顱內(nèi)壓的療效，將200例顱內(nèi)壓增高患者隨機(jī)分為2 組, 見表1。表1 試驗(yàn)組和對(duì)照組降低顱內(nèi)壓療效的比較組別有效無效有效率(%)合計(jì)試驗(yàn)組99595.20104對(duì)照組752178.1396合計(jì)1742687.00200 一、一般四格表（2 2 列聯(lián)表）資料組別有效無效 數(shù)據(jù)錄入： 打開SPSS；點(diǎn)擊Variable View 定義變量，變量 1 Name 為“group”，Type 為“String”；變量 2 Name 為“effect”，Ty

3、pe 為“String”；變量3 Name 為 “count”；Type 為 “Numeric”，Decimals 為“0”，其它為默認(rèn)設(shè)置；點(diǎn)擊Data View 輸入數(shù)據(jù)。 分析過程： 頻數(shù)加權(quán)（所有列聯(lián)表資料均需經(jīng)過頻數(shù)加權(quán), 以下例題分析中省略該過程）： Data Weight Cases Weight Cases by： Frequency Variable：count OK 數(shù)據(jù)錄入： 分析： Analyze Descriptive Statistics Crosstabs Rows：group Columns: effect Statistics： Chi- square：Con

4、tinue OK 分析： 主要結(jié)果與解釋：根據(jù)22表資料 c2 檢驗(yàn)的應(yīng)用條件：總例數(shù)（n） 40 且所有的理論值（T） 5 時(shí)，選用一般的 c2 檢驗(yàn)；P a 時(shí), 改用Fisher 確切概率法； n 40, 有1 T 5 時(shí),選用連續(xù)校正的檢驗(yàn)；或改用Fisher確切概率法; n 40，或T 1 時(shí)，用Fisher 確切概率法。 該資料n大于40，所有T 均大于5，可取Pearson c2 值和似然比（Likelihood ratio）c2值，二者 c2 值分別為12.123 和12.864，P 40，且所有 T 5，取Pearson c2 檢驗(yàn)結(jié)果，c2 = 9.981， P 0.01，

5、可認(rèn)為嬰兒腹瀉與喂養(yǎng)方式有關(guān)聯(lián)；Pearson 列聯(lián)系數(shù) C = 0.329，根據(jù) C 0.7為高度關(guān)聯(lián)、C 0.4 為中度關(guān)聯(lián)、C 0.4 為低度關(guān)聯(lián)的判斷原則，二者關(guān)聯(lián)強(qiáng)度較弱。 主要結(jié)果與解釋： 二、配對(duì)四格表（22 列聯(lián)表）資料 例3. 某研究者對(duì)同一批標(biāo)本進(jìn)行兩個(gè)指標(biāo)的檢測(cè)，結(jié)果見表3。問： 兩個(gè)指標(biāo)的分布有無差別? 兩個(gè)指標(biāo)有無關(guān)聯(lián)?表3 同一批標(biāo)本兩個(gè)指標(biāo)的檢測(cè)結(jié)果 甲指標(biāo)乙指標(biāo)+-+-2173423162 二、配對(duì)四格表（22 列聯(lián)表）資料甲指標(biāo)乙指標(biāo)+- 數(shù)據(jù)錄入： 定義變量，變量 1 Name 為 “A”，變量 2 Name 為“B”，變量 3 Name 為“count”，

6、Type 均為 “Numeric”，Decimals 均為 “0”，其它為默認(rèn)設(shè)置； 輸入數(shù)據(jù)，A= 1、B= 1：甲指標(biāo) “+”、乙指標(biāo)“+”；A= 1、B= 2：甲指標(biāo) “+”、乙指標(biāo)“-”； A= 2、B= 1：甲指標(biāo)“-”、乙指標(biāo)“+”；A= 2、B= 2：甲指標(biāo)“-”、乙指標(biāo)“-”。 分析過程： Analyze Descriptive Statistics Crosstabs Rows：A Columns：B Statistics： Chi square， Contingency coefficient， McNemar Continue OK 數(shù)據(jù)錄入： 主要結(jié)果與解釋： McNe

7、mar c2 檢驗(yàn)（SPSS McNemar 檢驗(yàn)法不給 c2 值），P = 0.185，說明兩個(gè)指標(biāo)分布的差別無統(tǒng)計(jì)學(xué)意義；一般c2 檢驗(yàn)，Pearson c2 = 235.821， P 0.01，說明兩個(gè)指標(biāo)存在關(guān)聯(lián)，Pearson 列聯(lián)系數(shù) C = 0.592，關(guān)聯(lián)程度中等。 主要結(jié)果與解釋： 例4. 某研究者用胃鏡和活檢對(duì)1 321 例患者進(jìn)行胃癌檢查，結(jié)果見表 4。問胃鏡和活檢的結(jié)果是否一致?表4 1 321 例患者胃鏡和活檢結(jié)果胃鏡活檢+-+-1434481 126 例4. 某研究者用胃鏡和活檢對(duì)1 321 例患者進(jìn) 數(shù)據(jù)錄入： 同例 3。 分析過程： Analyze Descr

8、iptive Statistics Crosstabs Rows: A Columns: B Statistics : Kappa, McNemar, Continue OK 數(shù)據(jù)錄入： 主要結(jié)果與解釋： McNemar 檢驗(yàn)，P 0.01，說明兩種檢查方法有差別；Kappa 檢驗(yàn)（SPSS 不給可信區(qū)間），Kappa = 0.824，P 0. 01，說明兩種檢查具有較好的一致性。 可以看出, 兩種檢驗(yàn)結(jié)果是矛盾的。為什么呢? McNemar 法一般用于樣本含量 n 不太大的資料，因只考慮結(jié)果不一致的情況，而未考慮樣本含量 n 和結(jié)果一致的情況，所以，當(dāng) n 很大且結(jié)果一致率高時(shí)，不一致的數(shù)值

9、相對(duì)較小，容易出現(xiàn)有統(tǒng)計(jì)學(xué)意義的檢驗(yàn)結(jié)果，但實(shí)際意義可能不大。本例即是如此，應(yīng)以一致性檢驗(yàn)結(jié)果為準(zhǔn)。實(shí)際上，對(duì)于兩種檢驗(yàn)（查）方法或診斷方法結(jié)果進(jìn)行分析時(shí)，主要分析的也就是一致性。根據(jù)Kappa 值判斷一致性強(qiáng)度的標(biāo)準(zhǔn)尚有爭(zhēng)議，一般認(rèn)為：Kappa 值 0.75 時(shí)，有極好的一致性。 主要結(jié)果與解釋： 三、R C 表資料 1. 雙向無序 R C 表資料 例5. 某醫(yī)生研究物理療法、藥物治療和外用膏藥三種療法治療周圍性面神經(jīng)麻痹的療效，資料見表5，問三種療法的有效率有無差別?表5 三種療法治療周圍性面神經(jīng)麻痹的療效療 法有效無效合計(jì)有效率（%）物理療法藥物療法外用膏藥19916411871826

10、26018214496.6090.1181.94合 計(jì)4815153290.41 三、R C 表資料療 法有效無效合計(jì)有效率（% 數(shù)據(jù)錄入與分析過程，見例1。 主要結(jié)果與解釋： 該資料的原因變量（療法）與結(jié)果變量（有效與無效）均無序，屬雙向無序 RC 表資料，選用一般 c2 檢驗(yàn)分析即可。RC 表資料c2 檢驗(yàn)時(shí)，要求所有T 1 且T 5，取Pearson c2 值或似然比c2 值，c2 = 21.038（或21.559），P 5，c2 = 213.162（或248.143）, P 0.01，認(rèn)為該地 ABO 血型和 MN 血型系統(tǒng)之間有關(guān)聯(lián)。列聯(lián)系數(shù)C = 0.188，關(guān)聯(lián)性很弱。 數(shù)據(jù)錄入

11、與分析過程: 2. 單向有序RC 表資料 例7 某醫(yī)生用 3 種藥物治療某病患者，療效如表7。問 3 種藥物的療效有無差別?表7 3 種藥物對(duì)某病患者的療效藥物種類治愈顯效好轉(zhuǎn)無效合計(jì)ABC154149915315045522241008585合 計(jì)207312651270 2. 單向有序RC 表資料 藥物種類治愈顯效好轉(zhuǎn)無 數(shù)據(jù)錄入： 定義變量，變量1 Name 為“medicine”，變量 2 Name 為“result”，變量 3 Name 為“count”，Type 均為“Numeric”，Decimals 均為“0”，其它為默認(rèn)設(shè)置； 輸入數(shù)據(jù)： medicine：1 為“A藥”，2

12、 為“B 藥”，3 為“C 藥”，result：1 為“治愈”，2 為“顯效”，3 為“好轉(zhuǎn)”，4 為“無效”。 數(shù)據(jù)錄入：分析過程： Analyze Nonparametric tests K Independent Samples Test Variable List: result Grouping Variable: medicine Define Range: Range for Grouping Variable: Minimum: 1; Maximum: 3; Continue Kruskal- Wallis H OK分析過程： 主要結(jié)果與解釋： 該資料原因變量“藥物種類”屬于無

13、序變量，而結(jié)果變量“療效”屬于有序變量；可以選用秩和檢驗(yàn)、Ridit分析或有序變量的Logistic回歸分析，本例選用秩和檢驗(yàn)分析。結(jié)果，c2 (Hc) = 61.146，P 0.05，故尚不能認(rèn)為高血壓患者的血壓控制情況與其肥胖程度有相關(guān)關(guān)系。Kendall 和 Spearman 等級(jí)相關(guān)分析結(jié)果：t = - 0.074，P = 0.401；rs = - 0.085，P = 0.381，也不支持兩者之間有相關(guān)關(guān)系。該資料若把“肥胖程度”視為無序的分組變量，則可按照“單向有序RC 表資料”進(jìn)行秩和檢驗(yàn)分析。 主要結(jié)果與解釋： 4. 雙向有序且屬性相同 RC 表資料 例9 兩位放射科醫(yī)生對(duì)一批矽

14、肺胸片獨(dú)自作出了矽肺分級(jí)診斷，請(qǐng)問他們的診斷結(jié)果是否一致。表9 兩位醫(yī)生對(duì)一批矽肺胸片的分級(jí)診斷結(jié)果甲醫(yī)生乙醫(yī)生合計(jì)1級(jí)2級(jí)3級(jí)1級(jí)2級(jí)3級(jí)32101554701245476752合計(jì)337657166 4. 雙向有序且屬性相同 RC 表資料 甲醫(yī)生乙 數(shù)據(jù)錄入與分析過程： 見例4。 主要結(jié)果與解釋： 該資料兩位醫(yī)生的診斷結(jié)果有序（等級(jí)）且屬性相同，選用一致性檢驗(yàn)。結(jié)果, Kappa= 0.676 0.75，P 0.01，認(rèn)為兩位醫(yī)生的診斷存在一致性，一致性中等。 數(shù)據(jù)錄入與分析過程： 四、多個(gè)樣本率間的兩兩比較 對(duì)于 K2 表或 2K 表資料，總體檢驗(yàn)有統(tǒng)計(jì)學(xué)意義時(shí)，只能說明各總體比率（分布

15、）之間總的說來有差別，但不能說明任兩個(gè)總體比率（分布）之間有差別。若要對(duì)每?jī)蓚€(gè)總體比率（分布）之間做出有無差別的推斷，需進(jìn)一步分析。 四、多個(gè)樣本率間的兩兩比較 1. 組間全面比較 Bonferroni 法 設(shè)總的檢驗(yàn)水準(zhǔn)為a，進(jìn)行比較的次數(shù)為m，令各次比較的檢驗(yàn)水準(zhǔn) a= a/ m, 并規(guī)定P a 時(shí)拒絕 H0，基于這樣的做法，就可以把 類錯(cuò)誤的累積概率控制在a。當(dāng)多組間比較次數(shù)不多時(shí)，該法的效果較好。但是，當(dāng)組間比較次數(shù)較多（如10 次以上）時(shí)，檢驗(yàn)水準(zhǔn)過低，導(dǎo)致結(jié)果過于保守，犯類錯(cuò)誤的概率增加。 以例5為例，分析三種療法任兩種療法之間療效有無差別。 1. 組間全面比較 Bonferron

16、i 法 求總的c2 值，P a 時(shí)，作兩兩比較。本例 Pearson c2 = 21.038，P 0.01； 求各四格表的 c2 值及 P 值； 按公式計(jì)算校正的 a 值，式中 k 為比較的組數(shù)，本例為 3； 各 P 值與校正的 a 值比較，判別有無統(tǒng)計(jì)學(xué)意義。表10 三種療法有效率的兩兩比較比較組c2值P值a值統(tǒng)計(jì)學(xué)意義物理療法與藥物治療物理療法與外用膏藥藥物治療與外用膏藥6.75621.3234.5910.0090.0000.0320.01250.01250.0125有有無 求總的c2 值，P a，表明尚不能認(rèn)為此兩種療法有效率有差別。 結(jié)果與解釋： 2. 多個(gè)實(shí)驗(yàn)組與同一對(duì)照組比較 Br

17、unden 法 若多個(gè)實(shí)驗(yàn)組與同一對(duì)照組進(jìn)行比較時(shí)，設(shè)總的檢驗(yàn)水準(zhǔn)為a，樣本率個(gè)數(shù)為k，則每個(gè)實(shí)驗(yàn)組與對(duì)照組比較的檢驗(yàn)水準(zhǔn)a，并規(guī)定 P a時(shí)拒絕 H0。仍以例 5 為例，以外用膏藥作為對(duì)照，分析物理療法和藥物治療兩種療法與物理療法之間療效有無差別。 2. 多個(gè)實(shí)驗(yàn)組與同一對(duì)照組比較 Brunden 求總的c2 值，P a時(shí)，作兩兩比較。本例Pearson c2 = 21.038，P 0.01； 求各實(shí)驗(yàn)組與對(duì)照組四格表的c2值及P 值； 按公式計(jì)算校正的a值； 各P值與校正的a值比較，判別有無統(tǒng)計(jì)學(xué)意義。表11 物理療法和藥物治療與外用膏藥法有效率的比較比較組c2值P 值a值統(tǒng)計(jì)學(xué)意義物理療

18、法與外用膏藥藥物治療與外用膏藥21.3234.5910.0000.0320.01250.0125有無 求總的c2 值，P a時(shí)，作兩兩 結(jié)果與解釋： 物理療法與外用膏藥比較，Pearson c2 = 21.323，P a，表明尚不能認(rèn)為此兩種療法的有效率有差別。 結(jié)果與解釋： 進(jìn)行多個(gè)樣本率的多重比較時(shí)，除上述介紹的兩種方法外，還有將RC 表分割成若干個(gè)四格表進(jìn)行分析的 c2 分割法、通過調(diào)整檢驗(yàn)界值進(jìn)行多個(gè)樣本率間兩兩比較的杜養(yǎng)志法和羅文海法以及通過計(jì)算兩率之差的可信區(qū)間來推斷比較組間有無差別的 Scheffe 法等可供選擇。 進(jìn)行多個(gè)樣本率的多重比較時(shí)，除上述介參考文獻(xiàn)： 李永紅. 實(shí)驗(yàn)研究列聯(lián)表資料的 SPSS 分析J .中國(guó)熱帶醫(yī)學(xué)， 2007， 7(10) : 19411958。 程琮，程瑋，范華. 行列相關(guān)的測(cè)度J. 中國(guó)衛(wèi)生統(tǒng)