Skip to main content

ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 11 August 2020
Sec. Dynamical Systems
This article is part of the Research Topic Mathematical modelling of the pandemic of 2019 novel coronavirus (COVID-19): Patterns, Dynamics, Prediction, and Control View all 34 articles

Mathematical Modeling of Business Reopening When Facing SARS-CoV-2 Pandemic: Protection, Cost, and Risk

\nHongyu Miao
Hongyu Miao1*Qianmiao GaoQianmiao Gao1Han FengHan Feng1Chengxue ZhongChengxue Zhong1Pengwei ZhuPengwei Zhu1Liang WuLiang Wu1Michael D. SwartzMichael D. Swartz1Xi LuoXi Luo1Stacia M. DeSantisStacia M. DeSantis1Dejian LaiDejian Lai1Cici BauerCici Bauer1Adriana PrezAdriana Pérez1Libin RongLibin Rong2David LairsonDavid Lairson3
  • 1Department of Biostatistics and Data Science School of Public Health, University of Texas Health Science Center at Houston, Houston, TX, United States
  • 2Department of Mathematics University of Florida Gainesville, FL, United States
  • 3Department of Management, Policy and Community Health School of Public Health, University of Texas Health Science Center at Houston, Houston, TX, United States

The sudden onset of the coronavirus (SARS-CoV-2) pandemic has resulted in tremendous loss of human life and economy in more than 210 countries and territories around the world. While self-protections such as wearing masks, sheltering in place, and quarantine policies and strategies are necessary for containing virus transmission, tens of millions of people in the U.S. have lost their jobs due to the shutdown of businesses. Therefore, how to reopen the economy safely while the virus is still circulating in population has become a problem of significant concern and importance to elected leaders and business executives. In this study, mathematical modeling is employed to quantify the profit generation and the infection risk simultaneously from a business entity's perspective. Specifically, an ordinary differential equation model was developed to characterize disease transmission and infection risk. An algebraic equation is proposed to determine the net profit that a business entity can generate after reopening and take into account the costs associated of several protection/quarantine guidelines. All model parameters were calibrated based on various data and information sources. Sensitivity analyses and case studies were performed to illustrate the use of the model in practice. The results show that with a combination of necessary infection protection measures implemented, a business entity may stand a good opportunity to generate a positive net profit while successfully controlling the within-business infection prevalence under that in the general population. The use of personal protective equipment (PPE) is also found of significant importance, especially at the early stage of business reopening.

Introduction

Severe Acute respiratory Syndrome Coronavirus 2 (SARS-COV-2) is an ongoing global health threat to every country in the world and has caused significant loss of almost every business entity [13]. In the United States alone, as of May 14, 2020, more than 1.3 million confirmed infections and more than 84,000 deaths had been reported since the COVID-19 outbreak, which is associated with a 4.8% drop in gross domestic product (GDP) in the first quarter of 2020. While the development of effective vaccination, treatment, and prevention strategies is underway, it is still unknown when such efforts will yield clinical and medical practices effective enough to allow a return to usual economic activity [4, 5]. Practicing a variety of stringent quarantine and shutdown policies is an effective way to protect our citizens' health and lives during a pandemic like SARS-CoV-2 [6]. However, there are costs that need to be considered. Besides the extraordinary physical and emotional stress and potentially significant medical expenses that COVID patients and their families must cope with, a large number of people have lost (or may lose) their jobs or businesses. This significant financial loss results in a pressure to reopen the economy prior to the availability of effective vaccination and treatment of COVID-19. Elected leaders and business executives as well as employees must address a critical question: how to reopen the economy in the midst of a SARS-CoV-2 pandemic safely?

In this study, we use mathematical modeling techniques to address particular challenges a business entity may face during reopening. In recognition that mathematical modeling results alone cannot stop pandemics, [7] the following behavioral and social guidelines are still strongly suggested for any business entities in planning to reopen:

1. Social distancing (including dinning in restaurants and manufacturing in factories), e.g., the number of workers (and clients or patrons) that a business environment can support while maintaining a 6ft distance between individuals;

2. Mask, glove, and goggle wearing while not alone;

3. Routine sanitization (e.g., entrance, exit, home, workplace, conference room, bathroom, public transportation, door knobs, shared electronic devices);

4. COVID-19 tests accessible to workers who have symptoms or recent exposure to confirmed infections;

5. Deployment of non-contact sensors (e.g., Kinsa smart thermometer) for real-time fever detection and location;

6. Case reporting and quarantine policy;

7. Determination of a maximum time duration of exposure to working environment;

8. Specific equipment (e.g., stronger ventilation system, UV purification system) for aerosol transmission prevention;

9. Employees in non-contact positions remain working from home;

The implementation of each guideline above can protect employees' health and lives but may add an economical cost. Indeed, the main purpose of reopening businesses is to prevent further deterioration of our economy by generating profits and provide incomes that many citizens desperately need. However, for employees in a contact-based position, the risk of COVID infection and transmission may remarkably increase if they return to work. When the infection rate of COVID-19 within a business entity becomes higher than the prevalence in the general population, it may trigger a domino effect and subsequent infections in a broader community originated from this business entity. Therefore, a delicate balance between profit/income generation and infection/transmission prevention must be the ideal. The focus of this modeling strategy is to provide a quantitative tool for investigating whether a business entity may reopen according to a simple definition of operational reopening feasibility: reopening is considered as feasible if a business entity can generate positive net profit after reopening while keeping the prevalence of COVID-19 infection within this entity less than or equal to the prevalence in the general population. The prevalence in the general population is assumed to remain under a certain threshold after reopening; otherwise, isolation and quarantine orders may be re-issued by government and reopening would halt, disregarding an individual business's performance.

Materials and Methods

Mathematical Model

To our knowledge, there exist only a few studies that have modeled the reopening of a country, a state, a city or a local community [811], and fewer studies have focused on the reopening of a business entity. This modeling work attempts to address the business-reopening problem by considering a transmission model together with a net profit equation.

Borrowing from the classical susceptible-infected-recovered (SIR) modeling framework as in several previous models [8, 12], the following definitions and notations are introduced. Let NT denote the total number of employees in a business entity, Nc denote the number of employees in a contact-based position, and NN denote the number of employees in a non-contact position. Employees in a non-contact position are expected to remain working from home (WFH) according to Guideline #9 and thus are excluded from the transmission model. These WFH employees still actively contribute to profit generation and are counted in the net profit equation. In addition, the total number of employees NT is assumed variable after the reopening due to various reasons, including infection-related death or other factors. Among the Nc employees in a contact-based position, let S denote the number of susceptible employees; C denote the number of asymptomatic, pre-symptomatic, and very-mildly-symptomatic (VMS) carriers (collectively known as silent spreaders); Q denote the number of infected employees on quarantine, under treatment (i.e., confirmed infection) or awaiting test result (i.e., unconfirmed infection); and D denote the number of deaths or resignations due to infection. Note that: (i) all infected employees confirmed by virus tests should be quarantined immediately, and any company fellow having a recent contact with the infected person is assumed to also be tested for the virus and initially be self-quarantined (see Guideline #4). The asymptomatic, pre-symptomatic, and very-mildly-symptomatic carriers are the three major categories of the so-called “silent spreaders.” While it remains unclear how quickly asymptomatic carriers can transmit the infection, some studies in China suggested that ~25% of those who tested positive without symptoms continued symptomless, and the remaining 75% turned out to be pre-symptomatic [13]. Other ongoing studies suggested that the proportion of asymptomatic infection could be as high as 50% [14]. People with very mild symptoms (e.g., occasional light cough or mild fever) may not recognize the infection but are also fully capable of spreading the disease like the presymptomatic carriers [15]. Finally, the WFH employees may still be infected and then practice at-home quarantine or receive treatment in hospital; for simplicity, instead of introducing another set of equations for characterizing the WFH transmission, the number of infected WFH employees NNQ is assumed to be collected by the business entity on a self-reported basis and directly accounted for in the net profit calculation equation.

