在公共衛(wèi)生監(jiān)測的漫長時間軸上，我們追蹤著疾病的發(fā)生與發(fā)展。一個核心問題始終縈繞：疾病的流行趨勢是鐵板一塊，還是分階段、有節(jié)奏地演變？一項重大干預(yù)措施的實施，是否真的在數(shù)據(jù)上留下了“印記”？

傳統(tǒng)回歸模型如同一個“一根筋”的畫家，只能用一支畫筆（一條直線或曲線）描繪整個趨勢。

而Joinpoint回歸則是一位“識時務(wù)”的大師，它能發(fā)現(xiàn)數(shù)據(jù)中的轉(zhuǎn)折，用多支畫筆（多條線段）更真實、更精細(xì)地還原出疾病流行的動態(tài)史詩。本文將從零開始，為你深度解讀這位“大師”的工作原理。

J oinpoint回歸，又名分段線性回歸。其基本思想非常直觀：通過引入一個或多個“連接點”，將整個時間序列分割成若干區(qū)間，并對每個區(qū)間分別進(jìn)行線性回歸擬合。

1998年Kim等首次提出Joinpoint回歸模型，該模型的核心思想是根據(jù)疾病分布的時間特征建立分段回歸，通過若干連接點將研究時間分割成不同區(qū)間，并對每個區(qū)間進(jìn)行趨勢擬合和優(yōu)化，進(jìn)而更詳細(xì)地評價全局時間范圍內(nèi)不同區(qū)間特異性的疾病變化特征。

Joinpoint回歸模型由美國國立癌癥研究所腫瘤控制與人口科學(xué)部開發(fā)，在腫瘤發(fā)病率和死亡率趨勢研究領(lǐng)域得到廣泛地應(yīng)用

如上圖所示，相比于一條直線（全局模型），由Joinpoint生成的折線顯然能更貼切地反映數(shù)據(jù)在不同時期的變化特征。

每一個拐點，都可能對應(yīng)著一個重要的公共衛(wèi)生事件或社會變遷。

Joinpoint回歸最常用的模型是對數(shù)線性模型，因為它能直接輸出具有明確公共衛(wèi)生意義的指標(biāo)——年度變化百分比。

對于一個包含 k個連接點的模型，其回歸方程可表示為：

讓我們來拆解這個看似復(fù)雜的方程：

y：因變量，通常是發(fā)病率、死亡率等指標(biāo)。

x：自變量，即時間變量（如年份）。

β?：截距項。

β?：第一個時間區(qū)間的斜率（基礎(chǔ)斜率）。

τ? … τ?：模型需要估計的連接點（即拐點發(fā)生的具體時間）。

δ? … δ?：斜率的變化量。在連接點 τ?之后，斜率將在 β?的基礎(chǔ)上增加 δ?。

（x - τ?）?：這是一個斷點函數(shù)。當(dāng)x > τ?時，其值為(x - τ?)；當(dāng)x ≤ τ?時，其值為 0。這個函數(shù)是實現(xiàn)分段擬合的關(guān)鍵。

通俗地講，這個模型就是在說：在第一個連接點 τ?之前，趨勢由 β?決定；過了 τ?點，趨勢就變成了 （β? + δ?）；過了 τ?點，趨勢進(jìn)一步變?yōu)?（β? + δ? + δ?），以此類推。

確定連接點的過程和數(shù)量，是Joinpoint回歸的統(tǒng)計精髓，絕非主觀臆斷。

1. 網(wǎng)格搜索法：尋找最優(yōu)拐點位置

對于可能成為連接點的每一個時間點，模型都會計算其擬合的誤差平方和。網(wǎng)格搜索法就像一張精密的網(wǎng)格，系統(tǒng)性地遍歷所有可能的分段方案，最終選擇那個擬合誤差最小的方案，從而確定連接點的最優(yōu)位置。

2. 置換檢驗：確定“幾個拐點”最合適——防止過度擬合

是不是連接點越多越好？絕非如此！過多的連接點會使模型過于復(fù)雜，甚至去擬合數(shù)據(jù)中的隨機噪聲（過度擬合）。那么，多少個連接點才是“恰到好處”的呢？

Joinpoint軟件默認(rèn)采用蒙特卡洛置換檢驗這一穩(wěn)健的方法來解決這個問題。它的基本流程是：

原假設(shè)（H?）：模型有 k個連接點就足夠了。

備擇假設(shè)（H?）：模型需要 k+1個連接點。

檢驗過程：通過數(shù)千次（默認(rèn)4500次）的數(shù)據(jù)隨機置換（模擬），計算出一個P值。如果P值小于顯著性水平（如0.05），則拒絕原假設(shè)，認(rèn)為增加一個連接點能顯著改善模型擬合效果，然后繼續(xù)測試增加下一個連接點的必要性。

