Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

On the optimal presence strategies for workplace during pandemics: A COVID-19 inspired probabilistic model

  • Mansoor Davoodi ,

    Roles Formal analysis, Investigation, Methodology, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    m.davoodi-monfared@hzdr.de

    Affiliation Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden Rossendorf (HZDR), Görlitz, Germany

  • Abhishek Senapati,

    Roles Writing – original draft, Writing – review & editing

    Affiliation Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden Rossendorf (HZDR), Görlitz, Germany

  • Adam Mertel,

    Roles Software, Writing – original draft, Writing – review & editing

    Affiliation Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden Rossendorf (HZDR), Görlitz, Germany

  • Weronika Schlechte-Welnicz,

    Roles Writing – original draft, Writing – review & editing

    Affiliation Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden Rossendorf (HZDR), Görlitz, Germany

  • Justin M. Calabrese

    Roles Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden Rossendorf (HZDR), Görlitz, Germany, Department of Ecological Modelling, Helmholtz Centre for Environmental Research – UFZ, Leipzig, Germany, Department of Biology, University of Maryland, College Park, MD, United States of America

Abstract

During pandemics like COVID-19, both the quality and quantity of services offered by businesses and organizations have been severely impacted. They often have applied a hybrid home office setup to overcome this problem, although in some situations, working from home lowers employee productivity. So, increasing the rate of presence in the office is frequently desired from the manager’s standpoint. On the other hand, as the virus spreads through interpersonal contact, the risk of infection increases when workplace occupancy rises. Motivated by this trade-off, in this paper, we model this problem as a bi-objective optimization problem and propose a practical approach to find the trade-off solutions. We present a new probabilistic framework to compute the expected number of infected employees for a setting of the influential parameters, such as the incidence level in the neighborhood of the company, transmission rate of the virus, number of employees, rate of vaccination, testing frequency, and rate of contacts among the employees. The results show a wide range of trade-offs between the expected number of infections and productivity, for example, from 1 to 6 weekly infections in 100 employees and a productivity level of 65% to 85%. This depends on the configuration of influential parameters and the occupancy level. We implement the model and the algorithm and perform several experiments with different settings of the parameters. Moreover, we developed an online application based on the result in this paper which can be used as a recommender for the optimal rate of occupancy in companies/workplaces.

Introduction

During the pandemics, particularly the very recent coronavirus (SARS-CoV-2, hereafter COVID-19) pandemic, most of the services provided by the companies, organizations and governments’ facilities have been affected significantly in terms of quality and quantity [16]. On one hand, most of the employees have been at risk of infection and on the other hand, the work-flows and functions have become slow. As COVID-19 spreads mostly through direct/close contacts between individuals, the risk of getting infected remains despite adopting some recommended personal hygiene measures in the workplace (e.g., wearing a mask, maintaining distance, etc.). To control and overcome these issues, most of the companies mainly followed three strategies: (i) reducing the number of on-site employees by teleworking, (ii) screening employees via periodically testing, and (iii) vaccinating the employees [79]. These strategies have worked quite well to decrease the number of infected cases among the employees, however, the efficiency and the productivity of the employees, particularly who use to work in collaborative environment or where the physical presence of the employee is necessary may get affected while working from home [1, 10]. Therefore, the managers mostly prefer to increase the availability of the employees in the offices to increase work and service productivity. It is very evident that there is a potential trade-off between productivity and the risk of infection, and an important question that arises here is “what would be the optimal presence rate in the workplace during the pandemic?”.

During the pandemic, many companies and institutions tried to adapt to the new situation by lowering the possibilities of getting infected at the workplace and keep the employees safe during the work time. Therefore, implementation of effective strategies in the workplace regarding the duration of working hour, presence rate of the employees, level of personal protection for avoiding infection is very important from the management perspective. An ample amount of mathematical models have been developed specifically to investigate several aspects of transmission mechanism and control strategies of COVID-19 [1115]. In addition, several mathematical and statistical model based investigation have also been carried out to gain better insights on the outcomes of the proposed strategies particularly tailored to different organizations. Most of these studies focused on implementing organizational strategies mainly in the workplaces like hospitals and health care facilities, nursing homes, offices, schools, etc. These studies investigated several aspects like effectiveness of control measures in preventing the disease transmission, scheduling work place arrangements to minimize the consequence of infection, strategies for reopening different activities in organizations, etc.

Regarding effectiveness of control measures and surveillance several attempts have been made by the previous studies with the help of mathematical models to quantify the optimal frequency of testing in healthcare environment [16, 17], nursing home [18], timeliness of contact tracing and importance of molecular in preventing disease transmission among health workers and other high risk groups [19], performance of combinatorial group testing [20, 21], importance of self-isolation and contact tracing measures [22], etc. Also, [23] proposed an optimal testing strategy that essentially minimizes the presence of pre-symptomatic and asymptomatic employees in the workplace.

In the context of scheduling the workplace arrangement, a desynchronization strategy has been proposed in which the workers in the healthcare system can be divided into two non-overlapping teams who will be working in alternating weeks to increase the workforce productivity in the facility [24]. Also incorporating the source of infection from inside and outside the hospital, [25] studied the effectiveness of regular testing and desynchronization protocol in preventing COVID-19 infection. The effect of cyclic 4-day work and 10-day lockdown strategy in workplace has been investigated by [26]. In nursing home and long-term care facilities, [27] constructed a mathematical model based on bipartite network consists of health care workers and the residents and investigated how the restructured interactions in the bipartite network can change the outcome of the epidemics in the facility by exploiting the strategy of shied immunity, where the recovered people increase their interactions with susceptible population to protect them from the infection. In order to support the decision-making process, [28] developed a strategy to predict the risk of infection among the health workers and other employees in the essential sectors by taking into account the working process, proximity between employees, type of activity, etc.

