WOLFRAM NOTEBOOK

Epidemiological Models for Influenza and COVID-19—part 4

Robert B. Nachbar
Original post: 11-Mar-2020
Original post: 11-Mar-2020
Original post: 11-Mar-2020
  • This set of notebooks uses a package which can be downloaded from https://www.wolframcloud.com/obj/rnachbar/Published/CompartmentalModeling--v1.wl
  • The package should be in the same directory as the notebooks, and is automatically loaded as part of the initialization.
  • Table of Contents

  • Use the subsection cells to navigate to the other notebooks
  • COVID-19

  • Models
  • SEIQRD Model with standard incidence
  • Hubei outbreak—SEIQRD
  • Outbreak start 20 Dec 2019
  • Fitting data
  • Manual fit
  • Optimization—
    500000
  • Optimization—
    100000
  • Optimization—
    50000
  • Outbreak start 30 Dec 2019
  • Fitting data
  • Manual fit
  • Optimization—
    500000
  • Optimization—
    100000
  • Optimization—
    50000
  • Comparisons
  • References

    Initialization

    COVID-19

    Models

    SEIQRD Model with standard incidence

  • These are the ODEs
  • forceOfInfectionSEIQRD=λ[t_]
    [t]
    -[t]-[t]
    ;{varsSEIQRD,odesSEIQRD}=Values@KineticCompartmentalModel
    βλ[t]
    ,
    ζ
    ,
    δ
    ,
    γ
    ,
    γ
    ,
    μ
    ,
    μ
    ,t{1,-1};Column[odesSEIQRD]
    Out[]=
    [t]μ[t]+μ[t]
    [t]-ζ[t]+β[t]λ[t]
    [t]ζ[t]-γ[t]-δ[t]-μ[t]
    [t]δ[t]-γ[t]-μ[t]
    [t]γ[t]+γ[t]
    [t]-β[t]λ[t]
    {susceptibleSEIQRD,exposedSEIQRD,infectedSEIQRD,quarantinedSEIQRD,recoveredSEIQRD,deadSEIQRD}={,,,,,}/.ParametricNDSolve[Join[odesSEIQRD/.forceOfInfectionSEIQRD,{[0]-I0,[0]0,[0]I0,[0]0,[0]0,[0]0}],varsSEIQRD,{t,0,500},{,I0,β,ζ,γ,δ,μ}];

    Hubei outbreak—SEIQRD

    Outbreak start 20 Dec 2019

    Fitting data

    In[]:=
    t0=DateObject["20 Dec 2019"]
    Out[]=
    Day: Fri 20 Dec 2019
    In[]:=
    fitData1=fitData=fitDataWDR
    Hubei, China
    ADMINISTRATIVE DIVISION
    ,t0,"makeCorrectionForHubei"True,"dateRange"{All,"26 Feb 2020"};Dimensions/@%
    Out[]=
    {{36,2},{36,2},{36,2}}
    In[]:=
    ListPlot[fitData,PlotLabel"Hubei",FrameLabel{"time (d)","# individuals"},PlotLegends{"confirmed"-("recovered"+"dead"),"recovered","dead"}]
    Out[]=
    confirmed-dead-recovered
    recovered
    dead

    Manual fit

    Out[]=
    effective population size
    50000
    100000
    500000
    1000000
    5000000
    β = 142. (infection)
    ζ = 0.85 (incubation)
    γ = 0.029 (recovery)
    δ = 124. (quarantine)
    μ = 0.0047 (death)
    t
    max
    y
    max
    0.5

    Optimization
    500000

    Optimization
    100000

    Optimization
    50000

    Outbreak start 30 Dec 2019

    Fitting data

  • Since we don’t know with certainty when the first infection occurred, let’s try a later date to see what effect that has
  • In[]:=
    t0=DateObject["30 Dec 2019"]
    Out[]=
    Day: Mon 30 Dec 2019
    In[]:=
    fitData2=fitData=fitDataWDR
    Hubei, China
    ADMINISTRATIVE DIVISION
    ,t0,"makeCorrectionForHubei"True,"dateRange"{All,"26 Feb 2020"};Dimensions/@%
    Out[]=
    {{36,2},{36,2},{36,2}}

    Comparisons

  • increasing the population size makes and rise more slowly, and decreases the epidemic size
  • shifting t0 to a later date (closer to the data) appears to make the peak narrower, and reverses the population size trends for β and δ
  • References
    [TG] T. Götz, “First attempts to model the dynamics of the coronavirus outbreak 2020”, https://arxiv.org/pdf/2002.03821.pdf
    [PYZ] L. Peng, W. Yang, D. Zhang, C. Zhuge, L. Hong “Epidemic analysis of COVID-19 in China by dynamical modeling”, https://www.medrxiv.org/content/10.1101/2020.02.16.20023465v1
    [ZCW] Y. Zhou, Z. Chen, X. Wu, Z. Tian, L. Cheng, L. Ye “The Outbreak Evaluation of COVID-19 in Wuhan District of China”, https://arxiv.org/pdf/2002.09640.pdf
    [JDL] J. Jia, J. Ding, S. Liu, G. Liao, J. Li, B. Duan, G. Wang, R. Zhang “Modeling the Control of COVID-19: Impact of
    Policy Interventions and Meteorological Factors”, https://arxiv.org/pdf/2003.02985.pdf
    [EGE] E. G. M E. “An SEIR like model that fits the coronavirus infection data”, https://community.wolfram.com/groups/-/m/t/1888335
    [AA] A. Antonov “Basic experiments workflow for simple epidemiological models”, https://community.wolfram.com/groups/-/m/t/1895686
    [AV] J. Arino, P. van den Driessche “Time delays in Epidemic Models; Modeling and Numerical Considerations” in “Delay Differential Equations and Applications”, O. Arino (ed.) Springer, 2006.
    [FB] F. Brauer “Reproduction numbers and final size relations”, https://www.fields.utoronto.ca/programs/scientific/10-11/drugresistance/emergence/fred1.pdf
    [BCR] M. Biggerstaff, S. Cauchemez, C. Reed, M. Gambhir, L. Finelli “Estimates of the reproduction number for seasonal, pandemic, and zoonotic influenza: a systematic review of the literature” BMC Infectious Diseases, 14, 480 (2014), http://www.biomedcentral.com/1471-2334/14/480
    [MM] M. Martcheva “An introduction to mathematical epidemiology” Springer, 2015.
    [HJR] H. J. Rose “The use of amantadine and influenza vaccine in a type A influenza epidemic in a boarding school”, Journal of Royal College of General Practitioners, 30, 619-621 (1980). PubMedCentral
    [FT] Z. Feng, H. R. Thieme “Recurrent Outbreaks of Childhood Diseases Revisited: The Impact of Isolation”, Math. Biosciences, 128, 93-130 (1995). https://doi.org/10.1016/0025-5564(94)00069-C
    [BK] S. Boseley, L. Kuo “Huge rise in coronavirus cases casts doubt over scale of epidemic”, The Guardian, 13 Feb 2020, https://www.theguardian.com/world/2020/feb/13/huge-rise-coronavirus-cases-raises-doubts-scale-epidemic-china
    [DWC] Z. Du, L. Wang, S. Cauchemex, X. Xu, X. Wang, B. J. Cowling, L. A. Meyers “Risk for Transportation of 2019 Novel Coronavirus (COVID-19) from Wuhan to Cities in China”, https://doi.org/10.1101/2020.01.28.20019299
    [CXL] J. Cai, J. Xu, D. Lin, Z. Yang, L. Xu, Z, Qu, Y. Zhang, H. Zhang, R. Jia, P. Liu, X. Wang, Y. Ge, A. Xia, H. Tian, H. Chang, C. Wang, J. Li, J. Wang, M. Zheng “A Case Series of children with 2019 novel coronavirus infection: clinical and epidemiological features”, Clinical Infectious Diseases, https://doi.org/10.1093/cid/ciaa198
    [CWB] B. J. Coburn, B. G. Wagner, S. Blower “Modeling influenza epidemics and pandemics: insights into the future of swine flu (H1N1)”, BMC Medicine, 7, (2009), http://www.biomedcentral.com/1741-7015/7/30
    Initialization

    General

    Fitting data

    Fitting error

    SEIQRD

    Wolfram Cloud

    You are using a browser not supported by the Wolfram Cloud

    Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.


    I understand and wish to continue anyway »

    You are using a browser not supported by the Wolfram Cloud. Supported browsers include recent versions of Chrome, Edge, Firefox and Safari.