SRAT without adjustment for additional covariates
srat0(Y, Z, n_ptb = 1000)
Y | outcome variable |
---|---|
Z | covariates for testing |
n_ptb | number of perturbations |
# Data from wilcox.test x <- c(0.80, 0.83, 1.89, 1.04, 1.45, 1.38, 1.91, 1.64, 0.73, 1.46) y <- c(1.15, 0.88, 0.90, 0.74, 1.21) wilcox.test(x, y)#> #> Wilcoxon rank sum test #> #> data: x and y #> W = 35, p-value = 0.2544 #> alternative hypothesis: true location shift is not equal to 0 #># SRAT p-value set.seed(123) Y <- c(x, y) Z <- rep(0:1, c(length(x), length(y))) srat0(Y, Z)$p_value#> [1] 0.2758489