Several recommendations have been proposed in reopening different activities in workplace. How enhanced level of testing, contact-tracing, adoption of facial mask, home quarantine can play a significant role in relaxing social distancing interventions, which is not always feasible to maintain in the facilities like hospital [29, 30]. It has been recommended that with the adoption of large-scale trace and test intervention along with moderate social distancing, the schools in Île-de-France region can be reopened with maximum attendance of 50% that could avoid second wave [31]. How designing an optimal employee screening strategy can reduce the on-site infection in a workplace and based on this strategy, the feasibility of return to work policy are investigated in [32].

From the above discussions, it is clear that most of the studies have focused on quantifying the impact of implementing different measures in work to reduce infection by developing mathematical and statistical models. Also, some of the studies are devoted to recommending strategies for scheduling workplace arrangements and reopening different activities. However, the aspect of how the influential factors related to preventive measures such as testing frequency, contact rate, and vaccination rate of employees play significant roles in quantifying the optimal presence rate as well as home-office productivity factor paid no attention in the previous studies. Motivated by this circumstance, in this paper, we develop a model of disease transmission within companies and organizations that allows us to identify the optimal balance of the mentioned components. To this end, we incorporate influential factors such as the local incidence level, transmission rate of COVID-19, testing frequency, rate of vaccination and rate of contact among the employees. These factors are the most important determinants of the optimal strategy for home-vs-office work. In addition to them, one crucial factor is the work efficiency at home compared to the work efficiency at the office. This factor varies from organization to organization depending on the type of organization and their inputs, processes and outputs. For the sake of simplicity, hereafter, we call the ratio of efficiency at home to efficiency at office by productivity.

The proposed model in this paper is based on a probabilistic framework that allows to model spreading of the COVID-19 under a setting of the mentioned influential factors and computes the expected number of infected employees when some infection arrives at the facility. As the advantages of this probabilistic framework, it is simple, fast and practical with the ability to be customized for a specific organization. Moreover, the proposed probabilistic model can be efficiently used for any size of the population. Based on this model, we developed an optimization application to find Pareto optimal solutions for the objectives of maximizing productivity and minimizing infection risk in organizations and companies. The model and the graphical interface of the application are implemented in such a way that it provides the users full ability to reassign all input parameters and renders the results on any change.

This paper is organized into four sections. In the next section, we introduce the influential parameters and state the problem. In addition, we propose the probabilistic analysis and algorithm to compute the risk of infection as well as the Pareto optimal solutions. In the third section, we perform several experiments and test the model and algorithm in deriving Pareto optimal solutions. Finally, we finished the paper in the fourth section by drawing a conclusion and future directions.

Modeling the office/workplace presence problem and solution approach

In the last two decades, teleworking and home office have rarely been applied for some companies, organizations, corporations, and government facilities (for the sake of simplicity, in the rest of this paper we use the company to refer to all of them) regarding the type of services they provide. There are comparison studies that discuss the advantages and disadvantages of teleworking [33, 34]. However, since the beginning of the COVID-19 pandemic, most companies have been forced to follow home office partially. This usually has decreased their productivity rate, but significantly reduces the effect of COVID-19 by preventing (or at least reducing) the spreading of the virus in the facilities. Therefore, there is a trade-off between the rate of presence at the offices and the rate of infection at the offices, and it is more challenging for the companies that their employee work at the offices with a higher rate of productivity compared to working at home.

In this section, we formulate this trade-off as a bi-objective optimization problem to find the optimal rate of presence by considering the risk of infection at the facility. It is notable if the productivity rate of employees of some company (because of the type of task and process it provides) is not reduced by following the home-office strategy, clearly, there is no trade-off between the risk of infection in the offices and productivity, and the optimal trivial strategy is minimum occupancy in the offices. Note that, following this strategy, there is still some possibility of infection for the employees working at home, however, there is no propagation of the disease in the office space. Thus, in this paper, we only consider the companies for which the total productivity is increased by the availability of their employees in the offices. It is notable that there will be still some natural risk for the employees working at home, but the companies do not have direct control over this.

Table 1 shows the key parameters and variables used in the model. The exact definition of them and the model’s assumptions are explained as follows.

thumbnail
Table 1. The parameters and variables used in the proposed model.

https://doi.org/10.1371/journal.pone.0285601.t001

