王朝陽(yáng)，劉泉泉，田正福，馬佳琪，張利萍

（新疆醫(yī)科大學(xué)醫(yī)學(xué)工程技術(shù)學(xué)院，新疆 烏魯木齊 830017）

自2019 年年底至今，在全球200 多個(gè)國(guó)家和地區(qū)已發(fā)現(xiàn)新型冠狀病毒感染，且其感染病例的數(shù)量還在逐步上升[1]。我國(guó)新型冠狀病毒感染情況已得到基本控制，僅有少數(shù)境外新型冠狀病毒感染輸入病例。針對(duì)新型冠狀病毒感染的大流行，很多學(xué)者從各個(gè)方面如病毒學(xué)[2]、流行病學(xué)[3]等進(jìn)行了研究和報(bào)道。該病毒傳染性強(qiáng)，并且在世界范圍內(nèi)迅速傳播。截至2021 年8 月31 日，全球累積超過2 億新型冠狀病毒感染者，我國(guó)累積報(bào)告確診新型冠狀病毒感染94 879 例，累計(jì)境外輸入型冠狀病毒感染者8321例。Joinpoint 模型是由美國(guó)國(guó)立癌癥研究所開發(fā)的用于腫瘤時(shí)間趨勢(shì)分析的統(tǒng)計(jì)方法。分析始于存在最小的“轉(zhuǎn)折點(diǎn)”數(shù)量，檢驗(yàn)是否有1 個(gè)或更多的“轉(zhuǎn)折點(diǎn)”，將年份發(fā)病率曲線分為若干有統(tǒng)計(jì)學(xué)意義的區(qū)段[4]。目前，我國(guó)未再出現(xiàn)全國(guó)范圍的新型冠狀病毒感染，但是國(guó)外疫情的防控形勢(shì)依然不容樂觀，每日新增的感染人數(shù)仍然維持在較高的水平[5]，我國(guó)“外防輸入”壓力依然較大。相關(guān)研究或報(bào)道均表明新型冠狀病毒感染防控工作常態(tài)化后我國(guó)幾起規(guī)模較大的本土疫情與境外輸入人或物有密切聯(lián)系[6]。本研究根據(jù)我國(guó)每日新增新型冠狀病毒病感染人數(shù)和境外日輸入新型冠狀病毒感染人數(shù)的時(shí)間趨勢(shì)分別擬合對(duì)數(shù)線性Joinpoint 模型，得到重要節(jié)點(diǎn)，再依據(jù)我國(guó)采取的政策分析節(jié)點(diǎn)前后感染情況的波動(dòng)原因。

1 資料與方法

1.1 資料來(lái)源 收集2021 年5 月1 日-8 月31 日我國(guó)新型冠狀病毒感染日新增感染的病例數(shù)及境外日輸入新型冠狀病毒感染的病例數(shù)，數(shù)據(jù)均來(lái)自中華人民共和國(guó)國(guó)家衛(wèi)生健康委員會(huì)官網(wǎng)。

1.2 方法

1.2.1 Joinpoint 回歸模型及計(jì)算原理 Joinpoint 回歸模型（JPR 模型）主要用于分析時(shí)間序列或其他有序數(shù)據(jù)的趨勢(shì)變化。Joinpoint 回歸模型又稱片段回歸，線段回歸或多階段回歸，是將一個(gè)長(zhǎng)期趨勢(shì)線分成若干段，每段用連續(xù)性的線性進(jìn)行描述的方法[7]。JPR 模型將輸入的數(shù)據(jù)擬合成一個(gè)長(zhǎng)期趨勢(shì)線，其中會(huì)包含若干個(gè)轉(zhuǎn)折點(diǎn)和連續(xù)的線段，并對(duì)這些點(diǎn)和線段進(jìn)行分析。通過日平均變化百分比（AAPC）反映平均變化速度和日變化百分比（APC）反映階段性變化速度來(lái)表示，從而揭示每日境外輸入和我國(guó)每日新增新冠感染例數(shù)趨勢(shì)變化的特征。