After taking silent spreaders into consideration, [16] the proposed model structure is given below:

dSdt=-βO·κ·S-βA·αCN·S-βP·(1-α)CN·S    (1)
dCdt=βO·κ·S+βA·αCN·S+βP·(1-α)CN·S           -τ·(1-α)C+ω·Q    (2)
dQdt=τ·(1-α)C-(γ+ω)·Q-δ·Q    (3)
dRdt=γ·Q    (4)
dDdt=δ·Q    (5)

and the net profit is defined as the following:

Pntot=ξ1·ξ7·ρ·((NN-NNQ)+S+C+R)                    -i1or 7ξi·(S+C+R)-w·Q-w·NNQ    (6)

where = S+C+R denotes the number of employees in contact-based positions who are working on site. The definitions of all model variables and parameters are listed in Table 1, and a diagram (Figure 1) is also given to illustrate the transmission model mechanisms and assumptions. As suggested in Figure 1 and Equation (1), susceptible employees can be infected and become virus carriers by contacting people outside of the business at a rate βO, or interacting with asymptomatic carriers within the business at a rate βA and with pre-symptomatic and VMS carriers within the business at a rate βP. Also, α denotes the proportion of asymptomatic carriers among all tested-positive infections, and κ denote the probability of transmission via an average number of contacts per day of one person with other persons. In Equation (2), the infected susceptible become asymptomatic, pre-symptomatic or VMS carriers, and those pre-symptomatic and VMS carriers further progress to symptomatic infections at a rate τ. In Equation (3), symptomatic infections may recover at a rate γ and become immune to virus infection due to, e.g., memory immunity, [17] regress to symptomless carriers (e.g., some patients can test positive and shed viruses after symptoms end [18]) at a rate ω, or die at a rate δ. Also, Equations (4) and (5) characterize the dynamics of recovered and death, respectively. In the net profit Equation (6), ρ denotes the net profit per day generated by healthy WFH, susceptible, and silent carrier employees, ξi (i≠1, 7) is the average cost per person per day in a contact-based position associated with the i-th protection Guideline (i = 1, …, 9), and w is the average wage per day per person for employees on quarantine or under treatment. It is assumed that employees on quarantine or under treatment still receive their wages given that the typical quarantine or treatment length is 14–28 days even considering them as sick days.

TABLE 1
www.frontiersin.org

Table 1. Variable and parameter definitions, values, and sources.

FIGURE 1
www.frontiersin.org

Figure 1. Transmission model diagram. Each rectangle represents a model variable, a directed arrow shows the transition from one to another model variable, and the associated model parameters are labeled next to an arrow.

Data Source and Parameter Values

The datasets generated from recently published studies and surveys are heterogeneous in terms of population demographics, geospatial characteristics, medical resource availability, treatment regimens, prevention polices and implementation, among others. Extraordinary efforts are needed to integrate such heterogeneous datasets and perform a variety of statistical analyses, which is beyond the scope of this study. Instead, the summary statistics or previously-calibrated parameter values from previous studies and surveys are used as the primary information sources for our model parameter calibration. Note that studies reported different accuracy (in terms of decimals), it is thus difficult to standardize the number of decimals without losing accuracy so we keep the original numbers.

To calibrate the transmission model parameters in Equations (1)–(5), we started with the case death rate and the recovery time. In the United States, the case fatality rate is currently 5.65% while Omer et al. [29] previously estimated the infection-related death rate as low as 1.7%. In other countries, the reported case death rates may range from 0.56% (Iceland) to 13.53% (Italy); [29]1,2. to preserve parameter uncertainty, this wide range of infection fatality rate (0.56~13.53%) is adopted in our calculations. To calculate the mean recovery time, we consider the following observations:3 (i) 8 out of 10 infected persons with symptoms have only mild illness; (ii) the average recovery time for mild cases is about 2 weeks; (iii) for severe cases, recovery may take up to 6 weeks; (iv) the overall recovery rate is between 97 and 99.75%. The conservative mean recovery time for all cases is thus 1/[(0.814+0.242)×0.97]=16.65 days among people with symptomatic infections, which is longer than the reported median hospitalization period of 12 days [33] among survivors. It was also reported that the recovery time for mild cases can be as short as 7 days, therefore, this study assumes that recovery time ranges from 7 to 42 days4.

In Equation (1), for the proportion of asymptomatic carriers (denoted by α), multiple studies have been conducted to estimate this parameter among different populations with different methods [14, 1922]. The range of the estimate of α^ is (13.8%, 75%), and about half of such studies reported a result around 50%. Sun et al. (26) analyzed 391 cases in Zhejiang Province, China from Jan. 20th, 2020 to Feb. 10th, 2020, and found 54 (13.8%) cases were asymptomatic. Nishiura et al. [22] estimated α in the evacuated Japanese citizens to be 30.8% (95% CI: 7.7–53.8%). Mizumoto et al. [14] obtained an estimate 17.9% (95% CI: 15.5–20.2%) using the data from the Diamond Princess cruise. Kimball et al. [21] analyzed the results of symptom assessment and SARS-CoV-2 testing in King County, Washington, and found that 56.5% of the tested positive was asymptomatic. Day [20] suggested that the proportion of asymptomatic infection was between 50 and 75% in northern Italy, and others5 suggested that α could be between 25 and 50% on Apr. 5th, 2020 [34].

Furthermore, κ is the probability of one employee having contacts with infected carriers outside a business entity, and βO denotes the associated transmission rate. A limited number of studies were found helpful for estimating these two parameters, we thus calibrated such parameter values primarily using the simulation results from Tang et al. [23]. In Tang's study, the daily average contact number was estimated to be 14.78 (SE 0.90) contacts per day per person, and the probability of successful transmission per contact was estimated to be 2.10 × 10−8 (SE 1.19 × 10−9). According to the definitions of κ and βO, we have βOκ=14.78×2.10×10-8= 3.10×10-7 day-1, and the range of βO κ, (2.42×10-7, 3.88×10-7), can be calculated using the standard errors (i.e., 2.42 × 10−7 = (14.78−2 × 0.90) × (2.10 × 10−8−2 × 1.19 × 10−9), and 3.88 × 10−7 = (14.78+2 × 0.90) × (2.10 × 10−8+2 × 1.19 × 10−9).

