各位學(xué)友大家好，醫(yī)院的數(shù)據(jù)科學(xué)死氣沉沉，2025春節(jié)假期不遠(yuǎn)了，很多外院的老師問我，時(shí)間依賴性AUC我會(huì)做了，怎么都是出來一堆值??？有沒有綜合的指標(biāo)，綜合的指標(biāo)是有的，叫integrated AUC，簡稱iAUC。

iAUC 的計(jì)算是基于時(shí)間依賴性 AUC 的結(jié)果。具體做法是對(duì)每個(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)平均值

# 對(duì)應(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)平均：我們對(duì)每個(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)。