save.seed <- .Random.seed; set.seed(1)
X <- runif(20000, min=0.5, max=4.5) #covariate
H <- rbinom(20000,1,0.1) #hypothesis true or false
Z <- rnorm(20000, H*X) #Z-score
.Random.seed <- save.seed
pvalue <- 1-pnorm(Z) #pvalue
ihw_res <- ihw(pvalue, X, .1)
plot(ihw_res)
Run the code above in your browser using DataLab