
出品 | CDA數(shù)據(jù)分析研究院
數(shù)據(jù)科學(xué)是一個(gè)在過(guò)去幾年里真正成長(zhǎng)的術(shù)語(yǔ),似乎每個(gè)人都想加入。最具吸引力的目標(biāo)之一是利用數(shù)據(jù)資產(chǎn)的力量來(lái)創(chuàng)建能夠預(yù)測(cè)各種結(jié)果的機(jī)器學(xué)習(xí)模型。通過(guò)定義明確的模型,可以確定能夠預(yù)測(cè)結(jié)果的最重要影響因素,為戰(zhàn)略假設(shè)開(kāi)發(fā)有價(jià)值的洞察力,甚至可以通過(guò)友好的用戶界面將模型的邏輯實(shí)現(xiàn)到軟件應(yīng)用程序中。
然而,在任何這種魔法發(fā)生之前,我們需要知道模型創(chuàng)建的預(yù)測(cè)是否有益!例如,如果我們的電子郵件程序的垃圾郵件分類器只能檢測(cè)到50%的不需要的電子郵件或請(qǐng)求,我們都會(huì)非常憤怒。在這篇文章中,我將討論如何使用每個(gè)數(shù)據(jù)科學(xué)家應(yīng)該熟悉的經(jīng)典工具來(lái)評(píng)估預(yù)測(cè)模型:接收器操作特性(ROC)曲線。
說(shuō)明性實(shí)例:預(yù)測(cè)冠狀動(dòng)脈疾病
我在這篇文章中使用的數(shù)據(jù)集來(lái)自Bom等人對(duì)“疑似冠狀動(dòng)脈疾病患者冠狀動(dòng)脈斑塊形態(tài)的靶向蛋白質(zhì)組學(xué)的預(yù)測(cè)價(jià)值”。
這項(xiàng)研究檢查了兩種不同的結(jié)果,我將重點(diǎn)關(guān)注這篇文章的是沒(méi)有冠狀動(dòng)脈疾?。–AD)。作者評(píng)估了蛋白質(zhì)組學(xué)生物標(biāo)志物在有癥狀患者中檢測(cè)CAD缺失的預(yù)測(cè)能力。確定一組可以區(qū)分無(wú)CAD患者和需要立即干預(yù)的蛋白質(zhì)的蛋白質(zhì)組將提供更準(zhǔn)確和成本有效的非侵入性診斷測(cè)試。
該數(shù)據(jù)集是測(cè)試多個(gè)數(shù)據(jù)科學(xué)主題的絕佳來(lái)源。通過(guò)少量觀察,它很容易使用,但它還包含大量變量以增加復(fù)雜性。結(jié)果(使用冠狀動(dòng)脈計(jì)算機(jī)斷層掃描血管造影術(shù),或CCTA)的可行“黃金標(biāo)準(zhǔn)”被提供來(lái)測(cè)試預(yù)測(cè)。對(duì)于這篇文章,我將主要關(guān)注從單一模型構(gòu)建ROC曲線,但可以在后面的帖子中使用此數(shù)據(jù)集深入研究更高級(jí)的主題。
預(yù)測(cè)精度測(cè)量概述
在開(kāi)始創(chuàng)建曲線之前,了解一些評(píng)估預(yù)測(cè)準(zhǔn)確性的常用指標(biāo)非常重要。
積極類:我將正類定義為我試圖檢測(cè)的結(jié)果類。在這種情況下,缺少 CAD。雖然我意識(shí)到這可能會(huì)導(dǎo)致診斷測(cè)試領(lǐng)域中的“積極”和“消極”這兩個(gè)術(shù)語(yǔ)混淆,但這種方式定義積極類更能推廣到其他情況。
在這種情況下處理有癥狀患者的默認(rèn)方法是讓他們接受進(jìn)一步的測(cè)試和程序。通過(guò)充分檢測(cè)沒(méi)有CAD的患者,我們將不再需要不必要的,更具侵入性的手術(shù)。
一般準(zhǔn)確性:簡(jiǎn)單地說(shuō),有多少科目正確分類?
敏感性:正確識(shí)別真陽(yáng)性的比例。在這種情況下,健康患者的比例由診斷工具正確識(shí)別。這有時(shí)被稱為“召回”。
特異性:正確識(shí)別真陰性的比例。在這種情況下,通過(guò)診斷工具正確識(shí)別CAD患者的比例。
積極預(yù)測(cè)價(jià)值:該工具報(bào)告的陽(yáng)性比例,實(shí)際上是積極的。對(duì)于診斷工具報(bào)告缺乏CAD的患者組,PPV是實(shí)際上沒(méi)有患病的患者的比例。這有時(shí)被稱為“精確度”。
負(fù)面預(yù)測(cè)值:該工具報(bào)告的負(fù)面影響的比例,實(shí)際上是負(fù)面的。對(duì)于診斷工具報(bào)告存在CAD的患者組,NPV是實(shí)際上沒(méi)有CAD的患者的比例。
不要擔(dān)心,如果你的大腦試圖理清所有這些指標(biāo)。這些術(shù)語(yǔ)很容易混淆,特別是在第一次學(xué)習(xí)它們時(shí)。我發(fā)現(xiàn)在2x2表中可視化內(nèi)容最容易。
二進(jìn)制類的混淆矩陣(正=沒(méi)有CAD,負(fù)= CAD)
曲線的解剖
ROC曲線的優(yōu)點(diǎn)在于它是一種易于理解的圖形工具,可應(yīng)用于您創(chuàng)建的任何預(yù)測(cè)模型。以下是曲線的基礎(chǔ)知識(shí):
軸:靈敏度和誤報(bào)率
首先,我們需要為情節(jié)創(chuàng)造空間。通過(guò)將靈敏度與1-特異性(或假陽(yáng)性率)作圖來(lái)建立ROC曲線。
預(yù)測(cè)概率
現(xiàn)在我們需要繪制一些東西?;叵胍幌拢A(yù)測(cè)模型會(huì)將每個(gè)觀察分配給最可能的類(在這種情況下,沒(méi)有CAD與CAD的存在)。模型實(shí)際上在做的是計(jì)算屬于特定類的概率。在0和1之間選擇截止值,并且如果計(jì)算的概率超過(guò)該閾值,則將觀察值分配給該類。
您可能會(huì)發(fā)現(xiàn)在大多數(shù)軟件包中,默認(rèn)截止值設(shè)置為0.5,邏輯是,在二進(jìn)制類中,可以將觀察值分配給最可能的類。但是,正如我們所看到的,最好在考慮靈敏度和特異性之間的權(quán)衡之后選擇截止值。
通過(guò)繪制所有可能的截止值來(lái)生成ROC曲線,所述截止值是分配給每個(gè)觀察的概率。選擇不同的截止值將改變預(yù)測(cè)工具的靈敏度和特異性,因此,可以使用相關(guān)的靈敏度和1-特異性作為坐標(biāo)在圖的空間中繪制每個(gè)截止概率。最靠近左上角的點(diǎn)(SN = 1,F(xiàn)PR = 0)提供精度度量之間的最大平衡。
曲線下面積(AUC)
AUC是一個(gè)類似于二元模型的一致性或c統(tǒng)計(jì)量的度量。這是具有正類的觀察將具有比在負(fù)類中的觀察更大的預(yù)測(cè)概率的概率。
如果AUC = 1,則表示模型有完美的預(yù)測(cè)。如果AUC = 0.5,則意味著模型無(wú)法區(qū)分類。它在邏輯回歸中的行為與r平方相似,因?yàn)樘砑痈囝A(yù)測(cè)因子會(huì)增加AUC。因此,在分析中包括對(duì)外部數(shù)據(jù)的交叉驗(yàn)證或驗(yàn)證非常重要。
AUC可用于評(píng)估不同的預(yù)測(cè)模型。
參考線
繪制參考線通常是個(gè)好主意,其中AOC = 0.5在ROC圖上。這提供了比較曲線的基線視覺(jué)效果。
生成ROC曲線
雖然可以使用幾個(gè)不同的程序來(lái)開(kāi)發(fā)預(yù)測(cè)模型和ROC曲線,但我已經(jīng)在R中實(shí)現(xiàn)了分析。
為了預(yù)測(cè)結(jié)果,我使用xgboost創(chuàng)建了一個(gè)分類模型,其中蛋白質(zhì)生物標(biāo)記物作為預(yù)測(cè)因子,并使用重新取代方法來(lái)預(yù)測(cè)數(shù)據(jù)集中的結(jié)果。使用預(yù)測(cè)的概率,我可以形成上面顯示的ROC曲線。
注意:該程序也設(shè)置為使用k-fold交叉驗(yàn)證,或者可以更新以利用外部驗(yàn)證,例如列車和測(cè)試集的70/30分割(盡管我想要更多數(shù)據(jù))。
解釋曲線
AUC為0.738,實(shí)際上低于作者用他們的機(jī)器學(xué)習(xí)方法所發(fā)現(xiàn)的(堆疊可能提高了分類器的準(zhǔn)確性)。但是,我更關(guān)心的是選擇合適的截止值并權(quán)衡工具的成本和收益。
模型的總體準(zhǔn)確性是不錯(cuò)的,但是當(dāng)發(fā)生假陽(yáng)性結(jié)果時(shí),會(huì)發(fā)生CAD患者而不是通過(guò)必要的程序?讓我們同意這不會(huì)是一件好事。
為了安全起見(jiàn),盡可能減少診斷工具的誤報(bào)率可能是個(gè)好主意。在生成混淆矩陣時(shí),很容易選擇保持低誤報(bào)率的概率截止值。這將作為誰(shuí)被標(biāo)記為無(wú)CAD的決定標(biāo)準(zhǔn)。
通過(guò)選擇保持90%特異性(10%假陽(yáng)性率)的概率截止值,分類工具能夠檢測(cè)到42%沒(méi)有CAD的患者(敏感性= 0.42)。在這種情況下,這實(shí)際上是一個(gè)非常好的結(jié)果,因?yàn)楝F(xiàn)狀是讓所有有癥狀的患者接受更具侵入性的診斷程序。盡管42%的靈敏度可能聽(tīng)起來(lái)很低,但在這種特定情況下,這是一個(gè)很有希望的結(jié)果。
在開(kāi)發(fā)用于預(yù)測(cè)有癥狀患者的CAD狀態(tài)的分類器時(shí),最重要的是在保持低誤報(bào)率的同時(shí)實(shí)現(xiàn)盡可能高的靈敏度。
ROC曲線的優(yōu)點(diǎn)在于您可以從單個(gè)圖像中可視化所有這些性能指標(biāo)。將曲線與競(jìng)爭(zhēng)模型進(jìn)行比較是一種快速簡(jiǎn)便的方法,可以選擇合適的分類或診斷工具。
我只研究了使用ROC曲線和二進(jìn)制分類工具,但它可以很容易地?cái)U(kuò)展到多類方案。在同一圖形上,可以繪制與每個(gè)特定類別中的概率相對(duì)應(yīng)的多條曲線。
在撰寫(xiě)本文時(shí),我詳細(xì)介紹了如何構(gòu)建ROC曲線,目的是增加對(duì)預(yù)測(cè)精度度量的理解,以及如何評(píng)估您可能正在構(gòu)建的這些數(shù)據(jù)科學(xué)分類工具。
標(biāo)題:Illustrating Predictive Models with the ROC Curve
數(shù)據(jù)分析咨詢請(qǐng)掃描二維碼
若不方便掃碼,搜微信號(hào):CDAshujufenxi
DSGE 模型中的 Et:理性預(yù)期算子的內(nèi)涵、作用與應(yīng)用解析 動(dòng)態(tài)隨機(jī)一般均衡(Dynamic Stochastic General Equilibrium, DSGE)模 ...
2025-09-17Python 提取 TIF 中地名的完整指南 一、先明確:TIF 中的地名有哪兩種存在形式? 在開(kāi)始提取前,需先判斷 TIF 文件的類型 —— ...
2025-09-17CDA 數(shù)據(jù)分析師:解鎖表結(jié)構(gòu)數(shù)據(jù)特征價(jià)值的專業(yè)核心 表結(jié)構(gòu)數(shù)據(jù)(以 “行 - 列” 規(guī)范存儲(chǔ)的結(jié)構(gòu)化數(shù)據(jù),如數(shù)據(jù)庫(kù)表、Excel 表、 ...
2025-09-17Excel 導(dǎo)入數(shù)據(jù)含缺失值?詳解 dropna 函數(shù)的功能與實(shí)戰(zhàn)應(yīng)用 在用 Python(如 pandas 庫(kù))處理 Excel 數(shù)據(jù)時(shí),“缺失值” 是高頻 ...
2025-09-16深入解析卡方檢驗(yàn)與 t 檢驗(yàn):差異、適用場(chǎng)景與實(shí)踐應(yīng)用 在數(shù)據(jù)分析與統(tǒng)計(jì)學(xué)領(lǐng)域,假設(shè)檢驗(yàn)是驗(yàn)證研究假設(shè)、判斷數(shù)據(jù)差異是否 “ ...
2025-09-16CDA 數(shù)據(jù)分析師:掌控表格結(jié)構(gòu)數(shù)據(jù)全功能周期的專業(yè)操盤(pán)手 表格結(jié)構(gòu)數(shù)據(jù)(以 “行 - 列” 存儲(chǔ)的結(jié)構(gòu)化數(shù)據(jù),如 Excel 表、數(shù)據(jù) ...
2025-09-16MySQL 執(zhí)行計(jì)劃中 rows 數(shù)量的準(zhǔn)確性解析:原理、影響因素與優(yōu)化 在 MySQL SQL 調(diào)優(yōu)中,EXPLAIN執(zhí)行計(jì)劃是核心工具,而其中的row ...
2025-09-15解析 Python 中 Response 對(duì)象的 text 與 content:區(qū)別、場(chǎng)景與實(shí)踐指南 在 Python 進(jìn)行 HTTP 網(wǎng)絡(luò)請(qǐng)求開(kāi)發(fā)時(shí)(如使用requests ...
2025-09-15CDA 數(shù)據(jù)分析師:激活表格結(jié)構(gòu)數(shù)據(jù)價(jià)值的核心操盤(pán)手 表格結(jié)構(gòu)數(shù)據(jù)(如 Excel 表格、數(shù)據(jù)庫(kù)表)是企業(yè)最基礎(chǔ)、最核心的數(shù)據(jù)形態(tài) ...
2025-09-15Python HTTP 請(qǐng)求工具對(duì)比:urllib.request 與 requests 的核心差異與選擇指南 在 Python 處理 HTTP 請(qǐng)求(如接口調(diào)用、數(shù)據(jù)爬取 ...
2025-09-12解決 pd.read_csv 讀取長(zhǎng)浮點(diǎn)數(shù)據(jù)的科學(xué)計(jì)數(shù)法問(wèn)題 為幫助 Python 數(shù)據(jù)從業(yè)者解決pd.read_csv讀取長(zhǎng)浮點(diǎn)數(shù)據(jù)時(shí)的科學(xué)計(jì)數(shù)法問(wèn)題 ...
2025-09-12CDA 數(shù)據(jù)分析師:業(yè)務(wù)數(shù)據(jù)分析步驟的落地者與價(jià)值優(yōu)化者 業(yè)務(wù)數(shù)據(jù)分析是企業(yè)解決日常運(yùn)營(yíng)問(wèn)題、提升執(zhí)行效率的核心手段,其價(jià)值 ...
2025-09-12用 SQL 驗(yàn)證業(yè)務(wù)邏輯:從規(guī)則拆解到數(shù)據(jù)把關(guān)的實(shí)戰(zhàn)指南 在業(yè)務(wù)系統(tǒng)落地過(guò)程中,“業(yè)務(wù)邏輯” 是連接 “需求設(shè)計(jì)” 與 “用戶體驗(yàn) ...
2025-09-11塔吉特百貨孕婦營(yíng)銷案例:數(shù)據(jù)驅(qū)動(dòng)下的精準(zhǔn)零售革命與啟示 在零售行業(yè) “流量紅利見(jiàn)頂” 的當(dāng)下,精準(zhǔn)營(yíng)銷成為企業(yè)突圍的核心方 ...
2025-09-11CDA 數(shù)據(jù)分析師與戰(zhàn)略 / 業(yè)務(wù)數(shù)據(jù)分析:概念辨析與協(xié)同價(jià)值 在數(shù)據(jù)驅(qū)動(dòng)決策的體系中,“戰(zhàn)略數(shù)據(jù)分析”“業(yè)務(wù)數(shù)據(jù)分析” 是企業(yè) ...
2025-09-11Excel 數(shù)據(jù)聚類分析:從操作實(shí)踐到業(yè)務(wù)價(jià)值挖掘 在數(shù)據(jù)分析場(chǎng)景中,聚類分析作為 “無(wú)監(jiān)督分組” 的核心工具,能從雜亂數(shù)據(jù)中挖 ...
2025-09-10統(tǒng)計(jì)模型的核心目的:從數(shù)據(jù)解讀到?jīng)Q策支撐的價(jià)值導(dǎo)向 統(tǒng)計(jì)模型作為數(shù)據(jù)分析的核心工具,并非簡(jiǎn)單的 “公式堆砌”,而是圍繞特定 ...
2025-09-10CDA 數(shù)據(jù)分析師:商業(yè)數(shù)據(jù)分析實(shí)踐的落地者與價(jià)值創(chuàng)造者 商業(yè)數(shù)據(jù)分析的價(jià)值,最終要在 “實(shí)踐” 中體現(xiàn) —— 脫離業(yè)務(wù)場(chǎng)景的分 ...
2025-09-10機(jī)器學(xué)習(xí)解決實(shí)際問(wèn)題的核心關(guān)鍵:從業(yè)務(wù)到落地的全流程解析 在人工智能技術(shù)落地的浪潮中,機(jī)器學(xué)習(xí)作為核心工具,已廣泛應(yīng)用于 ...
2025-09-09SPSS 編碼狀態(tài)區(qū)域中 Unicode 的功能與價(jià)值解析 在 SPSS(Statistical Product and Service Solutions,統(tǒng)計(jì)產(chǎn)品與服務(wù)解決方案 ...
2025-09-09