Message Boards

GROUPS:

Anton Antonov

SEI2HR model with quarantine scenarios

Anton Antonov, Accendo Data LLC

Posted 3 months ago

705 Views

0 Replies

3 Total Likes

Follow this post

SEI2HR model with quarantine scenarios

Version 0.8

Anton Antonov
MathematicaForPrediction at WordPress
SystemModeling at GitHub
March, April 2020

Introduction

The epidemiology compartmental model, [Wk1], presented in this notebook -- SEI2HR -- deals with the left-most and middle rectangles in this diagram:

ImageResize[Import["https://github.com/antononcube/SystemModeling/raw/master/Projects/Coronavirus-propagation-dynamics/Diagrams/Coronavirus-propagation-simple-dynamics.jpeg"],900]

In[]:=

Out[]=

“SEI2HR” stands for “Susceptible, Exposed, Infected two, Hospitalized, Recovered” (populations.)

In this notebook we also deal with quarantine scenarios.

Remark: We consider the contagious disease propagation models as instances of the more general System Dynamics (SD) models. We use SD terminology in this notebook.

The models

SEI2R

The model SEI2R are introduced and explained in the notebook [AA2]. SEI2R differs from the classical SEIR model , [Wk1, HH1], with the following elements:

1. Two separate infected populations: one is "severely symptomatic", the other is "normally symptomatic"

2. The monetary equivalent of lost productivity due to infected or died people is tracked.

SEI2HR

For the formulation of SEI2HR we use a system of Differential Algebraic Equations (DAE’s). The package [AAp1] allows the use of a formulation that has just Ordinary Differential Equations (ODE’s).

Here are the unique features of SEI2HR:

People stocks

◼

Two types of infected populations: normally symptomatic and severely symptomatic

◼

Hospitalized population

◼

Deceased from infection population

◼

Hospital beds

◼

Hospital beds are a limited resource that determines the number of hospitalized people

◼

Only severely symptomatic people are hospitalized according to the available hospital beds

◼

The hospital beds stock is not assumed constant, it has its own change rate.

◼

Money stocks

◼

The money from lost productivity are tracked

◼

The money for hospital services are tracked

◼

SEI2HR’s place a development plan

This graph shows the “big picture” of the model development plan undertaken in [AAr1] and SEI2HR (discussed in this notebook) is in that graph:

Out[]=

Notebook structure

The rest of notebook has the following sequence of sections:

Package load section

◼

SEI2HR structure in comparison of SEI2R

◼

Explanations of the equations of SEI2HR

◼

Quarantine scenario modeling preparation

◼

Parameters and initial conditions setup

◼

Populations, hospital beds, quarantine scenarios

◼

Parametric simulation solution

◼

Interactive interface

◼

Sensitivity analysis

◼

Load packages

The epidemiological models framework used in this notebook is implemented with the packages [AAp1, AAp2, AA3]; many of the plot functions are from the package [AAp4].

Import["https://raw.githubusercontent.com/antononcube/SystemModeling/master/Projects/Coronavirus-propagation-dynamics/WL/EpidemiologyModels.m"];Import["https://raw.githubusercontent.com/antononcube/SystemModeling/master/Projects/Coronavirus-propagation-dynamics/WL/EpidemiologyModelModifications.m"];Import["https://raw.githubusercontent.com/antononcube/SystemModeling/master/Projects/Coronavirus-propagation-dynamics/WL/EpidemiologyModelingVisualizationFunctions.m"];Import["https://raw.githubusercontent.com/antononcube/SystemModeling/master/WL/SystemDynamicsInteractiveInterfacesFunctions.m"];

In[]:=

SEI2HR extends SEI2R

The model SEI2HR is an extension of the model SEI2R, [AA2].

Here is SEI2R:

reprTP="AlgebraicEquation";lsModelOpts={"Tooltips"True,TooltipStyle{BackgroundYellow,CellFrameColorGray,FontSize20}};modelSEI2R=SEI2RModel[t,"InitialConditions"True,"RateRules"True,"TotalPopulationRepresentation"reprTP];ModelGridTableForm[modelSEI2R,lsModelOpts]

In[]:=

Stocks

#	Symbol	Description
1	TP[t]	Total Population
2	SP[t]	Susceptible Population
3	EP[t]	Exposed Population
4	INSP[t]	Infected Normally Symptomatic Population
5	ISSP[t]	Infected Severely Symptomatic Population
6	RP[t]	Recovered Population
7	MLP[t]	Money of Lost Productivity

,Rates