這個過程反復(fù)進(jìn)行，直到找到那個“增加連接點也不再帶來顯著改善”的模型為止。這種方法有效避免了人為選擇的主觀性，保證了結(jié)果的客觀性。

模型最終的輸出結(jié)果中，最具解讀價值的是兩個指標(biāo)：

1.第一是年度變化百分比（APC）：用于描述每個獨立時間段內(nèi)的趨勢快慢。其計算公式由對數(shù)線性模型的斜率 β?推導(dǎo)而來：

如何判斷趨勢是否顯著？看t檢驗及其P值或者看APC的95%置信區(qū)間。如果區(qū)間不包含0，則認(rèn)為該時間段內(nèi)的上升（APC>0）或下降（APC<0）趨勢具有統(tǒng)計學(xué)意義。

2.第二是平均年度變化百分比（AAPC）：當(dāng)整個研究期間存在多個趨勢段時，我們需要一個指標(biāo)來概括全局的平均變化速度。

AAPC應(yīng)運而生，它是以各時間段的長度（權(quán)重）對各個區(qū)間的APC進(jìn)行加權(quán)平均計算得出的：

AAPC提供了一個簡潔的“整體故事”，便于對不同人群、不同地區(qū)或不同時期的長期趨勢進(jìn)行比較。

1. 主要應(yīng)用領(lǐng)域

腫瘤流行病學(xué)：評估癌癥發(fā)病/死亡率的長期變化，分析篩查項目、新療法引入的效果。

傳染病監(jiān)測：追蹤疫苗普及后傳染病發(fā)病率的變化趨勢。

公共衛(wèi)生政策評估：驗證控?zé)?、限酒等政策是否在?shù)據(jù)上產(chǎn)生“拐點”效應(yīng)。

2. 重要注意事項

數(shù)據(jù)要求：時間序列不宜過短，通常需要10年以上的數(shù)據(jù)，且每個趨勢段內(nèi)至少應(yīng)有4-5個數(shù)據(jù)點，否則模型識別能力會減弱。

關(guān)聯(lián)≠因果：Joinpoint回歸能提示趨勢變化的可能時間點，并將其與歷史事件關(guān)聯(lián)，但不能證明因果關(guān)系。合理解釋必須結(jié)合流行病學(xué)背景知識。

模型假設(shè)：使用對數(shù)線性模型時，默認(rèn)率的對數(shù)與時間呈線性關(guān)系，在實際分析中需留意這一假設(shè)是否合理。

3.總結(jié)

Joinpoint回歸通過其獨特的分段擬合思想和嚴(yán)謹(jǐn)?shù)慕y(tǒng)計檢驗，將我們從對趨勢的籠統(tǒng)認(rèn)知，帶入到對疾病動態(tài)演變的精細(xì)化解讀階段。它幫助我們回答的不再僅僅是“總體是升是降”，而是“在什么時候、以何種速度、發(fā)生了怎樣的轉(zhuǎn)變”。

掌握了這些原理，你就拿到了理解Joinpoint回歸的鑰匙。在下一篇文章中，我們將手把手帶你實操美國國家癌癥研究所（NCI）的Joinpoint軟件，從數(shù)據(jù)準(zhǔn)備、參數(shù)設(shè)置到結(jié)果解讀，讓你真正將這一強大工具應(yīng)用于實際工作之中！敬請期待。

編輯：普通疾控人 | 審核：詩酒趁年華

文章來源 | 原創(chuàng)

說明 | 轉(zhuǎn)載只為分享，如有侵權(quán)聯(lián)系刪除

?版權(quán)聲明 | 部分信息和圖片來自公開網(wǎng)絡(luò)

轉(zhuǎn)載請注明

再次轉(zhuǎn)載請注明出處

科普健康 | 宣傳疾控

本號為多位疾控機構(gòu)從業(yè)者運營

重點關(guān)注國內(nèi)外健康事件

致力于疾控科普

在做好科普服務(wù)大眾的同時

做好疾控機構(gòu)的宣傳

讓更多的人了解疾控，擁抱健康

歡迎加「小編」微信(cdcjkr126com)

本文具體說明

本文為原創(chuàng)內(nèi)容，文章為個人理解所學(xué)，不涉及疫情信息及內(nèi)部保密數(shù)據(jù)，發(fā)表的目的為自我總結(jié)及給有需求的人士學(xué)習(xí)使用。如有不妥之處，歡迎聯(lián)系小編修改、刪除。

更多精彩視頻，盡在“CDC疾控人”視頻號