Next Article in Journal
On the Numerical Treatment of the Temporal Discontinuity Arising from a Time-Varying Point Mass Attachment on a Waveguide
Next Article in Special Issue
Ant-Balanced Multiple Traveling Salesmen: ACO-BmTSP
Previous Article in Journal
Fourier Neural Operator for Fluid Flow in Small-Shape 2D Simulated Porous Media Dataset
Previous Article in Special Issue
Dynamic SAFFRON: Disease Control Over Time via Group Testing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Infection Spread Model Based Group Testing

Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA
*
Author to whom correspondence should be addressed.
Algorithms 2023, 16(1), 25; https://doi.org/10.3390/a16010025
Submission received: 4 December 2022 / Revised: 27 December 2022 / Accepted: 30 December 2022 / Published: 2 January 2023

Abstract

:
Group testing idea is an efficient approach to detect prevalence of an infection in the test samples taken from a group of individuals. It is based on the idea of pooling the test samples and performing tests to the mixed samples. This approach results in possible reduction in the required number of tests to identify infections. Classical group testing works consider static settings where the infection statuses of the individuals do not change throughout the testing process. In our paper, we study a dynamic infection spread model, inspired by the discrete time SIR model, where infections are spread via non-isolated infected individuals, while infection keeps spreading over time, a limited capacity testing is performed at each time instance as well. In contrast to the classical, static group testing problem, the objective in our setup is not to find the minimum number of required tests to identify the infection status of every individual in the population, but to control the infection spread by detecting and isolating the infections over time by using the given, limited number of tests. In order to analyze the performance of the proposed algorithms, we focus on the average-case analysis of the number of individuals that remain non-infected throughout the process of controlling the infection. We propose two dynamic algorithms that both use given limited number of tests to identify and isolate the infections over time, while the infection spreads, while the first algorithm is a dynamic randomized individual testing algorithm, in the second algorithm we employ the group testing approach similar to the original work of Dorfman. By considering weak versions of our algorithms, we obtain lower bounds for the performance of our algorithms. Finally, we implement our algorithms and run simulations to gather numerical results and compare our algorithms and theoretical approximation results under different sets of system parameters.

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.,  t { 0 , 1 , 2 , } . 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 U i ( t ) 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 ( t = 0 ), we introduce the initial infections in the system, independently with probability p, where U i ( 0 ) is a Bernoulli random variable with parameter p. Random variables U i ( 0 ) are mutually independent for i [ n ] . Let α ( t ) denote the number of susceptible individuals at time t, λ ( t ) denote the number of non-isolated infected individuals at time t and γ ( t ) denote the number of isolated individuals at time t. Starting from t = 1 , 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 t = 1 , 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 t 1 is equal to
q P U j ( t 1 ) = 1 , U i ( t 1 ) = 0 f o r i , j [ n ] .
Testing Phase: At each time instance starting from t = 1 , 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 T × n binary test matrices, X ( t ) , which specify the pooling scheme for the tests at each time t. For each time instance t 1 , we have the test result vectors y ( t ) , which are equal to
y i ( t ) = j [ n ] X i j ( t ) 𝟙 { U j ( t ) = 1 } , i [ T ]
where y i ( t ) denotes the ith test result at time t, X i j ( t ) denotes the ith row, jth column of the test matrix X ( t ) .
Note that, since the previous test matrices and test results are available while designing these test matrices, X ( t ) can depend on the previous test results y ( t ) for t < t . We assume that when tests are performed at some time instance t , the test results y ( t ) will be available before the infection spread phase at time t + 1 . 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 t , then U i ( t ) = 2 . Recall that, after an infected individual is isolated at some time t , they cannot infect others at times greater than t and their infection status cannot change, i.e.,  U i ( t ) = 2 for t t .
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 t ¯ to be the time when U i ( t ¯ ) 1 for all individuals i [ n ] for the first time and we say that the infection is under control at t ¯ . Note that, after  t ¯ , infection status of the individuals cannot change and the steady state is achieved: They are either isolated ( U i ( t ) = 2 ) or non-infected ( U i ( t ) = 0 ). 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, t ¯ 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 t ¯ 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 t ¯ and less number of total infections at the time of infection control γ ( t ¯ ) are favored. Proposed algorithms may not simultaneously improve both metrics: One policy may bring the infection spread under control fast (i.e., low t ¯ ) but may result in high number of total infections (i.e., high γ ( t ¯ ) ) 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,
P ( U i ( t ) = k ) = P ( U j ( t ) = k ) , i , j [ n ] k { 0 , 1 , 2 } t 0
and convergence criterion,
lim t P ( U i ( t ) = 1 ) = o ( 1 / n ) , i [ n ]
Furthermore, we assume that the probability of an individual not being identified in the tests at time t, denoted by p ( t ) , only depends on the testing capacity T, α ( t ) , λ ( t ) and γ ( t ) . Note that, α ( t ) + λ ( t ) + γ ( t ) = n for all time instances t.
Infection Spread Probability: We consider q = o ( 1 / n ) 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 E [ α ( t ) ] , E [ λ ( t ) ] and E [ γ ( t ) ] when t is large enough.
Lemma 1.
When a symmetric and converging dynamic testing algorithm is implemented,
lim t E [ λ ( t ) ] = o ( 1 )
and thus, the system approaches to steady state, in the average-case.
Proof. 
Note that all three system functions α ( t ) , λ ( t ) and γ ( t ) can be written as the summation of n indicator functions
E [ λ ( t ) ] = E i = 1 n 𝟙 { U i ( t ) = 1 }
= i = 1 n E 𝟙 { U i ( t ) = 1 }
= i = 1 n P ( U i ( t ) = 1 )
= n P ( U i ( t ) = 1 )
which results in lim t E [ λ ( t ) ] = n o ( 1 / n ) , due to converging algorithm assumption (4), which is equal to o ( 1 ) .    □
Note that, when the system reaches a state where λ ( t ) = 0 , 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 lim t P ( λ ( t ) > ϵ ) = o ( 1 ) for arbitrarily small, constant, ϵ R .
Proof. 
Since λ ( t ) 0 for all t 0 , we can apply Markov’s inequality,
lim t P ( λ ( t ) > ϵ ) lim t E [ λ ( t ) ] ϵ
= o ( 1 )
where (10) follows from the fact that P ( λ ( t ) > ϵ ) E [ λ ( t ) ] ϵ for all t 0 , 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 E [ α ( t ) ] , we have to analyze the long term behavior of P ( U i ( t ) = 0 ) . A direct calculation of this probability is not analytically tractable, however, by conditioning on λ ( t 1 ) , we give a recursive asymptotic calculation. Before stating the recursive relation, we first prove a lemma that will be useful.
Lemma 3.
For q = o ( 1 / n ) and for all t 0 , we have
c o v P ( U i ( t ) = 0 | λ ( t ) ) , ( 1 q ) λ ( t ) 0
Proof. 
For the proof, we use the covariance inequality, i.e.,  | c o v ( X , Y ) | v a r ( X ) v a r ( Y ) which is a direct application of the Cauchy-Schwarz inequality, applied to the random variables X E [ X ] and Y E [ Y ] . Using the covariance inequality, we have
| c o v P ( U i ( t ) = 0 | λ ( t ) ) , ( 1 q ) λ ( t ) |
v a r ( P ( U i ( t ) = 0 | λ ( t ) ) ) v a r ( ( 1 q ) λ ( t ) )
v a r ( ( 1 q ) λ ( t ) )
= E [ ( 1 q ) 2 λ ( t ) ] ( E [ ( 1 q ) λ ( t ) ] ) 2
( 1 q ) E [ 2 λ ( t ) ] ( ( 1 q ) E [ λ ( t ) ] ) 2
= 0
where (14) follows from the fact that the random variable P ( U i ( t ) = 0 | λ ( t ) ) is bounded above by 1 and below by 0, and (16) follows from the linear approximation of the function ( 1 q ) x for small q = o ( 1 / n ) and λ ( t ) that is bounded above by n.    □
Lemma 4.
When a symmetric and converging dynamic testing algorithm is implemented, we have
P ( U i ( t ) = 0 ) ( 1 p ) ( 1 q ) n j = 0 t 1 P ( U 1 ( j ) = 1 )
Proof. 
Conditioned on λ ( t 1 ) , we have the following recursive relation for P ( U i ( t ) = 0 )
P ( U i ( t ) = 0 ) = E [ P ( U i ( t ) = 0 | λ ( t 1 ) ) ]
= E [ P ( U i ( t 1 ) = 0 | λ ( t 1 ) ) ( 1 q ) λ ( t 1 ) ]
E [ P ( U i ( t 1 ) = 0 | λ ( t 1 ) ) ] E [ ( 1 q ) λ ( t 1 ) ]
= P ( U i ( t 1 ) = 0 ) E [ ( 1 q ) λ ( t 1 ) ]
P ( U i ( t 1 ) = 0 ) ( 1 q ) E [ λ ( t 1 ) ]
= P ( U i ( t 1 ) = 0 ) ( 1 q ) j = 1 n P ( U j ( t 1 ) = 1 )
= P ( U i ( t 1 ) = 0 ) ( 1 q ) n P ( U 1 ( t 1 ) = 1 )
where (21) follows from Lemma 3, (23) follows from the linear approximation of the function ( 1 q ) x 1 q x , and (25) follows from the symmetry criterion of the implemented algorithm. Recursively using the result in (25) and the initial value P ( U i ( 0 ) = 0 ) = ( 1 p ) 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 P ( U i ( t ) = 1 ) . Similar to the case of P ( U i ( t ) = 0 ) , a direct calculation is not analytically tractable, however, we have a recursive relation when conditioned on λ ( t 1 ) , and we obtain the following lemma.
Lemma 5.
When a symmetric and converging dynamic testing algorithm is implemented and c o v P ( U i ( t ) = 0 | λ ( t ) ) , p λ ( t ) ( t + 1 ) and c o v P ( U i ( t ) = 1 | λ ( t ) ) , p λ ( t ) ( t + 1 ) are arbitrarily small for all t 0 , we have
P ( U i ( t ) = 1 ) p ( 1 + n q ( 1 p ) ) t j = 1 t p ( j )
where the conditional probability of an individual not being identified in the tests at time t given λ ( t 1 ) is denoted by p λ ( t 1 ) .
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 t 0 , we have
E [ α ( t ) ] n ( 1 p ) ( 1 q ) n p i = 0 t 1 ( 1 + n q ( 1 p ) ) i j = 1 i p ( j )
Proof. 
Expressing α ( t ) in terms of the corresponding indicator random variables and using the symmetry criterion and results of Lemmas 4 and 5 yields
E [ α ( t ) ] = E i = 1 n 𝟙 { U i ( t ) = 0 }
= i = 1 n E [ 𝟙 { U i ( t ) = 0 } ]
= i = 1 n P ( U i ( t ) = 0 )
= n P ( U i ( t ) = 0 )
n ( 1 p ) ( 1 q ) n p i = 0 t 1 ( 1 + n q ( 1 p ) ) i j = 1 i p ( j )
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 t 1 , 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 T < n , since when T n , 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
lim t P ( U i ( t ) = 1 ) = 0 , i [ n ]
Proof. 
First, the probability that an infected individual is detected at a time instance t, denoted by 1 p ( t ) is
1 p ( t ) = E [ 1 p γ ( t 1 ) ( t ) ]
= E T n γ ( t 1 )
T n
where p γ ( t 1 ) ( t ) denotes the probability of the conditional event that an infected individual is not detected at the time instance t given γ ( t 1 ) . 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
i = 1 ( 1 p ( i ) ) i = 1 T n
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 D t denote the event that the individual i is identified at time t, then
P ( lim sup t D t ) = 1
which yields the desired result of lim t P ( U i ( t ) = 1 ) = 0 .    □
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 n γ ( t ) 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
1 p ( t ) = T n , t > 0
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
lim t E [ α o r i g ( t ) ] lim t E [ α w e a k ( t ) ]
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 p ( t ) 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 ( 1 T n ) ( 1 + n q ( 1 p ) ) < 1 , we have
lim t E [ α w e a k ( t ) ] n ( 1 p ) ( 1 q ) n p 1 ( 1 T n ) ( 1 + n q ( 1 p ) )
which is a lower bound for lim t E [ α o r i g ( t ) ] , 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 p ( t ) = 1 T / n . On the other hand, in the case of ( 1 T n ) ( 1 + n q ( 1 p ) ) 1 , we have lim t E [ α w e a k ( t ) ] 0 .    □

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 t 1 .
At each time instance, dynamic Dorfman type group testing algorithm uniformly randomly partitions the set of all non-isolated individuals to equal sized T / 2 disjoint sets (with possible 1 unequal sized set if the total number of non-isolated individuals is not divisible by T / 2 ). Then, test samples of the individuals are mixed with others in the same group: T / 2 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 T / 2 , depending on the system parameters) is uniformly randomly selected to be individually tested. T / 2 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
lim t P ( U i ( t ) = 1 ) = 0 , i [ n ]
Proof. 
The probability that an individual is identified at a time instance t, which is 1 p ( t ) , satisfies the following
1 p ( t ) T 2 n
since T / 2 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 T / 2 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
i = 1 ( 1 p ( t ) ) i = 1 T 2 n
grows to infinity, we conclude that lim t P ( U i ( t ) = 1 ) = 0 , 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 T / 2 group tests are discarded and it is basically equivalent to the uniformly random individual testing of T / 2 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 1 p ( t ) , is always less than the original dynamic Dorfman type group testing algorithm, due to the discarded T / 2 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 p λ ( t 1 ) ( t ) , 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 ( 1 T 2 n ) ( 1 + n q ( 1 p ) ) < 1 , we have
lim t E [ α w e a k ( t ) ] n ( 1 p ) ( 1 q ) n p 1 ( 1 T 2 n ) ( 1 + n q ( 1 p ) )
which is a lower bound for lim t E [ α o r i g ( t ) ] , 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 T / 2 ), 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 α ˜ ( t ) susceptible and λ ˜ ( t ) 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, T λ ˜ ( t ) α ˜ ( t ) + λ ˜ ( t ) infections are detected and isolated at time t.
Proof. 
When T individuals from α ˜ ( t ) + λ ˜ ( t ) individuals are uniformly randomly selected, we have
E i = 1 T 𝟙 U ˜ i ( t ) = 1 = T P ( U ˜ i ( t ) = 1 )
= T λ ˜ ( t ) α ˜ ( t ) + λ ˜ ( t )
where U ˜ i ( t ) 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, T / 2 individuals to be individually tested are chosen from a set of individuals of size 2 ( α ˜ ( t ) + λ ˜ ( t ) ) T that is guaranteed to have at least one infected individual, in the case of α ˜ ( t ) + λ ˜ ( t ) T 2 / 4 . When α ˜ ( t ) + λ ˜ ( t ) < T 2 / 4 , T / 2 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 α ˜ ( t ) + λ ˜ ( t ) T 2 / 4 , which, in general, holds for practical applications with low testing capacity. Moreover, the following result is also a lower bound for the case of α ˜ ( t ) + λ ˜ ( t ) < T 2 / 4 , where T / 2 individuals to be individually tested are selected from multiple positive groups.
Lemma 9.
When there are α ˜ ( t ) susceptible and λ ˜ ( t ) non-isolated infected individuals in a system after the infection spread phase and just before the testing phase at time instance t, with  α ˜ ( t ) + λ ˜ ( t ) T 2 / 4 , and dynamic Dorfman type group testing algorithm is being used, if  α ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T , on average,
T λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) 1 α ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T α ˜ ( t ) + λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T 1
infections are detected and isolated at time t. If  α ˜ ( t ) < 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T , then, on average, T λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) infections are detected and isolated at time t. In the case of α ˜ ( t ) + λ ˜ ( t ) < T 2 / 4 , (48) is a lower bound for the average number of detected and isolated individuals at time t.
Proof. 
When T / 2 individuals are uniformly randomly selected from a set of individuals that are guaranteed to have at least one infection, with size 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T , we have
E i = 1 T / 2 𝟙 U ˜ i ( t ) = 1 | i = 1 C 𝟙 U ˜ i ( t ) = 1 1
= E i = 1 T / 2 𝟙 U ˜ i ( t ) = 1 P i = 1 C 𝟙 U ˜ i ( t ) = 1 1
= T λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) 1 α ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T α ˜ ( t ) + λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T 1
where U ˜ i ( t ) represents the infection status of the ith selected individual for testing, at the time of testing phase and C = 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T .
In the case of α ˜ ( t ) + λ ˜ ( t ) < T 2 / 4 , T / 2 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.,  P i = 1 C 𝟙 U ˜ i ( t ) = 1 1 , 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 α ˜ ( t ) + λ ˜ ( t ) < T 2 / 4 .    □
For a given state of the system at the time of the testing phase, i.e.,  α ˜ ( t ) and λ ˜ ( t ) , 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 α ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T and
1 / 2 < α ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T α ˜ ( t ) + λ ˜ ( t ) 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T
= i = 0 C ( α ˜ ( t ) i ) i = 0 C ( α ˜ ( t ) + λ ˜ ( t ) i )
where C = 2 ( α ˜ ( t ) + λ ˜ ( t ) ) / T .
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 t = 0 , 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 t = 500 , and obtain the sample paths of the random processes α ( t ) , λ ( t ) and γ ( t ) . 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 α ( t ) , λ ( t ) and γ ( t ) , numerically. In Figure 1, Figure 2 and Figure 3, we plot these expected values of α ( t ) , λ ( t ) and γ ( t ) 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 p ( t ) 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 α ( t ) approximation curve. We calculate and plot the average of the α ( t ) 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 α ( t ) . Pseudo-code of the simulations that we perform is given below in Algorithm 1.
Algorithm 1 Simulations to Obtain Numerical Results.
i 1 ▹ Simulations are repeated to obtain average sample paths.
while i i m a x do
     t 0
     U 0 n x 1
    for individual j in [ n ]  do
        if random roll in [ 0 , 1 ] p  then
            U j 1
        end if
    end for
     t 1 ▹ Initial infections are realized
    while  t t m a x  do
        for individual j in infections do
           for individual k in susceptibles do
               if random roll in [ 0 , 1 ] q  then
                    U k 1
               end if
           end for
        end for▹ End of the infection spread phase
        for test τ in [ T ]  do
            y τ j [ n ] X τ j 𝟙 { U j = 1 }
        end for▹ Testing is done
        for individual j in [n] do
           if j has infection and detected then
                U j 2
           end if
        end for▹ Recovery is done
         t t + 1      ▹ α ( t ) , λ ( t ) , γ ( t ) are saved
    end while
