sas卡方檢驗(yàn)編程語句

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

2、2 2 檢驗(yàn)，可用基本檢驗(yàn)，可用基本 公式或?qū)Ｓ霉焦交驅(qū)Ｓ霉?a、b、c、d為實(shí)際頻數(shù)為實(shí)際頻數(shù) 2 2 (ad bc) n (a b)(c d)(a c)(b d) 四格表四格表 2 2 值的校正公式值的校正公式 當(dāng)當(dāng)1T 5，且且n40時，需校正時，需校正 2 2 檢驗(yàn)檢驗(yàn) a .基本公式的校正：基本公式的校正： b.專用公式的校正：專用公式的校正： 2 2 0.5AT T 2 2 2ad bc nn a b c d a c b d 3、行、行列表資料的列表資料的 2 2 檢驗(yàn)檢驗(yàn) 公式： 2 2 RC A n1 n n RC表資料的表資料的 2 2檢驗(yàn)適用檢驗(yàn)適用條件條件 1.理論

3、數(shù)不能小于理論數(shù)不能小于1 2.理論數(shù)理論數(shù)140n40且各格子的期望值均大于且各格子的期望值均大于5 5， 因而選用因而選用Chi-SquareChi-Square的的 2 2統(tǒng)計量及其顯統(tǒng)計量及其顯 著性水平，即著性水平，即 2 239.92739.927，P=0.0001P=0.0001， 拒絕拒絕H H0 0，認(rèn)為內(nèi)科療法對兩種類型胃潰，認(rèn)為內(nèi)科療法對兩種類型胃潰 瘍的治愈率差別有統(tǒng)計學(xué)意義，一般類瘍的治愈率差別有統(tǒng)計學(xué)意義，一般類 型的治愈率高于特殊型。型的治愈率高于特殊型。 例例2 某省三地區(qū)花生黃曲霉素某省三地區(qū)花生黃曲霉素B1污染率比較污染率比較 地區(qū)地區(qū)未污染未污染污染污染合

4、計合計污染率（污染率（%） 甲甲 乙乙 6 30 8 23 14 3 29 44 11 79.3 31.8 27.3丙丙 合計合計44408447.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)輸入法：循環(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/

5、chisq expected nopercent nocol; weight count; run; 結(jié)果 本例各格子期望值均大于本例各格子期望值均大于5 5，選用，選用Chi-Chi- SquareSquare的的 2 2統(tǒng)計量及其顯著水平，即統(tǒng)計量及其顯著水平，即 2 2=17.907=17.907， P=0.0001P=0.0001，按按 =0.05=0.05的檢驗(yàn)水準(zhǔn)拒受的檢驗(yàn)水準(zhǔn)拒受H H0 0，認(rèn)為，認(rèn)為 三地花生黃曲霉素三地花生黃曲霉素B B1 1污染率有差別。污染率有差別。 配對設(shè)計的 2 2檢驗(yàn)（檢驗(yàn)（SASSAS程序）程序） 200名已確診的血吸蟲患者，治名已確診的血吸蟲患

6、者，治 療前經(jīng)皮試法及糞檢法檢查，結(jié)果療前經(jīng)皮試法及糞檢法檢查，結(jié)果 如下表，問兩種檢查方法的結(jié)果有如下表，問兩種檢查方法的結(jié)果有 無差別？無差別？ 糞檢法糞檢法合計合計 +- 皮試法皮試法+11240152 -123648 合計合計12476200 表表2 血吸蟲患者兩種檢查方法的結(jié)果血吸蟲患者兩種檢查方法的結(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

7、; weight f; run; 此處用此處用Agree選項(xiàng)實(shí)現(xiàn)選項(xiàng)實(shí)現(xiàn)SAS的的McNemar檢驗(yàn)，檢驗(yàn)， 但要注意，但要注意，SAS的的McNemar檢驗(yàn)不考慮檢驗(yàn)不考慮 校正問題，只適合于校正問題，只適合于較大樣本較大樣本（b+c40）， 所得的結(jié)果與所得的結(jié)果與 公式的計算結(jié)果公式的計算結(jié)果 相同。相同。 )( )( 2 2 cb cb 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+cKappa0.4兩者兩者 一致性一般；一致性一般；Kappa40b+c40，故選用，故選用McnemarMcnemar檢驗(yàn)的檢驗(yàn)的 2 2 統(tǒng)計量及其顯著水平，即統(tǒng)計量及其顯著水平，即 2 2= =15.076915.0769，P P 0.00010.0001，按按 =0.05=0.05的檢驗(yàn)水準(zhǔn)拒受的檢驗(yàn)水準(zhǔn)拒受H H0 0， 認(rèn)為兩種檢查方法有差別。認(rèn)為兩種檢查方法有差別。 公式法得到的結(jié)