The transmission rate (βP) can be calculated from R0recovery time, where R0 denotes the reproduction number (i.e., the average number of new infections by one carrier) that has been frequently investigated in many SARS-CoV-2 studies [13, 23, 3543]. For instance, according to the study of Liu et al. [26], the estimated R0 of SARS-CoV-2 in China at the early stage of the pandemic ranges from 1.4 to 6.49 with a mean 3.28 and a median 2.79. Note that the population in Liu et al. [26] was mostly not aware of the pandemic and used little protections such as mask wearing, hand washing, and social distancing. Thus, for the unprotected general population, the estimated transmission rate is 3.28/16.65= 0.197 day−1. The corresponding range is between 0.0333 (=1.4/42) and 0.927 (=6.49/7) per day. Note that in the early stage of a pandemic, it can be assumed that most of the cases are pre-symptomatic or mildly-symptomatic such that the parameter value calculated above may be close to the true value of βP. Another study by Li et al. (43) suggested that the transmission rate βP could be 1.12 (95% CI: 1.07–1.18), 0.52 (95% CI: 0.42–0.72), or 0.35 (95% CI: 0.28–0.45) per day by fitting data of different pandemic stages in China. For these reasons, we assumed a range for βP of (0.0333, 1.18) per day.

For the transmission rate of asymptomatic infection (βA), very little useful information was found in the scientific literature. Here we borrowed the idea of Li ‘s work, [25] in which the transmission rate of asymptomatic infection was calculated by multiplying a reduction factor μ with the transmission rate of symptomatic infection. The reduction factor μ was estimated as 0.55 (95% CI: 0.46–0.62), 0.50 (95% CI: 0.37–0.69) and 0.43 (95% CI: 0.31, 0.61), corresponding to different stages of SARS-CoV-2 outbreak in China. Using the same reduction factor μ in our model, βA can be estimated by μβp, leading to a mean of 0.197 × 0.5 = 0.099 per day and a range from 0.0103 (= 0.333 × 0.31) to 0.814 (= 1.18 × 0.69) per day.

In Equation (2), τ is the rate of silent carriers progressing to symptomatic infections. The mean incubation period was previously reported to be 5.2 days with a 95% confidence interval (4.1, 7.0) days, [13, 24] or 5.1 days with a 95% confidence interval (2.2, 11.5) days [44]. Based on these results, we estimated τ as the reciprocal of the incubation period with a mean 1/5.2=0.192 per day and a 95% confidence interval (0.143, 0.244) per day.

In Equation (3), for the clearance rate of symptomatic infection γ and the death rate due to COVID infection δ, multiple studies developed various methods to estimate them [6, 23, 45]. In Piguillem's method, (6) the calculations were mainly based on the case mortality rate and the recovery time. After substituting 16.65 days for the mean recovery time and 5.65% for the mean death rate aforementioned, we obtain the estimate of γ as (1−0.0565)2/16.65=0.0535 per day, and the death rate δ as (1-0.0565)×0.056516.65=0.0032 per day. The range of γ is found to be (0.0178, 0.141) per day, and the range of δ is (1.3 × 10−4, 1.67 × 10−2) per day.

In Equations (2) and (3), we also introduced ω, the rate of regression to carriage, considering the fact that 3.23% of the patients recovered from SARS-CoV-2 infection were tested positive after hospital discharge [27, 28]. However, the potential infectivity of these carriers remains unclear. Here ω is estimated to be 3.23% ×γ^=0.0323 × 0.0535= 0.00173 per day, and its range is between 0.000575 (=0.0178 × 3.23%) and 0.00455 (=0.141 × 3.23%) per day.

Now consider the additional parameters in the net-profit Equation (6). The average net profit per capita ρ in the U.S. is found to be $400.73 per person per day, which is calculated by dividing Jan 2020 U.S. corporate profit $1908.02 Billion US dollars6. by 158,714,000 (the number of employed persons in the United States)7 and then by 30 days. The average weekly wage of one employee in the U.S. is $1,093 per person in the third quarter of Year 2019,8 so the average daily wage of employee w is $218.60 per day (dividing $1,109 by the number of weekdays 5). In the previous study of Thumstrom and Newbold (2020), [12] GDP loss was considered as one cost of social distancing; i.e., an immediate GDP decline associated with practicing social distancing alone (i.e., house-hold quarantine) in the United States was predicted to be 13.7−6.49 = 7.21 trillion US dollars. The projected GDP of Fiscal Year 2020 is 22.11 trillion US dollars according to the United States Congressional Budget Office (CBO)9. Therefore, for Guideline #1, the cost of social distancing is the loss of productivity by 7.21 ÷ 22.11 = 32.6%, and thus ξ1 = 1−32.6% = 67.4%. For Guideline #2, the cost of personal protective equipment (PPE), [30] including surgical mask, gloves, goggle wearing, hand sanitization, and soap, is calculated as ξ2=$3.60 per person per day [31] under the assumption that each person will consume two surgical masks per day, two pairs of gloves per day, and one goggles. According to the Personal Protective Equipment (PPE) guideline from Perdue University,10 goggles can last for years if kept clean by using mild soap and water, and if stored in a protected, dry, and temperate storage location. So the use life of goggles is much longer than our typical setting for simulation time window length (i.e., 100 days) and the cost of one goggles per day is $5 ÷ 100 = $0.05 per day for simplicity. The detailed costs of each PPE item as well as hand sanitization and soap can be found from online resources11,12. For Guideline #3, the cost of routine sanitization (e.g., cleaning the workplace, bathroom, and shared electronic devices) is calculated as ξ3= $10.45 per environmental service staff per day [32]. Specifically, in the work of White (2019), 11 hospitals with a total of 1,700 environmental service staffs, they found to spend $11,308 per week to maintain the hygiene by consuming disinfection products. We thus calculate ξ3 using the following equation:

ξ3=cost of disinfection products(Total number of enviromental services staffsNumber of hospitals)÷7 days.

For Guideline #4, the current policy dictates that the test of COVID is free to the business entity (ξ4 = $0 per person per day). For Guideline #5, as the current market prices of a non-contact sensor may range from $50 to $100, and one non-contact sensor is required for each worksite. With the assumption of having 100 employees per worksite, the cost of deploying non-contact sensor ξ5 is $0.005 to $0.01 per person per day for a 100-day time horizon. We assume that Guideline #6 does not cost any money (ξ6 = $0 per person per day). For Guideline #7, the cost is proportional to the current maximum working hours (i.e., ξ7=current maximum working hoursregular working hours). In our model, the current maximum working hours is assumed to be 70% of the regular working hours. For Guideline #8, referring to the work of Chen (2013), [31] the cost of deploying specific equipment for reducing aerosol transmission such as ultraviolet germicidal irradiation (UVGI) and high-efficiency particular air (HEPA) filtration are $182.37 and $136.78 per person per year, respectively, ξ8 is thus equal to $0.874 per person per day via dividing the total cost of the aforementioned equipment by 365 days. Finally, we assume that implementing Guideline #9 does not incur any cost (ξ9 = $0 per person per day).