Parameters, variables, and assumptions of model

  • n: The total number of employees who regularly work in the facility and the productivity in the facility would be affected if they got infected.
  • nv: The number of employees who are (fully) vaccinated with one of the available COVID-19 vaccines. We count partially vaccinated (e.g., one dose) individuals as non-vaccinated
  • nI: The expected number of employees who are infected and can spread the disease in the facility. As soon as they do a test or the disease’s symptoms appear, they should be in quarantine, so the facility does not count them as an active employee.
  • βu and βv: The probability of the virus transmitting from one infected individual to one unvaccinated (βu) or one vaccinated (βv) individual per one potential contact.
  • prod: The productivity factor is the productivity of working at home compared to working at the office. This factor varies from facility to facility and from employee to employee, based on their mission and the type of services they provide. Managers (or any decision-maker) can determine this factor based on their experience in working under the COVID-19 pandemic and comparing it with normal situations. Precisely, . So, prod = 1 means the manager believes that there is no difference between productivity for working at home and office, and prod = 0 means the manager strictly prefers the employees work at the office.
  • τ: On average, how often do the employees conduct a test to detect COVID-19? Precisely, it is the average time interval between two sequential tests and significantly helps to detect the infection before it highly spreads in the facility. Since the symptoms of the COVID-19 disease appear within two weeks [35, 36], we assume τ varies from one day to 14 days (incubation period). It is worth mentioning, for the sake of simplicity, in this paper, we assume the tests are perfect and no sensitivity exists. Otherwise, as an advantage of the proposed model, we can easily apply the test’s error as a coefficient in computing the expected time to detect the infection.
  • ρ: The probability that some infections arrive at the facility per day. This can be computed by considering the number of employees at the facility and the incidence level in the neighborhood of the facility. For example, the number of infections for the last 7 days is reported daily (e.g., see https://www.coronavirus.sachsen.de). So, by utilizing it, we can compute such probability. However, to achieve a more realistic value, we apply the ratio of the local vaccination rate to the employees’ vaccination rate.
  • κ: The average number of contacts among the employees that have the potential to transmit the disease. We assume the employees and clients mostly wear masks and keep a distance of at least 1.5 meters when they contact each other.
  • occup: The rate (percentage) of employees present at the office. This is the only decision variable of the model which should be determined optimally.

As mentioned, all the variables and parameters are given except the decision variable occup which is the output of the model. Indeed, the employees are divided into three groups: (i) who work at the office, (ii) who work from home, and (iii) who are infected. The productivity factor of the first group is one, and that of the second and the third groups are prod and zero, respectively. So, nnI is the number of healthy and active employees and out of them occup percentage of employees who works at the offices with full productivity equals one, and the remaining 1 − occup percentage of them work from home with productivity prod. Therefore, the total productivity of a facility can be determined by the following formula. Note that it is possible to customize this formula and define it based on the input/fellows/output in a company. In this paper, we skip such details and only focus on this basic and general definition of total productivity. (1)

Therefore, one objective function of the model is maximizing the total productivity, and the other one is minimizing the expected number of infections, . There is a clear conflict between these two objectives. Note that, prod is the productivity factor determined by the decision-maker. Without loss of generality, we assume it is between zero and one, and if prod ≥ 1, there is a trivial optimal solution occup = 0, that is working from home strictly is preferred to working at the office, and it results in minimizing the number of infected employees. The most important part of the objective functions is computing the expected number of infected employees, , by arriving some infections at the facility. Indeed, we need to compute it as a function over time and based on the influential parameters mentioned in Table 1 which are explained in detail in the next subsections.

Probabilistic framework for computing the number of infected employees

We start by a full contact network of unvaccinated employees and present the approach. Then we extend it to two groups of vaccinated and unvaccinated employees. Employees can contact each other with the same probability, however, the average number of contacts per day is bounded by κ. We suppose an infection arrives at the facility at the time (day) t = 0, and compute the probability of an arbitrary individual being infected after Δt days when the transmission rate of the virus is βu and the number of contacts is κ per day. Let denote this value with the probability function PIt, n, βu, κ). By having this probability function, the expected number of infected individuals can be easily computed by (2)

For the sake of simplicity, in the rest of the paper, let denote the probability function PIt, n, βu, κ) by PIt) for known variables of n and κ. Thus, computing the probability of infection per employee results in achieving the expected number of infected employees. Let denote the source of infection at t = 0 by s. Let u be an employee who stays healthy till Δt − 1. There are two ways for infecting u at day Δt; via some direct contact with s, or via contact with one of the n − 2 other employees. The transmission probability for each of the contacts is βu, and the probability of infection for the source is one, while the probability of infection for the other employees is PIt − 1). See Fig 1. That means the probability of an arbitrary employee like u staying healthy after c contacts with s is (1 − 1 × βu)c. While the probability of u stays healthy after c contacts with the other n − 2 employees in day Δt is (1 − PIt − 1) × βu)c. Therefore, the infection probability can be presented as a recursive function. Since we assume κ contacts for each employee per day and all of them have an equal chance to occur, this function can be written as (3) where is the expected number of contacts between employee u and the infected source employee s, and is the expected number of contacts between employee u and the other employees except s. Note that, (1 − PIt − 1)) is the probability of the employee staying healthy till day Δt − 1.

thumbnail
Fig 1. Infection possibilities of n employees which one of them (s) is infected.

https://doi.org/10.1371/journal.pone.0285601.g001

Eq 3 computes the probability of an arbitrary employee getting infected after t = Δt days when the employees have κ contacts per day and present at the office with occup = 1, that is, full occupancy. For 0 ≤ occup ≤ 1 as the rate of presence at the office, all the employees (include s) stay (healthy) at home with the expected probability 1 − occup. So the probability of having one contact between u and s is occup2. Therefore, the infection probability formula for known n, κ and occup can be extended as follows (4)

Also, we set PI(0) = 0 as the scenario supposed here. After computing the probability of infection, the expected number of infected employees can be computed by applying Eq 2 easily.

Now, let’s extend the above computation to two different groups, say vaccinated and unvaccinated employees. We consider a probability of the disease transmission βv (βvβu) for the vaccinated group. Without loss of generality, we consider the source of infection s as an unvaccinated employee and compute the probability of infection in two cases (i) for an unvaccinated employee u, and (ii) for a vaccinated employee v (see Fig 2). There are three ways for infecting u; via direct contact with s, via contacting with one of the nv vaccinated employees, and via contact with one of the nnv − 2 unvaccinated employees. The transmission probability for each of such contacts is βu because u is an unvaccinated employee, however, the transmission probability for the vaccinated group is βv. Let and be the probability of infection for an unvaccinated employee and a vaccinated employee, respectively. Thus, the recursive equation for them can be written below (5) and (6) where . One advantage of these recursive formulas is that and can be efficiently computed in linear time to Δt using a simple bottom-up approach. Finally, the expected number of infected employees after Δt days of introducing an infection can be computed as follows (7)

