各位學(xué)友大家好,醫(yī)院的數(shù)據(jù)科學(xué)死氣沉沉,2025春節(jié)假期不遠(yuǎn)了,很多外院的老師問我,時(shí)間依賴性AUC我會(huì)做了,怎么都是出來一堆值???有沒有綜合的指標(biāo),綜合的指標(biāo)是有的,叫integrated AUC,簡稱iAUC。 iAUC 的計(jì)算是基于時(shí)間依賴性 AUC 的結(jié)果。具體做法是對每個(gè)時(shí)間段的 AUC 進(jìn)行加權(quán)平均,權(quán)重是相鄰時(shí)間點(diǎn)之間的時(shí)間差。 # 獲取AUC值和時(shí)間點(diǎn) auc_values <- time_roc_result$AUC # 獲取AUC time_points <- times # 獲取時(shí)間點(diǎn) # 計(jì)算時(shí)間間隔 time_intervals <- diff(time_points) # 計(jì)算每個(gè)時(shí)間段的AUC加權(quán)平均值 # 對應(yīng)時(shí)間段的AUC為兩個(gè)相鄰時(shí)間點(diǎn)的AUC平均值 weighted_auc <- sapply(1:(length(auc_values) - 1), function(i) { (auc_values[i] + auc_values[i + 1]) / 2 * time_intervals[i] }) # 計(jì)算iAUC iAUC <- sum(weighted_auc) / (max(time_points) - min(time_points)) print(iAUC) 解釋: AUC 值:我們從 timeROC 的輸出中提取每個(gè)時(shí)間點(diǎn)的 AUC 值(time_roc_result$AUC)。 時(shí)間間隔:通過 diff() 函數(shù)計(jì)算相鄰時(shí)間點(diǎn)之間的時(shí)間差。 加權(quán)平均:我們對每個(gè)時(shí)間段的 AUC 取相鄰時(shí)間點(diǎn)的 AUC 平均值,并乘以時(shí)間間隔,得到該時(shí)間段的加權(quán) AUC。 iAUC:將所有時(shí)間段的加權(quán) AUC 加和,再除以整個(gè)時(shí)間段的長度,得到 iAUC。 這種方法可以幫助你總結(jié)模型在整個(gè)時(shí)間范圍內(nèi)的表現(xiàn)。 |
|