hMeanChiSq {ReplicationSuccess} | R Documentation |
The p-value from the harmonic mean chi-squared test is computed based on study-specific z-values.
hMeanChiSq(z, w=rep(1, length(z)), alternative="greater", bound=TRUE)
z |
A vector of z-values. |
w |
A vector of weights. |
alternative |
Either |
bound |
Determines whether p-value that cannot be computed are
reported as "> bound" ( |
The p-value from the harmonic mean chi-squared test
Leonhard Held
Held, L. (2020). The harmonic mean chi-squared test to substantiate scientific findings. Journal of the Royal Statistical Society: Series C (Applied Statistics), 69, 697-708. https://doi.org/10.1111/rssc.12410
## Example from Fisher (1999) as discussed in Held (2020) library(ReplicationSuccess) pvalues <- c(0.0245, 0.1305, 0.00025, 0.2575, 0.128) lower <- c(0.04, 0.21, 0.12, 0.07, 0.41) upper <- c(1.14, 1.54, 0.60, 3.75, 1.27) se <- ci2se(lower, upper, ratio=TRUE) estimate <- ci2estimate(lower, upper, ratio=TRUE) hMeanChiSq(p2z(pvalues, alternative="less"), alternative="less") hMeanChiSq(p2z(pvalues, alternative="less"), alternative="two.sided") hMeanChiSq(p2z(pvalues, alternative="less"), alternative="none") hMeanChiSq(p2z(pvalues, alternative="less"), w=1/se^2, alternative="less") hMeanChiSq(p2z(pvalues, alternative="less"), w=1/se^2, alternative="two.sided") hMeanChiSq(p2z(pvalues, alternative="less"), w=1/se^2, alternative="none")