Description Usage Arguments Value Author(s) Examples
View source: R/extra_functions.R
Estimate the true effect at the causal variant using estimated effect sizes and their standard errors
1 | est_mu_bhat(bhat, V, N0, N1, p1 = 1e-04, W = 0.2)
|
bhat |
Vector of estimated effect sizes |
V |
Prior variance for estimated effect sizes |
N0 |
Number of controls |
N1 |
Number of cases |
p1 |
Prior probability a SNP is associated with the trait, default 1e-4 |
W |
Prior for the standard deviation of the effect size parameter, beta |
Estimate of the true effect at the causal variant
Anna Hutchinson
1 2 3 4 5 6 7 8 9 10 11 | 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.