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

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 ;

接著玩玩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方式估算的信賴區間。

​​

註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.