#	Symbol	Description
1	μ[TP]	Population death rate
2	μ[INSP]	Infected Normally Symptomatic Population death rate
3	μ[ISSP]	Infected Severely Symptomatic Population death rate
4	sspf[SP]	Severely Symptomatic Population Fraction
5	β[INSP]	Contact rate for the normally symptomatic population
6	β[ISSP]	Contact rate for the severely symptomatic population
7	aip	Average infectious period
8	aincp	Average incubation period
9	lpcr[ISSP,INSP]	Lost productivity cost rate (per person per day)

,Equations

#	Equation
1	′ SP t- INSP[t] SP[t] β[INSP] TP[t] - ISSP[t] SP[t] β[ISSP] TP[t] - SPt μTP
2	′ EP t INSP[t] SP[t] β[INSP] TP[t] + ISSP[t] SP[t] β[ISSP] TP[t] - EPt 1 aincp + μTP
3	′ INSP t- INSP[t] aip + EP[t] (1- sspf[SP] ) aincp - INSPt μINSP
4	′ ISSP t- ISSP[t] aip + EP[t] sspf[SP] aincp - ISSPt μISSP
5	′ RP t INSP[t] + ISSP[t] aip - RPt μTP
6	′ MLP [t] lpcr[ISSP,INSP] (- RP[t] - SP[t] + TP[t] )
7	TP[t] Max[0, EP[t] + INSP[t] + ISSP[t] + RP[t] + SP[t] ]

,RateRules

#	Symbol	Value
1	μTP	1 45625
2	μISSP	0.035 aip
3	μINSP	0.01 aip
4	β[ISSP]	0.15
5	β[INSP]	0.15
6	aip	26
7	aincp	6
8	sspf[SP]	0.2
9	lpcr[ISSP,INSP]	1

,InitialConditions

#	Equation
1	SP[0]99998
2	EP[0]0
3	ISSP[0]1
4	INSP[0]1
5	RP[0]0
6	MLP[0]0
7	TP[0]100000



Out[]=

Here is SEI2HR:

modelSEI2HR=SEI2HRModel[t,"InitialConditions"True,"RateRules"True,"TotalPopulationRepresentation"reprTP];ModelGridTableForm[modelSEI2HR,lsModelOpts]

In[]:=

Stocks

#	Symbol	Description
1	TP[t]	Total Population
2	SP[t]	Susceptible Population
3	EP[t]	Exposed Population
4	INSP[t]	Infected Normally Symptomatic Population
5	ISSP[t]	Infected Severely Symptomatic Population
6	RP[t]	Recovered Population
7	MLP[t]	Money of Lost Productivity
8	HP[t]	Hospitalized Population
9	DIP[t]	Deceased Infected Population
10	HB[t]	Hospital Beds
11	MHS[t]	Money for Hospital Services

,Rates

#	Symbol	Description
1	μ[TP]	Population death rate
2	μ[INSP]	Infected Normally Symptomatic Population death rate
3	μ[ISSP]	Infected Severely Symptomatic Population death rate
4	sspf[SP]	Severely Symptomatic Population Fraction
5	β[INSP]	Contact rate for the normally symptomatic population
6	β[ISSP]	Contact rate for the severely symptomatic population
7	aip	Average infectious period
8	aincp	Average incubation period
9	lpcr[ISSP,INSP]	Lost productivity cost rate (per person per day)
10	μ[HP]	Hospitalized Population death rate
11	β[HP]	Contact rate for the hospitalized population
12	nhbr[TP]	Number of hospital beds rate
13	hscr[ISSP,INSP]	Hospital services cost rate (per bed per day)
14	nhbcr[ISSP,INSP]	Number of hospital beds change rate (per day)

,Equations

Equation

TP[t]

Max[0,

EP[t]

INSP[t]

ISSP[t]

RP[t]

SP[t]

]

′

t-

HP[t]

SP[t]

β[HP]

TP[t]

INSP[t]

SP[t]

β[INSP]

TP[t]

Max[0,-

HP[t]

ISSP[t]

]

SP[t]

β[ISSP]

TP[t]

SPt

μTP

′

t

HP[t]

SP[t]

β[HP]

TP[t]

INSP[t]

SP[t]

β[INSP]

TP[t]

Max[0,-

HP[t]

ISSP[t]

]

SP[t]

β[ISSP]

TP[t]

EPt

aincp

μTP

′

INSP

t-

INSP[t]

aip

EP[t]

(1-

sspf[SP]

)

aincp

INSPt