thumbnail
Fig 2. Infection possibilities of n employees which one of them (s) is infected and there are two groups of vaccinated and unvaccinated employees.

The left subfigure shows the scenario for an unvaccinated employee u and the right one shows the scenario for a vaccinated employee v.

https://doi.org/10.1371/journal.pone.0285601.g002

Thus, Eq 7 provides a recursive formula to compute the number of infections over time. An implementation of this equation and comparison with simulation results is presented in the Appendix.

Computing optimal presence rate

As explained, the expected number of infected employees after arriving some infection to the facility can be computed using Eq 7 as a function over time and based on the mentioned influential parameters. For the sake of simplicity, we denote this expected number by . Let ρ be the probability of infection arriving at the facility. This probability can be determined using two straightforward approaches. First, using the (recent) historical data of incidence in the company, and second, using the local incidences and applying the ratio of the vaccination rate in the neighborhood of the facility to the vaccination rate in the company. Therefore, the cumulative number of infected employees for a time interval of T days can be computed as follows (8) where Z(T) is defined as below (9)

Note that Z(T) is a linear recursive formula to compute the cumulative number of infections for t = 0, 1, …, T by removing the expected overlapped infected employees over time.

Eq 8 computes the number of all infected in a time interval T days where the probability of arriving an infection is ρ per day. If the employees do a test every τ days on average, and they uniformly distribute in the test interval (having tests per day on average), the expected days to detect the infection can be computed using the following equation (10) where is the probability of infection per employee after t days of arriving some infection, so, 1 − (1 − Pr(t))k is the probability of detecting at least one infected employee by testing a subgroup of k after t days an infection arrives the facility. For detecting the infection after t > 1 days of arriving, we need to compute the probability of it has been not detected in days before t, i.e. t′ < t, which is computed using .

in Eq 10 shows the expected time to detect an infection if the employees do a test every τ days on average. Therefore, the expected number of infected employees before detecting an infection is obtained by . So, the total productivity (see Eq 1) of a facility can be determined for any given occup in [0, 1], and based on the other introduced influential parameters.

A practical multi-objective solution approach

As described before, the optimal office presence strategy is formulated in the framework of a bi-objective optimization problem as follows (11)

The first objective is minimizing the expected number of infected employees presented in Eq (8), and the second objective is maximizing the total productivity of the company presented in Eq (1). The only constraint which we considered in the model is related to minimum possible occupancy in the companies, Occupancy Threshold. That is a limit on the minimum number of employees that have to be present at the facility in order to process some tasks for which the physical presence of the employees is required. Further, it is possible to add more lower- or upper-bound constraints on the occupancy variable. The model has two conflicting objectives. Clearly, by increasing the occupancy in a company, the expected number of infected employees will increase as well. However, this is not the case for the second objective, because by increasing the occupancy, the number of infected employees will also increase, and consequently, since the productivity of the infected employees is zero, the total productivity will decrease. Fig 3 shows the effect of occupancy on the first and second objective separately for two different scenarios of the test interval, one test per week (blue curve) and one test per two weeks (red curve). As expected, the number of infections is a strictly increasing function of the rate of presence, while, the productivity function is increasing for the low rate of presence and after reaching the maximum productivity (e.g., for occup ≈ 88% in the first scenario and for occup ≈ 65% in the second one) it decreases by increasing presence rate.

thumbnail
Fig 3. Effect of occupancy on the expected number of infections (left panel) and on the total productivity (right panel) for two intervals, one test per week and one test per two weeks.

In these simulations, we assumed a company with 100 employees which half of them are vaccinated, with 15% contact rate, and the home productivity rate 0.90.

https://doi.org/10.1371/journal.pone.0285601.g003

The outcome of the optimization problem presented in model 11 is a set of trade-off solutions– called Pareto optimal– a solution that is not improved for an objective unless it sacrifices the other objective [37]. There are different approaches to solving multi-objective optimization problems, such as the weighted-sum method [38], the lexicographic method [39], the ϵconstraint approach [40], the goal programming [41], evolutionary algorithms [37, 42], etc. All of these approaches have their own advantages and disadvantages. Some of them have high time complexity, some of them work only for differentiable functions or convex models, some of them have no guarantee to find all the Pareto-optimal solutions. Thus, we propose a practical quick method to compute all Pareto optimal solutions of the model 11 efficiently.

The method utilizes some observations to solve the model. The, occup is a continuous variable, but in reality, there are at most n possible choices for the number of employees present in the facility. So, at first, for any which satisfy the occupancy threshold, we compute the expected detection time based on Eq (10), and the expected number of infected employees and its corresponding total productivity based on Eqs (1) and (7), respectively. All of these steps can be handled in O() time using the presented linear recursive equations. Let denote this solution set by S. Now, we can find the non-dominated solutions of S and report them as the Pareto-optimal optimal solutions of the problem. A solution is a non-dominated solution, if there is no other solution in S such that it is better than in both objectives, i.e., its expected number of infections is less than the expected number of infections of , and simultaneously, its productivity is more than the productivity of . Since the solutions of S are constructed one by one, from the minimum expected infections to the maximum one, the Pareto-optimal solutions of the problem can be computed in linear time using a sweep-line approach [43]. The pseudocode of the proposed algorithm is presented as follows.

Algorithm 1 Computing Pareto Optimal Office Presence Strategies

Input: Company’s Parameters (the notation mentioned in Table 1)

Output: All Pareto Optimal Strategies

