# NOT RUN {
nsnps <- 100
z_scores <- rnorm(nsnps, 0, 3) # simulate a vector of Z-scores
N0 <- 5000 # number of controls
N1 <- 5000 # number of cases
maf <- runif(nsnps, 0.05, 0.5)
est_mu(z = z_scores, f = maf, N0 = N0, N1 = N1)
# }
Run the code above in your browser using DataLab