The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.
The rema (rare event meta-analysis) package implements a
permutation-based approach for binary meta-analyses of 2x2 tables,
founded on conditional logistic regression, that provides more reliable
statistical tests when heterogeneity is observed in rare event data
(Zabriskie et al. 2021). To adjust for the effect of heterogeneity, this
method conditions on the sufficient statistic of a proxy for the
heterogeneity effect as opposed to estimating the heterogeneity
variance. While this results in the model not strictly falling under the
random-effects framework, it is akin to a random-effects approach in
that it assumes differences in variability due to treatment. Further,
this method does not rely on large-sample approximations or continuity
corrections for rare event data.
This method uses the permutational distribution of the test statistic instead of asymptotic approximations for inference. The number of observed events drives the computation complexity for creating this permutational distribution. Accordingly, for this method to be computationally feasible, it should only be applied to meta-analyses with a relatively low number of observed events. To create this permutational distribution, a network algorithm, based on the work of Mehta et al. (1992) and Corcoran et al. (2001), is employed using C++ and integrated into the package.
You can install the released version of rema from CRAN with
install.packages("rema")and the development version from GitHub through the devtools
package with
# install.packages("devtools")
devtools::install_github("ZabStatLab/rema")Here is an example using rema for a rare event
meta-analysis on a small example data set. For more detailed analyses,
please refer to the package vignette.
library(rema)
trt.events <- c(2, 4, 6, 7, 7, 11)
trt.total <- c(39, 44, 107, 103, 110, 154)
ctrl.events <- c(1, 4, 4, 5, 3, 4)
ctrl.total <- c(43, 44, 110, 100, 106, 146)
rema(trt.events, trt.total, ctrl.events, ctrl.total)
#> Call:
#> rema(trt.events = trt.events, trt.total = trt.total, ctrl.events = ctrl.events,
#> ctrl.total = ctrl.total)
#>
#> OR 95%-CI p-value
#> 0.6457 [0.1512; 3.2015] 0.2423
#>
#> Details on meta-analytical method:
#> - Rare event, heterogeneous meta-analysis method
#> - Two-sided p-value returned (mid.p = TRUE)
#> - Conditional Maximum Likelihood Estimate (CMLE) used when computing the odds ratioCorcoran C, Ryan L, Senchaudhuri P, Mehta C, Patel N, Molenberghs G (2001). “An Exact Trend Test for Correlated Binary Data.” Biometrics, 57, 941–948, doi:10.1111/j.0006-341x.2001.00941.x.
Mehta CR, Patel N, Senchaudhuri P (1992). “Exact Stratified Linear Rank Tests for Ordered Categorical and Binary Data.” Journal of Computational and Graphical Statistics, 1(1), 21–40, doi:10.2307/1390598.
Zabriskie BN, Corcoran C, Senchaudhuri P (2021). “A Permutation-Based Approach for Heterogeneous Meta-Analyses of Rare Events.” Statistics in Medicine, 40(25), 5587-5604, doi:10.1002/sim.9142.
These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.