max_prod ← −1

i ← 0

while in do

if occupOccupancy_Threshold then

  Compute using Eq (8)

  Compute expected number of infections using Eq (7)

  Compute the total productivity using Eq (1) and denote it by TP

  if TP > max_prod then

   Report the current strategy as a Pareto optimal strategy

   max_prodTP

  end if

end if

ii + 1

end while

Simulation results and discussions

In this section, we show some implementation results of the proposed model and method for computing Pareto optimal solutions of the objectives, minimizing the expected number of infections and maximizing productivity. Since there are several influential input parameters in the proposed model, there are so many possible combinations of them as well. Therefore, we shortly display some outputs of the model, and instead, we provided an online optimization tool based on the model to give any possible combinations of the influential parameters as the input and observe the output. The basic version of this tool is available at www.where2test.de.

To run the simulations, we consider a middle-size company with n = 100 employees and suppose the 7-day incidence rate per 100,000 population is 500 individuals. Thus, the average probability of infection per employee in a week will be , and consequently, the probability of infection arrives at the company is . In this study, we specifically consider two scenarios based on the disease transmission rate for unvaccinated employees (βu). First, we choose the baseline value as βu = 0.04 from a possible range of values reported in the previous study [44]. This baseline scenario can be regarded as the situation where the different variants of COVID-19 like Alpha and Delta are the dominant variants and cause infections. On the other hand, the growing evidences [4547] suggest that the recently identified Omicron variant is more transmissible than the previous ones. Therefore, as the second choice, we consider the situation where the transmissibility increases 2.5-folds, i.e. βu = 0.1. Also, the transmission rate of vaccinated employees to βv = (1 − 0.80)βu. That means 80% immunity for vaccinated individuals is assumed which is in accordance with the efficacy of the available COVID-19 vaccines [4851].

Further, we assume two different values for the other input parameters. The home productivity factors, prod = 0.6 and prod = 0.9, and two vaccination rates of employees 0.4 and 0.8. For the test interval, τ, we consider the simulations for two test intervals, one test every week (i.e., τ = 7 days), and one test every two weeks (i.e., τ = 14 days), and for the number of contacts, κ, we consider low and high rate of contacts. The low-case scenario assumes every employee who presented at the office has κ = 5 + 0.10 × occup × n contacts with the other employees, where occup is the rate of presence at the offices. In fact, we assume 5 contacts (can be considered as unavoidable contacts) per employee as a constant number, plus 10 percent of the employees who are present at the offices. Likewise, the high rate scenario assumes the number of contacts is κ = 5 + 0.20 × occup × n. Note that, this is just a sample setting of the influential parameters and any other setting may replace regarding the availability of the information. Fig 4 displays the result of 15 scenarios for a different settings of input parameters, whereas Table 2 shows explains these 15 input settings.

thumbnail
Fig 4. Pareto optimal solutions for 15 different settings of input parameters.

See Table 2. In each subfigure, the horizontal axis shows the possible occupancy, the blue diagram and the left vertical axis illustrate the obtained total productivity, and the red diagram and the right vertical axis illustrate the expected number of infected employees per week.

https://doi.org/10.1371/journal.pone.0285601.g004

thumbnail
Table 2. 15 different settings of input parameters.

The corresponding Pareto optimal solutions are depicted in Fig 4.

https://doi.org/10.1371/journal.pone.0285601.t002

Each subfigure in Fig 4 shows Pareto optimal solutions and their objective values. The horizontal axis shows the possible occupancy values. The blue diagram and the left vertical axis illustrate the obtained total productivity, and the red diagram and the right vertical axis illustrate the expected number of infected employees per week. Note that, the Pareto optimal solutions are displayed only for the occupancy values which result in such solutions, and the dominated solutions are not displayed. These results together provide useful information for a decision-maker who may consider a threshold as a maximum risk of infection and try to find the maximum productivity by some changes in the vaccination rate of the employees, test interval, or even their number of contacts.