All the parameter definitions, values and ranges are summarized in Table 1. However, it should be stressed that the parameter values calibrated above are for heterogeneous populations. More importantly, at the early stage of the pandemic, the estimates of certain transmission model parameters (e.g., the reproduction number and transmission rates) are expected to be larger due to the absence of protection policies and self-protection awareness; and the transmission rates are expected to have a notable drop after the implementation of various protection and quarantine strategies (PQS). Such a hypothesis is supported by several recent studies, which showed that the overall transmission rates may decrease by 2.1–3.2 folds after implementing PQS. Also, the study of Seto et al. [46] quantified the odds ratios of SARS infection as 13 (95% CI: 3, 60), 2 (95% CI: 0.6, 7), or 5 (95% CI: 1, 19), corresponding to the use of masks, gloves or hand-washing, respectively. Thus, the values of certain model parameters (e.g., transmission rates) will be different from (e.g., smaller than) the values calibrated in this section after implementing Guidelines #1-9, which will be elaborated in the result Section.

Implementation and Computing Configuration

All the computing codes were implemented in MATLAB® (MathWorks, Natick, MA), and the ordinary differential equation (ODE) solver ode15s was employed for solving the transmission model numerically. The relative error tolerance was set at 10−7, the absolute error tolerance at 10−7, and the maximum step size at 10−2.

Results

Sensitivity Analysis

We evaluated the local sensitivity of the transmission model in Equations (1–5). Specifically, we evaluated the changes in the model outcome variables (i.e. S, C, Q, R, D) corresponding to a 1% increase in one parameter value, with the other seven parameter values fixed at their default values as in Table 1. Initial values for the sensitivity analysis were set as 299 susceptible and 1 silent carrier which approximated the proportion of infections in the U.S. general population as reported by May 5, 2020. The initial numbers of quarantined, recovered, and death were all set as 0 for simplicity. The results were visualized in Figure 2, showing that the transmission model was most sensitive to the transmission rate of pre-symptomatic and VMS infections (βP) and the rate of progression to symptomatic infection (τ), and least sensitive to the product βO·κ. For the most sensitive parameters, a 1% increase in parameters βP or τ resulted in a change <1.5% in the outcomes throughout a period of 200 days. At the end of the 200 days, all changes in the outcomes reached a plateau or tended to diminish. The results suggested that the transmission model was not locally sensitive to parameter value changes and could make robust predictions. Additional sensitivity analyses were then performed to understand the outcome changes over the entire range of parameter values. As shown in Figure 2, S was one of the most sensitive outcomes to parameter value changes; therefore, we plotted S against time for each of the eight parameter value ranges in Figure 3. For each parameter range, the minimum, the 1st quantile, the mean, the 3rd quantile, and the maximum parameter values were used to generate the trajectories of S. The results suggested that the dynamic model did not crash into chaos and the outcome S varied smoothly and predictably within the entire parameter range.

FIGURE 2
www.frontiersin.org

Figure 2. Local sensitivity analysis of the transmission model, showing the changes in the model outcome variables corresponding to 1% increase in one parameter value (with the other seven parameter values fixed at their default values as in Table 1).

FIGURE 3
www.frontiersin.org

Figure 3. The trajectories of S over the entire range of each parameter. For each parameter, the trajectories of S are generated corresponding to the minimum, the 1st quantile, the mean, the 3rd quantile, and the maximum parameter values.

To evaluate the global sensitivity of this model, the Partial Spearman Rank Correlation Coefficient (PRCC) method [47] was employed. The model parameters were randomly sampled 100 times using the Latin hypercube sampling technique over uniform distribution. The range of the uniform distribution for the global sensitivity for each model parameter was reported in Table 1. Figure 4 showed the plot of the PRCC values (next to the y axis) between model parameters and ODE outputs against time. As suggested in many previous studies, an absolute PRCC value >0.4 was deemed as practically significant [48]. For the susceptible population (S), the corresponding subfigure showed that the transmission rates βA and βP were strongly and negatively (< −0.7) correlated with S. This result was expected as the higher the transmission rates, the smaller the number of the susceptible persons. Please note that the percentage of asymptomatic carriers (α) was positively correlated with S initially (> 0.4) and then became negatively correlated with S at the end (< −0.6). The initial positive correlation between S and α should not be interpreted as that S will increase when α increases, but that S will decrease less when α increases. This will happen when βA·α is less than βP·(1−α). Around the end of the time window of the simulation, a larger α corresponded to a smaller N = S+C+R due to the increase in death, and thus S will decrease more given a smaller N in the denominator at the righthand side of Equation (1). From all the subfigures, the transmission rates βA and βP were always found strongly correlated with the outcomes while it was not the case for βO·κ, the value of which was too small to have a substantial impact on the outcomes. In addition, parameters α, γ and τ were also found strongly correlated with the outcomes, and interventions like more effective drug treatment or vaccination would affect these parameter values.

FIGURE 4
www.frontiersin.org

Figure 4. Global sensitivity analysis (PRCC) of the transmission model.

Case Study

Business executives are strongly recommended to follow the guidelines listed in the Background Section to prevent the reopening from causing any exacerbation of the ongoing pandemic. However, practical difficulties may arise due to, e.g., insufficiency of budget or medical recourses such that only some of the guidelines will be implemented by employers. Through this case study, we illustrated the use of the transmission model and the cost equation in different scenarios to evaluate the feasibility of reopening. Note that our simulation results were obtained under many assumptions and subject to both model structure and parameter value uncertainty; therefore, decisions of business executives should not be made solely based on the results presented in this section.

In this case study, we focused on four scenarios: (i) none of the nine Guidelines was constantly implemented (baseline scenario); (ii) all Guidelines 1–9 were constantly implemented; (iii) PPE and all other indispensable guidelines (1, 2, 3, 4, 6, 9) were constantly implemented; (IV) Indispensable guidelines (1, 3, 4, 6, 9) were constantly implemented but without the use of PPE. Scenarios I was the baseline scenario, corresponding to a complete devaluing of infection risk. Across scenarios II-IV, Guidelines 1, 3, 4, 6, 9 were assumed to be always implemented given the necessity and indispensability of these five Guidelines to business reopening. Scenario III was designed to be less restrictive than Scenario II, considering that some business entities might not have the budget to deploy non-contact sensors for real-time fever detection, reduce the number of working hours, or acquire specific equipment for aerosol transmission prevention. In Scenario IV, the use of PPE was intentionally dropped to account for the possible shortage of such materials on the market and to understand the importance of the use of PPE. The implementation of behavioral and social practice guidelines in Scenarios II-IV led to changes in the values of these three, βO·κ, βA, and βP, parameters in our transmission model. The other model parameters were assumed not directly affected by Guidelines 1-9. For instance, the recovery rate γ primarily depends on subject-specific immunity, the availability, affordability of effective medical intervention and health care, instead of behavioral and social practice patterns. For these reasons, we adjusted the values of βO·κ, βA, and βP for scenarios II–IV, respectively. The guidelines followed in each scenario and the corresponding adjusted parameter values were listed in Table 2. For simplicity, we assumed that the effects of different guidelines were independent and remained constant throughout the entire simulation time window. We also assumed that after reopening, every employee rigorously followed the guidelines and immediately reported their symptoms or infections once identified.

TABLE 2
www.frontiersin.org

Table 2. Parameters value adjustment in different scenarios.

