卡方檢驗(yàn)(FREQ過程)

《卡方檢驗(yàn)(FREQ過程)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、實(shí)驗(yàn)五χ2檢驗(yàn)?zāi)康囊?1.掌握分類變量資料的統(tǒng)計(jì)推斷的方法和計(jì)算過程。2.掌握procfreq語句的格式、用途，且能對(duì)結(jié)果作全面分析。2021/10/41一、復(fù)習(xí)有關(guān)內(nèi)容1.應(yīng)用：（1）分類變量資料；（2）推斷兩個(gè)或兩個(gè)以上的樣本率或構(gòu)成比之間有無差異；（3）檢驗(yàn)頻數(shù)分布的擬合優(yōu)度。2021/10/422.四格表資料的χ2檢驗(yàn)（1）四格表的形式（2）χ2檢驗(yàn)的基本思想（3）基本公式自由度?=(行數(shù)-1)(列數(shù)-1)2021/10/43①當(dāng)Ti≥5，且n≥40時(shí)使用普通?2檢驗(yàn)，可用基本公式或?qū)Ｓ霉絘、b、c、d為實(shí)際頻數(shù)2021/10/44②四格表?2值的校正公式當(dāng)1

2、≤T<5，且n≥40時(shí)，需校正?2檢驗(yàn)a.基本公式的校正：b.專用公式的校正：2021/10/453、行×列表資料的?2檢驗(yàn)公式：2021/10/46R×C表資料的?2檢驗(yàn)適用條件1.理論數(shù)不能小于12.理論數(shù)1

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

4、型。2021/10/413例2某省三地區(qū)花生黃曲霉素B1污染率比較地區(qū)未污染污染合計(jì)污染率（%）甲乙63082314329441179.331.827.3丙合計(jì)44408447.62021/10/414程序1：一般輸入方法dataex2;inputrccount@@;cards;116122321302214318323;2021/10/415程序2：循環(huán)輸入法dataex2;dor=1to3;doc=1to2;inputcount@@;output;end;end;cards;623301483;2021/10/416procfreq;tablesr*c/chisqexp

5、ectednopercentnocol;weightcount;run;2021/10/417結(jié)果2021/10/418本例各格子期望值均大于5，選用Chi-Square的?2統(tǒng)計(jì)量及其顯著水平，即?2=17.907，P=0.0001，按?=0.05的檢驗(yàn)水準(zhǔn)拒受H0，認(rèn)為三地花生黃曲霉素B1污染率有差別。2021/10/419配對(duì)設(shè)計(jì)的?2檢驗(yàn)（SAS程序）200名已確診的血吸蟲患者，治療前經(jīng)皮試法及糞檢法檢查，結(jié)果如下表，問兩種檢查方法的結(jié)果有無差別？2021/10/420糞檢法合計(jì)+-皮試法+11240152-123648合計(jì)12476200表2血吸蟲患者兩種檢查方法

6、的結(jié)果2021/10/421dataex3;dor=1to2;doc=1to2;inputf@@;output;end;end;cards;112401236;2021/10/422procfreq;tablesr*c/chisqexpectednopercentnocolagree;weightf;run;此處用Agree選項(xiàng)實(shí)現(xiàn)SAS的McNemar檢驗(yàn)，但要注意，SAS的McNemar檢驗(yàn)不考慮校正問題，只適合于較大樣本（b+c>40），所得的結(jié)果與公式的計(jì)算結(jié)果相同。2021/10/423Dataex4;Inputf11f12f21f22;Iff12+f21<40t

7、henChisq=(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的時(shí)候，使用公式編輯法計(jì)算?2統(tǒng)計(jì)量及對(duì)應(yīng)的P值。2021/10/4242021/10/425此部分結(jié)果是普通四格表?2檢驗(yàn)的結(jié)果，不適于配對(duì)?2檢驗(yàn)使用。2021/10/4261960年Cohen等提出用Kappa值作為評(píng)價(jià)判斷的一致性程度的