metahet.hybrid: Alternative Tests and Measures for Between-Study...

View source: R/metahet.hybrid.R

metahet.hybridR Documentation

Alternative Tests and Measures for Between-Study Inconsistency in Meta-Analysis

Description

Performs several alternative tests (including a hybrid test) and calculates the associated measures for assessing between-study inconsistency in meta-analysis.

Usage

metahet.hybrid(y, s2, data, gam, testinf = TRUE, iter.resam = 200)

Arguments

y

Numeric vector of observed effect estimates, or the name of a column in data.

s2

Numeric vector of within-study variances corresponding to y, or the name of a column in data.

data

Optional data frame containing the variables y and s2.

gam

Numeric vector of power parameters \gamma, e.g., 1:8.

testinf

Logical; if TRUE, includes the Q_{\infty} statistic (maximum standardized deviation).

iter.resam

Integer; number of bootstrap resamples used to estimate expected values and p-values.

Details

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.

Value

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.

Author(s)

Zhiyuan Yu, Xing Xing, Lifeng Lin

References

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")}>

See Also

metahet

Examples

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)

altmeta documentation built on April 30, 2026, 9:08 a.m.