1.2.2 轉(zhuǎn)折點(diǎn)及分析方法 JPR 模型通過GSM 法（GSM：網(wǎng)絡(luò)搜索法是Joinpoint 默認(rèn)的建模方法）建立所有可能存在的區(qū)間分段函數(shù)連接點(diǎn)，其中GSM將參數(shù)所處的空間劃分成網(wǎng)格，計(jì)算各種情況下的誤差平方和（SSE）和均方差（MSE），選擇MSE 最小的網(wǎng)格位點(diǎn)為最優(yōu)轉(zhuǎn)折點(diǎn)[8]。根據(jù)我國(guó)在此次疫情中采取的措施結(jié)合拐點(diǎn)分析APC 值下降或上升來(lái)分析每日境外輸入對(duì)我國(guó)日新增新冠感染例數(shù)的影響。將中國(guó)新型冠狀病毒感染日新增感染的病例數(shù)及境外日輸入新型冠狀病毒感染的病例數(shù)建立Excel 數(shù)據(jù)表，并轉(zhuǎn)為CSV 數(shù)據(jù)文件，將日期2021 年5月1 日至8 月31 日編碼為1-123 替代。采用JPR軟件中泊松分布模型對(duì)我國(guó)每日新型冠狀病毒感染新增確診病例和境外輸入病例時(shí)間趨勢(shì)分析，最大轉(zhuǎn)折點(diǎn)設(shè)置為2。

1.3 統(tǒng)計(jì)學(xué)方法 采用SPSS 20.0 軟件進(jìn)行描述性統(tǒng)計(jì)分析和正態(tài)分布檢驗(yàn)；采用Joinpoint Regression Program 4.6.0.0 分析疫情階段趨勢(shì)特征；JPR 軟件在癌癥研究中常以年度為觀察單位，本研究資料以日為觀察單位。

2 結(jié)果

2.1 描述性分析結(jié)果

2.1.1 基本情況 2021 年5 月-8 月，我國(guó)共有24 個(gè)?。ㄗ灾螀^(qū)、直轄市）報(bào)告境外輸入新型冠狀病毒感染病例2623 例，云南省境外輸入感染人數(shù)高達(dá)607例，見表1。

表1 我國(guó)2021 年5 月-8 月各省境外輸入新型冠狀病毒感染病例數(shù)和本土病例數(shù)（n）

2.1.2 時(shí)間分布 2021 年7 月-8 月，境外輸入新型冠狀病毒感染人數(shù)較多，2021 年7 月20 日-8 月20日，國(guó)內(nèi)每日新增病例速度較快，最高值的到達(dá)時(shí)間分別為7 月19 日和8 月9 日，最高值分別為42 例和143 例，其余時(shí)間總體呈現(xiàn)平穩(wěn)波動(dòng)變化，見圖1。

圖1 我國(guó)每日新型冠狀病毒感染確診病例和境外輸入病例報(bào)告日期分布

2.2 Joinpoint 回歸模型分析 對(duì)2021 年5 月1 日-2021 年8 月31 日我國(guó)新型冠狀病毒日新增感染人數(shù)及境外日輸入新型冠狀病毒感染人數(shù)的時(shí)間趨勢(shì)分別擬合對(duì)數(shù)線性JPR 模型，兩者的MSE 分別為2.95 和2.08，兩者的回歸方程不平行，表明該模型的擬合效果較好，兩組數(shù)據(jù)擬合的回歸方程及其趨勢(shì)變化特征指標(biāo)均存在差異。其中，境外日輸入新型冠狀病毒感染人數(shù)的時(shí)間趨勢(shì)出現(xiàn)2 個(gè)轉(zhuǎn)折點(diǎn)，分別在2021 年7 月11 日和7 月19 日，5 月1 日至7 月11 日為緩慢上升期（APC=0.894%，t=4.219，P＜0.05）；7 月11 日-7 月19 日為快速上升期（APC=5.995%，t=1.129，P＞0.05）；7 月19 日-8 月31 日為緩慢下降期（APC=-0.656%，t=-2.023，P＜0.05）。國(guó)內(nèi)每日新增新型冠狀病毒感染時(shí)間趨勢(shì)共出現(xiàn)2 個(gè)轉(zhuǎn)折點(diǎn)，分別在2021 年7 月12 日和8 月8 日，其中5月1 日-7 月12 日為上升期（APC=0.970%，t=4.083，P＜0.05）；7 月13 日-8 月8 日為快速上升期（APC=5.408%，t=8.585，P＜0.05）；8 月9 日-31 日為快速下降期（APC=-8.249%，t=-10.976，P＜0.05），見圖2。

圖2 我國(guó)2021 年5 月-8 月每日新增新型冠狀病毒確診病例和境外輸入病例階段趨勢(shì)