According to the study of Koo et al. [49], the practice of social distancing together with disease testing, reporting and quarantine policy (Guidelines 1, 4, and 6) could reduce the number of infections by 78.20% (IQR: 59.0–94.4%) compared with the baseline scenario when Ro was 2.5. That is, the transmission rates might drop to 1 – 78.2% = 21.80% of their baseline values after implementing Guidelines 1, 4, and 6. To quantify the effects of wearing mask, glove, goggle and hand washing (Guideline 2) on infection transmission, we adopted the results from Seto et al. [46] and Yin et al. [50] for severe acute respiratory syndrome (SARS). In their results, the odds ratio for mask wearing was 13 (95% CI: 3–60, attack rate 1183=13.25%), for glove wearing was 2 (95% CI: 0.6–7, attack rate 9133=6.77%), for hand-washing was 5 (95% CI: 1–19, attack rate 317=17.65%), and for goggles wearing was 10.2=5 (95% CI: 2.44−10, attack rate 61.50%). Correspondingly, the transmission rates might drop to 8.83% (95% CI: 1.92–36.56%) for mask wearing, 48.85% (95% CI: 15.17–100%) for glove wearing, 39.37% (95% CI: 22.40–64.33%) for goggle wearing, and 23.91% (95% CI: 6.32–100%) for hand washing. Thus, due to the implementation of Guideline 2, the transmission rates might drop to 8.83% × 48.85% × 39.37% × 23.91% = 0.41% of their baseline values. Zhang et al. [51] investigated influenza A transmission in a student office setting, which showed that the infection risk could reduce by 2.14% with implementation of routine sanitization. Routine sanitation is expected to reduce the transmission rates to 1−2.14% = 97.86% of their baseline values after implementing Guideline 3. The implementation of Guideline 5 would improve the implementation of Guidelines 4 and 6 so its effect on infection transmission was not explicitly quantified in this study. For Guideline 7, we assumed a linear relationship between the number of working hours and the infection risk. The average working hours per day after reopening was assumed to be 70% of the regular working hours. Correspondingly, the transmission rates might drop to 70% (range 50–100%) with Guideline 7. For Guideline 8, according to the report of Mendell et al. [52], a 10–14% reduction in communicable respiratory infections might result from improved work environments. Taking 12% as the median, we assumed that the transmission rates might drop to 100% − 12% = 88% (range 86–90%) with implementing Guideline 8. Furthermore, the transmission rates among employees in contact-based positions were not affected by Guideline 9. Finally, in scenarios II, III, and IV, the transmission rates dropped to21.80% × 0.41% × 97.86% × 70% × 88% = 0.05%, 21.80% × 0.41% × 97.86% = 0.09% and 21.80% × 97.86% = 21.33% of their baseline values, respectively. See Table 2 for the adjusted transmission rates.

All the simulation results in this section were generated using the same set of initial values for simplicity. Taking the DELLTM center at Austin, Texas as an example, the total number of active workers after reopening could be around 14,000. The ratio between the on-site workers and the WFH workers was assumed as 2:1 (i.e., 9,333 on-site workers and 4,667 WFH workers). On the first day of reopening (day 0), the proportion of silent carriers among all the 14,000 employees was approximated by the proportion of infections in the general population of U.S. estimated using the number of reported cases as of May 5, 2020. The proportion of workers that had recovered from the infection and thus acquired immunity (referred as “recovered” in our model) were also estimated from the reported number of recoveries of the general U.S. population. The numbers of quarantined employees and death on day 0 were assumed to be 0. We conducted the simulations for 200 days for Scenario I and 100 days for Scenarios II–IV to verify the short-term feasibility of reopening. The outcome trajectories in the four scenarios were shown in Figures 5, 6 showed the reopening feasibility by plotting the net profit and the prevalence of infections in the workplace as well as in the general population (the left column), and the corresponding values of cost, profit, and net profit (the right column).

FIGURE 5
www.frontiersin.org

Figure 5. Outcome trajectories for the four scenarios. For the same scenario, the five outcome variables of different orders of magnitude are visualized by three subfigures in one row. The scale of y-axis is adjusted to enable a clear view of each single trajectory over time. (A,D,G,J) Susceptible (S); (B,E,H,K) Recovered (R); (C,F,I,L) Carrier (C), Quarantined (Q), and Death (D).

FIGURE 6
www.frontiersin.org

Figure 6. Cost, profit, net profit, and infection prevalence for the four scenarios (each row for one scenario and Scenarios I–IV in the order of the top row to the bottom). For each scenario, the left subplot shows the reopening feasibility based on both net profit and infection prevalence in the workplace, with the prevalence in the population presented for comparison; and the right subplot shows the cost, profit and net profit. (A,C,E,G) Reopening feasibility, (B,D,F,H) cost and profit.

As shown in Figures 5A–C, without implementing any of the guidelines, reopening merely led to a large number of infections and deaths (230 deaths by the last day, 649 quarantined and 422 silent carriers at their peaks). The prevalence of infections within the business entity exceeded that in the general population throughout the 200-days window, and peaked at 114 infections per 1,000 people, which was 30 folds higher than that of the general population. Even by the end of the 200 days, the within-entity prevalence was 8 folds higher (32.9 cases per 1,000 persons) than the population average. Also, even after we extended the simulation time window to 200 days, the transmission model still did not reach its steady states, suggesting a less predictable risk of disease transmission. In short, despite that the net profit remained positive (Figure 6B), since the prevalence in workplace was constantly (much) higher than that in the general population, reopening turned out to be infeasible in Scenario I.

If a business entity strictly followed all guidelines 1–9 as in Scenario II, the number of infections and deaths reduced remarkably, as shown in Figures 5D–F (2 deaths among 9,333 onsite workers, at most 16 quarantined and 35 silent carriers). Figure 5 also showed that the outcomes of the transmission model reached a plateau toward the end of the 100-days window. In other words, the infection could be contained in this scenario, with the number of carriers and the number of the quarantined reaching and staying at 0 by day 48 and day 92, respectively. According to Figure 6C, under the assumption that the prevalence of infections within the business entity was the same as in the general population at the beginning of reopening, the within-business prevalence dropped under the general population prevalence immediately (on day 2), and continuously decreased down to 0.035 infections per 1,000 people by the last day of simulation. This within-business prevalence was 104 folds lower than the general population prevalence reported in other studies. The business also attained higher and more stable net profit than in Scenario I (Figure 6D).

In Scenario III, guidelines 5, 7, 8 were skipped and the transmission rate values changed accordingly (see Table 2); however, the simulation results were surprisingly similar to those in Scenario II during the 100-day time window (Figures 5G–I, 6E,F). Specifically, the numbers of deaths, the quarantined and silent carriers were nearly the same as those in Scenario II; also, it took approximately the same amount of time for the numbers of carriers and the quarantined to drop to 0 (Figures 5G–I). The nearly same predicted population trajectories resulted in a nearly same prediction on net profit (Figures 6E,F). To confirm such results, additional local sensitivity analysis was conducted at the parameter values in Scenarios II and III, suggesting that the simulation outputs were not sensitive to parameter value changes (e.g., a change of <0.1 in all the five output variables corresponding to 1% parameter value change).

