1、sas 2 檢驗(yàn)編程語(yǔ)句,一、復(fù)習(xí)有關(guān)內(nèi)容,1.應(yīng)用： （1）分類(lèi)變量資料； （2）推斷兩個(gè)或兩個(gè)以上的樣本率或構(gòu)成比之間有無(wú)差異； （3）檢驗(yàn)頻數(shù)分布的擬合優(yōu)度。,2. 四格表資料的2 檢驗(yàn),（1）四格表的形式 （2）2 檢驗(yàn)的基本思想 （3） 基本公式 自由度 =(行數(shù)-1)(列數(shù)-1), 當(dāng)Ti5，且n 40時(shí) 使用普通2 檢驗(yàn)，可用基本 公式或?qū)Ｓ霉?a、b、c、d為實(shí)際頻數(shù),四格表2 值的校正公式,當(dāng)1T 5，且n40時(shí)，需校正2 檢驗(yàn) a .基本公式的校正： b.專(zhuān)用公式的校正：,3、行列表資料的2 檢驗(yàn),公式：,RC表資料的2檢驗(yàn)適用條件,1.理論數(shù)不能小于1 2.理論數(shù)1T，

2、且小于5的格子數(shù)不超過(guò)總格子數(shù)的1/5。 若條件不適合，需作如下處理： A 增大樣本例數(shù) B 刪除理論數(shù)太小的行或列 C 合并（性質(zhì)相同）,二、完全隨機(jī)設(shè)計(jì)的2檢驗(yàn)（SAS）,例1 某醫(yī)院用內(nèi)科療法治療一般類(lèi)型胃潰瘍患者80例，治愈63例；治療特殊類(lèi)型胃潰瘍患者99例，治愈31例。 問(wèn)內(nèi)科療法對(duì)兩種類(lèi)型胃潰瘍的治愈率差別有無(wú)顯著意義？,表 兩種類(lèi)型胃潰瘍內(nèi)科療法治療結(jié)果,data ex1; input r c count ; cards; 1 1 63 1 2 17 2 1 31 2 2 68 3 1 322 3 2 ; proc freq; tables r*c/chisq expected

3、 nopercent nocol; weight count; run;,結(jié)果,結(jié)果解釋,本例n40且各格子的期望值均大于5，因而選用Chi-Square的2統(tǒng)計(jì)量及其顯著性水平，即239.927，P=0.0001，拒絕H0，認(rèn)為內(nèi)科療法對(duì)兩種類(lèi)型胃潰瘍的治愈率差別有統(tǒng)計(jì)學(xué)意義，一般類(lèi)型的治愈率高于特殊型。,例2 某省三地區(qū)花生黃曲霉素B1污染率比較,地區(qū),未污染,污染,合計(jì),污染率（%）,甲 乙,6 30 8,23 14 3,29 44 11,79.3 31.8 27.3,丙,合計(jì),44,40,84,47.6,程序1：一般輸入方法,data ex2; input r c count; car

4、ds; 1 1 6 1 2 23 2 1 30 2 2 14 3 1 8 3 2 3 ;,程序2：循環(huán)輸入法,data ex2; do r=1 to 3; do c=1 to 2; input count; output;end;end; cards; 6 23 30 14 8 3 ;,proc freq; tables r*c/chisq expected nopercent nocol; weight count; run;,結(jié)果,本例各格子期望值均大于5，選用Chi-Square的2統(tǒng)計(jì)量及其顯著水平，即2=17.907，P=0.0001，按=0.05的檢驗(yàn)水準(zhǔn)拒受H0，認(rèn)為三地花生黃曲

5、霉素B1污染率有差別。,配對(duì)設(shè)計(jì)的2檢驗(yàn)（SAS程序）,200名已確診的血吸蟲(chóng)患者，治療前經(jīng)皮試法及糞檢法檢查，結(jié)果如下表，問(wèn)兩種檢查方法的結(jié)果有無(wú)差別？,表2 血吸蟲(chóng)患者兩種檢查方法的結(jié)果,data ex3; do r=1 to 2; do c=1 to 2; input f; output;end;end; cards; 112 40 12 36 ;,proc freq; tables r*c/chisq expected nopercent nocol agree; weight f; run;,此處用Agree選項(xiàng)實(shí)現(xiàn)SAS的McNemar檢驗(yàn)，但要注意，SAS的McNemar檢驗(yàn)不考

6、慮 校正問(wèn)題，只適合于較大樣本（b+c40），所得的結(jié)果與 公式的計(jì)算結(jié)果相同。,Data ex4; Input f11 f12 f21 f22; If f12+f2140 then Chisq=(abs(f12-f21)-1)*2/(f12+f21); Else chisq=abs(f12-f21)*2/(f12+f21); P=1-probchi(chisq,1); Cards; 112 40 12 36 ;proc print; Run;,此法適用于b+c=40的時(shí)候， 使用公式編輯法計(jì)算2統(tǒng)計(jì)量及對(duì)應(yīng)的P值。,此部分結(jié)果是普通四格表2檢驗(yàn)的結(jié)果， 不適于配對(duì)2檢驗(yàn)使用。,本例b+c40，故選用Mcnemar檢驗(yàn)的2統(tǒng)計(jì)量及其顯著水平，即2=15.0769，P0.0001