2021 年5 月1 日-7 月22 日，境內(nèi)新型冠狀病毒感染病例累計(jì)量以較慢的速度增加（APC=0.436%，t=137.080，P＜0.05）。7 月22 日-月14 日，新型冠狀病毒感染累計(jì)病例開始以超過前一時(shí)間段的速度增加（APC=1.336%，t=68.655，P＜0.05），最后從8 月14 日-31 日增加的速度急劇放緩（APC=0.360%，t=13.081，P＜0.05），7 月22 日-8 月14 日從5800 例增加到7820 例，在這種嚴(yán)峻的形勢(shì)中，8 月8 日，新型冠狀病毒感染新增患者為119 例，達(dá)到日確診人數(shù)的最高峰。8 月14 日-8 月31 日增長(zhǎng)速度下降（APC=0.36%，t=13.081，P＜0.05）。2021 年5 月1 日-8 月31 日，每日新增新型冠狀病毒感染病例數(shù)一共出現(xiàn)2 個(gè)轉(zhuǎn)折點(diǎn)，為2021 年7 月12 日和8月8 日。伴隨著境外輸入新型冠狀病毒感染人數(shù)不斷增加，7 月20 日南京祿口機(jī)場(chǎng)首次檢測(cè)出9 例陽(yáng)性病例后，疫情蔓延多省多市，這與7 月12 日轉(zhuǎn)折點(diǎn)的出現(xiàn)相符。7 月27 號(hào)南京市新型冠狀病毒疫情防控工作指揮部發(fā)布“關(guān)于進(jìn)一步嚴(yán)格疫情防控工作的通告（第6 號(hào)）”。通過進(jìn)一步落實(shí)疫情防控措施，快速有效阻斷疫情擴(kuò)散蔓延，嚴(yán)格執(zhí)行疫情防控工作通告的十項(xiàng)要求，所以出現(xiàn)了8 月2 日的轉(zhuǎn)折點(diǎn)。

圖3 我國(guó)2021 年5 月-8 月新型冠狀病毒感染累計(jì)病例時(shí)間趨勢(shì)圖

3 討論

2021 年5 月1 日-8 月31 日，新型冠狀病毒借助全球化的人口流動(dòng)大潮，跨國(guó)界肆虐全球，我國(guó)的疫情防控進(jìn)入了防境外輸入的新階段[5]。已有研究報(bào)道[9，10]，無(wú)癥狀的新型冠狀病毒感染者在潛伏期亦可有效傳播病毒，無(wú)癥狀感染者已成為疫情傳播的主要風(fēng)險(xiǎn)之一。從5 月1 日到8 月31 日，境外輸入病例主要集中于云南、廣東、上海、四川、福建5個(gè)省市。云南輸入境外新型冠狀病毒感染人數(shù)最多，主要原因是瑞麗市旅游業(yè)所帶來(lái)的國(guó)外游客和回國(guó)的游客，導(dǎo)致云南省境外輸入人數(shù)居高不下，3 月17 日云南省采取緊急措施控制境外輸入疫情的輸入和擴(kuò)散[11]。江蘇省在這期間本土感染病例最高共821 例，感染源是國(guó)外疫情的輸入，導(dǎo)致江蘇省感染防控失控的原因主要分為4 點(diǎn)：病毒毒株變異、傳染速度極快、感染的場(chǎng)所比較密集、新冠防控意識(shí)降低。國(guó)內(nèi)多起本土暴發(fā)疫情確診病例占比較境外輸入疫情高[12]。

2021 年7 月11 日-19 日境外日輸入新型冠狀病毒感染人數(shù)快速上升，這進(jìn)一步導(dǎo)致我國(guó)2021 年7 月13 日-8 月8 日新型冠狀病毒感染人數(shù)快速上升[13-15]。局部地區(qū)的疫情爆發(fā)和失控，一定程度上是由境外新型冠狀病毒感染患者的輸入和該地區(qū)存在防疫松懈導(dǎo)致疫情二次暴發(fā)，最終導(dǎo)致本土疫情的增長(zhǎng)[16]。本次研究從全國(guó)每日新增新型冠狀病毒感染人數(shù)的下降時(shí)間段得出，通過積極隔離、集中觀察和檢查密切接觸者、限制人員流動(dòng)、防護(hù)措施宣傳教育及社區(qū)網(wǎng)格化管理、加強(qiáng)境外輸入的管控等措施，可有效遏制新型冠狀病毒感染的蔓延[17]。

當(dāng)前，隨著新型冠狀病毒變異、疫情變化、疫苗接種普及和防控經(jīng)驗(yàn)增加，我國(guó)新型冠狀病毒感染防控面臨著新形勢(shì)，防控工作進(jìn)入新的階段。為更好地防控新型冠狀病毒感染，應(yīng)進(jìn)一步提高老年人新型冠狀病毒疫苗接種率，加快提高疫苗加強(qiáng)免疫接種覆蓋率；強(qiáng)化對(duì)新型冠狀病毒感染的監(jiān)測(cè)與應(yīng)對(duì)，及時(shí)跟蹤研判并采取針對(duì)性措施；廣泛宣傳倡導(dǎo)“每個(gè)人都是自己健康第一責(zé)任人”的理念，倡導(dǎo)堅(jiān)持個(gè)人防護(hù)措施。