bgtvs {binGroup}R Documentation

Binomial group testing with variable group sizes

Description

Calculates confidence intervals for binomial group testing if groups of different size are evaluated

Usage

bgtvs(n, s, Y, conf.level = 0.95, alternative = "two.sided",
 maxiter = 100)

Arguments

n numeric vector, specifying the number of groups of the different sizes
s numeric vector, specifying the group sizes, must be of same length as n
Y numeric vector, specifying the number of positive groups among n
conf.level a single numeric value, the confidence level of the interval
alternative a character string, with options "two.sided", "less", "greater"
maxiter maximal number steps in iteration of confidence limits

Details

Hepworth (1996) describes methods for constructing confidence intervals in binomial group testing, if groups of different size are used. Currently, only the exact method (Hepworth, 1996, equation5, Table.5) is implemented. Note, that the algorithm becomes computationally very extensive if the number of different groups becomes larger than 3. See Hepworth (1996) for application of the method

Value

A list containing

conf.int a numerioc vector, the lower and upper limits of the confidence interval
estimate the point estimate

moreover, some of the input arguments.

References

Hepworth, G (1996): Exact confidence intervals for proportions estimated by group testing. Biometrics 52, 1134-1146.

Examples

# Consider a very simple example,
# given in Hepworth (1996), table 5:
# 2 groups each containing 5 units,
# and 3 groups, each containing 2 units
# evaluated for containing at least one 
# positive unit:

# For the result Y=1 positive group among the groups
# of size 5, and Y=2 positive groups among the groups
# of size 2 we yield the interval:

bgtvs(n=c(2,3), s=c(5,2), Y=c(1,2)) 

###############################################

# The examples of Hepworth (1996), table 5 are:

 bgtvs(n=c(2,3), s=c(5,2), Y=c(0,0))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(0,1))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(0,2))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(0,3))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(1,0))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(1,1))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(1,2))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(1,3))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(2,0))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(2,1))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(2,2))

 bgtvs(n=c(2,3), s=c(5,2), Y=c(2,3))


[Package binGroup version 0.3-2 Index]