sas卡方檢驗編程語句_第1頁
sas卡方檢驗編程語句_第2頁
sas卡方檢驗編程語句_第3頁
sas卡方檢驗編程語句_第4頁
sas卡方檢驗編程語句_第5頁
已閱讀5頁,還剩26頁未讀, 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領

文檔簡介

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

2、2 2 檢驗,可用基本檢驗,可用基本 公式或?qū)S霉焦交驅(qū)S霉?a、b、c、d為實際頻數(shù)為實際頻數(shù) 2 2 (ad bc) n (a b)(c d)(a c)(b d) 四格表四格表 2 2 值的校正公式值的校正公式 當當1T 5,且且n40時,需校正時,需校正 2 2 檢驗檢驗 a .基本公式的校正:基本公式的校正: b.專用公式的校正:專用公式的校正: 2 2 0.5AT T 2 2 2ad bc nn a b c d a c b d 3、行、行列表資料的列表資料的 2 2 檢驗檢驗 公式: 2 2 RC A n1 n n RC表資料的表資料的 2 2檢驗適用檢驗適用條件條件 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,認為內(nèi)科療法對兩種類型胃潰,認為內(nèi)科療法對兩種類型胃潰 瘍的治愈率差別有統(tǒng)計學意義,一般類瘍的治愈率差別有統(tǒng)計學意義,一般類 型的治愈率高于特殊型。型的治愈率高于特殊型。 例例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的檢驗水準拒受的檢驗水準拒受H H0 0,認為,認為 三地花生黃曲霉素三地花生黃曲霉素B B1 1污染率有差別。污染率有差別。 配對設計的 2 2檢驗(檢驗(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àn)選項實現(xiàn)SAS的的McNemar檢驗,檢驗, 但要注意,但要注意,SAS的的McNemar檢驗不考慮檢驗不考慮 校正問題,只適合于校正問題,只適合于較大樣本較大樣本(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檢驗的檢驗的 2 2 統(tǒng)計量及其顯著水平,即統(tǒng)計量及其顯著水平,即 2 2= =15.076915.0769,P P 0.00010.0001,按按 =0.05=0.05的檢驗水準拒受的檢驗水準拒受H H0 0, 認為兩種檢查方法有差別。認為兩種檢查方法有差別。 公式法得到的結(jié)

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論