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

上單元我們已實例示範RCS在存活分析及邏輯斯迴歸分析的應用，本單位我們繼續介紹RCS在線性迴歸分析的應用

範例三、分析醫師的執刀經驗(累積執刀數目)情形與病患住院日數的相關性

1. 屬於線性迴歸分析，R軟體使用ols function。
2. 自變項為執刀經驗(累積執刀數目)數值：cum_op。
3. 依變項為病患住院日數：hospital_stays。
4. 使用之非線性模擬模式為rcs function，knots設定3（建議在3-5間）。

語法如下

線性迴歸通常不需指定參考點（參考點設定適用於呈現相對的風險時例如存活分析的HR或邏輯斯迴歸分析的OR等），呈現自變項（原始值）與依變項（原始值）的相關性。可以進行下一步鍵入繪圖指令輸出下圖

RCS分析完成後，我們再繼續分享圖形的製作與美化，步驟如下

(1)匯出RCS分析結果為DTA報表，匯出檔案語法如下，以此指定匯出的路徑、檔名及檔案類型（PS. dta形式的檔案是可以用SPSS軟體開啟的喔）。

(2)DTA報表後續轉檔為EXCEL檔案進行後續繪圖。我們以範例一（Creatinine上升情形與病患一年死亡）的結果為例，匯出之EXCEL檔案報表如下，報表會依照Creatinine數值range切分成200組，列出Creatinine及其對應的HR及95%上下限。

簡易繪圖版本為繪製3條折線圖

進階繪圖版本可按右鍵點選「變更圖表類型」，修改上下限的圖表類型為區域圖，上限之填滿色彩為灰色；下限之填滿色彩為白色。

按右鍵點選「選取資料」，呈現順序要先上限(upper)再下限(lower)，即為以下限覆蓋上限之意

即可呈現美化之相關性圖型如下

以上的說明與實例示範希望對大家分析有所助益，預祝大家分析順利，報表順利產出 !

參考資源

RMS package: https://www.rdocumentation.org/packages/rms/versions/6.0-1