
SPSS操作:多項(xiàng)測(cè)量指標(biāo)的ROC曲線分析
在前面幾講中,我們已向大家介紹過如何依據(jù)一項(xiàng)測(cè)量指標(biāo)設(shè)計(jì)和評(píng)價(jià)診斷試驗(yàn)。但在實(shí)際臨床工作中,我們往往是根據(jù)多項(xiàng)指標(biāo)綜合判斷病情的。比如,在診斷高血壓時(shí),我們會(huì)根據(jù)年齡、性別、并發(fā)癥等多個(gè)因素同時(shí)評(píng)價(jià)受試者的患病情況。那么在這種情況下,我們應(yīng)如何判斷診斷結(jié)果的真實(shí)性呢?
一、問題與數(shù)據(jù)
某呼吸內(nèi)科醫(yī)生擬通過性別、年齡、BMI、COPD病史和是否吸煙等因素預(yù)測(cè)受試者的肺癌患病情況。他招募了85名肺癌患者,259名非肺癌患者,并通過查閱病歷、問卷調(diào)查的方式收集了上述信息。變量的賦值和部分原始數(shù)據(jù)見表1和表2。如果該醫(yī)生依據(jù)這幾項(xiàng)因素預(yù)測(cè)受試者是否患肺癌,那么應(yīng)如何預(yù)測(cè),準(zhǔn)確性又如何呢?
表1 肺癌危險(xiǎn)因素分析研究的變量與賦值
表2 部分原始數(shù)據(jù)
從本質(zhì)上講,該研究也是結(jié)局變量為二分類的診斷試驗(yàn)。但是該診斷試驗(yàn)的測(cè)量指標(biāo)很多,應(yīng)該如何預(yù)測(cè)每一位受試者是否患肺癌呢?
我們可以通過二分類Logistic回歸模型,用性別、年齡、BMI、COPD病史和是否吸煙等因素,計(jì)算受試者患肺癌的預(yù)測(cè)概率。
二、SPSS分析方法
1. 數(shù)據(jù)錄入SPSS
2. Logistic回歸分析(關(guān)于Logistic回歸每一步設(shè)置的意義,可參考“SPSS實(shí)例教程:二分類Logistic回歸”。)
選擇Analyze→Regression→Binary Logistic
(1)主對(duì)話框設(shè)置
將因變量cancer送入Dependent框中,將納入模型的自變量sex, age, BMI和COPD變量Covariates中,選擇Forward: LR的自變量篩選方法(Method對(duì)話框)。
(2)Categorical設(shè)置
本研究中,COPD是多分類變量,我們指定“無COPD病史”的研究對(duì)象為參照組,分別比較“輕/中度”和“重度”組相對(duì)于參照組患肺癌的風(fēng)險(xiǎn)。
點(diǎn)擊Categorical→將左側(cè)Covariates中的COPD變量送入右側(cè)Categorical Covariates中。在Reference Category的右側(cè)選擇First(表示選擇變量COPD中,賦值最小的,即“0”作為參照。)→點(diǎn)擊Change→點(diǎn)擊Continue。
(3)Save設(shè)置
點(diǎn)擊Save→選擇Probabilities→點(diǎn)擊Continue。
三、Logistic回歸結(jié)果
1. 納入Logistic回歸模型的變量
最終模型納入了性別(sex)、COPD病史(COPD)和吸煙(smoke)三個(gè)變量。也就是說,這該Logistic回歸模型認(rèn)為,這三個(gè)變量可以預(yù)測(cè)是否患肺癌,而年齡和BMI并沒有預(yù)測(cè)意義。
2. 個(gè)體患肺癌的概率
根據(jù)上述Logistic回歸的結(jié)果,我們可以寫出每個(gè)受試者根據(jù)性別、COPD病史和是否吸煙三個(gè)因素,預(yù)測(cè)是否患肺癌的危險(xiǎn)得分Logit(P):
Logit(P)= -3.062 + 0.836*sex(男=1;女=0) + 0.454*COPD(輕度) + 1.281COPD(中/重度) + 1.237*smoke(無=0;曾吸/現(xiàn)吸=1)
并可以按照以下公式計(jì)算得到每一個(gè)受試者患肺癌的預(yù)測(cè)概率:
實(shí)際上,當(dāng)點(diǎn)選了上述2.5的操作,運(yùn)行該回歸分析后,SPSS會(huì)自動(dòng)生成每一位受試者的預(yù)測(cè)概率(PRE_1),而不需要上述的手工計(jì)算。
至此,我們就可以根據(jù)受試者的真實(shí)患病情況和預(yù)測(cè)概率,評(píng)估根據(jù)性別、COPD病史和是否吸煙三個(gè)因素,預(yù)測(cè)個(gè)體是否患肺癌的準(zhǔn)確性了。
四、ROC曲線的繪制
1. 選擇Analyze→ROC Curve
2. 主對(duì)話框設(shè)置
將已知的疾病情況cancer送入State Variable框中,預(yù)測(cè)概率Predicted probability送入Test Variable中,并在Value of State Variable框中填1→OK。
五、結(jié)果解讀
SPSS的ROC曲線結(jié)果會(huì)給出ROC曲線和曲線下面積。
根據(jù)結(jié)果,我們可以知道該診斷試驗(yàn)的ROC曲線下面積是0.718,判斷其準(zhǔn)確性,并用于與其他診斷試驗(yàn)的比較。至于評(píng)價(jià)診斷試驗(yàn)的其它指標(biāo),需要我們根據(jù)預(yù)測(cè)概率(PRE_1)確定診斷截點(diǎn)(cut-off值)后再計(jì)算,有興趣的小伙伴可以自己嘗試計(jì)算哦。
數(shù)據(jù)分析咨詢請(qǐng)掃描二維碼
若不方便掃碼,搜微信號(hào):CDAshujufenxi
SQL Server 中 CONVERT 函數(shù)的日期轉(zhuǎn)換:從基礎(chǔ)用法到實(shí)戰(zhàn)優(yōu)化 在 SQL Server 的數(shù)據(jù)處理中,日期格式轉(zhuǎn)換是高頻需求 —— 無論 ...
2025-09-18MySQL 大表拆分與關(guān)聯(lián)查詢效率:打破 “拆分必慢” 的認(rèn)知誤區(qū) 在 MySQL 數(shù)據(jù)庫管理中,“大表” 始終是性能優(yōu)化繞不開的話題。 ...
2025-09-18CDA 數(shù)據(jù)分析師:表結(jié)構(gòu)數(shù)據(jù) “獲取 - 加工 - 使用” 全流程的賦能者 表結(jié)構(gòu)數(shù)據(jù)(如數(shù)據(jù)庫表、Excel 表、CSV 文件)是企業(yè)數(shù)字 ...
2025-09-18DSGE 模型中的 Et:理性預(yù)期算子的內(nèi)涵、作用與應(yīng)用解析 動(dòng)態(tài)隨機(jī)一般均衡(Dynamic Stochastic General Equilibrium, DSGE)模 ...
2025-09-17Python 提取 TIF 中地名的完整指南 一、先明確:TIF 中的地名有哪兩種存在形式? 在開始提取前,需先判斷 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ù)庫表、Excel 表、 ...
2025-09-17Excel 導(dǎo)入數(shù)據(jù)含缺失值?詳解 dropna 函數(shù)的功能與實(shí)戰(zhàn)應(yīng)用 在用 Python(如 pandas 庫)處理 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è)操盤手 表格結(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)求開發(fā)時(shí)(如使用requests ...
2025-09-15CDA 數(shù)據(jù)分析師:激活表格結(jié)構(gòu)數(shù)據(jù)價(jià)值的核心操盤手 表格結(jié)構(gòu)數(shù)據(jù)(如 Excel 表格、數(shù)據(jù)庫表)是企業(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ù)法問題 為幫助 Python 數(shù)據(jù)從業(yè)者解決pd.read_csv讀取長(zhǎng)浮點(diǎn)數(shù)據(jù)時(shí)的科學(xué)計(jì)數(shù)法問題 ...
2025-09-12CDA 數(shù)據(jù)分析師:業(yè)務(wù)數(shù)據(jù)分析步驟的落地者與價(jià)值優(yōu)化者 業(yè)務(wù)數(shù)據(jù)分析是企業(yè)解決日常運(yùn)營(yíng)問題、提升執(zhí)行效率的核心手段,其價(jià)值 ...
2025-09-12用 SQL 驗(yàn)證業(yè)務(wù)邏輯:從規(guī)則拆解到數(shù)據(jù)把關(guān)的實(shí)戰(zhàn)指南 在業(yè)務(wù)系統(tǒng)落地過程中,“業(yè)務(wù)邏輯” 是連接 “需求設(shè)計(jì)” 與 “用戶體驗(yàn) ...
2025-09-11塔吉特百貨孕婦營(yíng)銷案例:數(shù)據(jù)驅(qū)動(dòng)下的精準(zhǔn)零售革命與啟示 在零售行業(yè) “流量紅利見頂” 的當(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)景中,聚類分析作為 “無監(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-10