若要在羅吉斯迴歸中計算C統計量(C Statistic),可以使用R軟體進行,此處的C統計量即為ROC曲線下的面積。
羅吉斯迴歸分析(Logistic regression, logit model)相關說明可參考:
http://www.yongxi-stat.com/logistic-regression/
ROC曲線法(ROC curve)相關說明可參考:
http://www.yongxi-stat.com/receiver-operating-characteristic-roc/
一、使用R計算C統計量的步驟
(一) 首先大家可自行安裝R軟體,建議可以安裝操作介面較為友善的R studio免費版。
(二) 接著點擊右下角Packages中的Install,並安裝「DescTools」套件。
(三) 安裝完成後,以「Library(“DescTools”)」指令呼叫上述套件。
(四) 接著匯入所欲分析的檔案。
(五) 使用「glm()」函數進行羅吉斯迴歸:r.glm <- glm(plate_exposure ~ Postop_radiation + ALT_M + Chimeric, data=final_data1029, family=binomial)
(六) 使用剛剛安裝好的套件計算C統計量:Cstat(r.glm)即可得出結果。
(七)若是上述函數有報錯情形,也可使用另一個函數:Cstat(x = predict(r.glm, method=”response”), resp = model.response(model.frame(r.glm))) 亦可得出相同結果。
【參考文獻】
https://rdrr.io/cran/DescTools/man/Cstat.html