μINSP

′

ISSP

t-

ISSP[t]

aip

EP[t]

sspf[SP]

aincp

HPt

μHP

HPt

ISSPt

μISSP

′

t

Min HBt - HPt , EP[t] sspf[SP] aincp 	HPt < HBt
0	True

HP[t]

aip

HPt

μHP

′

t

INSP[t]

ISSP[t]

aip

RPt

μTP

′

DIP

[t]

HP[t]

μ[HP]

INSP[t]

μ[INSP]

+(-

HP[t]

ISSP[t]

)

μ[ISSP]

′

[t]

HB[t]

nhbcr[ISSP,INSP]

′

MHS

[t]

HP[t]

hscr[ISSP,INSP]

′

MLP

[t]

lpcr[ISSP,INSP]

(

INSP[t]

ISSP[t]

HP[t]

μ[HP]

INSP[t]

μ[INSP]

+(-

HP[t]

ISSP[t]

)

μ[ISSP]

)

,RateRules

#	Symbol	Value
1	μTP	1 45625
2	μISSP	0.035 aip
3	μINSP	0.01 aip
4	β[ISSP]	0.15
5	β[INSP]	0.15
6	aip	26
7	aincp	6
8	sspf[SP]	0.2
9	lpcr[ISSP,INSP]	1
10	μ[HP]	0.25 μ[ISSP]
11	β[HP]	0.1 β[ISSP]
12	nhbr[TP]	0.0029
13	nhbcr[ISSP,INSP]	0
14	hscr[ISSP,INSP]	1

,InitialConditions

#	Equation
1	SP[0]99998
2	EP[0]0
3	ISSP[0]1
4	INSP[0]1
5	RP[0]0
6	MLP[0]0
7	TP[0]100000
8	HP[0]0
9	DIP[0]0
10	HB[0] nhbr[TP] TP[0]
11	MHS[0]0



Out[]=

Here are the “differences” between the two models:

