Compare group testing results

Compare group testing results from objects of class "opchar" returned by operatingCharacteristics1 (opChar1) or operatingCharacteristics2 (opChar2).

CompareConfig(object1, object2)

Arguments

object1

An object of class "opChar" containing group testing results.

object2

A second object of class "opChar" containing group testing results.

Value

A data frame with the expected percent reduction in tests (PercentReductionTests) and the expected increase in testing capacity (PercentIncreaseTestCap) when using the second testing configuration rather than the first testing configuration. Positive values for these quantities indicate that the second testing configuration is more efficient than the first.

Details

The CompareConfig function compares group testing results from two objects of class "opChar". The function creates a data frame with these comparisons.

Author

Brianna D. Hitt and Christopher R. Bilder

Examples

config.mat1 <- matrix(data = c(rep(1, 10), rep(1:2, each = 5), 1:10),
                      nrow = 3, ncol = 10, byrow = TRUE)
res1 <- opChar1(algorithm = "D3", p = 0.05, Se = 0.99, Sp = 0.99,
                hier.config = config.mat1)
#> 
#>  Number of minutes running:  0 
#>  
config.mat2 <- matrix(data = c(rep(1, 10), 1:10),
                      nrow = 2, ncol = 10, byrow = TRUE)
res2 <- opChar1(algorithm = "D2", p = 0.05, Se = 0.99, Sp = 0.99,
        hier.config = config.mat2)
#> 
#>  Number of minutes running:  0 
#>  
CompareConfig(res2, res1)
#> Testing configurations compared to individual testing: 
#>   ExpTests ExpTestsPerIndividual PercentReductionTests PercentIncreaseTestCap
#> 1   5.0324                0.5032                 49.68                  98.71
#> 2   4.0416                0.4042                 59.58                 147.43
#> 
#> Percent reduction in tests when using the second testing
#>      configuration rather than the first: 19.69 
#>  
#> Percent increase in testing capacity when using the second testing
#>      configuration rather than the first: 24.52
#>  

config.mat3 <- matrix(data = c(rep(1, 10), rep(1, 5),
                               rep(2, 4), 3, 1:9, NA),
                      nrow = 3, ncol = 10, byrow = TRUE)
Se <- matrix(data = rep(0.95, 6), nrow = 2, ncol = 3,
             dimnames = list(Infection = 1:2, Stage = 1:3))
Sp <- matrix(data = rep(0.99, 6), nrow = 2, ncol = 3,
             dimnames = list(Infection = 1:2, Stage = 1:3))
res3 <- opChar2(algorithm = "D3", p.vec = c(0.95, 0.02, 0.02, 0.01),
                Se = Se, Sp = Sp, hier.config = config.mat3)
#> 
#>  Number of minutes running:  0 
#>  
config.mat4 <- matrix(data = c(rep(1, 12), rep(1, 6), rep(2, 6),
                               rep(1, 4), rep(2, 2), rep(3, 3),
                               rep(4, 3), 1:12),
                    nrow = 4, ncol = 12, byrow = TRUE)
Se <- matrix(data = rep(0.95, 8), nrow = 2, ncol = 4,
             dimnames = list(Infection = 1:2, Stage = 1:4))
Sp <- matrix(data = rep(0.99, 8), nrow = 2, ncol = 4,
             dimnames = list(Infection = 1:2, Stage = 1:4))
res4 <- opChar2(algorithm = "D4", p.vec = c(0.92, 0.05, 0.02, 0.01),
                Se = Se, Sp = Sp, hier.config = config.mat4)
#> 
#>  Number of minutes running:  0 
#>  
CompareConfig(res4, res3)
#> Testing configurations compared to individual testing: 
#>   ExpTests ExpTestsPerIndividual PercentReductionTests PercentIncreaseTestCap
#> 1   6.2334                0.5195                 48.05                  92.51
#> 2   3.9765                0.3977                 60.23                 151.48
#> 
#> Percent reduction in tests when using the second testing
#>      configuration rather than the first: 23.45 
#>  
#> Percent increase in testing capacity when using the second testing
#>      configuration rather than the first: 30.63
#>