SAS 2017/7/28出版的文件【SAS 9.4 新功能】中,介紹許多SAS 9.4M1到SAS 9.4M4的增強功能及調整。其中,在Proc Freq程序中,增強估計勝算比 (odds ratio, OR)的信賴區間。以下將透過簡單的範例玩玩這個功能。

首先,先以proc logistic的程序估計勝算比及勝算比的信賴區間,程式碼如下所示(Mortality=1為死亡):

 

proc logistic data=a.ex descending ;

class Gender(ref="1") ;

model mortality=Gender / rl ;

run ;

 

1.jpg

 

接著玩玩SAS 9.4之下,Proc Freq程序新功能估計勝算比,程式碼如下所示:

 

proc freq data=a.ex ; tables Gender*mortality

/ nopercent norow chisq BINOMIAL(CL=ALL) RELRISK OR(CL=SCORE) RISKDIFF(COMMON) ; run ;

 

SAS 9.4的Freq程序擴充了三種新功能,分別是BINOMIAL(CL=)、RELRISK(CL=SCORE)及RISKDIFF(COMMON)。

  1. BINOMIAL估計的信賴區間,可以依照使用者需求,選擇許多估算形式,例如:AC, AGRESTICOULL, ALL, CLOPPERPEARSON, CP, EXACT, J, JEFFREYS, SCORE, W, WALD, WALDC, WILSON, WILSONC (此範例中選擇ALL)。
  2. RELRISK是要求計算相關風險的指令,其後方需指名估算相關風險的形式,支援的形式有BIN, BINOMIAL, ODDSRATIO, OR, RISKDIFF (此範例中選擇OR)。
  3. RISKDIFF(COMMON)可以輸出分層的風險差異值以及Newcombe方式估算的信賴區間。

 

2.jpg

 

註1:勝算比跟信賴區間估計的小數位數多了一位,勝算比的數值位在表標題列處,加上報表小小一張,若是不仔細看,挺容易錯過的。

註2:在此僅貼出勝算比估計值跟信賴區間的報表,以利跟上方proc logistic程序的報表數值對照,其它報表本篇文章省略,有興趣的同好可以自行玩玩喔。

延伸閱讀

SAS Institute Inc. 2013. SAS 9.4 的新功能. Cary, NC: SAS Institute Inc.

What’s New in SAS/STAT 13.2 (SAS官方網站: http://support.sas.com/documentation/cdl/en/statug/67523/HTML/default/viewer.htm#statug_whatsnew_toc.htm)

 

關鍵字:SAS 9.4, proc freq, proc logistic, odds ratio, OR, 信賴區間, CL, confidence interval, confidence limits.

創作者介紹

晨晰統計部落格新站(統計、SPSS、BIG DATA討論園地)

晨晰部落格新站 發表在 痞客邦 留言(0) 人氣()