After a business entity further dropped guideline 2 (Mask, glove, and goggle wearing while not alone) in Scenario IV, our model predicted that the spread of infections did not become uncontrollable within 100 days but associated with a higher cost (Figures 5J–L, 6G,H). While the initial within-business prevalence was the same as the general population prevalence, the workplace quickly became a “hot spot” of infection spreading, and a prevalence much higher than the population average was reached (Figure 6G). The within-business prevalence continued to stay above the general population average for 11 days, and peaked at a level of 4 infections per 1,000 people on day 5. The results suggested that reopening should stop to prevent this entity from developing into a source of infection and posing significant risk on its workers as well as their close social contacts. Note that the within-business prevalence dropped to 0.09 infections per 1,000 people by the end of the 100-day time window (Figure 6G), and the numbers of deaths and the quarantined were controlled under 3 and 20 among 9,333 onsite workers, respectively. Business executives should not rely on such optimistic predictions and underestimate the infection risk for two reasons. First, constant parameter values were used in our simulations, which were not capable of capturing every possible time-varying characteristic of disease transmission over time (i.e., parameter values could be time-varying instead of being constant). Second, our simulation was performed under the assumption that all workers strictly followed the selected guidelines. In reality, it is unlikely that every single employee would faithfully stick to such rules and guidelines daily and constantly over time. Moreover, Figure 6G suggested that the use of PPE was very important at the early stage of reopening. In Scenarios II and III, with the use of PPE, we did not observe the rapid increase of prevalence since the beginning of reopening; however, in Scenario IV without the protection of PPE, the within-business prevalence started to increase and exceed the population average on day one. This result was consistent with the recent study by Kai et al. [53], which demonstrated the significant effect of universal use of facial masks (e.g., at least 80% population wear masks) on impeding the spread of infections.

Discussion

SARS-CoV-2 struck the whole world since 2019 and caused significant loss of human lives and economy. While various lockdown, quarantine, and isolation rules and polices worked effectively to control COVID infection transmission, numerous businesses were closed and tens of millions of people lost their jobs. As of May 20, 2020, all the states in the U.S. strategically moved toward gradual reopening to save the economy. However, given that neither effective drug treatments nor vaccines were available, the risk of infection spreading within a business entity became a key issue that any business executive had to consider. In this study, an ODE-based transmission model was developed together with a net profit equation to quantitatively evaluate the trade-off among profit, cost, and infection risk within a business entity after reopening. Model parameter values were calibrated from heterogeneous sources to enable computer simulations. Both local and global sensitivity analyses were conducted to understand our model behavior and result robustness. Finally, a case study assessed scenarios, in which different combinations of behavioral and social practice guidelines were implemented. The simulation results suggested that infection transmission was controllable within a business entity and a positive net profit could be generated after reopening only if a combination of selected guidelines were implemented. Also, our results suggested that the use of PPE could be significantly important at the early stage of reopening to prevent infection spreading.

Modeling business reopening is a novel and important problem during the COVID-19 pandemic or a similar situation. This study adopted a business entity's perspective instead of a social perspective (e.g., one can further consider social cost of death due to infection and model the interactions out of a workplace). However, a study from a business entity's point of view has its own added value. Note that reopening decision is not just the call of local governments but also business executives. For instance, even if a state or local government issues the reopening orders, a business entity may decide by itself whether it should reopen immediately or postpone the reopening after taking into consideration of both business-specific parameters and benefit-cost-risk trade-off. With that said, future researches from different perspectives are more than welcome to shed novel insights into the business reopening problem.

There are several issues worth of further discussions. First, while the proposed model did not explicitly introduce a term for the incompliance of employees on the use of PPE, we have investigated both Scenario III and Scenario IV in the Case Study Section. The only difference between the two scenarios is whether PPE is used; therefore, the results corresponding to partial or full incompliance by employees will fall between the results of Scenario III and those of Scenario IV. Second, our model implicitly accounted for the fact that a virus testing process might not be instant and completely error free. Recall that the outcome variable for “carriers” (C) in our model is the total number of asymptomatic, pre-symptomatic, and very-mildly-symptomatic (VMS) carriers (collectively known as silent spreaders), as defined in the Mathematical Model Section. It is reasonable to expect that any person with notable symptoms should practice self-quarantine or receive medical treatment even when waiting for virus test results or tested false negative; therefore, cases that are not immediately identified or quarantined are accounted for in C. Third, βo × κ was found not to play a very significant role in this study. Note that after reopening, employees will spend the majority of their ‘gathering' time in the workplace so the workplace transmission would and should play a major role. Consequently, the more time employees spend on work-related activities, the less number of interactions they could have with others outside the workplace, suggesting a small value for βo × κ. Finally, while mathematically it could be interesting to examine all the 29 = 512 possible combinations of the nine guidelines to identify the optimal combination(s), it may not serve the practical purpose of this study. For instance, we can hardly imagine a worker who would strictly stick to the use of PPE but not practice any handwash; or a workplace gives up on the cost-effective measures such as social distancing and routine sanitization, while spending lots of budget on equipment for aerosol transmission prevention. More importantly, the infection protection measures considered in Scenario III are those strongly recommended by CDC due to their efficacy and cost-effectiveness. That is also why the optional protection measures were only considered in Scenario II to quantify the additional benefits such measures may provide.

We also recognize a number of limitations of this study. First, the mathematical model and the net profit equation were developed under multiple assumptions. While many of those assumptions are commonly adopted by related professional communities (e.g., epidemiology or mathematical modeling), some assumptions were introduced due to the lack of accurate and/or complete information (e.g., asymptomatic infection transmission rate). With additional efforts invested in future SARS-CoV-2 related research, we expect that informative and high-quality data will become available such that less assumptions are needed and more accurate results can be generated by our approach. Second, while we have compiled information from many different sources to calibrate the possible ranges of model parameter values, the parameter uncertainty may not be completely characterized by such parameter ranges due to, e.g., the heterogeneity in population demographics, health conditions, behavioral patterns, and social networks. Third, although the costs associated with guidelines #1-9 were included in the net profit calculation, our cost estimation is subject to market fluctuations and many other (unpredictable) factors. It is thus suggested for users of our model to fine tune parameter values upon available business-specific data or information. Finally, in this study, herd immunity and vaccination strategies were not considered because they were not available as of the current moment. Considering the active research on SARS-CoV-2 vaccine and drug development, the availability of effective vaccines and medications can be expected and it should be taken into consideration in our modeling work at some point.

In summary, this modeling work provides a quantitative tool for decision-makers to explore and evaluate the business reopening option in the midst of COVID pandemics, and it can be extended to similar scenarios (e.g., outbreak of unknown or new virus) by re-calibrating related parameter values. We expect further research efforts in this direction to better prepare for possible strikes of infectious diseases in the future.

Data Availability Statement

All datasets generated for this study are included in the article/supplementary material.

Author Contributions

HM designed and supervised this study, implemented computing code, and wrote the manuscript. QG, HF, CZ, PZ, and LW calibrated parameter values, implemented computing code, and contributed to manuscript writing. MS, XL, SD, DLai, CB, AP, LR, and DLair co-developed the modeling idea, provided parameter values, and contributed to manuscript editing.

Funding

