# NOT RUN {
nsnps <- 100
N0 <- 5000 # number of controls
N1 <- 5000 # number of cases
maf <- runif(nsnps, 0.05, 0.3)
varbeta <- Var.data.cc(f = maf, N = N0 + N1, s = N1/(N0+N1))
bhats = rnorm(nsnps,0,0.2) # log(OR)
est_mu_bhat(bhat = bhats, V = varbeta, N0 = N0, N1 = N1)
# }
Run the code above in your browser using DataLab