卡方檢驗FREQ過程

卡方檢驗FREQ過程第1頁/共31頁一、復(fù)習(xí)有關(guān)內(nèi)容1.應(yīng)用：（1）分類變量資料；（2）推斷兩個或兩個以上的樣本率或構(gòu)成比之間有無差異；（3）檢驗頻數(shù)分布的擬合優(yōu)度。第2頁/共31頁2.四格表資料的χ2檢驗（1）四格表的形式（2）χ2檢驗的基本思想（3）基本公式

自由度=(行數(shù)-1)(列數(shù)-1)第3頁/共31頁①當(dāng)Ti≥5，且n≥40時

使用普通2檢驗，可用基本

公式或?qū)Ｓ霉?/p>

a、b、c、d為實際頻數(shù)第4頁/共31頁②四格表2值的校正公式當(dāng)1≤T<5，且n≥40時，需校正2檢驗a.基本公式的校正：b.專用公式的校正：第5頁/共31頁3、行×列表資料的2

檢驗公式：第6頁/共31頁R×C表資料的2檢驗適用條件1.理論數(shù)不能小于12.理論數(shù)1<T，且小于5的格子數(shù)不超過總格子數(shù)的1/5。若條件不適合，需作如下處理：

A

增大樣本例數(shù)

B

刪除理論數(shù)太小的行或列

C

合并（性質(zhì)相同）第7頁/共31頁二、完全隨機設(shè)計的2檢驗（SAS）例1某醫(yī)院用內(nèi)科療法治療一般類型胃潰瘍患者80例，治愈63例；治療特殊類型胃潰瘍患者99例，治愈31例。問內(nèi)科療法對兩種類型胃潰瘍的治愈率差別有無顯著意義？第8頁/共31頁

表兩種類型胃潰瘍內(nèi)科療法治療結(jié)果第9頁/共31頁dataex1;inputrccount@@;cards;1163121721312268;procfreq;tablesr*c/chisqexpectednopercentnocol;weightcount;run;第10頁/共31頁結(jié)果實際頻數(shù)理論頻數(shù)第11頁/共31頁第12頁/共31頁結(jié)果解釋本例n>40且各格子的期望值均大于5，因而選用Chi-Square的2統(tǒng)計量及其顯著性水平，即2＝39.927，P=0.0001，拒絕H0，認(rèn)為內(nèi)科療法對兩種類型胃潰瘍的治愈率差別有統(tǒng)計學(xué)意義，一般類型的治愈率高于特殊型。第13頁/共31頁

例2某省三地區(qū)花生黃曲霉素B1污染率比較

地區(qū)未污染污染合計污染率（%）甲乙630823143

29441179.331.827.3丙合計44408447.6第14頁/共31頁程序1：一般輸入方法dataex2;inputrccount@@;cards;116122321302214318323;第15頁/共31頁程序2：循環(huán)輸入法dataex2;dor=1

to

3;doc=1

to

2;inputcount@@;output;end;end;cards;623301483;第16頁/共31頁proc

freq;tablesr*c/chisq

expected

nopercent

nocol;weightcount;run;

第17頁/共31頁結(jié)果第18頁/共31頁本例各格子期望值均大于5，選用Chi-Square的2統(tǒng)計量及其顯著水平，即2=17.907，P=0.0001，按=0.05的檢驗水準(zhǔn)拒受H0，認(rèn)為三地花生黃曲霉素B1污染率有差別。

第19頁/共31頁配對設(shè)計的2檢驗（SAS程序）200名已確診的血吸蟲患者，治療前經(jīng)皮試法及糞檢法檢查，結(jié)果如下表，問兩種檢查方法的結(jié)果有無差別？第20頁/共31頁糞檢法合計+-皮試法+11240152-123648合計12476200表2血吸蟲患者兩種檢查方法的結(jié)果第21頁/共31頁dataex3;dor=1to2;doc=1to2;inputf@@;output;end;end;cards;112401236;第22頁/共31頁procfreq;tablesr*c/chisqexpectednopercent

nocolagree;weightf;run;此處用Agree選項實現(xiàn)SAS的McNemar檢驗，但要注意，SAS的McNemar檢驗不考慮校正問題，只適合于較大樣本（b+c>40），所得的結(jié)果與公式的計算結(jié)果相同。第23頁/共31頁Dataex4;Inputf11f12f21f22;Iff12+f21<40thenChisq=(abs(f12-f21)-1)**2/(f12+f21);Elsechisq=abs(f12-f21)**2/(f12+f21);P=1-probchi(chisq,1);Cards;112401236;procprint;Run;此法適用于b+c<=40的時候，使用公式編輯法計算2統(tǒng)計量及對應(yīng)的P值。第24頁/共31頁第25頁/共31頁此部分結(jié)果是普通四格表2檢驗的結(jié)果，不適于配對2檢驗使用。第26頁/共31頁1960年Cohen等提出用Kappa值作為評價判斷的一致性程度的指標(biāo)。當(dāng)兩個診斷完全一致時，Kappa值為1。當(dāng)觀測一致率大于期望一致率時，Kappa值為正數(shù)，且Kappa值越大，說明一致性越好。當(dāng)觀察一致率小于期望一致率時，Kappa值為負(fù)數(shù)，這種情況一般來說比較少見。根據(jù)邊緣概率的計算，Kappa值的范圍值應(yīng)在-1~1之間。Kappa≥0.75兩者一致性較好；0.75>Kappa≥0.4兩者一致性一般；Kappa<0.4兩者一致性較差。第27頁/共31頁本例b+c>40，故選用Mcnemar檢驗的2統(tǒng)計量及其顯著水平，即2=15.0769，P＝0.0001，按=0.05的檢驗水準(zhǔn)拒受H0，認(rèn)為兩種檢查方