j09188 / SIR

The source code of the article "針對台灣公衛政策的SIR模型擴增:新冠肺炎疫情的擬合與應用"
1 stars 0 forks source link

你好,請問可以咨詢這個R0的計算方法嗎,想學習一下,沒有看懂,謝謝您 #4

Open HiJackk opened 1 year ago

HiJackk commented 1 year ago

請問可以咨詢這個R0的計算方法嗎,Next Generation Matrix,沒有看懂這個NGM方法,請問可以解答一下或者提供一些文獻嗎?謝謝您!

2023-03-21 00 35 32
j09188 commented 1 year ago

這裡參考的論文是(The construction of next-generation matrices for compartmental epidemic models)。 這裡對應的方程式(義大利隔離模型)如圖所示,其中T(transmission matrix)代表從未感染S到感染{C,I,Q_C,Q_I}的可能途徑,因為我們只假設S有α的機率變成是C,因此在T matrix中只有一個entry=α。 Σ(transition matrix)代表{C,I,Q_C,Q_I}之間彼此轉換的關係,像是第一列第一行的 -(ω+ρ) 就代表是d/dt C這個式子中關於C的係數、第一列第二行的 -αβτ(1-ρ) 則是代表d/dt C中關於I的係數,其他就依此類推。關於用T & Σ求出R_0值的部分,我也只是按照參考論文的方法做,詳細的證明還是要看本來的論文。 義大利隔離模型