tests/spjotvoll_weights_tests.R

In dobriban/pweight: P-Value Weighting

# # Unit tests for Spjotvoll weights
#
# #source("./R/bayes_weights.R")
# #source("./R/spjotvoll_weights.R")
#
# # Plots  ----------------
# # generate means
# J <- 2000
# mu <- -abs(rnorm(J))
# sigma <- 0.1 * rep(1, J)
# alpha <- 1 / J
#
# # find weights
# w_1 <- bayes_weights(mu, sigma, alpha)
# plot(mu, w_1$w)
# w_2 <- spjotvoll_weights(mu, alpha)
# points(mu, w_2,col="red", pch=2)