QT2: Calculating p-value for bootstrap-adjusted efficient score... In nshi-stat/netiim3: Improved Inference Methods for Network Meta-Analysis

Description

Calculating p-value for adjusted efficient score test using parametric bootstrap.

Usage

 `1` ```QT2(x, x0) ```

Arguments

 `x` A vector of bootstrap samples of efficient score statistic. `x0` The value of the ordinary efficient score statistic.

Details

Please see Noma et al. (2017) for details.

Value

Caluculated p-value for the adjusted test.

References

Noma, H., Nagashima, K., Maruo, K., Gosho, M., Furukawa, T. A. (2017). Bartlett-type corrections and bootstrap adjustments of likelihood-based inference methods for network meta-analysis. ISM Research Memorandum 1205.

Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28``` ```dae <- data.aug.edit(smoking) y <- dae\$y S <- dae\$S beta1e <- 0.80 C <- 0.95 alpha <- 1 - C ml1 <- ML(y, S) a1 <- ml1\$Coefficients[, 1] a2 <- (ml1\$`Between-studies_SD`)^2 a3 <- a2*(ml1\$`Between-studies_COR`) a4 <- c(a1, a2, a3) mu13 <- ml1\$Coefficients[1, 1] ci3 <- ml1\$Coefficients[1, 3:4] R <- qchisq(C, df=1) beta1 <- log(beta1e) PBS2 <- PBS0.ES(y, S, ml0 = a4, mu0 = beta1, B = 400) ES0 <- ES(y, S, ml0 = a4, mu0 = beta1) QT2(PBS2, ES0) # adjusted p-value ```

nshi-stat/netiim3 documentation built on Aug. 1, 2017, 12:11 a.m.