This work was partially supported by NSF grant DMS-1620957 (HM), Weatherhead Foundation (DLai), and NSF grant DMS-1758290 (LR).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Also, thank Dr. Eric Boerwinkle for motivating this study and sharing constructive comments and thank Drs. Wenyaw Chan and Hulin Wu for insightful discussions and feedbacks.

Footnotes

References

1. Atkeson A. What Will Be the Economic Impact of COVID-19 in the US? Rough Estimates of Disease Scenarios. National Bureau of Economic Research Working Paper Series No. 26867 (2020). doi: 10.3386/w26867

CrossRef Full Text | Google Scholar

2. Anderson RM, Heesterbeek H, Klinkenberg D, Hollingsworth TD. How will country-based mitigation measures influence the course of the COVID-19 epidemic? Lancet. (2020) 395:931–4. doi: 10.1016/S0140-6736(20)30567-5

PubMed Abstract | CrossRef Full Text | Google Scholar

3. McKibbin WJ, Fernando R. The Global Macroeconomic Impacts of COVID-19: Seven Scenarios. CAMA Working Paper No 19/2020 (2020). Available online at SSRN: https://ssrn.com/abstract=3547729.

Google Scholar

4. Wang Y, Zhang D, Du G, Du R, Zhao J, Jin Y, et al. Remdesivir in adults with severe COVID-19: a randomised, double-blind, placebo-controlled, multicentre trial. Lancet. (2020) 365:1569–78. doi: 10.1016/S0140-6736(20)31022-9

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Williamson BN, Feldmann F, Schwarz B, Meade-White K, Porter DP, Schulz J, et al. Clinical benefit of remdesivir in rhesus macaques infected with SARS-CoV-2. bioRxiv. (2020). doi: 10.1101/2020.04.15.043166. [Epub ahead of print].

CrossRef Full Text | Google Scholar

6. Piguillem F, SHI L. Optimal Covid-19 quarantine and testing policies. CEPR Discussion Paper No DP14613. (2020). Available online at SSRN: https://ssrn.com/abstract=3594243.

Google Scholar

7. Huremović D. Social Distancing, Quarantine, and Isolation. In: Huremović D, editor. Psychiatry of Pandemics: A Mental Health Response to Infection Outbreak. Cham: Springer International Publishing (2019). p. 85–94. doi: 10.1007/978-3-030-15346-5_8

CrossRef Full Text | Google Scholar

8. Song B, Hei X. Models and strategies on reopening lockdown societies due to COVID-19. OSF Preprints. (2020). doi: 10.31219/osf.io/umtvh

CrossRef Full Text | Google Scholar

9. Bhatia R, Klausner J. A Step-by-Step Plan to Reopen California. (2020). Available online at: https://thebolditaliccom/covid-19-next-steps-for-california-c01632c8e8b6 (accessed January 5, 2020).

Google Scholar

10. Aboelkassem Y. COVID-19 pandemic: a Hill type mathematical model predicts the US death Number and the reopening date. medRxiv. (2020). doi: 10.1101/2020.04.12.20062893

CrossRef Full Text | Google Scholar

11. Wang X, Tang S, Chen Y, Feng X, Xiao Y, Zongben X. When will be the resumption of work in Wuhan and its surrounding areas during COVID-19 epidemic? A data-driven network modeling. Sci Sinica Math. (2020) 50:969–78. doi: 10.1360/SSM-2020-0037

CrossRef Full Text | Google Scholar

12. Thunstrom L, Newbold S, Finnoff D, Ashworth M, Shogren JF. The benefits and costs of using social distancing to flatten the curve for COVID-19. Forthcoming J Benefit Cost Anal. (2020) 11:179–95. doi: 10.1017/bca.2020.12

CrossRef Full Text | Google Scholar

13. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y, et al. Early transmission dynamics in wuhan, china, of novel coronavirus–infected pneumonia. N Engl J Med. (2020) 382:1199–207. doi: 10.1056/NEJMoa2001316

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Eurosurveillance. (2020) 25:2000180. doi: 10.2807/1560-7917.ES.2020.25.10.2000180

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Bai Y, Yao L, Wei T, Tian F, Jin D-Y, Chen L, et al. Presumed asymptomatic carrier transmission of COVID-19. JAMA. (2020) 323:1406–7. doi: 10.1001/jama.2020.2565

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Chisholm RH, Campbell PT, Wu Y, Tong SYC, McVernon J, Geard N. Implications of asymptomatic carriers for infectious disease transmission and control. Royal Soc Open Sci. (2018) 5:172341. doi: 10.1098/rsos.172341

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Dijkstra J, Hashimoto K. Expected immune recognition of COVID-19 virus by memory from earlier infections with common coronaviruses in a large part of the world population [version 1; peer review: awaiting peer review]. F1000Research. (2020) 9:285. doi: 10.12688/f1000research.23458.2

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wölfel R, Corman VM, Guggemos W, Seilmaier M, Zange S, Müller MA, et al. Virological assessment of hospitalized patients with COVID-2019. Nature. (2020) 581:465–9. doi: 10.1038/s41586-020-2196-x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sun WW, Ling F, Pan JR, Cai J, Miao ZP, Liu SL, et al. Epidemiological characteristics of 2019 novel coronavirus family clustering in Zhejiang Province. Zhonghua Yu Fang Yi Xue Za Zhi. (2020) 54:E027. doi: 10.3760/cma.j.cn112150-20200227-00199

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Day M. Covid-19: identifying and isolating asymptomatic people helped eliminate virus in Italian village. BMJ. (2020) 368:m1165. doi: 10.1136/bmj.m1165

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Kimball A, Hatfield KM, Arons M. Asymptomatic and presymptomatic SARS-CoV-2 infections in residents of a long-term care skilled nursing facility — King County, Washington. MMWR Morb Mortal Wkly Rep. (2020) 69:377–81. doi: 10.15585/mmwr.mm6913e1

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Nishiura H, Kobayashi T, Miyama T, Suzuki A, Jung S, Hayashi K, et al. Estimation of the asymptomatic ratio of novel coronavirus infections (COVID-19). Int J Infect Dis. (2020) 94:154–5. doi: 10.1016/j.ijid.2020.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Tang B, Wang X, Li Q, Bragazzi NL, Tang S, Xiao Y, et al. Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions. J Clin Med. (2020) 9:462. doi: 10.3390/jcm9020462

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Linton NM, Kobayashi T, Yang Y, Hayashi K, Akhmetzhanov AR, Jung S, et al. Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data. J Clin Med. (2020) 9:538. doi: 10.3390/jcm9020538

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li R, Pei S, Chen B, Song Y, Zhang T, Yang W, et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science. (2020) 368:489–93. doi: 10.1126/science.abb3221

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Liu Y, Gayle AA, Wilder-Smith A, Rocklöv J. The reproductive number of COVID-19 is higher compared to SARS coronavirus. J Travel Med. (2020) 27:taaa021. doi: 10.1093/jtm/taaa021

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Lan L, Xu D, Ye G, Xia C, Wang S, Li Y, et al. Positive RT-PCR test results in patients recovered from cOVID-19. JAMA. (2020) 323:1502–3. doi: 10.1001/jama.2020.2783

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Xing Y, Mo P, Xiao Y, Zhao O, Zhang Y, Wang F. Post-discharge surveillance and positive virus detection in two medical staff recovered from coronavirus disease 2019 (COVID-19), China, January to February 2020. Eurosurveillance. (2020) 25:2000191. doi: 10.2807/1560-7917.ES.2020.25.10.2000191

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Omer SB, Malani P, del Rio C. The COVID-19 Pandemic in the US: a clinical update. JAMA. (2020) 323:1767–8. doi: 10.1001/jama.2020.5788

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Wang Q, Shi N, Huang J, Cui T, Yang L, Ai J, et al. Effectiveness and cost-effectiveness of public health measures to control COVID-19: a modelling study. medRxiv. (2020). doi: 10.1101/2020.03.20.20039644

