poolbinom.lrt {binomSamSize} R Documentation

Calculate LRT based confidence interval for binomial proportion for pooled samples

Description

Calculate LRT based confidence interval for the Bernoulli proportion of k*n individuals, which are pooled into n pools each of size k. Observed is the number of positive pools x.

Usage

```poolbinom.lrt(x, k, n, conf.level=0.95, bayes=FALSE, conf.adj=FALSE)
```

Arguments

 `x` Number of positive pools (can be a vector). `k` Pool size (can be a vector). `n` Number of pools (can be a vector). `conf.level` The level of confidence to be used in the confidence interval `bayes` See `binom.bayes` `conf.adj` See `binom.bayes`

Details

Compute LRT based intervals for the binomial response X \sim Bin(n, θ), where θ = 1 - (1-π)^k. As a consequence,

π = g(θ) = 1 - (1-π)^{1/k}

.

One then knows that the borders for π are just transformations of the borders of theta using the above g(θ) function.

For further details about the pooling setup see `poolbinom.logit`.

Value

A data.frame containing the observed proportions and the lower and upper bounds of the confidence interval. The output is similar to the `binom.confint` function of the `binom` package

M. HÃ¶hle

Examples

```binom.lrt(x=0:34,n=34)
poolbinom.lrt(x=0:34,k=1,n=34)
poolbinom.lrt(x=0:34,k=10,n=34)
```

[Package binomSamSize version 0.1-5 Index]