zj_pp: Simulate posterior probabilities of causality from joint...

In corrcoverage: Correcting the Coverage of Credible Sets from Bayesian Genetic Fine Mapping

Description

Simulate nrep marginal Z-scores from joint Z-scores and convert these to posterior probabilities of causality

Usage

zj_pp(Zj, V, nrep = 1000, W = 0.2, Sigma)

Arguments

Zj	Vector of joint Z-scores (0s except at CV)
V	Variance of the estimated effect size (can be obtained using Var.beta.cc function)
nrep	Number of posterior probability systems to simulate (default 1000)
W	Prior for the standard deviation of the effect size parameter, beta (default 0.2)
Sigma	SNP correlation matrix

Details

Does not include posterior probabilities for null model

Value

Matrix of simulated posterior probabilties, one simulation per row

Author(s)

Anna Hutchinson

Examples

set.seed(1)
nsnps <- 100
Zj <- rep(0, nsnps)
iCV <- 4 # index of CV
mu <- 5 # true effect at CV
Zj[iCV] <- mu

## generate example LD matrix and MAFs
library(mvtnorm)
nsamples = 1000

simx <- function(nsnps, nsamples, S, maf=0.1) {
    mu <- rep(0,nsnps)
    rawvars <- rmvnorm(n=nsamples, mean=mu, sigma=S)
    pvars <- pnorm(rawvars)
    x <- qbinom(1-pvars, 1, maf)
}

S <- (1 - (abs(outer(1:nsnps,1:nsnps,`-`))/nsnps))^4
X <- simx(nsnps,nsamples,S)
LD <- cor2(X)
maf <- colMeans(X)

## generate V (variance of estimated effect sizes)
varbeta <- Var.data.cc(f = maf, N = 5000, s = 0.5)

res <- zj_pp(Zj, V = varbeta, nrep = 5, W = 0.2, Sigma = LD)

res[c(1:5), c(1:5)]

corrcoverage documentation built on Dec. 7, 2019, 1:07 a.m.