View source: R/metahet.hybrid.R
| metahet.hybrid | R Documentation |
Performs several alternative tests (including a hybrid test) and calculates the associated measures for assessing between-study inconsistency in meta-analysis.
metahet.hybrid(y, s2, data, gam, testinf = TRUE, iter.resam = 200)
y |
Numeric vector of observed effect estimates, or the name of a column in |
s2 |
Numeric vector of within-study variances corresponding to |
data |
Optional data frame containing the variables |
gam |
Numeric vector of power parameters |
testinf |
Logical; if |
iter.resam |
Integer; number of bootstrap resamples used to estimate expected values and p-values. |
The function generalizes the adaptive sum of powered score (aSPU) framework by combining standardized powered deviations across studies. Specifically, for each power parameter gamma, it calculates a statistic that measures the average standardized deviation of individual study effects from the overall mean, with greater gamma values placing more weight on extreme deviations.
For each gamma, the expected value of this statistic is obtained through a parametric bootstrap under the null hypothesis of homogeneity. The observed statistic is then standardized by subtracting its expected value and dividing by the observed total deviation, resulting in the expected gamma (E-gamma) measure.
In addition, a hybrid test is computed by combining the p-values across all gamma values using the minimum-p combination rule. It provides an overall summary of between-study inconsistency that integrates evidence from multiple power levels.
A list named out containing:
E_gamma: Vector of E-gamma statistics for each \gamma, including E_{\infty} and the hybrid E_{Lhyb} term.
p_gamma: Bootstrap p-values corresponding to each \gamma.
p_hyb: Hybrid p-value combining all \gamma values.
gam: Vector of power parameters used in the tests.
table: Combined summary table for reporting.
Zhiyuan Yu, Xing Xing, Lifeng Lin
Pan W, Kim J, Zhang Y, Shen X, Wei P (2014). "A powerful and adaptive association test for rare variants." Genetics, 197(4), 1081–1095. <\Sexpr[results=rd]{tools:::Rd_expr_doi("10.1534/genetics.114.165035")}>
Xu G, Lin L, Wei P, Pan W (2016). "An adaptive two-sample test for high-dimensional means." Biometrika, 103(3), 609–624. <\Sexpr[results=rd]{tools:::Rd_expr_doi("10.1093/biomet/asw029")}>
Yu Z, Xiao M, Xing X, Lin L (2025). "Alternative tests and measures for between-study inconsistency in meta-analysis." BMC Medical Research Methodology, 25, 261. <\Sexpr[results=rd]{tools:::Rd_expr_doi("10.1186/s12874-025-02719-7")}>
metahet
data("dat.hughes")
set.seed(123)
## increase iter.resam to reduce Monte Carlo error due to resampling
metahet.hybrid(y, s2, data = dat.hughes, gam = 1:8,
testinf = TRUE, iter.resam = 200)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.