Graphpad繪制ROC曲線
在生物標志物模型分析中姆蘸,還有一種常見的分析為受試者工作特征(ROC)曲線分析怨喘。如下圖所示厅贪,ROC曲線也是這類文獻中常見的分析圖灶平,結(jié)合臨床數(shù)據(jù)伺通,可以驗證某個基因或模型作為疾病診斷和預(yù)后標志物。
那么如何使用Graphpad繪制ROC曲線呢逢享?我們就以下圖為例罐监,來一步步演示吧。
1? 數(shù)據(jù)準備
繪制ROC曲線需要生存狀態(tài)瞒爬、基因表達量弓柱,將數(shù)據(jù)整理成表格。示例文件示例數(shù)據(jù).xlsx所示疮鲫。
(1)patient:患者編號吆你;(2)status:生存狀態(tài),其中0表示存活俊犯,1表示死亡;(3)expression:基因的表達量(或模型的風(fēng)險評分)伤哺。
2? 使用Graphpad進行ROC分析
接下來燕侠,打開Graphpad,選擇Column立莉,點擊Create绢彤。
將Excel中ROC中的數(shù)據(jù)復(fù)制到Graphpad中。
點擊Analysis中的Analyze蜓耻,在Column analysis中選擇ROC Curve茫舶,點擊OK。
在Parameters:ROC Curve中選中Line of identity(是否繪制參考線刹淌,AUC=0.5)饶氏,點擊OK讥耗。
點擊Results中ROC of Data 1中的Area,其中Area under the ROC curve下的Area值就是我們需要的AUC值疹启。
3? 繪制ROC曲線
點擊Graphs中的ROC curve: ROC of Data 1古程。
雙擊坐標軸,修改相關(guān)參數(shù)喊崖。點擊Frame and Origin挣磨,在Size中選擇Square(長寬相等),同時修改Width和Height荤懂。在Frame style中選擇Plain Frame(加邊框)茁裙。
點擊X axis,取消選擇Automatically determine the range and interval节仿。修改Minimum(x軸最小值)晤锥、Maximum(x軸最大值)和Major ticks interval(坐標間隔)。點擊Left Y axis粟耻,修改相同參數(shù)查近。
雙擊繪圖區(qū)域,在Data set中選擇ROC of Data 1:ROC curve: Sensitivity%挤忙,取消選擇Show symbols(顯示關(guān)鍵點)霜威、Color(顏色)、Style(類型改成折線)册烈、Show legend(圖例)等參數(shù)戈泼;將Data set修改成ROC of Data 1:ROC curve:Indentity%,修改相同的參數(shù)赏僧。
修改x軸和y軸標題大猛,并將圖例調(diào)整到合適的位置,ROC圖就繪制完成了淀零。
結(jié)果還是非常漂亮的對不對挽绩?