ModelGridTableForm@Merge[{modelSEI2HR,modelSEI2R},If[AssociationQ[#〚1〛],KeyComplement[#],Complement@@#]&]

In[]:=

Stocks

#	Symbol	Description
1	HP[t]	Hospitalized Population
2	DIP[t]	Deceased Infected Population
3	HB[t]	Hospital Beds
4	MHS[t]	Money for Hospital Services

,Rates

#	Symbol	Description
1	μ[HP]	Hospitalized Population death rate
2	β[HP]	Contact rate for the hospitalized population
3	nhbr[TP]	Number of hospital beds rate
4	hscr[ISSP,INSP]	Hospital services cost rate (per bed per day)
5	nhbcr[ISSP,INSP]	Number of hospital beds change rate (per day)

,Equations

Equation

′

DIP

[t]

HP[t]

μ[HP]

+INSP[t]μ[INSP]+(ISSP[t]-

HP[t]

)μ[ISSP]

′

t

SP[t]

HP[t]

β[HP]

TP[t]

INSP[t]SP[t]β[INSP]

TP[t]

Max[0,ISSP[t]-

HP[t]

]SP[t]β[ISSP]

TP[t]

-EPt

aincp

+μTP

′

[t]

HB[t]

nhbcr[ISSP,INSP]

′

t

Min EP[t]sspf[SP] aincp , HBt - HPt 	HPt < HBt
0	True

HP[t]

aip

HPt

μHP

′

ISSP

t-

ISSP[t]

aip

EP[t]sspf[SP]

aincp

HPt

μHP

-ISSPt-

HPt

μISSP

′

MHS

[t]

HP[t]

hscr[ISSP,INSP]

′

MLP

[t]lpcr[ISSP,INSP](INSP[t]+ISSP[t]+

HP[t]

μ[HP]

+INSP[t]μ[INSP]+(ISSP[t]-

HP[t]

)μ[ISSP])

′

t-

SP[t]

HP[t]

β[HP]

TP[t]

INSP[t]SP[t]β[INSP]

TP[t]

Max[0,ISSP[t]-

HP[t]

]SP[t]β[ISSP]

TP[t]

-SPtμTP

,RateRules

#	Symbol	Value
1	μ[HP]	0.25μ[ISSP]
2	β[HP]	0.1β[ISSP]
3	nhbr[TP]	0.0029
4	nhbcr[ISSP,INSP]	0
5	hscr[ISSP,INSP]	1

,InitialConditions

#	Equation
1	DIP[0]0
2	HB[0] nhbr[TP] TP[0]
3	HP[0]0
4	MHS[0]0



Out[]=

Equations explanations

In this section we provide rationale of the model equations of SEI2HR.

The equations for Susceptible, Exposed, Infected, Recovered populations of SEI2R are "standard" and explanations about them are found in [WK1, HH1]. For SEI2HR those equations change because of the stocks Hospitalized Population and Hospital Beds.

The equations time unit is one day. The time horizon is one year. Since we target COVID-19, [Wk2, AA1], we do not consider births.

Remark: For convenient reading the equations in this section have tooltips for the involved stocks and rates.

Verbalization description of the model

We start with one infected (normally symptomatic) person, the rest of the people are susceptible. The infected people meet other people directly or get in contact with them indirectly. (Say, susceptible people touch things touched by infected.) For each susceptible person there is a probability to get the decease. The decease has an incubation period: before becoming infected the susceptible are (merely) exposed. The infected recover after a certain average infection period or die. A certain fraction of the infected become severely symptomatic. If there are enough hospital beds the severely symptomatic infected are hospitalized. The hospitalized severely infected have different death rate than the non-hospitalized ones. The number of hospital beds might change: hospitals are extended, new hospitals are build, or there are not enough medical personnel or supplies. The deaths from infection are tracked (accumulated.) Money for hospital services and money from lost productivity are tracked (accumulated.)

The equations below give mathematical interpretation of the model description above.

Code for the equations

Each equation in this section are derived with code like this:

ModelGridTableForm[modelSEI2HR,lsModelOpts]["Equations"][[1,EquationPosition[modelSEI2HR,RP]+1,2]]

and then the output cell is edited to be “DisplayFormula” and have CellLabel value corresponding to the stock of interest.

The infected and hospitalized populations

SEI2HR has two types of infected populations: a normally symptomatic one and a severely symptomatic one. A major assumption for SEI2HR is that only the severely symptomatic people are hospitalized. (That assumption is also reflected in the diagram in the introduction.)

Each of those three populations have their own contact rates and mortality rates.

Here are the contact rates from the SEI2HR dictionary

ColumnForm@Cases[Normal@modelSEI2HR["Rates"],HoldPattern[β[_]_],∞]

In[]:=

β[INSP]Contact rate for the normally symptomatic population

β[ISSP]Contact rate for the severely symptomatic population

β[HP]Contact rate for the hospitalized population

Out[]=

Here are the mortality rates from the SEI2HR dictionary

ColumnForm@Cases[Normal@modelSEI2HR["Rates"],HoldPattern[μ[_]_],∞]

In[]:=

μ[TP]Population death rate

μ[INSP]Infected Normally Symptomatic Population death rate

μ[ISSP]Infected Severely Symptomatic Population death rate

μ[HP]Hospitalized Population death rate

Out[]=

Remark: Below with “Infected Population” we mean both stocks Infected Normally Symptomatic Population (INSP) and Infected Severely Symptomatic Population (ISSP).

Total Population

In this notebook we consider a DAE’s formulation of SEI2HR. The stock Total Population has the following (obvious) algebraic equation:

TP[t]

Max[0,

EP[t]

INSP[t]

ISSP[t]

RP[t]

SP[t]

]

(TP)

Note that with Max we specified that the total population cannot be less than

Remark: As mentioned in the introduction, the package [AAp1] allows for the use of non-algebraic formulation, without an equation for TP.

Susceptible Population

The stock Susceptible Population (SP) is decreased by (1) infections derived from stocks Infected Populations and Hospitalized Population (HP), and (2) morality cases derived with the typical mortality rate.

′

[t]-

HP[t]

SP[t]

β[HP]

TP[t]

INSP[t]

SP[t]

β[INSP]

TP[t]

Max[0,-

HP[t]

ISSP[t]

]

SP[t]

β[ISSP]

TP[t]

SP[t]

μ[TP]

(SP)

Because we hospitalize the severely infected people only instead of the term

ISSP[t]SP[t]β[ISSP]

TP[t]

we have the terms

HP[t]

SP[t]

β[HP]

TP[t]

Max[0,-

HP[t]

ISSP[t]

]

SP[t]

β[ISSP]

TP[t]

The first term is for the infections derived from the hospitalized population. The second term for the infections derived from people who are infected severely symptomatic and not hospitalized.

Births term

Note that we do not consider in this notebook births, but the births term can be included in SP’s equation:

Block[{m=SEI2HRModel[t,"BirthsTerm"True]},ModelGridTableForm[m]["Equations"]〚1,EquationPosition