Updates to SAS® Power and Sample Size Software in SAS/STAT® 9.2

Link: http://www2.sas.com/proceedings/forum2008/368-2008.pdf

PROC POWER程序從SAS 9.1版發表出來後，讓許多power和樣本計算的過程不需要靠使用者自行寫程式才能完成。但一開始的功能並不完備，不過在SAS 9.2版釋出後，PROC POWER 內建的功能就逐漸被補齊了。Wayne Watson 於 SAS Global Forum 2008 發表了一份技術文件，整理了 SAS 9.2版的 PROC POWER 一些新增的功能。


(1) LOGISTIC REGRESSION
羅吉斯迴歸的 power analysis 在此已經可以進行了。使用前要先寫 logistic 告知 SAS 要跑這種模型的 power analysis，並把要檢定的預測變數、勝算比以及反應變數 Y=1 的機率都寫上去，重點是要把檢定的預測變數的分配設定好（藍色部分）。

proc power;
    logistic
        alpha = 0.05
        vardist(’Duration’) = normal(4, 1.5)
        testpredictor = ’Duration’
        testoddsratio = 1.7
        responseprob = 0.65
        ntotal = 50 60 70
        power = . ;
run;

(2) CONFIDENCE INTERVAL FOR ONE PROPORTION
二項變數的信賴區間也可以算了，但這結果其實跟用 PROC FREQ 算出來的結果是一樣的，所以如果你已經會用 PROC FREQ 來算二項變數的信賴區間的話，就不用多學這一段了。

proc power;
    onesamplefreq ci = Wilson
        alpha = 0.05
        proportion = 0.3
        halfwidth = 0.1
        ntotal = 70
        probwidth = .;
run;

(3) EQUIVALENCE TESTS FOR ONE PROPORTION
假設檢定中的比例檢定的 power analysis 也可以做了。重點就是要把紅色那段程式碼寫上去，其餘的參數設定跟以往都差不多。

proc power;
    onesamplefreq test = equiv_exact
        alpha = 0.05
        proportion = 0.35
        lower = 0.2
        upper = 0.4
        ntotal = 500
        power = .;
run;

(4) WILCOXON MANN-WHITNEY TEST FOR TWO INDEPENDENT GROUPS
無母數 Wilcoxon Mann-Whitney 檢定的 power analysis 可以用 twosamplewilcoxon 呼叫出。比較麻煩的是類別變數的分配要自己將每一層的 coding 和機率值打上去（藍色部分），注意後面機率值的總和要等於一。

proc power;
    twosamplewilcoxon
        alpha = 0.05
        vardist(’lidocaine’) =ordinal( (-3 -2 -1 0 1 2 3): (.01 .04 .20 .50 .20 .04 .01) )
        vardist(’mironel plus lidocaine’) =ordinal( (-3 -2 -1 0 1 2 3): (.01 .03 .15 .35 .30 .10 .06) )
        variables = ’lidocaine’ | ’mironel plus lidocaine’
        sides = u
        ntotal = .
        power = 0.85;
run;

最後原文有介紹一個SAS專門拿來做 Power analysis 的副產品Power and Sample Size 3.1，簡稱 PSS。不過這個軟體並沒有附在 SAS/STAT 裡面，需要另外安裝。我自己的 SAS 9.2 光碟片裡面好像也沒有這個附加軟體，研判應該是需要另外付費才有。我想應該也不會有很多人有 PSS，所以在此略過不提。把 PROC POWER 熟用應該就可以應付大部分的 Power analysis 以及相關繪圖。

CONTACT INFORMATION
Wayne Watson
Building S, Room 3040
SAS Institute, Inc.
SAS Campus Drive
Cary, NC 27513
Work Phone: 919-531-6770
E-mail: wayne.watson@sas.com