For example, when the 7-day incidences are 500 individuals of 100,000 population (this is almost the average number of incidences in Saxony, Germany from January first till the middle of February 2022, e.g., see https://www.where2test.de/saxony), it means the natural risk (or say the background risk) of infection per employee is per week, and for the whole of the company it is . The background risk can be interpreted as the risk of infection if the employees do not come to the office and stay at home for work and maintain the normal social contacts by e.g., going to restaurants and shopping. While the number of incidences is not so high or so low, the decision-makers and managers can use it as a mental guidance like a reference point in determining a proper threshold for the risk of infection in the company. The horizontal green line in each of the subfigure displays such background risk.

For example, if a decision-maker (say a manager in a company) would like to follow a presence strategy whose result never exceeds the natural risk, he/she can consider the occupancy corresponding to the intersection point between the background risk and the computed infection risk, i.e., the green dashed line and the red diagram. For instance, such occupancy in the first scenario in Fig 4(a) (τ = 7, prod = 0.6, vaccination rate = 0.5, Low contact rate and βu = 0.04) corresponds with occupancy 64%, and it will result in 85% of productivity, while in the second scenario, Fig 4(o)(τ = 14, prod = 0.9, vaccination rate = 0.8, High contact rate and βu = 0.10), it is 46% with more than 94% of productivity.

In a reverse usage of the Pareto optimal solutions, a decision-maker may be interested to know what is the risk of infection (e.g., compared to the natural risk) if he/she would like to achieve a particular level of productivity in the company. For instance, in the second scenario, Fig 4(b)(τ = 7, prod = 0.6, vaccination rate = 0.8, Low contact rate and βu = 0.04), a level of 70% of the productivity can be achieved with only half of the background risk, i.e., 0.2, while only 25% of employees presence at offices. To compute this, we first find the intersection point of productivity 70% with the productivity diagram (the blue diagram). In the second scenario, it is almost 25% occupancy. Then find the infection risk on the red diagram which is corresponding with such occupancy, which is almost 0.2.

Finally, decision-makers can use several diagrams together to find a sense of the amount of obtained productivity or infection risk by changing some influential parameters, such as test interval and/or vaccination rate among the employees. For instance, if the current situation of a company is similar to the 13th scenario (see Fig 4(l)), and occupancy 54% will result in a risk of infection is the same as the background risk and the productivity 81%. Now, if the test interval among the employees rises to τ = 14 days, (see 15-th scenario, Fig 4(l)), and the decision-maker still would like to the risk of infection in the company does not exceed the natural risk, he/she have to decrease occupancy to almost 46% and so achieve a productivity level of 78%. A similar analysis can be extracted for investigating the effect of other input parameters. As aforementioned, because of the so many possible combinations of the influential input parameters, there is no way to discuss all sensitivity analyses and observe the changes in the objectives by changing the parameters. Instead, we suggest using the provided online application for this purpose.

Conclusion and future directions

In this study, we developed a model to compute the optimal trade-off solutions for reducing the risk of infection and increasing the total productivity in companies during the pandemics such as COVID-19. By including fundamentally important parameters such as the local incidence level, the number of contacts among the employees, their average test interval, and vaccination rate, we proposed a probabilistic analysis approach to compute the expected number of infected employees over the course of time. This probabilistic approach can be used to compute the number of infected individuals and probability of infection over time for different groups of people in terms of the probability of infection and spreading the virus.

We assumed two groups of employees with different infection probabilities; vaccinated and unvaccinated, and presented a practical approach to compute the Pareto optimal presence rate of employees to reach the maximum productivity and minimum infected risk. In addition to the incidence level, three critical parameters to obtain the maximum productivity are the home productivity rate of the employees, the contact rate between the employees and the average test interval among the employees. We designed the model as simple as possible to be able to cover and interpret the effect of all influential parameters. The presented approach is linear in terms of time complexity, and it can be simply extended to consider the sensitivity of COVID-19’s tests.

The manager of the business and organization can benefit from the outcome of the model in decision processes such as self-testing regulations of the employees and rate of remote working. The Pareto-optimal solutions illustrated a trade-off range between infection risk and productivity, which considerable managers in their decision-making. Further, the model can be extended to more than two groups of employees, e.g, different age groups. A basic implementation of this model and algorithms are available online at https://test.where2test.de. Future extensions of this study may consider the following subjects

  • In addition to employees, visitors (i.e., clients and customers) who have direct contact with the employees can be considered for computing the probability of an infection arriving at the facility.
  • Heterogeneous groups of employees with different contact rates, physical networks, productivity factors and test frequencies can be considered.
  • We defined the total productivity of the company based on individual productivity, that was, the summation of all employees’ productivity. However, it is possible for some companies the total productivity is defined based on the type of fellows and tasks in different sections of the company. To this end, the physical graph of contacts and more details on the outcomes may need as given.

Appendix: Evaluation of the proposed probabilistic approach

However, we explained the theory behind the probabilistic approach, in this section, we implement the proposed probabilistic approach for computing the expected number of infected employees as well. To this end, we assume different settings of input parameters, and compute the expected number of infections along the time. Indeed, we suppose one infection at time t = 0, and compute the expected infections for one month, i.e., t < 30 using Eq 7. Fig 5 illustrates the results for different values of βu, κ, occup and nv. In these results, we assumed n = 150 employees and βv = 0.15βu. In each subfigure, three cases of transmission rates, βu = 0.05, βu = 0.10 and βu = 0.15, are shown. Also, for comparison, we show simulation results (the dots), that is, simulating the companies using a set of n = 150 agents (employees) based on the parameter’s setting. To reduce the effect of random number generators, we ran the simulations for 100 times and report the average number of infected agents. Table 3 shows the mean absolute percentage error for each figure and the average of them. As it can be seen, the average (the last row of the Table) percentage of difference between the results is 5.6%, 3.8% and 2.7% for transmission rates βu = 0.05, βu = 0.10 and βu = 0.15, respectively. The results show the probabilistic analysis estimates the number of infected individuals with high accuracy, while it is fast and flexible approach to apply for heterogeneous group of people with different transmission rate and number of contacts.

thumbnail
Fig 5. The simulation results and the probabilistic analysis of the expected number of infected employees for a facility with n = 150 employees along time, Δt = 1, 2, …, 29 days, after the infection arrives.

Six different settings of the influential parameters are shown. The computations are based on Eq 7.

https://doi.org/10.1371/journal.pone.0285601.g005

thumbnail
Table 3. Mean absolute percentage error between the simulation results (average of 100 iterations) and the proposed probabilistic analysis for estimating the number of infected employees for six different settings of the influential parameters.

https://doi.org/10.1371/journal.pone.0285601.t003

References

  1. 1. Yang L., Holtz D., Jaffe S., Suri S., Sinha S., Weston J., et al., The effects of remote work on collaboration among information workers, Nature Human Behaviour (2021) 1–12. pmid:34504299
  2. 2. Majer J., Udoh K., Beleke A., Dugisye A., Deepak K., Aimee S., et al., Operational challenges and considerations for COVID-19 research in humanitarian settings: a qualitative study of a project in eastern Democratic Republic of the Congo and South Sudan. Plos one 17, no. 6 (2022): e0267822. pmid:35771799
  3. 3. Martin L., Hauret L., Fuhrer C., Digitally transformed home office impacts on job satisfaction, job stress and job productivity. COVID-19 findings. Plos one, 17(3) (2020), p.e0265131.
  4. 4. Arntz M., Ben Yahmed S., Berlingieri F., Working from home and COVID-19: The chances and risks for gender gaps, Intereconomics 55 (6) (2020) 381–386. pmid:33281218
  5. 5. N. Bloom, P. Bunn, P. Mizen, P. Smietanka, G. Thwaites, The impact of COVID-19 on productivity, Tech. rep., National Bureau of Economic Research (2020).
  6. 6. Cifuentes-Faura Javier, Circular economy and sustainability as a basis for economic recovery post-COVID-19, Circular Economy and Sustainability, Springer, 2 (2) 1–7 (2022). pmid:34888561
  7. 7. WHO, Coronavirus disease (COVID-19): Health and safety in the workplace, https://www.who.int/newsroom/qa-detail/coronavirus-disease-covid-19-health-and-safetyin-the-workplace, [Accessed: December 26, 2021] (2020).
  8. 8. Ho K.-f., Ho K., Wong S., Cheung A., Yeoh E., Workplace safety and coronavirus disease (COVID-19) pandemic: survey of employees, Bulletin Of The World Health Organization.
  9. 9. Cirrincione L., Plescia F., Ledda C., Rapisarda V., Martorana D., Moldovan R. E., et al., COVID-19 pandemic: Prevention and protection measures to be adopted at the workplace, Sustainability 12 (9) (2020) 3603.
  10. 10. Tavares F., Santos E., Diogo A., Ratten V., Teleworking in Portuguese communities during the covid-19 pandemic, Journal of Enterprising Communities: people and places in the global economy.
  11. 11. Kucharski A. J., Russell T. W., Diamond C., Liu Y., Edmunds J., Funk S., et al., Early dynamics of transmission and control of COVID-19: a mathematical modelling study, The Lancet Infectious Diseases 20 (5) (2020) 553–558. pmid:32171059
  12. 12. Vespignani A., Tian H., Dye C., Lloyd-Smith J. O., Eggo R. M., Shrestha M., et al., Modelling COVID-19, Nature Reviews Physics 2 (6) (2020) 279–281. pmid:33728401
  13. 13. Zhu H., Li Y., Jin X., Huang J., Liu X., Qian Y., et al., Transmission dynamics and control methodology of COVID-19: a modeling study, Applied Mathematical Modelling 89 (2021) 1983–1998. pmid:32982019
  14. 14. Biala T., Afolabi Y., Khaliq A., How efficient is contact tracing in mitigating the spread of COVID-19? a mathematical modeling approach, Applied Mathematical Modelling 103 (2022) 714–730. pmid:34815616
  15. 15. Furati K., Sarumi I., Khaliq A., Fractional model for the spread of COVID-19 subject to government intervention and public perception, Applied mathematical modelling 95 (2021) 89–105. pmid:33619419
  16. 16. Chin E. T., Huynh B. Q., Chapman L. A., Murrill M., Basu S., Lo N. C., Frequency of routine testing for COVID-19 in high-risk healthcare environments to reduce outbreaks, MedRxiv.
  17. 17. Hoover C. M., Skaff N. K., Blumberg S., Fukunaga R. Aligning staff schedules, testing, and isolation reduces the risk of COVID-19 outbreaks in carceral and other congregate settings: A simulation study. PLOS Global Public Health, 3 (1), (2023) pmid:36962883
  18. 18. See I., Paul P., Slayton R. B., Steele M. K., Stuckey M. J., Duca L., et al., Modeling effectiveness of testing strategies to prevent coronavirus disease 2019 (COVID-19) in nursing homes—United States, 2020, Clinical Infectious Diseases 73 (3) (2021) e792–e798.
  19. 19. Grassly N. C., Pons-Salort M., Parker E. P., White P. J., Ferguson N. M., Ainslie K., et al., Comparison of molecular testing strategies for COVID-19 control: a mathematical modelling study, The Lancet Infectious Diseases 20 (12) (2020) 1381–1389. pmid:32822577
  20. 20. McDermott J. H., Stoddard D., Woolf P. J., Ellingford J. M., Gokhale D., Taylor A., et al., A nonadaptive combinatorial group testing strategy to facilitate health care worker screening during the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak, The Journal of Molecular Diagnostics 23 (5) (2021) 532–540. pmid:33549858
  21. 21. Arasli B., Ulukus S. Group testing with a graph infection spread model. Information, 14 (1), 48, (2023)
  22. 22. Kucharski A. J., Klepac P., Conlan A. J., Kissler S. M., Tang M. L., Fry H., et al., Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: a mathematical modelling study, The Lancet Infectious Diseases 20 (10) (2020) 1151–1160. pmid:32559451
  23. 23. Hernandez X., Valentinotti S., On an optimal testing strategy for workplace settings operating during the COVID-19 pandemic, medRxiv.
  24. 24. Sánchez-Taltavull D., Castelo-Szekely V., Candinas D., Roldán E., Beldi G., Modelling strategies to organize healthcare workforce during pandemics: application to COVID-19, Journal of Theoretical Biology 523 (2021) 110718. pmid:33862091
  25. 25. Sanchez-Taltavull D., Castelo-Szekely V., Murugan S., Rollenske T., Ganal-Vonarburg S. C., Buchi I., et al., Regular testing of asymptomatic healthcare workers identifies cost-efficient SARS-CoV-2 preventive measures, medRxiv.
  26. 26. Karin O., Bar-On Y. M., Milo T., Katzir I., Mayo A., Korem Y., et al., Adaptive cyclic exit strategies from lockdown to suppress COVID-19 and allow economic activity, MedRxiv.
  27. 27. Lucia-Sanz A., Magalie A., Rodriguez-Gonzalez R., Leung C., Weitz J. S., Modeling shield immunity to reduce COVID-19 transmission in long-term care facilities, medRxiv.
  28. 28. Iavicoli S., Boccuni F., Buresti G., Gagliardi D., Persechino B., Valenti A., et al., Risk assessment at work and prevention strategies on COVID-19 in Italy, PloS One 16 (3) (2021) e0248874. pmid:33740016
  29. 29. Aleta A., Martin-Corral D., Piontti A. P. y, Ajelli M., Litvinova M., Chinazzi M., et al., Modelling the impact of testing, contact tracing and household quarantine on second waves of COVID-19, Nature Human Behaviour 4 (9) (2020) 964–971. pmid:32759985
  30. 30. M. D’Orazio, Bernardini G., Quagliarini E., How to restart? an agent-based simulation model towards the definition of strategies for COVID-19 “second phase” in public buildings, arXiv preprint arXiv:2004.12927.
  31. 31. Di Domenico L., Pullano G., Sabbatini C. E., Boëlle P., Colizza V., Expected impact of reopening schools after lockdown on COVID-19 epidemic in île-de-France, MedRxiv.
  32. 32. Meier K., Curnow K. J., Vavrek D., Moon J., Farh K., Chian M., et al., A working model to inform risk-based back to work strategies, medRxiv.
  33. 33. Tavares A. I., Telework and health effects review, International Journal of Healthcare 3 (2) (2017) 30.
  34. 34. Nakrošienė A., Bučiūnienė I., Goštautaitė B., Working from home: characteristics and outcomes of telework, International Journal of Manpower.
  35. 35. Paul S., Lorin E., Distribution of incubation period of COVID-19 in the Canadian context: Modeling and computational study, medRxiv.
  36. 36. Qin J., You C., Lin Q., Hu T., Yu S., Zhou X.-H., Estimation of incubation period distribution of COVID-19 using disease onset forward time: a novel cross-sectional and forward follow-up study, Science Advances 6 (33) (2020) eabc1202. pmid:32851189
  37. 37. Coello C. A. C., Lamont G. B., Van Veldhuizen D. A., et al., Evolutionary algorithms for solving multi-objective problems, Vol. 5, Springer, 2007.
  38. 38. Kim I. Y., De Weck O. L., Adaptive weighted-sum method for bi-objective optimization: Pareto front generation, Structural and Multidisciplinary Optimization 29 (2) (2005) 149–158.
  39. 39. Rentmeesters M. J., Tsai W. K., K.-J. Lin, A theory of lexicographic multi-criteria optimization, in: Proceedings of ICECCS’96: 2nd IEEE International Conference on Engineering of Complex Computer Systems (held jointly with 6th CSESAW and 4th IEEE RTAW), IEEE, 1996, pp. 76–79.
  40. 40. Miettinen K., Nonlinear Multiobjective Optimization, Vol. 12, Springer Science & Business Media, 2012.
  41. 41. Chang C.-T., Multi-choice goal programming, Omega 35 (4) (2007) 389–396.
  42. 42. Deb K., Multi-objective optimisation using evolutionary algorithms: an introduction, in: Multi-objective evolutionary optimisation for product design and manufacturing, Springer, 2011, pp. 3–34. https://doi.org/10.1007/978-0-85729-652-8_1
  43. 43. Jensen M. T., Reducing the run-time complexity of multiobjective eas: The NSGA-II and other algorithms, IEEE Transactions on Evolutionary Computation 7 (5) (2003) 503–515.
  44. 44. Lelieveld J., Helleis F., Borrmann S., Cheng Y., Drewnick F., Haug G., et al., Model calculations of aerosol transmission and infection risk of COVID-19 in indoor environments, International Journal of Environmental Research and Public Health 17 (21) (2020) 8114. pmid:33153155
  45. 45. Yang W., Shaman J., SARS-CoV-2 transmission dynamics in South Africa and epidemiological characteristics of the Omicron variant, medRxiv.
  46. 46. WHO, Enhancing response to Omicron SARS-CoV-2 variant: Technical brief and priority actions for member states (2022).
  47. 47. Lyngse F. P., Mortensen L. H., Denwood M. J., Christiansen L. E., Møller C. H., Skov R. L., et al., SARS-CoV-2 Omicron VOC transmission in Danish households, medRxiv.
  48. 48. Polack F. P., Thomas S. J., Kitchin N., Absalon J., Gurtman A., Lockhart S., et al., Safety and efficacy of the bnt162b2 mrna covid-19 vaccine, New England Journal of Medicine.
  49. 49. Baden L. R., El Sahly H. M., Essink B., Kotloff K., Frey S., Novak R., et al., Efficacy and safety of the mrna-1273 sars-cov-2 vaccine, New England Journal of Medicine.
  50. 50. Emary K. R., Golubchik T., Aley P. K., Ariani C. V., Angus B., Bibi S., et al., Efficacy of chadox1 ncov-19 (azd1222) vaccine against sars-cov-2 variant of concern 202012/01 (b. 1.1. 7): an exploratory analysis of a randomised controlled trial, The Lancet 397 (10282) (2021) 1351–1362.
  51. 51. Fiolet T., Kherabi Y., MacDonald C.-J., Ghosn J., Peiffer-Smadja N., Comparing covid-19 vaccines for their characteristics, efficacy and effectiveness against sars-cov-2 and variants of concern: A narrative review, Clinical Microbiology and Infection. pmid:34715347