Abstract
A mathematical model of COVID-19 is presented where the decision to increase or decrease social distancing is modelled dynamically as a function of the measured active and total cases as well as the perceived cost of isolating. Along with the cost of isolation, we define a healthcare cost and a total cost. We explore these costs by adjusting parameters that could change with policy decisions. We observe that minimum costs are not always associated with increased spending and increased vigilance which is due to the desire for people to not distance and the fatigue they experience when they do. We demonstrate that an increased in the number of lock-downs, each of shorter duration can lead to minimal costs. Our results are compared to case data in Ontario, Canada from March to August 2020 and details of extracting the results to other regions is presented.
1. Introduction
As of early October 2020 there have been over 33 million cases of COVID-19 worldwide, over 150 000 cases in Canada, and over 50 000 cases in the province of Ontario. The early stages of the outbreak focused on mathematical modelling of disease dynamics such as transmission and the basic reproduction number [1,2]. It quickly became clear that asymptomatic spreading was important and that undetected infections were important to consider in models [14].
This caused a global policy shift towards travel restrictions, community closures, and social distancing implementations. The impacts of mathematical modelling on policy are documented in [20].
The implementation of non-pharmaceutical intervention (NPI) such as social distancing quickly became an important mathematical modelling task (cf. [16,17,19,31]). The majority of these models focus on fixed policy implementations such as reducing contacts on a given date and reinstating them on another. There are two main issues with this, the first is that it requires knowledge of the implementation and relaxation times. While this can be explored in model simulations and optimized for best results, its independence from the model itself can make it hard to adapt to other diseases, strains, or important factors. A second problematic issue is that it assumes a instantaneous policy compliance, i.e., that people will immediately reduce contacts upon implementation and stop upon relaxation. While this can be impacted by an adherence parameter, it does not allow for a dynamic response which is more realistic of human choice. Therefore a dynamics social distancing model that reacts to the disease dynamics is more realistic.
A dynamic intervention strategy where intervention was turned on and off based on the state of the epidemic was considered in [31] where a decrease in both total infections and social distancing duration was observed compared to a fixed-duration intervention which they also considered. However, modelling the dynamics of intervention entirely on the disease progression assumes that people will immediately distance or relax at some threshold. This suggests that a periodic solution will emerge centered around the critical disease threshold and this appears to happen in [31]. While it is quite realistic that disease dynamics drive people into isolation, it is a separate mechanism, namely the cost of staying home, that people consider when relaxing their isolation habits. Cost is seldom considered in models, with most of the focus on larger economic influence [4,22]. These economic factors certainly play a role in individual cost but psychological factors such as loneliness and habit displacement are important as well
For this paper we propose a differential equation model for the spread of COVID-19 with separate dynamics for isolation and relaxation dependent on disease progression and relaxation cost respectively. The disease progression information typically comes from media reports and has been investigated in the context of infectious diseases such as influenza (cf. [7,8,30,33]) and is usually used to reduce the susceptibility of individuals who are positively influenced by media. The relaxation cost is less often considered and its inclusion recognizes that repeated lock-downs would have diminishing returns as the cost to stay home becomes too overwhelming. A dynamic response model allows for more realistic policy strategies for disease mitigation and mortality prevention. Our model focuses on the spread of the disease in Ontario, Canada, but could be adapted with other parameters to other regions.
Our study is outlined as follows. In section 2 we introduce the model and the dynamic response functions for social distancing and relaxation. We also introduce the parameters including those which we fit to data from [25]. We define health, economic, and total costs of the pandemic. The health cost is based on overloading existing healthcare resources while the economic cost is the personal or societal cost of social distancing. We show the excellence of fit to our data in section 3 and present a series of results based on different scenarios where policy parameters that control distancing and relaxing are varied. We consider scenarios where both health and relaxation costs are equally weighted or where health cost is much more strongly influencing the total cost. We consider a modification to the relaxation rate so that it depends on both cost and cases and see that multiple outbreak peaks can occur. We discuss the implications and conclusions of our work in section 4.
2. Model
We consider a mathematical model for COVID19 consisting of classes of people with various exposure to the disease. These classes are listed in Table 1where we note that the removed groups include people who have died from the virus which we do not separately consider. For each of the population classes we assume there are three levels of social distancing indicated by a variable subscript zero, one, or two. If the subscript is 0 then there is no social distancing, subscript 1 indicates that there is social distancing which reduces the contact probability by some percentage while for subscript 2, the contact probability is zero, i.e. full isolation. We introduce a further subscript, M which represents the mitigation on spread due to individuals who have tested positive and are isolated. We assume that only P, IS, and IA populations can test positive and that these people will immediately and completely isolate effectively placing them in the social distance 2 category for the duration of their disease.
We follow the usual SEIR model framework (cf. [13,15]) which we illustrate in Figure 1 with equations detailed in Appendix A.
For the model we assume that vital statistics are not important on the time scales we consider so we take a fixed population N. We also normalize the model by another population Ncrit which is the amount of people needing healthcare resources that puts the system at full capacity. Figure 1 also shows the various model parameters which are summarized in Table 1. Each transition parameter has a subscript (i, j) with i the originating class and j the terminal class. The exception to this is βi,j where i is the class of the susceptible person and j the class of the infected contact. Whenever there are multiple paths for a class, a parameter qi ∈ [0, 1] is used to partition them into each final compartment.
We make the following assumptions about the model parameters:
the parameters δ, σ, ϕ, γ, and qi are constant and the same for each social distancing class as the disease progression characteristics are unaffected by social distancing.
a constant which incorporates both contact and disease transmission probability. We assume that people not showing symptoms shed a lower viral load and hence reduce transmission by a constant factor α and that those in social distancing class 1 reduce their contacts by a constant factor of δ, effectively also reducing their transmission. For example and .
people in the infected symptomatic class choose to isolate at a constant rate µI with qI going into and (1 − qI) going into . They stay in the social distancing class until they have recovered from the disease, i.e., νI = 0. Furthermore, this means that someone already social distancing in state 1 or 2 who becomes symptomatic remains in that social distance class. We note that individuals in , and know that they are sick but have not tested positive for the disease. If they test positive, they transition to and are completely quarantined.
there are two testing rates ρA and ρS for asymptomatic (including presymptomatic) and symptomatic individuals respectively with ρS > ρA as we assume that symptomatic people are more likely to seek out a test as they would have symptoms. Asymptomatic people are likely to only seek a test out of curiosity or if they have been in contact with someone who has tested positive. Despite the fact that testing numbers fluctuate with the progression of the disease we take the testing rates to be constant which makes them an effective average testing rate. This is consistent with studies estimating global infections per symptomatic test case (cf. [18,26]). To help restrict the model we take ρS = 4ρA which is a similar value as observed in [18] who compared data to Germany, South Korea, and the USA.
we assume that only people in the P, IS or IA classes will test positive if a test is administered. Therefore, we explicitly assume people in the E class do not have a high enough viral load to shed.
people who have tested positive are isolated (effectively put in social class 2) until recovery and cannot transmit the disease. This ignores infections to family members living in a household with an isolated member or infections to healthcare workers who are conducting tests or treating COVID-19 patients. See [34] for considerations of a model with household structure included.
(a) Social Distancing and Testing
Since people without symptoms are unaware of whether they have the virus or not, we assume that both social distancing and relaxing rates are independent of the disease class they are in. We therefore define µ as the rate of social distancing from state 0 with proportion q0 going to state 1 and (1 − q0) going to state 2. We similarly define ν as the rate of decreasing social distance from state 2 with proportion q2 going to state 1 and (1 − q2) going to state 0. We define those social distancing from state 1 to state 2 as µ/2 to account for the fact that anybody in state 1 has already undergone one transition and so they should be slower at making a secondary transition. For a similar reason we define the social distancing relaxation from state 1 to state 0 as ν/2.
Testing provides two important quantities reported by the media that can help inform social distancing, the total number of cases M and the active cases MA (each also scaled by Ncrit) which are defined by
If the disease is in the exponential phase of spread then the doubling rate can be deduced from the cumulative case information, M, to yield
We assume that this can approximate the doubling rate at all times and is what locally drives social distancing. However, while the disease growth rate is important, it should be weighed against the number of active cases as well and therefore we propose a social distance transition function where µmax is the maximal rate of social distancing, [·]+ is defined such that,
Social distancing is not something that people want to do and the parameters kc and Mc represent critical doubling rates and active case numbers respectively below which people will not social distance which is the role of the maximum function. These parameters can be thought of as policy parameters since implementing lock-downs, closing businesses, and halting social gatherings will impact these values. The parameters k0 and M0 represent doubling rates and case numbers where social distancing reaches its half-maximum.
We assume that social relaxation, ν, is proportional to the cost of social distancing, C. This may be economic cost in the form of people staying home from their jobs, but also the psychological cost of being isolated for a long period of time. We model the cost as, where we only include susceptible and exposed classes because we assume there is a greater benefit to having transmitting classes (Pi, Ii) stay home. Arguably, exposed people who will soon become infectious should stay home too, but as they would test negative, they would think they are healthy and therefore we assume they contribute to the cost. For simplicity we ignore the cost of recovered people social distancing which will be invalid if many people have recovered but policy prevents them from returning to their workplaces etc. As written, the cost accumulates with time. We could remedy this by including a decay factor −µCC in (2.4) but we assume that the time scale of recovery is much longer than that of the pandemic. Since the cost does not factor in day-to-day economic costs in non-pandemic times, it is normalized so that zero cost represents the cost of society pre-pandemic. Similarly then the maximum additional costs come from those isolating completely in social distance state 2. Those in social distance state 1 effect their transmissibility by a factor δ and we assume this comes at a reciprocal burden cost of (1 − δ). For example, if δ = 1 then the 1 and 2 states are both fully isolated and contribute an equal maximal cost. We note that since populations are non-dimensional in the model, the cost is denoted in days where one day represents the cost of every person being in social distance state 2.
Defining the cost, we propose ν be modelled by where νmax is the maximal rate at which social distancing can be relaxed, C0 is the cost which triggers the half-maximal rate, and Cc is the cost required to trigger social relaxation. Cc is also a policy parameter as mental health promotion, economic stimulus, and wage subsidy programs can influence the cost people can endure before social relaxation.
To understand the true cost of the pandemic we must balance the relaxation cost, C, with the healthcare cost, H, which we define as where t0 is the time where active cases begin to strain the healthcare system and t1 is the time they stop. This measures the burden on the healthcare system, and likely on mortality rates, if active cases get too large. We have defined Ncrit as the value above which healthcare resources are exceeded and we assume the burden on the system begins when Ncrit/2 is reached. The choice of defining H balances intensity of the outbreak along with duration. The sum allows for multiple outbreaks where the hospital resources are exceeded. The reason we measure active cases is that we assume all COVID-19 cases entering the hospital will be tested. Realistically a portion of the untested symptomatic cases will also exasperate the healthcare system and therefore this can be considered an underestimated cost. Having defined the healthcare cost we can then define the total cost as where H∞ is the healthcare cost with no social distancing intervention (µmax = 0) and C∞ is the largest isolation cost allowable. We define ω as a weighting factor between the two cost contributions. In the absence of actual dollar amounts, the total cost as defined in (2.7) is a relative cost.
(b) Parameter Determination
We first consider a variant of the full model that does not include social distancing or testing (green-dash rectangle of Figure 1) which is given by,
This reduced model represents the disease transmission dynamics prior to widespread knowledge of COVID-19. Following [10] we assume that people without symptoms are half as infectious as those with symptoms and therefore take α = 1/2.
The disease free state is [S, E, P, IS, IA, RS, RA] = [N/Ncrit, 0, 0, 0, 0, 0, 0] and we identify the basic reproduction number R0 as the non-zero eigenvalue of the next generation matrix produced from (2.8) (cf. [9,12]),
Taking R0 from measurements such as the studies in [10,24] which estimate R0 = 2.4 we can rearrange (2.9) to determine a value for β.
We assume that this base transmission rate between the susceptible and symptomatic populations is the same as that between S0 and in the social distancing model (A 1), i.e. . The parameters considered for the base model are presented in Table 1 and we comment on some of the assumptions made.
As of July 16, 2020, the hospitalization rate of COVID-19 in Ontario Canada was 12.3% and there were approximately 10000 hospital beds available for people which together define Ncrit = 81301. We choose µmax = νmax = 1 under the assumption that people generally need at least one day to change their routines. We arbitrarily assume that C0 = 2Cc and that k0 = 4kc to help constrain the model. This means that the relaxation cost needed to initiate the half-maximum rate is twice as many days as the onset of social relaxation while the disease needs to double twice for the half-maximal social distancing rate to occur. We chose Cc = 50 based on Ontario imposing a stage-one lock-down that lasted almost 100 days with the fact that it did not impact the entire province.
We predicted the values of kc, Mc, and ρA (recalling that ρS = 4ρA) by fitting our model to active and total case data from [25] between March 10, 2020 August 18, 2020 inclusive. We used a non-linear least squares method for the fitting, the details of which are in Appendix B. Using the Ncrit scale we can convert the values of Mc and M0 from Table 1 to 2090 and 4180 people respectively. The values of µI, q0, q2, and qI are arbitrarily chosen. However, as is seen in Appendix B where a sensitivity analysis is performed, these parameters are not very influential on model results so their values are not likely very consequential on model results. The most influencing parameter is µI. Considering that at the half-maximal rates, the social distancing rate is µ = 1/4 and that a recent study from [3] suggested that up to 90% of Americans go into work sick then a further 90% reduction would yield µI ≈ 0.025 which is the same order of magnitude to the chosen value.
We took a median value for the symptomatic rate of 69% following a variety of studies (see [6,21,23,27–29]).
3. Results
We simulated (A 1) using MATLAB 2020a with parameters in Table 1. We took March 10, 2020 as the initial time with an initial condition that 0.02% of the population was infected with symptoms and placed the remaining 99.98% of the population in the susceptible class.
We demonstrate the results from data-fitting the parameters kc, Mc, and ρA in Figure 2. Comparing data to simulation, we observe a difference in the early peak-time of 4.7 days and a difference in peak active cases of 359 people.
We observe the impact of the disease on total cost (2.7) by simulating the full model (A 1) and varying the critical threshold at which people social distance (Mc) and the critical cost before social relaxation begins (Cc). We consider 1/4, 1/2, 2, and 4 times the base values given in Table 1. The value for kc from data fitting is already quite extreme and we do not vary this. We plot heat maps for H, C, and CT given by (2.6), (2.4), and (2.7) respectively in Figure 3 where they are coloured relative to the maximal and minimal costs. We compute H∞ by simulating the model with parameters in Table 1 and taking µmax = 0. We take C∞ as the highest cost that emerges from all of the simulations. Intuitively this occurs at and as evidenced in Figure 3b. We consider an even weighting between H and C (ω = 1/2) in Figure 3c as well as a 90:10 H: C weighting (ω = 0.9) in Figure 3d.
The overall trend, particularly with even weighting, is unsurprising, namely that increased vigilance (larger Mc) and increased spending (larger Cc) contribute to a larger cost. However, there are some important things to note. Firstly, this trend is not always true. For example in Figure 3a at , the health cost increases when transitioning from to , i.e. doubling the critical cost allowing people to stay home longer, causes a worse health outcome. We plot the active cases for these two scenarios in Figure 4 where we see that indeed spending more adds a delay to the peak time and decreases the peak value, but the duration above the critical health threshold of Ncrit/2 is three days longer. This added time impact leads to a larger overall health cost.
A second interesting feature of note in the health cost is that the decision of whether to increase vigilance or the critical cost is not always clear. For example if the vigilance is at and the cost at in Figure 3a then it is more beneficial (less relative health cost) to double vigilance than to double critical cost. However starting at and it is more beneficial to double critical cost than double vigilance. An extreme example is at and where doubling either critical cost or vigilance each lead to higher healthcare costs. This is important to note from a policy standpoint because it can look like both policy initiatives are failures, however it is only stating that the improvement measures were not sufficient enough and more needs to be done. For example quadrupling the critical cost does reduce the healthcare cost. This same non-uniform decision process is not true for social relaxation cost, C, as spending more always increases the cost and increasing vigilance always decreases the cost. This trend persists for the total cost when both are equally weighted, however in the 90/10 weighting, a weighting that may be more consistent with prioritizing human lives, a richer structure emerges. For example, in Figure 3d starting at and it is more costly to both half and double spending indicating that an optimal spending level is met at that vigilance. As well, at and it is more costly to both double and half vigilance indicating a optimal vigilance balance there. The global minimum in the parameter space considered is at and where decreasing vigilance as well as increasing or decreasing critical cost leads to overall more costly outcomes. Non-uniform trends and local minima in costs are important considerations that should be made in policy decisions.
A third important feature to recognize in Figure 3 is that even when intuitive trends of cost increases are followed, the relative increase in cost can be quite different. For example in Figure 3d starting at and doubling the critical cost is a 1.3% increase in total cost but doubling the critical cost again is a 3.6% increase in total cost. For Figure 3c starting at a high vigilance of with and doubling critical cost results in a 19.7% increase in total cost, but doubling it again results in a 33.8% increase in total cost. A lower relative increase in cost indicates that the strategy was more effective compared to others.
Plots of all of the active cases and costs for the scenarios in Figure 3 can be found in the supplementary material. The cumulative number of symptomatic-infected people for each scenario is plotted in Figure 5.
It can be seen that decreasing vigilance causes a more uniform accumulation of cases while increasing relaxation cost delays the peak of infection. In all cases, the same number of total people are infected as in the baseline case where testing occurs but no social distancing happens. This is expected as complacency and fatigue from non-pharmaceutical interventions eventually force the cost of social distancing to be too high for people to remain away from others. However, these delays can provide time for vaccination and other medical efforts to minimize the impact of the disease.
(a) Multiple Secondary Waves
The model as derived only allows for one large secondary wave following the peak in Ontario around early May 2020. Since the relaxation rate ν is solely a function of relaxation cost (2.4) which is always increasing, isolation fatigue becomes too overwhelming that there is resistance for entering secondary lock-downs. This model is likely appropriate for regions that have a strong aversion to social distancing. For other regions, it is likely that relaxation will be a function of cost and active cases as people will prioritize their health in a sustained outbreak and thus not want to relax if case numbers are sufficiently large. If we refer to the rate in (2.5) as ν0 then we propose modifying ν to where η is the factor of critical cases Mc that stops social relaxation regardless of cost. Implementing this change allows for secondary infection peaks as evidenced in Figure 6 where we use parameters in Table 1 and arbitrarily take η = 1/2 for Figure 6a and η = 1/5 for Figure 6b. It is important to note that these changes do not impact the initial peak fitted to data in Figure 2 and only alter future projections. For this reason, it is difficult to estimate η as several peaks will need to have occurred.
We repeat the cost analysis as in Figure 3 for the modified relaxation cost (3.1), however we fix , the value from Table 1 and instead modify η. The results are presented in Figure 7.
Figure 7b shows a different result compared to the case of Figure 3b when (2.4) was used for the relaxation. In the latter case, there was a general trend upward in cost that had little difference between the value of Mc, however in Figure 7b there is a strong dependence in η. This is because the multi-secondary outbreaks caused by reducing relaxation with high active cases extends the duration of the epidemic which only increases the cumulative cost. The high impact of this is noticed as well with C∞ = 230 days for Figure 3 while C∞ = 614 days for Figure 7. H∞ remains the same in both cases since that is calculated with no social distancing at all (and therefore no social relaxation).
The introduction of the modified relaxation cost (3.1) has an impact on the health cost as seen in Figure 7a where a seemingly strange result is that increasing spending leads to worse health outcomes. This is particularly interesting for a case such as η = 1/4 and as there is no healthcare cost, i.e. the tested active cases never surpasses the critical threshold. However, spending four times as much to encourage people to stay home leads to a significantly higher healthcare cost. The rationale for this is similar to what was observed in [31] where keeping people isolated for a longer duration increases their fatigue and resistance to staying isolated in future instances. Figures 7c and 7d demonstrate the impact of this fatigue since regardless of weighting to health cost, the cheapest total cost options are with minimal encouragement to isolate. However, for a given cost, a higher reactivity (smaller η) to active cases does not lead to a minimum total cost. This is because there is a critical value of η below which no additional healthcare savings occur but increasing expenses occur for relaxation. These results suggest a careful policy direction with more lock-downs of shorter duration.
The plots of cases and total infections for each of the scenarios in Figure 7 can be found in the supplementary material. We plot the cumulative symptomatic-infected proportions for each scenario in Figure 8.
Unlike the scenarios in Figure 5 associated with the relaxation cost (2.4), basing relaxation on active cases as well can impact the terminal number of cumulative infections. Continually reducing η decreases the total number of infected people which provides more evidence that increasing the relaxation cost threshold Cc can cause more people to become infected. However similar to Figure 5, increasing Cc leads to a longer delay before significant infection numbers occur.
4. Conclusion
We have presented a model for COVID-19 that allows for dynamic social distancing and relaxation based on the measured active cases and individual cost of isolating. The aim of this approach is that it more accurately reflects human behaviour and psychology unlike the modelling approach where behaviours are turned on and off at pre-determined times. Understanding how people will react to a change in policy surrounding lock-downs or bans on social gatherings is essential in gauging the impact that COVID-19 and mitigation strategies will have on infections and mortality. Improving this modelling aspect can make sure that policies are put into place at the right time so people will react accordingly.
By modelling behaviour dynamically we were able to produce unintuitive results regarding the relative total cost of the disease, namely that increasing vigilance and relaxation cost as well as decreasing reactivity to growing case loads does not necessarily lead to a decrease in total cost. This is because of the desire for people to socialize and fatigue for being intrinsically part of the model. We have demonstrated that in certain circumstances the additional healthcare cost can be eliminated entirely, but at a relatively large cost to keeping people home.
An advantage of the dynamic framework used in this model is that it is not restricted to Ontario nor is it even restricted to COVID-19. Changing the disease and behaviour parameters will allow this model to adapt to other scenarios. For COVID-19, policy makers would be advised to use data such as that used here from [25] in the relatively early stages of a lock-down to fit behaviour parameters. Earlier time point data helps reduce the likelihood that the relaxation cost threshold has been exceeded so that the behaviour parameters are more accurate. Otherwise, and become stronger functions of the choice for . The limited data and types of data available should discourage too much parameter fitting. Having determined the parameters to a given set of data, cost analysis using (2.7) can be done leading to results similar to Figures 3 and 7. Understanding the influence of tangible actions such as forced closures, wage subsidies, etc. on parameters such as Mc, Cc, and η likely requires surveys and other followup studies.
It is important to acknowledge that this model does not take into account vaccination or other pharmaceutical interventions. These have an important role in not only limiting the healthcare impact, but also in outbreak peak time and duration. The introduction of the modified cost in (3.1) causes significant delays between peaks at increased economic and psychological cost. The duration of the pandemic can then be several years longer than when fatigue is strong resulting in a single outbreak peak. The duration may be significantly longer than the time for an effective vaccine to be developed and deployed and this needs to be considered in future work. Another important consideration is that social distancing is not truly discrete in that people do not suddenly reduce their contacts. In reality it is a spectrum with fluid contact rates and this needs to be further explored. Finally, taking the testing rates constant should be relaxed and reflect that both the testing capacity and willingness for individuals to test is a function of disease progression.
Data Availability
All data and code is available upon request
Funding
IRM acknowledges funding from an NSERC Discovery grant (2019-06337). JMH acknowledges funding from an NSERC Discovery grant and Discovery Accelerator Supplement.
A. Differential Equation Model
The differential equation model visualized in Figure 1 is given by where is the force of infection,
In deriving the model we have normalized the population by Ncrit which represents the population that causes the healthcare system to be at capacity.
B. Parameter Fitting and Model Sensitivity
Many of the parameters associated with the natural progression of the disease are unknown as are several of the intervention parameters such as social distancing and testing. We use data from [25] for the 161 days between March 10 and August 18, 2020 on active and total cases to elucidate some key parameters. To do this we solve our model (A 1) with the parameters in Table 1 excluding kc, k0, Mc, M0, ρA and ρS which we fit to the data. We run our model for 161 days using an initial condition that 0.02% of the population was initially infected with the remaining 99.98% being symptomatic. The actual number of people with COVID-19 is a matter of speculation and the arbitrary choice of the initial value will affect the fitting parameters, particularly the testing rates which are intimately linked. We do not partition any of the initial infected population into symptomatic or asymptomatic also due to the lack of clarity on true numbers. To help constrain the model we take k0 = 4kc so that the doubling rate has to double twice to trigger the half-maximal social distancing. We also take M0 = 2Mc so that the number of active cases need to double to trigger the half-maximal social distancing rate. Finally, since we assume that symptomatic people are more likely to get tested than asymptomatic people we take ρS = 4ρA. These constraints should not be too restricting, likely impacting the fitted values of the remaining free variables kc, Mc, and ρA.
We use a non-linear least squares iterative procedure to identify the parameters. This leads to the values kc=1/16.24 d−1, Mc = 2.57 × 10−2, and ρA = 8.7 × 10−3 d−1 as in Table 1 with a residual norm of 9.5 × 10−3. A graphical representation of the fit between data and model is in Figure 2.
(a) Sensitivity Analysis
Having fit data to the model we then performed a sensitivity analysis on the parameters ρS, ρA, qI, µI, q2, q0, q, γ, ϕ, σ, and α. We used Table 1 for the fixed parameters and as mean values for the varying parameters, each distributed uniformly between maximum and minimum values. We used the latin hypercube sampling technique with 10000 iterations and a Spearman partial rank correlation coefficient to measure monotonicity. We tested the sensitivity to the cumulative infected (symptomatic and asymptomatic), susceptibles, peak time for the outbreak, and the value at the peak of the outbreak. We plot the results in Figure 9
The most significant parameter is the recovery rate γ (assumed the same for both classes) which seemingly has inverse behaviour to what is expected. That is an increase in the recovery rate seems to cause more people to become sick. This is because the basic reproduction number is fixed at 2.4 and therefore changing γ effectively changes the transmission β making a higher recovery rate lead to a more transmissible disease. Aside from this, the most significant parameters unsurprisingly are the testing rates ρS and ρA as well as the symptomatic proportion q. This supports the importance of testing and social distancing. Interestingly, there is not much sensitivity to the peak time and value of the outbreak confirming the need for long-term planning regarding vaccination and hospital resource management.
We did not perform a sensitivity analysis on kc or Mc as these are policy parameters. Their sensitivity is effectively measured by comparing costs in Figures 3 and 7.
Acknowledgements
The authors are grateful to members from the Centre for Disease Modelling for feedback on the model. We thank Suzan Sardroodi for a thorough review of the manuscript.