系列簡介

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

今天是第10節(jié)，如果上一節(jié)的內(nèi)容你理解了，那么這一節(jié)就是水到渠來的事情。

我們上一節(jié)其實(shí)已經(jīng)講到了CUSUM的計算公式，當(dāng)前的CUSUM值 = max[0, 上一期的CUSUM值 + 當(dāng)前觀察值與目標(biāo)值的偏差 - 參考值]。

今天我們詳細(xì)講解累積 和控制圖詳細(xì)操作步驟。

第1步：人為設(shè)定標(biāo)準(zhǔn)差系數(shù)k=2，預(yù)警閾值h=2。

第2步：計算每日病例數(shù)的7日平均移動值和7日移動標(biāo)準(zhǔn)差值。

第3步：根據(jù)敏感性的不同，分為了C1、C2和C3三種模式，首先講C1的計算。

第一天的初始CUSUM值=0。

最重要的點(diǎn)來了，第二天的CUSUM值=MAX{【（日病例數(shù)-7日移動均值）/7日移動標(biāo)準(zhǔn)差】-標(biāo)準(zhǔn)差系數(shù)K+上一天CUSUM值}。

第三天的和第二天的一樣，所以算出了第二天的之后只需要拖墜公式即可。

第4步：計算C1預(yù)警方式是否會報警：用算出來的本天的CUSUM值減去預(yù)警閾值h,如果＞0就代表預(yù)警。

就這么簡單，C1以及其是否預(yù)警就計算出來了，我們再來理解一下

CUSUM值的計算公式，其實(shí)就是當(dāng)前的CUSUM值 = max[0, 當(dāng)前觀察值與目標(biāo)值的偏差- 參考值標(biāo)準(zhǔn)差系數(shù)+上一期的CUSUM值 +]，偏差怎么衡量，就是用（日病例數(shù)-7日移動均值）/7日移動標(biāo)準(zhǔn)差，就是看這個差值是幾倍的標(biāo)準(zhǔn)差。

第5步：接下來計算C2

C2第一天的初始CUSUM值=0。

C2第二天的CUSUM值=MAX{【（日病例數(shù)-兩天前的7日移動均值）/兩天前7日移動標(biāo)準(zhǔn)差】-標(biāo)準(zhǔn)差系數(shù)K+上一天CUSUM值}。

第三天的和第二天的一樣，所以算出了第二天的之后只需要拖墜公式即可。

第6步：計算C2預(yù)警方式是否會報警，計算方式和C1一樣，其實(shí)就是用算出來的本天的CUSUM值減去預(yù)警閾值h,如果＞0就代表預(yù)警。

可以看出，C2的計算其實(shí)是和C1類似的，只不過用的移動均值和移動標(biāo)準(zhǔn)差不是本日的，而是2天前的，目的就是使其更加敏感。

第7步：計算C3，C3的計算就簡單了，C3當(dāng)前CUSUM值=當(dāng)天及之前2天的共3個C2相加，是否預(yù)警也是一樣，用算出來的本天的CUSUM值減去預(yù)警閾值h,如果＞0就代表預(yù)警。

所有的原理和關(guān)鍵公式見下圖

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

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

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

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

轉(zhuǎn)載請注明

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

科普健康 | 宣傳疾控

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

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

致力于疾控科普

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

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

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

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

本文具體說明

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

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