end while▹ Averages of α ( t ) , λ ( t ) , γ ( t ) are calculated over iterations i
In Figure 1, we present numerical results for the system with the parameters n = 1000 , T = 80 , q = 0.00003 and p = 0.2 . 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 α ( t ) . 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 α ( t ) , as expected. The difference of the average α ( t ) curves between the original and weak versions of the dynamic Dorfman type group testing algorithm is higher than the difference of the average α ( t ) 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 T / 2 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 α ( t ) curves in both dynamic Dorfman type group testing and dynamic individual testing algorithms. Similarly, the average α ( t ) 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 n = 1000 , T = 80 , q = 0.0001 and p = 0.01 . 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 α ( t ) curve that is a loose lower bound for the average α ( t ) 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 α ( t ) 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.

Author Contributions

Conceptualization, B.A. and S.U.; methodology, B.A. and S.U.; software, B.A. and S.U.; validation, B.A. and S.U.; formal analysis, B.A. and S.U.; investigation, B.A. and S.U.; resources, B.A. and S.U.; data curation, B.A. and S.U.; writing—original draft preparation, B.A. and S.U.; writing—review and editing, B.A. and S.U.; visualization, B.A. and S.U.; supervision, S.U.; project administration, S.U.; funding acquisition, S.U. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Lemma A1.
When a symmetric and converging dynamic testing algorithm is implemented and c o v P ( U i ( t ) = 0 | λ ( t ) ) , p λ ( t ) ( t + 1 ) and c o v P ( U i ( t ) = 1 | λ ( t ) ) , p λ ( t ) ( t + 1 ) are arbitrarily small for all t 0 , we have
P ( U i ( t ) = 1 ) p ( ( 1 + n q ( 1 p ) ) ) t j = 1 t p ( j )
where the conditional probability of an individual not being identified in the tests at time t given λ ( t 1 ) is denoted by p λ ( t 1 ) .
Proof. 
Conditioned on λ ( t 1 ) , we have the following recursive relation
P ( U i ( t ) = 1 )
= E [ P ( U i ( t ) = 1 | λ ( t 1 ) ) ] = E [ P ( U i ( t 1 ) = 0 | λ ( t 1 ) ) ( 1 ( 1 q ) λ ( t 1 ) ) p λ ( t 1 ) ( t )
+ P ( U i ( t 1 ) = 1 | λ ( t 1 ) ) p λ ( t 1 ) ( t ) ] E [ 1 ( 1 q ) λ ( t 1 ) ] E [ P ( U i ( t 1 ) = 0 | λ ( t 1 ) ) p λ ( t 1 ) ( t ) ]
+ E [ P ( U i ( t 1 ) = 1 | λ ( t 1 ) ) p λ ( t 1 ) ( t ) ] E [ p λ ( t 1 ) ( t ) ] E [ 1 ( 1 q ) λ ( t 1 ) ] E [ P ( U i ( t 1 ) = 0 | λ ( t 1 ) ) ]
+ E [ p λ ( t 1 ) ( t ) ] E [ P ( U i ( t 1 ) = 1 | λ ( t 1 ) ) ] = p ( t ) ( P ( U i ( t 1 ) = 0 ) ( 1 E [ ( 1 q ) λ ( t 1 ) ] )
+ P ( U i ( t 1 ) = 1 ) ) p ( t ) ( ( 1 p ) ( 1 q ) n j = 0 t 2 P ( U i ( j ) = 1 ) ( 1 ( 1 q ) E [ λ ( t 1 ) ] )
+ P ( U i ( t 1 ) = 1 ) ) = p ( t ) ( ( 1 p ) ( 1 q ) n j = 0 t 2 P ( U i ( j ) = 1 ) ( 1 ( 1 q ) n P ( U i ( t 1 ) = 1 ) )
+ P ( U i ( t 1 ) = 1 ) ) = p ( t ) ( ( 1 p ) ( ( 1 q ) n j = 0 t 2 P ( U i ( j ) = 1 ) ( 1 q ) n j = 0 t 1 P ( U i ( j ) = 1 ) )
+ P ( U i ( t 1 ) = 1 ) ) p ( t ) ( ( 1 p ) ( q n j = 0 t 1 P ( U i ( j ) = 1 ) q n j = 0 t 2 P ( U i ( j ) = 1 ) )
+ P ( U i ( t 1 ) = 1 ) )
= p ( t ) n q ( 1 p ) P ( U i ( t 1 ) = 1 ) + P ( U i ( t 1 ) = 1 )
= p ( t ) 1 + n q ( 1 p ) P ( U i ( t 1 ) = 1 )
where (A4) follows from the arbitrarily small variance of ( 1 q ) λ ( t ) similar to the proof of Lemma 3, (A5) follows from the given vanishing covariance assumptions in the statement of the lemma, (A7) follows from Lemma 4, and (A10) follows from the linear approximation ( 1 q ) x 1 q x . Recursively applying (A12) yields the desired result. □

