Estimate the true effect at the causal variant using estimated effect sizes and their standard errors

Estimate the true effect at the causal variant using estimated effect sizes and their standard errors

est_mu_bhat(bhat, V, N0, N1, p1 = 1e-04, W = 0.2)

Arguments

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

Value

Estimate of the true effect at the causal variant

Examples


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)
#> [1] 13.12425