ratesci (version 0.2-0)

moverci: MOVER confidence intervals for comparisons of independent binomial or Poisson

Description

Confidence intervals applying the MOVER method ("Method of Variance Estimates Recovery", developed from the Newcombe method for binomial RD) across different contrasts (RD, RR, OR) and distributions (binomial, Poisson) using equal-tailed Jeffreys intervals instead of the Wilson score method for the event rates. Also allows more general Beta and Gamma priors for an approximate Bayesian confidence interval incorporating prior beliefs about the group event rates.

Usage

moverci(x1, n1, x2 = NULL, n2 = NULL, a1 = 0.5, b1 = 0.5, a2 = 0.5,
  b2 = 0.5, cc = 0, level = 0.95, distrib = "bin", contrast = "RD",
  type = "jeff", adj = FALSE, ...)

Arguments

x1, x2

Numeric vectors of numbers of events in group 1 & group 2 respectively.

n1, n2

Numeric vectors of sample sizes (for binomial rates) or exposure times (for Poisson rates) in each group.

a1, b1, a2, b2

Numbers defining the Beta(ai,bi) prior distributions for each group (default ai = bi = 0.5 for Jeffreys method). Gamma priors for Poisson rates require only a1, a2.

cc

Number or logical specifying (amount of) continuity correction (default 0).

level

Number specifying confidence level (between 0 and 1, default 0.95).

distrib

Character string indicating distribution assumed for the input data: "bin" = binomial (default), "poi" = Poisson.

contrast

Character string indicating the contrast required: "RD" (default), "RR", or "OR". "p" gives an interval for the single proportion x1/n1.

type

Character string indicating the method used for the intervals for the individual group rates. "jeff" = Jeffreys equal-tailed intervals (default), "exact" = Clopper-Pearson exact intervals (also obtained using type = "jeff" with cc = 0.5), "wilson" = Wilson score intervals (as per Newcombe 1998). NB: "wilson" option is included only for legacy validation against previous published method by Newcombe. It is not recommended, as type="jeff" achieves much better coverage properties.

adj

Logical (default FALSE) indicating whether to apply the boundary adjustment for Jeffreys intervals recommended on p108 of Brown et al. (type = "jeff" only: set to FALSE if using informative priors)

...

Additional arguments.

Value

A matrix containing the confidence interval for the requested contrast

References

Laud PJ. Equal-tailed confidence intervals for comparison of rates. Pharmaceutical Statistics [in press]. Newcombe RG. Interval estimation for the difference between independent proportions: comparison of eleven methods. Statistics in Medicine 1998; Donner A, Zou G. Closed-form confidence intervals for functions of the normal mean and standard deviation. Statistical Methods in Medical Research Fagerland MW, Newcombe RG. Confidence intervals for odds ratio and relative risk based on the inverse hyperbolic sine transformation. Statistics in Medicine 2013; 32(16):2823-2836. Li HQ, Tang ML, Wong WK. Confidence intervals for ratio of two poisson rates using the method of variance estimates recovery. Computational Statistics 2014; 29(3-4):869-889.

Examples

 
  #Binomial RD, MOVER-J method:
  moverci(x1 = 5, n1 = 56, x2 = 0, n2 = 29)
  
  #Binomial RD, Newcombe method:
  moverci(x1 = 5, n1 = 56, x2 = 0, n2 = 29, type = "wilson")