系列簡(jiǎn)介

這是我們一系列原創(chuàng)技術(shù)貼，從易到難，每天學(xué)習(xí)一點(diǎn)。所有內(nèi)容均為疾控?cái)?shù)據(jù)分析、科研論文相關(guān)，或者說(shuō)很多和現(xiàn)在的熱門監(jiān)測(cè)預(yù)警相關(guān)，所以我們這個(gè)系列就叫“監(jiān)測(cè)預(yù)警基礎(chǔ)”。

今天是第26節(jié)，也是MEM的最后一節(jié)。

上一節(jié)我們已經(jīng)講了很多了，可以采用函數(shù)計(jì)算流行期了，但是這并不是MEM的真正目的，MEM的目的是對(duì)目標(biāo)時(shí)間進(jìn)行預(yù)警。

預(yù)警就需要閾值，我們這一節(jié)就講閾值的設(shè)置。

MEM模型通過(guò)歷史數(shù)據(jù)劃分流行期后，下一步是利用劃分結(jié)果計(jì)算閾值：

流行閾值（Epidemic Threshold）：基于流行前期和流行后期（非流行期）數(shù)據(jù)，確定流行開(kāi)始和結(jié)束的臨界值。

強(qiáng)度閾值（Intensity Threshold）：基于流行期數(shù)據(jù)，定義中、高、極高流行強(qiáng)度的分級(jí)標(biāo)準(zhǔn)。

這些閾值用于監(jiān)測(cè)當(dāng)前數(shù)據(jù)：當(dāng)監(jiān)測(cè)指標(biāo)超過(guò)流行閾值時(shí)，認(rèn)為流行開(kāi)始；超過(guò)強(qiáng)度閾值時(shí)，判定流行強(qiáng)度等級(jí)。

我們先計(jì)算流行閾值，也就是流行開(kāi)始和流行結(jié)束閾值。

流行開(kāi)始閾值是指當(dāng)監(jiān)測(cè)指標(biāo)（如發(fā)病率）超過(guò)此值時(shí)，認(rèn)為流行期開(kāi)始。流行結(jié)束閾值是指當(dāng)監(jiān)測(cè)指標(biāo)低于此值時(shí)，認(rèn)為流行期結(jié)束。這兩個(gè)閾值共同定義了流行期的邊界，確保預(yù)警信號(hào)及時(shí)且準(zhǔn)確。

1.計(jì)算原理 利用已經(jīng)劃分的流行前期和流行后期的監(jiān)測(cè)指標(biāo)值計(jì)算流行開(kāi)始和結(jié)束閾值。選擇各流行前期中監(jiān)測(cè)指標(biāo)最大的n個(gè)值納人分析(n=30/N，N為納入分析的流行季數(shù)，取值范圍為5~10),計(jì)算這些值的算術(shù)均數(shù)的單側(cè)95%置信區(qū)間、將置信區(qū)間的上限作為流行開(kāi)始閾值。同理，利用各流行后期數(shù)據(jù)計(jì)算流行結(jié)束閾值。

2.數(shù)據(jù)來(lái)源

使用歷史流行季的流行前期和流行后期數(shù)據(jù)（即非流行期數(shù)據(jù)），例如我們上一節(jié)的數(shù)據(jù)流行前期是第1-18周，流行后期是第33-52周。

篩選數(shù)據(jù)點(diǎn)：例如，如果N=9（即9個(gè)流行季），則n = 30 / 9 ≈ 3.33，通常取整為4，然后從每個(gè)歷史流行季的流行前期中，選擇監(jiān)測(cè)指標(biāo)值最大的4個(gè)值組成一個(gè)數(shù)據(jù)集。

3.計(jì)算這些值的算術(shù)均數(shù)

計(jì)算算術(shù)均數(shù)的單側(cè)95%置信區(qū)間（上限），作為流行閾值。流行開(kāi)始閾值使用流行前期的數(shù)據(jù)計(jì)算。流行結(jié)束閾值使用流行后期的數(shù)據(jù)計(jì)算，方法相同。

4.為什么這樣做？ 使用非流行期的高值數(shù)據(jù)，可以捕捉到“基線波動(dòng)”的上限，避免將正常波動(dòng)誤判為流行。95%置信區(qū)間提供了統(tǒng)計(jì)顯著性，確保閾值可靠。
1.計(jì)算目的

在流行期內(nèi)，進(jìn)一步將流行強(qiáng)度分為中、高、極高等級(jí)，以便采取差異化防控措施。例如，超過(guò)“高流行閾值”時(shí)，可能需要啟動(dòng)應(yīng)急響應(yīng)。

2.計(jì)算原理

利用已經(jīng)劃分的流行期的監(jiān)測(cè)指標(biāo)值計(jì)算中、高和極高流行強(qiáng)度閾值。選擇各流行期監(jiān)測(cè)指標(biāo)最大的n個(gè)值納人分析(n=30/N),計(jì)算這些值的幾何均數(shù)的單側(cè)50%、90%和95%(或40%、90%和97.5%)置信區(qū)間，將此3個(gè)置信區(qū)間的上限依次作為中、高和極高流行強(qiáng)度閾值。

3.數(shù)據(jù)來(lái)源

使用歷史流行季的流行期數(shù)據(jù)（即流行開(kāi)始到結(jié)束的周次數(shù)據(jù)），例如之前的例子中，流行期是第19-32周。

4.篩選數(shù)據(jù)點(diǎn)

從每個(gè)歷史流行季的流行期中，選擇監(jiān)測(cè)指標(biāo)值最大的n個(gè)值（n = 30 / N，與流行閾值相同）。

5.統(tǒng)計(jì)計(jì)算

先計(jì)算數(shù)據(jù)對(duì)數(shù)的均數(shù)和標(biāo)準(zhǔn)誤，然后反推回原始尺度。

中流行強(qiáng)度閾值：?jiǎn)蝹?cè)50%置信區(qū)間上限（即中位數(shù)水平）。

高流行強(qiáng)度閾值：?jiǎn)蝹?cè)90%置信區(qū)間上限。

極高流行強(qiáng)度閾值：?jiǎn)蝹?cè)95%或97.5%置信區(qū)間上限。

將所有選出的n個(gè)值合并為一個(gè)數(shù)據(jù)集，計(jì)算這些值的幾何均數(shù)（而非算術(shù)均數(shù)），因?yàn)閭魅静?shù)據(jù)常呈偏態(tài)分布，幾何均數(shù)更能代表集中趨勢(shì)，計(jì)算幾何均數(shù)的單側(cè)置信區(qū)間上限，作為強(qiáng)度閾值：公式示意（幾何均數(shù)的置信區(qū)間需對(duì)數(shù)轉(zhuǎn)換）。

到此為止，MEM就基本講完了。

它的具體實(shí)現(xiàn)有兩種方式，一是R軟件加載運(yùn)行mem包，二是通過(guò)網(wǎng)頁(yè)版MEM Dashboard實(shí)現(xiàn)，（http://memapp.iecscyl.com:8080/），我們后面再詳細(xì)講述！

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

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

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

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

轉(zhuǎn)載請(qǐng)注明

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

科普健康 | 宣傳疾控

本號(hào)為多位疾控機(jī)構(gòu)從業(yè)者運(yùn)營(yíng)

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

致力于疾控科普

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

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

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

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

本文具體說(shuō)明

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

更多精彩視頻，盡在“CDC疾控人”視頻號(hào)