CrossRef Full Text | Google Scholar

31. Chen SC, Liao CM. Cost-effectiveness of influenza control measures: a dynamic transmission model-based analysis. Epidemiol Infect. (2013) 141:2581–94. doi: 10.1017/S0950268813000423

PubMed Abstract | CrossRef Full Text | Google Scholar

32. White NM, Barnett AG, Hall L, Mitchell BG, Farrington A, Halton K, et al. Cost-effectiveness of an Environmental Cleaning Bundle for Reducing Healthcare-associated Infections. Clin Infect Dis. (2019) 70:2461–3. doi: 10.1093/cid/ciz717

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. (2020) 382:1708–20. doi: 10.1056/NEJMoa2002032

CrossRef Full Text | Google Scholar

34. Gudbjartsson DF, Helgason A, Jonsson H, Magnusson OT, Melsted P, Norddahl GL, et al. Spread of SARS-CoV-2 in the icelandic population. N Engl J Med. (2020) 382:2302–15. doi: 10.1101/2020.03.26.20044446

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet. (2020) 395:689–97. doi: 10.1016/S0140-6736(20)30260-9

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Shen M, Peng Z, Xiao Y, Zhang L. Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China. bioRxiv. (2020). doi: 10.1101/2020.01.23.916726

CrossRef Full Text | Google Scholar

37. Liu T, Hu J, Kang M, Lin L, Zhong H, Xiao J, et al. Transmission dynamics of 2019 novel coronavirus (2019-nCoV). bioRxiv. (2020). doi: 10.2139/ssrn.3526307

CrossRef Full Text | Google Scholar

38. Read JM, Bridgen JR, Cummings DA, Ho A, Jewell CP. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. medRxiv. (2020). doi: 10.1101/2020.01.23.20018549

CrossRef Full Text | Google Scholar

39. Majumder M, Mandl KD. Early Transmissibility Assessment of a Novel Coronavirus in Wuhan, China. Available online at SSRN: https://ssrncom/abstract=3524675 (accessed January 5, 2020).

PubMed Abstract | Google Scholar

40. Cao Z, Zhang Q, Lu X, Pfeiffer D, Jia Z, Song H, et al. Estimating the effective reproduction number of the 2019-nCoV in China. medRxiv. (2020). doi: 10.1101/2020.01.27.20018952

CrossRef Full Text | Google Scholar

41. Zhao S, Lin Q, Ran J, Musa SS, Yang G, Wang W, et al. Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: a data-driven analysis in the early phase of the outbreak. Int J Infect Dis. (2020) 92:214–7. doi: 10.1101/2020.01.23.916395

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Imai N, Cori A, Dorigatti I, Baguelin M, Donnelly CA, Riley S, et al. Report 3: Transmissibility of 2019-nCov. COVID 19 Resources, Imperial College London COVID-19 Reports (2020).

Google Scholar

43. Riou J, Althaus CL. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020. Eurosurveillance. (2020) 25:2000058. doi: 10.2807/1560-7917.ES.2020.25.4.2000058

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, et al. The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Anna Int Med. (2020) 172:577–82. doi: 10.7326/M20-0504

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet. (2020) 395:1054–62. doi: 10.1016/S0140-6736(20)30566-3

CrossRef Full Text | Google Scholar

46. Seto WH, Tsang D, Yung RWH, Ching TY, Ng TK, Ho M, et al. Effectiveness of precautions against droplets and contact in prevention of nosocomial transmission of severe acute respiratory syndrome (SARS). Lancet. (2003) 361:1519–20. doi: 10.1016/S0140-6736(03)13168-6

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Marino S, Hogue IB, Ray CJ, Kirschner DE. A methodology for performing global uncertainty and sensitivity analysis in systems biology. J Theoret Biol. (2008) 254:178–96. doi: 10.1016/j.jtbi.2008.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Overholser BR, Sowinski KM. Biostatistics Primer: Part 2. Nutrit Clin Pract. (2008) 23:76–84. doi: 10.1177/011542650802300176

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Koo JR, Cook AR, Park M, Sun Y, Sun H, Lim JT, et al. Interventions to mitigate early spread of SARS-CoV-2 in Singapore: a modelling study. Lancet Infect Dis. (2020) 20:678–88. doi: 10.1016/S1473-3099(20)30162-6

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Yin WW, Gao LD, Lin WS, Du L, Zhang XC, Zou Q, et al. Effectiveness of personal protective measures in prevention of nosocomial transmission of severe acute respiratory syndrome. Zhonghua Liu XIng Bing Xue Za Zhi. (2004) 25:18–22.

PubMed Abstract | Google Scholar

51. Zhang N, Li Y. Transmission of Influenza A in a student office based on realistic person-to-person contact and surface touch behaviour. Int J Environ Res Public Health. (2018) 15:1699. doi: 10.3390/ijerph15081699

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Mendell MJ, Fisk WJ, Kreiss K, Levin H, Alexander D, Cain WS, et al. Improving the health of workers in indoor environments: priority research needs for a national occupational research Agenda. Am J Public Health. (2002) 92:1430–40. doi: 10.2105/AJPH.92.9.1430

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kai D, Goldstein GP, Morgunov A, Nangalia V, Rotkirch A. Universal masking is urgent in the covid-19 pandemic: Seir and agent based models, empirical validation, policy recommendations. arXiv [Preprint] arXiv:200413553 (2020). doi: 10.13140/RG.2.2.21662.08001

CrossRef Full Text | Google Scholar

Keywords: SARS-CoV-2 pandemics, mathematical modeling, reopening business, parameter calibration, benefit-cost-risk trade-off

Citation: Miao H, Gao Q, Feng H, Zhong C, Zhu P, Wu L, Swartz MD, Luo X, DeSantis SM, Lai D, Bauer C, Pérez A, Rong L and Lairson D (2020) Mathematical Modeling of Business Reopening When Facing SARS-CoV-2 Pandemic: Protection, Cost, and Risk. Front. Appl. Math. Stat. 6:35. doi: 10.3389/fams.2020.00035

Received: 31 May 2020; Accepted: 13 July 2020;
Published: 11 August 2020.

Edited by:

Sen Pei, Columbia University, United States

Reviewed by:

Xin Wang, Beihang University, China
Steve Newbold, University of Wyoming, United States
Andreas Eilersen, University of Copenhagen, Denmark

Copyright © 2020 Miao, Gao, Feng, Zhong, Zhu, Wu, Swartz, Luo, DeSantis, Lai, Bauer, Pérez, Rong and Lairson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hongyu Miao, hongyu.miao@uth.tmc.edu

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.