Note: Please see pdf for full abstract with equations.
It is well and widely known that sample pooling could provide an effective and efficient way for fast coronavirus testing among massive asymptomatic individuals. The method of multi-level acceleration for asymptomatic COVID-19 screening has been introduced, and for one and two levels, the optimal group sizes have been obtained. However, there are still multiple challenges. First, it is not clear how to find the optimal group sizes for three or more levels. Second, there is lack of closed-form expressions for the optimal group sizes for two or more levels. Third, it is not clear how to determine the optimal number of levels. And last, it is not known what the maximum achievable speedup is. The motivation of this paper is to address all the above challenges. The optimization of a hierarchical pooling strategy includes its number of levels and the group size of each level. In this paper, based on multi-variable optimization and Taylor approximation, we are able to derive closed-form expressions for the optimal number of levels d* = ln(1/ln(1/q0)) − 1, the optimal group sizes m1* = ed* = 1/(ep0), m2* = ed*−1 = 1/(e2p0), ..., md** = e = 1/(ed*p0), and the maximum possible speedup of a hierarchical pooling strategy of 1/(ep0ln(1/p0)), where p0 is the fraction of infected people. The above speedup is nearly a linear function of the reciprocal of p0, in the sense that it is asymptotically greater than any sub-linear function (1/p0)1−ε of the reciprocal of p0 for any small ε > 0. Using the results in this paper, we can quickly and easily predict the performance of an optimal hierarchical pooling strategy. For instance, if the fraction of infected people is 0.0001, an 8-level hierarchical pooling strategy can achieve speedup of nearly 400.