1. Introduction
The group testing idea, introduced by Dorfman in their seminal work [
1], is an efficient approach to the detection of the prevalence of a certain infection in the test samples of a group of individuals. The group testing approach is based on the idea of dividing the individuals into groups, mixing the collected test samples within each group and testing those mixed samples. This way, a negative test result implies that every test sample included in that mixed sample is negative, while a positive test result implies that there is at least one positive sample in the mixed sample: each test acts like a logical OR operation of the test samples included in the mixed sample.
Dorfman’s original algorithm divides the population into disjoint groups and performs group tests by mixing the samples within each group. Subsequently, depending on the test results, positive groups are further tested individually to identify the status of every individual in the population. After Dorfman’s work in [
1], various adaptive (tests performed in multiple steps) and non-adaptive (tests performed in a single step) group testing algorithms have been proposed, capacity of the group testing problem has been studied for a variety of system models and family of algorithms, and extended analyses have been conducted for different regimes for the total number of infections in the population [
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23].
Most of the classical group testing problems have considered either combinatorial or i.i.d. probabilistic settings. Out of
n individuals, in the combinatorial setting, some fixed number of individuals,
k, are assumed to be infected and the problem is to identify the infected set that is uniformly randomly selected from the set of all size
k subsets of
n individuals with the minimum number of tests, while in the i.i.d. probabilistic group testing models, each individual is assumed to be randomly infected, with some fixed probability, independently. In the classical system models, under various assumptions, it is proven that the advantage of group testing over individual testing is considerable mostly in the scenarios where the infection prevalence rate in the population is not high [
5,
15,
16] and the advantage gets diminished as the infection prevalence rate gets higher within the population. More recently, there has been an increasing focus on modified system models, where practical considerations on the system models have significantly improved the performance of the group testing systems. Ref. [
24] considers a probabilistic model where each individual is randomly infected with an unequal probability, independently. Ref. [
25] studies non-identical infection probability distributions and proposes a novel adaptive group testing algorithm. Refs. [
26,
27] consider correlation between infection status of individuals, and Refs. [
28,
29,
30,
31] further model the correlations by considering community structures among the individuals.
Although these extended system models are practical and resemble real-life scenarios, in reality, the testing and infection identification processes are dynamic in nature, especially for contagious diseases. Instead of a static, single-shot identification procedure as in the classical group testing setup, testing and identification need to be done over a long time period while the infection keeps spreading and the infection status of the individuals are dynamically changing. Furthermore, rather than minimizing the required number of tests for perfect (or near perfect with vanishing error) identification of the infection status of a population as in the classical group testing setups, in most practical scenarios, the limit on the testing capacity might be fixed and the objective might be to identify as many infections as possible with the limited testing capacity.
Motivated by the above points, in this paper, we consider dynamic testing algorithms over discrete time for a dynamic infection spread model with fixed, limited testing capacity at each time instance, where a full identification is not possible. In our system, test results are available immediately, and thus, the disease spread is not due to the delay between applying tests and receiving test results, but rather due to the limited testing capacity at each time instance. We follow a dynamic infection spread model, which is inspired by the well-known SIR model where the individuals are divided into three groups: susceptible individuals (S), non-isolated infections (I) and isolated infections (R), i.e., recovered individuals in classical SIR model. We do not assume a community structure in our system. We initialize our system by introducing the initial infections, and after that, at each time instance, infection is spread by infected non-isolated individuals to the susceptible individuals. Meanwhile, at each time instance, after the infection spread phase, the testing phase is performed, where a limited number of T tests are performed to detect a number of infections in the system. In our system, the objective is not to minimize the number of required tests to identify everyone at each time instance, but to control the infection spread either as soon as possible or with minimum number of people that got infected throughout the process, by using the given, limited, testing capacity T at each time instance.
In this paper, we analyze the average-case performance of our system, i.e., the expected values of the number of susceptible individuals, non-isolated and isolated infections over time, which are random processes. For symmetric and converging algorithms, we state a general analytical result for the expected number of susceptible individuals in the system when the infection is brought under control, which is the time when there is no non-isolated infection left in the system. We present two dynamic algorithms: dynamic individual testing and dynamic Dorfman type group testing algorithm. We provide weak versions of these two algorithms and use our general result to obtain a lower bound on the expected number of susceptible individuals when the infection is under control. Finally, we run simulations to get numerical results of our proposed algorithms for different sets of parameters. The paper is organized as follows: in
Section 2 we discuss the papers that are closely related to this paper, in
Section 3 we present our system model, in
Section 4 we present our proposed dynamic algorithms and our theoretical analysis, in
Section 5 we present our simulation results, and finally in
Section 6 we present the conclusions.
2. Related Work
While the standard group testing models consider the one-shot infection status identification problem and do not consider dynamic infection status changes in the system over time, following works are related to our proposed model closer than the classical static group testing works. Ref. [
32] considers a limited identification scenario in a classical setting, addressing the limited test capacity with partial identification cases for real-life scenarios. In a continuous setting and without focusing on explicit testing algorithms, but approaching the problem as a control system with controlled variable testing rate, Ref. [
33] considers a dynamic testing infection spread model based on the well-known SIR model [
34]. However, the first works that focus on the testing of individuals in a dynamic setup with an emphasis on group testing are [
35,
36], where the authors consider a delay in test results and minimize the number of required tests to identify the infection status of the population at each time instance, stating an equivalence to the classical static group testing problem. Despite the fact that the tests are performed to identify the infection status of everyone in the population, new infections are introduced into the system, due to the assumed delay in test results, i.e., during the time between performing the tests and getting the test results. In a more recent work, Ref. [
37] considers a scenario where the tests are done in a dynamic manner, with a focus on two-step Dorfman testing with delay and disease spread between the two steps, in a discrete SIR model with a community structure; the authors in [
37] consider quarantining the possible infections, and analyze the trade-off between quarantine and test costs.
In this paper, we consider a dynamic disease spread model which is motivated by the discrete-time SIR model, similar to the disease spread model presented in [
35]. Distinct from the works that we discuss in this section, we propose a system model where the testing capacity is limited at each time instance and only partial identification can be done at each time instance, while disease keeps spreading over time.
3. System Model
We consider a population of n individuals whose infection statuses change over time. The time dimension t is discrete in our system, i.e., . Similar to the classical discrete SIR model, the population consists of three distinct subgroups: susceptible individuals who are not infected but can get infected by infected individuals (S), infected individuals who can infect the susceptible individuals (I), and isolated individuals who were infected, have been detected via performed tests and isolated indefinitely (R) (These are called recovered (R) individuals in the SIR model; we call them isolated individuals. As they are isolated indefinitely, they are recovered eventually). Let denote the infection status of individual i at time t, where 1 represents being infected, 0 represents not being infected and 2 represents being isolated. At the beginning (), we introduce the initial infections in the system, independently with probability p, where is a Bernoulli random variable with parameter p. Random variables are mutually independent for . Let denote the number of susceptible individuals at time t, denote the number of non-isolated infected individuals at time t and denote the number of isolated individuals at time t. Starting from , each time instance consists of two phases: infection spread phase and testing phase, in the respective order.
Infection Spread Phase: Infected individuals spread the infection to the susceptible members of the population. At each time instance, starting from
, the infection spreads independently across the individuals: Each infected individual can infect each susceptible individual with probability
q, independent across both infected individuals and susceptible individuals. Isolated individuals cannot infect others and their infection status cannot change after they are isolated. Thus, probability of the event that individual
i gets infected by another individual
j at time
is equal to
Testing Phase: At each time instance starting from
,
T tests can be performed to the individuals. Note that the testing capacity
T is a given parameter and thus, in contrast to the classical group testing systems, we do not seek to minimize the number of performed tests for full identification of the infection status of the population but aim to efficiently perform
T tests at each time instance to identify and isolate as many infections as possible to control the infection spread. Here, performed tests can be group tests, and we define the
binary test matrices,
, which specify the pooling scheme for the tests at each time
t. For each time instance
, we have the test result vectors
, which are equal to
where
denotes the
ith test result at time
t,
denotes the
ith row,
jth column of the test matrix
.
Note that, since the previous test matrices and test results are available while designing these test matrices, can depend on the previous test results for . We assume that when tests are performed at some time instance , the test results will be available before the infection spread phase at time . Thus, after the test results are available, detected infections are isolated immediately, i.e., if the ith individual is detected to be infected during the testing phase at time , then . Recall that, after an infected individual is isolated at some time , they cannot infect others at times greater than and their infection status cannot change, i.e., for .
Testing Policy: A testing policy is an algorithm that specifies how to allocate the given testing capacity T for each time instance, until the infection is under control. We define to be the time when for all individuals for the first time and we say that the infection is under control at . Note that, after , infection status of the individuals cannot change and the steady state is achieved: They are either isolated () or non-infected (). Since we do not consider re-entries of recoveries to the population, the infection spread is under control when all infections in the system are isolated. Otherwise, infection may keep spreading to the susceptible individuals by the non-detected infections.
Performance Metrics: The main objective is to bring the infection spread under control by detecting and isolating each infected individual by performing at most T tests at each time instance. Note that, meanwhile, infection keeps spreading and thus, detecting the infection status of an individual to be negative does not imply that they are identified for the rest of the process; they can get infected in the later time instances. As defined, is the time that the infection is under control and when the system has reached that state, further testing of the individuals is unnecessary. Therefore, there are two metrics to measure the performance of a testing policy : The time when the infection is brought under control and the total number of isolated individuals when the infection is under control, while comparing the performances of the testing policies, earlier infection control time and less number of total infections at the time of infection control are favored. Proposed algorithms may not simultaneously improve both metrics: One policy may bring the infection spread under control fast (i.e., low ) but may result in high number of total infections (i.e., high ) while another policy may bring the infection spread under control slowly but with lower number of total infections.
4. Proposed Algorithms and Analysis
In this section, we propose two algorithms and analyze their performances. The first algorithm does not utilize the group testing approach and it is based on the idea of dynamically and individually testing the population. The second algorithm consists of a group testing approach at each time instance, similar to the original idea of Dorfman [
1] in a dynamic setting. Before stating these two algorithms and further analyzing their performances individually, we first state general results.
Symmetric and Converging Dynamic Testing Algorithms: In our analysis, we focus on symmetric and converging dynamic testing algorithms, which satisfy the symmetry criterion,
and convergence criterion,
Furthermore, we assume that the probability of an individual not being identified in the tests at time t, denoted by , only depends on the testing capacity T, , and . Note that, for all time instances t.
Infection Spread Probability: We consider for the infection spread probability q. This is a practical assumption since q is the probability of the event of infection spread that is realized independently for every element of the set product of the infected individuals and susceptible individuals, at each time instance.
We analyze the long term behavior of the system in the average-case, i.e., we focus on the terms , and when t is large enough.
Lemma 1. When a symmetric and converging dynamic testing algorithm is implemented,and thus, the system approaches to steady state, in the average-case. Proof. Note that all three system functions
,
and
can be written as the summation of
n indicator functions
which results in
, due to converging algorithm assumption (
4), which is equal to
. □
Note that, when the system reaches a state where , then there will not be further change of the infection status of the individuals, i.e., the infection will be under control. The following lemma is useful for the justification of the average-case analysis of our system.
Lemma 2. When a symmetric and converging dynamic testing algorithm is implemented, we have for arbitrarily small, constant, .
Proof. Since
for all
, we can apply Markov’s inequality,
where (
10) follows from the fact that
for all
, and (
11) follows from the result of Lemma 1. □
The focus of our analysis is to give a lower bound for the number of susceptible individuals (who have never got infected throughout the process) when the infection is brought under control, in the average-case. In order to analyze the long term behavior of , we have to analyze the long term behavior of . A direct calculation of this probability is not analytically tractable, however, by conditioning on , we give a recursive asymptotic calculation. Before stating the recursive relation, we first prove a lemma that will be useful.
Lemma 3. For and for all , we have Proof. For the proof, we use the covariance inequality, i.e.,
which is a direct application of the Cauchy-Schwarz inequality, applied to the random variables
and
. Using the covariance inequality, we have
where (
14) follows from the fact that the random variable
is bounded above by 1 and below by 0, and (
16) follows from the linear approximation of the function
for small
and
that is bounded above by
n. □
Lemma 4. When a symmetric and converging dynamic testing algorithm is implemented, we have Proof. Conditioned on
, we have the following recursive relation for
where (
21) follows from Lemma 3, (
23) follows from the linear approximation of the function
, and (
25) follows from the symmetry criterion of the implemented algorithm. Recursively using the result in (
25) and the initial value
yields the desired result. □
In order to complete our analysis and give a lower bound for the expected number of susceptible individuals when the infection is under control, we further need to focus on . Similar to the case of , a direct calculation is not analytically tractable, however, we have a recursive relation when conditioned on , and we obtain the following lemma.
Lemma 5. When a symmetric and converging dynamic testing algorithm is implemented and and are arbitrarily small for all , we havewhere the conditional probability of an individual not being identified in the tests at time t given is denoted by . See
Appendix A for the proof of Lemma 5. Combining the results of Lemmas 4 and 5, we have the following result.
Theorem 1. When a symmetric and converging dynamic testing algorithm is implemented and vanishing covariance constraints in Lemma 5 are satisfied for all , we have Proof. Expressing
in terms of the corresponding indicator random variables and using the symmetry criterion and results of Lemmas 4 and 5 yields
which is the desired result. □
Our main result Theorem 1 is a general result and holds for the symmetric and converging dynamic testing algorithms as long as they satisfy the vanishing covariance conditions that we state in Lemma 5. In the remainder of this section, we propose and describe two dynamic testing algorithms and analyze their performance.
4.1. Dynamic Individual Testing Algorithm
In dynamic individual testing algorithm, we do not utilize the group testing approach, and uniformly randomly select T individuals to individually test at each time instance , from the non-isolated individuals.
In order to analyze the performance of our dynamic individual testing algorithm, we use the general result of Theorem 1. First, we show that the dynamic individual testing algorithm satisfies the symmetry and convergence criteria in (
3) and (
4).
Since the process of selection of individuals to be tested is repeated at each time instance with uniformly random selections, as well as the infection spread process, dynamic individual testing algorithm is symmetric. We show that dynamic individual testing algorithm also satisfies the convergence criterion (
4) in the following lemma. For the range of the testing capacity
T, we focus on the case of
, since when
, at one time instance, everyone can be tested individually and the infection will be under control trivially.
Lemma 6. For constant T and n, dynamic individual testing algorithm satisfies the convergence criterion Proof. First, the probability that an infected individual is detected at a time instance
t, denoted by
is
where
denotes the probability of the conditional event that an infected individual is not detected at the time instance
t given
. Now, since the conditional events of detection given that the individual is infected are independent across time due to the uniform random selection of tested individuals at each time instance, and the fact that
since the right hand side of (
37) grows to infinity, from the second Borel-Cantelli lemma, the conditional detection event occurs infinitely often, i.e., let
denote the event that the individual
i is identified at time
t, then
which yields the desired result of
. □
Next, we consider a weak version of our algorithm, where at each time instance, during the testing phase, instead of selecting
T individuals to test from
non-isolated individuals, we select
T individuals from
n individuals, including the isolated ones, whose test results will be negative. For the weak version of the dynamic individual testing algorithm, we have
which is the identification probability of an individual at time
t. Moreover, since it is an upper bound for the identification probability of an individual for the original dynamic individual testing algorithm, we have
Since the weak dynamic individual testing algorithm is a symmetric and converging algorithm (note that the result of Lemma 6 still holds) and due to the fact that is constant in the weak dynamic individual testing algorithm, we can directly use the result of Lemma 5, due to the fact that the vanishing covariance criteria are already satisfied. Now, using Theorem 1, we have the following result for the weak dynamic individual testing algorithm.
Theorem 2. When weak dynamic individual testing algorithm is used and , we havewhich is a lower bound for , i.e., the limit of the expected number of susceptible individuals for the dynamic individual testing algorithm. Proof. The weak dynamic individual testing algorithm satisfies the constraints for using Theorem 1. Thus, we can use Theorem 1 directly to derive the long term behavior of the expected number of the susceptible individuals by considering the limit of (
27) for constant
. On the other hand, in the case of
, we have
. □
4.2. Dynamic Dorfman Type Group Testing Algorithm
In dynamic Dorfman type group testing algorithm, we utilize the group testing idea while designing the test matrices at each time instance .
At each time instance, dynamic Dorfman type group testing algorithm uniformly randomly partitions the set of all non-isolated individuals to equal sized disjoint sets (with possible 1 unequal sized set if the total number of non-isolated individuals is not divisible by ). Then, test samples of the individuals are mixed with others in the same group: group tests are performed, and positive and negative groups are determined. Then, among the positive groups, one group (or multiple groups if the sizes of the groups are less than , depending on the system parameters) is uniformly randomly selected to be individually tested. individuals from the selected group are uniformly randomly selected and individually tested; here depending on the parameters, some individuals from the selected group may not be tested, as well as individuals from multiple positive groups may be selected. Detected infections are isolated and at the next time instance, the whole process is repeated with uniform random selections.
Since the partition selection and individuals within group selection are uniformly random at each time instance, dynamic Dorfman type group testing algorithm is symmetric. Similar to
Section 4.1, we proceed by showing that the dynamic Dorfman type group testing algorithm satisfies the convergence criterion in (
4) as well.
Lemma 7. For constant T and n, dynamic Dorfman testing algorithm satisfies the convergence criterion Proof. The probability that an individual is identified at a time instance
t, which is
, satisfies the following
since
individuals are individually tested at each time instance and due to the symmetry of the infection status in the system and the fact that the individuals are selected from a positive group (or from multiple positive groups), the probability of detection for the dynamic Dorfman type group testing algorithm, at each time instance, must be higher than uniformly randomized testing of
individuals. Now, since the events of identification of individuals are independent across time due to the uniform random selection of tested individuals at each time instance, and the fact that
grows to infinity, we conclude that
, from the second Borel-Cantelli lemma as in Lemma 6. □
Similar to the dynamic individual testing case, we focus on a weak version of the dynamic Dorfman type group testing algorithm to provide a lower bound for the expected number of susceptible individuals in the system at the steady state.
In the weak version of the dynamic Dorfman type group testing algorithm, the results from the
group tests are discarded and it is basically equivalent to the uniformly random individual testing of
individuals. Furthermore, the isolated individuals are also included in the testing procedure:
n individuals are divided into groups and then tested at each time instance, rather than only non-isolated individuals as in the original dynamic Dorfman type group testing algorithm. The probability of identification at time
t for the weak dynamic Dorfman type group testing algorithm, given by
, is always less than the original dynamic Dorfman type group testing algorithm, due to the discarded
group tests and included isolated individuals to the tests. Note that the weak dynamic Dorfman type group testing algorithm is also symmetric and satisfies the convergent criterion (
4), since Lemma 7 still holds; the lower bound in (
44) is the detection probability of the weak algorithm. Moreover, since the weak algorithm has constant value for
, it satisfies the vanishing covariance constraints given in the statement of Lemma 5. Using the general result of Theorem 1, we have the following result for dynamic Dorfman type group testing algorithm by following similar steps to those in Theorem 2.
Theorem 3. When weak dynamic Dorfman type group testing algorithm is used and , we havewhich is a lower bound for , i.e., the expected number of susceptible individuals for the dynamic Dorfman type group testing algorithm. Note that, this result of weak dynamic Dorfman type group testing algorithm is a loose lower bound for the performance of the algorithm, which is only significant because it shows that, the weak dynamic Dorfman type group testing algorithm performs in a similar manner with weak dynamic individual testing algorithm, order-wise (T replaced with ), which is a performance lower bound for the dynamic Dorfman type group testing algorithm.
4.3. Comparison of Dynamic Individual and Dorfman Algorithms
To compare the average number of detected infections at a given time instance for the dynamic individual testing and dynamic Dorfman type group testing algorithms, we obtain the following results stated in the following lemmas.
Lemma 8. When there are susceptible and non-isolated infected individuals in a system after the infection spread phase and just before the testing phase at time instance t, and dynamic individual testing algorithm is being used, on average, infections are detected and isolated at time t.
Proof. When
T individuals from
individuals are uniformly randomly selected, we have
where
represents the infection status of the
ith selected individual for testing at the time of testing phase. □
On the other hand, when dynamic Dorfman type group testing algorithm is used, individuals to be individually tested are chosen from a set of individuals of size that is guaranteed to have at least one infected individual, in the case of . When , individuals to be tested individually are chosen from multiple groups, each having at least one infected individual. The following lemma gives an average number of detected and isolated infections at each time instance for the case of , which, in general, holds for practical applications with low testing capacity. Moreover, the following result is also a lower bound for the case of , where individuals to be individually tested are selected from multiple positive groups.
Lemma 9. When there are susceptible and non-isolated infected individuals in a system after the infection spread phase and just before the testing phase at time instance t, with , and dynamic Dorfman type group testing algorithm is being used, if , on average,infections are detected and isolated at time t. If , then, on average, infections are detected and isolated at time t. In the case of , (48) is a lower bound for the average number of detected and isolated individuals at time t. Proof. When
individuals are uniformly randomly selected from a set of individuals that are guaranteed to have at least one infection, with size
, we have
where
represents the infection status of the
ith selected individual for testing, at the time of testing phase and
.
In the case of
,
individuals to be tested individually are chosen from multiple groups where each of them is guaranteed to have at least one infected individual. Therefore, the term in the denominator of the right hand side of (
49), i.e.,
, is replaced by the probability of the event that multiple subsets of size
C having at least one non-isolated infected member, which is a subset of the event that only one subset of individuals of size
C having at least one non-isolated infected member and thus, having lower probability. Therefore, (
49) is also a lower bound for the average number of detected and isolated infections at time instance
t, for the case of
. □
For a given state of the system at the time of the testing phase, i.e.,
and
, as we show in Lemmas 8 and 9, using dynamic Dorfman type group testing algorithm becomes advantageous with respect to the dynamic individual testing algorithm when
and
where
.
In the next section, we present numerical results of our two proposed dynamic algorithms, as well as their weak versions, under various sets of system parameters.
5. Numerical Results
In our numerical results, we implement the algorithms that we proposed: dynamic Dorfman type group testing algorithm, dynamic individual testing algorithm and the weak versions of these algorithms. In all of our simulations, we start with
n individuals with all of them susceptible. Then, at time
, we realize the initial infections in the system uniformly randomly with probability
p. At each time instance that follows, for the infection spread phase, we simulate the random infection spread from the non-isolated infections to the susceptible individuals. For the testing phase, we simulate the random selection of individuals to be tested and perform the tests. Depending on the test results, we simulate the isolation of the detected infections. We repeat these phases at each time instance until time
, and obtain the sample paths of the random processes
,
and
. We iterate this whole process 1000 times to obtain 1000 sample paths of the random processes, and then we calculate the average of the sample paths to obtain the expected values of
,
and
, numerically. In
Figure 1,
Figure 2 and
Figure 3, we plot these expected values of
,
and
for the algorithms that we propose. In our simulations, we also consider the value of the theoretical approximation result that we obtained in Theorem 1. For each sample path, at each time instance, we numerically calculate the values of
for both dynamic individual testing and dynamic Dorfman type group testing algorithms, and then use the expression that we obtained in Theorem 1 to calculate the
approximation curve. We calculate and plot the average of the
approximation curve. For the weak versions of the proposed algorithms, we use the results of Theorems 2 and 3 to directly calculate and plot the steady state approximations of
. Pseudo-code of the simulations that we perform is given below in Algorithm 1.
Algorithm 1 Simulations to Obtain Numerical Results. | |
| ▹ Simulations are repeated to obtain average sample paths. |
while do | |
| |
| |
for individual j in do | |
if random roll in then | |
| |
end if | |
end for | |
| ▹ Initial infections are realized |
while do | |
for individual j in infections do | |
for individual k in susceptibles do | |
if random roll in then | |
| |
end if | |
end for | |
end for | ▹ End of the infection spread phase |
for test in do | |
| |
end for | ▹ Testing is done |
for individual j in [n] do | |
if j has infection and detected then | |
| |
end if | |
end for | ▹ Recovery is done |
▹ are saved | |
end while | |
end while | ▹ Averages of are calculated over iterations i |
In
Figure 1, we present numerical results for the system with the parameters
,
,
and
. Due to the relatively high number of initial infections in the system, we observe that the dynamic individual testing algorithm performs better than the dynamic Dorfman type group testing algorithm in terms of the average steady state
. In the weak versions of the algorithms, we observe that their performance is strictly worse than their respective original algorithms, at each time instance, in terms of the average
, as expected. The difference of the average
curves between the original and weak versions of the dynamic Dorfman type group testing algorithm is higher than the difference of the average
curves between the original and weak versions of the dynamic individual testing algorithm. This is due to the fact that, in the weak dynamic individual testing algorithm, we still utilize
T tests at each time instance but can sample the isolated individuals to test while in the weak dynamic Dorfman type group testing algorithm, we ignore the group tests and only consider
individual tests. However, since the advantage of the group test is not effective for this set of parameters, as we present in
Figure 1, even the weak dynamic Dorfman type group testing algorithm provides a reasonable lower bound for its original version. Finally, we observe that our approximation results in Theorem 1 match with the average
curves in both dynamic Dorfman type group testing and dynamic individual testing algorithms. Similarly, the average
curves that we obtain from the weak versions of the proposed algorithms are also closely approximated by the results that we obtain in Theorems 2 and 3.
In
Figure 2, we run the same simulations as in
Figure 1, for the parameters
,
,
and
. Now, relative to the first set of parameters, the number of initial infections is lower but the infection spread probability is higher. Because of the targeted individual testing to the positive groups in the dynamic Dorfman type group testing algorithm, it outperforms the dynamic individual testing algorithm for this set of parameters, as we present in
Figure 2. Since the advantage of the group testing is more prevalent for this set of parameters, the weak version of the dynamic Dorfman type group testing algorithm results in an average
curve that is a loose lower bound for the average
curve of the original version, while the weak dynamic individual testing algorithm results in a proper lower bound for the original version. Furthermore, for this set of parameters, despite the fact that the Theorem 1 approximation matches the average
curves for both of the original versions of the proposed algorithms, the resulting Theorem 2 and Theorem 3 approximations cannot be used due to the non-convergent exponents in the expressions.
In our third and final set of parameters, we consider lower number of test capacity, T, than the first two sets of parameters, high number of initial infections, p, and high infection spread probability, q. As expected, for this set of parameters, for both of the algorithms, system reaches the steady state when almost everyone in the population gets infected. Due to the high number of infections at each time instance in the system, dynamic individual testing algorithm performs slightly better than the dynamic Dorfman type group testing algorithm, even though it still fails to control the infection spread in an effective manner.
6. Conclusions
In this paper, we considered a dynamic infection spread model over discrete time, inspired by the SIR model, widely used in the modelling of contagious infections in populations. Instead of recovered individuals in the system, we considered isolated infections, where infected individuals can be identified and isolated via testing. In our system model, infection status of the individuals are random processes, rather than random variables such as the infection status of the individuals in the classical group testing problems. In parallel with the dynamic configuration of our system, we considered dynamic group testing algorithms: At each time instance, after the infection is spread by infected individuals to the susceptible individuals randomly, a given limited number of (possibly group) tests are performed to identify and isolate infected individuals. This dynamic infection spread and identification system is more challenging than the classical group testing problem setup, since negative identifications are not finalized and can change over time while only the positive identifications are isolated for the rest of the process. We analyzed the performance of dynamic testing algorithms by providing approximation results for the expected number of susceptible individuals (that have never gotten infected) when the infection is brought under control where all infections are identified and isolated, for symmetric and converging algorithms. Then, we proposed two dynamic algorithms: dynamic individual testing algorithm and dynamic Dorfman type group testing algorithm. We considered the weak versions of these algorithms and used our general result to provide lower bounds on the expected number of susceptible individuals for these two algorithms. We compared the average identification performance of these two algorithms by deriving conditions when one algorithm outperforms the other. In our simulations, we implemented both the original and weak versions of the proposed algorithms and also simulated and compared the theoretical approximation results that we derived, for three different sets of parameters and we demonstrated various possible scenarios. Our work is unique in that the disease spread in our dynamic system is due to limited testing capacity as opposed to delay in obtaining (unlimited) test results in the existing literature.
Our numerical results suggest that the performance of the proposed algorithms strongly depends on the initial state of the system as well as the other system parameters. Each algorithm can perform better than the other one depending on the system parameters. In this sense, system model assumptions and system parameters play critical roles in the algorithm performances, as expected. This observation leads to the following future directions:
Further dynamic group testing algorithms can be proposed and their performance under varying system parameters can be analyzed.
Variations of the proposed dynamic system model can be studied.
Varying performance metrics can be proposed and studied for varying objectives, such as faster disease spread control time, rather than the number of susceptible individuals at the steady state.