References

  1. Dorfman, R. The Detection of Defective Members of Large Populations. Ann. Math. Stat. 1943, 14, 436–440. [Google Scholar] [CrossRef]
  2. Sobel, M.; Groll, P.A. Group Testing To Eliminate Efficiently All Defectives in a Binomial Sample. Bell Syst. Tech. J. 1959, 38, 1179–1252. [Google Scholar] [CrossRef]
  3. Hwang, F.K. A Method for Detecting All Defective Members in a Population by Group Testing. J. Am. Stat. Assoc. 1972, 67, 605–608. [Google Scholar] [CrossRef]
  4. Chen, H.B.; Hwang, F.K. Exploring the Missing Link Among d-Separable, d¯-Separable and d-Disjunct Matrices. Discret. Appl. Math. 2007, 155, 662–664. [Google Scholar] [CrossRef] [Green Version]
  5. Ruszinko, M. On the Upper Bound of the Size of the R-Cover-Free Families. J. Comb. Theory Ser. A 1994, 66, 302–310. [Google Scholar] [CrossRef] [Green Version]
  6. Chan, C.L.; Jaggi, S.; Saligrama, V.; Agnihotri, S. Non-Adaptive Group Testing: Explicit Bounds and Novel Algorithms. IEEE Trans. Inf. Theory 2014, 60, 3019–3035. [Google Scholar] [CrossRef] [Green Version]
  7. Baldassini, L.; Johnson, O.; Aldridge, M. The Capacity of Adaptive Group Testing. In Proceedings of the IEEE ISIT, Istanbul, Turkey, 7–12 July 2013. [Google Scholar]
  8. Mazumdar, A. Nonadaptive Group Testing With Random Set of Defectives. IEEE Trans. Inf. Theory 2016, 62, 7522–7531. [Google Scholar] [CrossRef]
  9. Scarlett, J.; Johnson, O. Noisy Non-Adaptive Group Testing: A (Near-)Definite Defectives Approach. IEEE Trans. Inf. Theory 2020, 66, 3775–3797. [Google Scholar] [CrossRef] [Green Version]
  10. Wu, S.; Wei, S.; Wang, Y.; Vaidyanathan, R.; Yuan, J. Partition Information and its Transmission Over Boolean Multi-Access Channels. IEEE Trans. Inf. Theory 2015, 61, 1010–1027. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, C.; Zhao, Q.; Chuah, C. Optimal Nested Test Plan for Combinatorial Quantitative Group Testing. IEEE Trans. Signal Proc. 2018, 66, 992–1006. [Google Scholar] [CrossRef]
  12. Zhu, D.Z.; Hwang, F.K. Combinatorial Group Testing and Its Applications, 2nd ed.; World Scientific: Singapore, 1999. [Google Scholar]
  13. Wolf, J. Born Again Group Testing: Multiaccess Communications. IEEE Trans. Inf. Theory 1985, 31, 185–191. [Google Scholar] [CrossRef]
  14. Aldridge, M.; Baldassini, L.; Johnson, O. Group Testing Algorithms: Bounds and Simulations. IEEE Trans. Inf. Theory 2014, 60, 3671–3687. [Google Scholar] [CrossRef] [Green Version]
  15. Riccio, L.; Colbourn, C.J. Sharper Bounds in Adaptive Group Testing. Taiwan. J. Math. 2000, 4, 669–673. [Google Scholar] [CrossRef]
  16. Aldridge, M.; Johnson, O.; Scarlett, J. Group Testing: An Information Theory Perspective; Now Foundations and Trends: Hanover, MA, USA, 2019. [Google Scholar]
  17. Cai, S.; Jahangoshahi, M.; Bakshi, M.; Jaggi, S. Efficient Algorithms for Noisy Group Testing. IEEE Trans. Inf. Theory 2017, 63, 2113–2136. [Google Scholar] [CrossRef]
  18. Karimi, E.; Kazemi, F.; Heidarzadeh, A.; Narayanan, K.R.; Sprintson, A. Non-adaptive Quantitative Group Testing Using Irregular Sparse Graph Codes. In Proceedings of the Allerton Conference, Monticello, IL, USA, 24–27 September 2019. [Google Scholar]
  19. Inan, H.A.; Kairouz, P.; Wootters, M.; Ozgur, A. On the Optimality of the Kautz-Singleton Construction in Probabilistic Group Testing. In Proceedings of the Allerton Conference, Monticello, IL, USA, 2–5 October 2018. [Google Scholar]
  20. Johnson, O.; Aldridge, M.; Scarlett, J. Performance of Group Testing Algorithms With Near-Constant Tests Per Item. IEEE Trans. Inf. Theory 2019, 65, 707–723. [Google Scholar] [CrossRef] [Green Version]
  21. Allemann, A. An Efficient Algorithm for Combinatorial Group Testing. In Information Theory, Combinatorics, and Search Theory: In Memory of Rudolf Ahlswede; Springer: Berlin/Heidelberg, Germany, 2013; pp. 569–596. [Google Scholar]
  22. Scarlett, J. Noisy Adaptive Group Testing: Bounds and Algorithms. IEEE Trans. Inf. Theory 2019, 65, 3646–3661. [Google Scholar] [CrossRef] [Green Version]
  23. Atia, G.K.; Saligrama, V. Boolean Compressed Sensing and Noisy Group Testing. IEEE Trans. Inf. Theory 2012, 58, 1880–1901. [Google Scholar] [CrossRef] [Green Version]
  24. Li, T.; Chan, C.L.; Huang, W.; Kaced, T.; Jaggi, S. Group Testing with Prior Statistics. In Proceedings of the IEEE ISIT, Honolulu, HI, USA, 29 June–4 July 2014. [Google Scholar]
  25. Doger, M.; Ulukus, S. Group Testing with Non-identical Infection Probabilities. In Proceedings of the 2021 XVII International Symposium “Problems of Redundancy in Information and Control Systems” (REDUNDANCY), Moscow, Russian, 25–29 October 2021. [Google Scholar]
  26. Lendle, S.D.; Hudgens, M.G.; Qaqish, B.F. Group Testing for Case Identification with Correlated Responses. Biometrics 2012, 68, 532–540. [Google Scholar] [CrossRef]
  27. Lin, Y.J.; Yu, C.H.; Liu, T.H.; Chang, C.S.; Chen, W.T. Positively Correlated Samples Save Pooled Testing Costs. IEEE Trans. Netw. Sci. Eng. 2021, 8, 2170–2182. [Google Scholar] [CrossRef]
  28. Nikolopoulos, P.; Srinivasavaradhan, S.R.; Guo, T.; Fragouli, C.; Diggavi, S. Group Testing for Connected Communities. In Proceedings of the AISTATS, Online, 13–15 April 2021. [Google Scholar]
  29. Nikolopoulos, P.; Srinivasavaradhan, S.R.; Guo, T.; Fragouli, C.; Diggavi, S. Group Testing for Overlapping Communities. In Proceedings of the IEEE ICC, Montreal, QC, Canada, 14–23 June 2021. [Google Scholar]
  30. Ahn, S.; Chen, W.N.; Ozgur, A. Adaptive Group Testing on Networks with Community Structure. In Proceedings of the IEEE ISIT, Melbourne, Australia, 12–20 July 2021. [Google Scholar]
  31. Arasli, B.; Ulukus, S. Graph and Cluster Formation Based Group Testing. In Proceedings of the IEEE ISIT, Melbourne, Australia, 12–20 July 2021. [Google Scholar]
  32. Aldridge, M. Pooled Testing to Isolate Infected Individuals. In Proceedings of the CISS, Baltimore, MD, USA, 24–26 March 2021. [Google Scholar]
  33. Acemoglu, D.; Fallah, A.; Giometto, A.; Huttenlocher, D.; Ozdaglar, A.; Parise, F.; Pattathil, S. Optimal Adaptive Testing for Epidemic Control: Combining Molecular and Serology Tests. arXiv 2021, arXiv:2101.00773. [Google Scholar]
  34. Wacker, B.; Schluter, J. Time-Continuous and Time-Discrete SIR Models Revisited: Theory and Applications. Adv. Differ. Equ. 2020, 2020, 556. [Google Scholar] [CrossRef] [PubMed]
  35. Srinivasavaradhan, S.R.; Nikolopoulos, P.; Fragouli, C.; Diggavi, S. Dynamic Group Testing to Control and Monitor Disease Progression in a Population. arXiv 2021, arXiv:2106.10765. [Google Scholar]
  36. Srinivasavaradhan, S.R.; Nikolopoulos, P.; Fragouli, C.; Diggavi, S. An Entropy Reduction Approach to Continual Testing. In Proceedings of the IEEE ISIT, Melbourne, Australia, 12–20 July 2021. [Google Scholar]
  37. Doger, M.; Ulukus, S. Dynamical Dorfman Testing with Quarantine. In Proceedings of the CISS, Baltimore, MD, USA, 24–26 March 2022. [Google Scholar]
