sas卡方檢驗編程語句

1、sas χ2 檢驗編程語句,一、復(fù)習(xí)有關(guān)內(nèi)容,1.應(yīng)用：（1）分類變量資料；（2）推斷兩個或兩個以上的樣本率或構(gòu)成比之間有無差異；（3）檢驗頻數(shù)分布的擬合優(yōu)度。,2. 四格表資料的χ2 檢驗,（1）四格表的形式（2）χ2 檢驗的基本思想（3） 基本公式　　　　　　　　 自由度 ?=(行數(shù)-1)(列數(shù)-1),① 當Ti≥5，且n ≥40時 使用普通?2 檢驗，可用基本 公式

2、或?qū)Ｓ霉?a、b、c、d為實際頻數(shù),②　四格表?2 值的校正公式,當1≤T < 5，且n≥40時，需校正?2 檢驗a .基本公式的校正：b.專用公式的校正：,3、行×列表資料的?2 檢驗,公式：,R×C表資料的?2檢驗適用條件,1.理論數(shù)不能小于12.理論數(shù)1<T，且小于5的格子數(shù)不超過總格子數(shù)的1/5。若條件不適合，需作如下處理： A 增大樣本例數(shù) B 刪除理論數(shù)太小的行或

3、列 C 合并（性質(zhì)相同）,,二、完全隨機設(shè)計的?2檢驗（SAS）,例1 某醫(yī)院用內(nèi)科療法治療一般類型胃潰瘍患者80例，治愈63例；治療特殊類型胃潰瘍患者99例，治愈31例。 問內(nè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;

4、proc freq;tables r*c/chisq expected nopercent nocol;weight count;run;,結(jié)果,結(jié)果解釋,本例n>40且各格子的期望值均大于5，因而選用Chi-Square的?2統(tǒng)計量及其顯著性水平，即?2＝39.927，P=0.0001，拒絕H0，認為內(nèi)科療法對兩種類型胃潰瘍的治愈率差別有統(tǒng)計學(xué)意義，一般類型的治愈率高于特殊型。,例2 某省三地區(qū)花生黃曲霉素B1污染率比較

5、,,,,地區(qū),未污染,污染,合計,污染率（%）,甲乙,6308,23 143,29 4411,79.331.827.3,丙,,合計,44,40,84,47.6,,程序1：一般輸入方法,data ex2;input r c count@@;cards;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

6、 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)計量及其顯著水平，即?2=17.907，P=0.0001，按?=0.05的檢驗水準拒受H0，認為三地花

7、生黃曲霉素B1污染率有差別。,,配對設(shè)計的?2檢驗（SAS程序）,200名已確診的血吸蟲患者，治療前經(jīng)皮試法及糞檢法檢查，結(jié)果如下表，問兩種檢查方法的結(jié)果有無差別？,表2 血吸蟲患者兩種檢查方法的結(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 exp

8、ected nopercent　 nocol agree;weight f;run;,此處用Agree選項實現(xiàn)SAS的McNemar檢驗，但要注意，SAS的McNemar檢驗不考慮校正問題，只適合于較大樣本（b+c>40），所得的結(jié)果與 公式的計算結(jié)果相同。,Data ex4; Input f11 f12 f21 f22; If f12+f21<40 then

9、 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的時候，使用公式編輯法計算?2統(tǒng)計量及對應(yīng)的P值。,此部分結(jié)果是普通四格表?2檢驗的結(jié)果，不適于配對?2檢