绘制ROC曲线.PDFVIP

SAS實用巨集簡介：繪製 ROC曲線和比較相關 ROC曲線下面積 中國醫藥大學 生物統計中心 統計分析師 林保萱 本篇文章主要介紹SAS官方網頁中關於 接受器操作特性曲線(receiver operating characteristic curve, ROC curve的巨集) (Macro)程式 ，程式 分別為「繪製ROC曲線」 和「以無母數方式比較相關 ROC曲線下面積」。「繪製 ROC曲線 」主要介紹如何可以 快速繪製出比 SAS說明檔中美觀且實用的 ROC曲線。 「以無母數方式比較相關 ROC 曲線下面積」主要介紹如何檢定相關 ROC曲線下面積有無達統計上顯著的差異。關於 ROC曲線的介紹請見生統中心生統 e報第 001 期的「統計學方法在醫學診斷上之應用 (1) 」。 本文中所提及之程式，僅為原作者的部分語法，詳細語法請見原文說明，兩篇原文 及 SAS巨集程式連結分別如下： Sample 25018: Plot ROC curve with labelled points for a binary-response model (/kb/25/018.html) rocplot.sas程式：/kb/25/addl/fusion25018_4_rocplot.sas.txt Sample 25017: Nonparametric comparison of areas under correlated ROC curves (/kb/25/017.html) roc.sas程式 ：/kb/25/addl/fusion_25017_6_roc.sas.txt 請先將程式下載至本機硬碟中，建議將 Sample 25017的程式取名為 roc.sas ，而 Sample 25018的程式取名為 rocplot.sas 。必須要注意的是，如果是想比較相關 ROC 曲線下面積 ，也就是roc.sas這個巨集程式，必須安裝 SAS/IML 。 繪製 ROC曲線 在 SAS程式中已經有介紹如何繪製出簡易的 ROC曲線，但在 我們諮詢的個案中比 較常被問到的問題為， ROC曲線 上的某個點其「原始數值為何」，利用rocplot.sas 便 可輕鬆完成這個工作。 假設有一組資料，研究者想以 年齡和某一個 BW指標來判斷罹患某疾病的機會。 問題：繪製以年齡作為判斷是否罹患此疾病的 ROC曲線。 第一步 ：先將資料讀入SAS中 data Age_data; input disease age bw@@; datalines; 0 50 65 0 39 61 0 21 70 0 61 67 0 30 55 0 35 63 0 25 72 0 41 66 0 43 52 0 36 54 0 37 76 0 25 61 1 2008/7 0 41 53 0 62 55 0 28 70 0 33 68 1 52 45 1 49 61 1 47 42 1 62 31 1 55 67 1 70 61 1 75 55 1 77 52 1 81 58 1 64 53 1 62 41 1 39 57 1 61 51 1 61 55 1 57 49 1 79 47 ; run; 第二步 ：將rocplot.sas程式預先讀入 SAS程式中 %include c:\temp\ROC\rocplot.sas; /* %include “程式所放的位置 ”; */ 透過%include這個巨集指令，可將 rocplot.sas呼叫至 SAS程式中。 第三步 ：執行proc logistic程式，並將結果輸出， 以供 rocplot.sas程式使用 proc logistic data=age_data; model disease (event=1)=age / outroc=roc1 roceps=0; output out=outp p=phat; run; 在這個程式中，我們將