Figure 1. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorems 1–3 when n = 1000 , T = 80 , q = 0.00003 , p = 0.2 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm, (c) weak dynamic Dorfman type group testing algorithm, (d) weak dynamic individual testing algorithm.
Figure 1. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorems 1–3 when n = 1000 , T = 80 , q = 0.00003 , p = 0.2 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm, (c) weak dynamic Dorfman type group testing algorithm, (d) weak dynamic individual testing algorithm.
Algorithms 16 00025 g001aAlgorithms 16 00025 g001b
Figure 2. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorem 1 when n = 1000 , T = 80 , q = 0.0001 , p = 0.01 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm, (c) weak dynamic Dorfman type group testing algorithm, (d) weak dynamic individual testing algorithm.
Figure 2. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorem 1 when n = 1000 , T = 80 , q = 0.0001 , p = 0.01 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm, (c) weak dynamic Dorfman type group testing algorithm, (d) weak dynamic individual testing algorithm.
Algorithms 16 00025 g002
Figure 3. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorem 1 when n = 1000 , T = 40 , q = 0.0002 , p = 0.2 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm.
Figure 3. Average values of the random processes α ( t ) , λ ( t ) and γ ( t ) , with obtained theoretical approximations given in Theorem 1 when n = 1000 , T = 40 , q = 0.0002 , p = 0.2 , for (a) dynamic Dorfman type group testing algorithm, (b) dynamic individual testing algorithm.
Algorithms 16 00025 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Arasli, B.; Ulukus, S. Dynamic Infection Spread Model Based Group Testing. Algorithms 2023, 16, 25. https://doi.org/10.3390/a16010025

AMA Style

Arasli B, Ulukus S. Dynamic Infection Spread Model Based Group Testing. Algorithms. 2023; 16(1):25. https://doi.org/10.3390/a16010025

Chicago/Turabian Style

Arasli, Batuhan, and Sennur Ulukus. 2023. "Dynamic Infection Spread Model Based Group Testing" Algorithms 16, no. 1: 25. https://doi.org/10.3390/a16010025

APA Style

Arasli, B., & Ulukus, S. (2023). Dynamic Infection Spread Model Based Group Testing. Algorithms, 16(1), 25. https://doi.org/10.3390/a16010025

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop