wakefield_pp: compute posterior probabilities using Wakefield's approximate...

Description Usage Arguments Value Author(s)

View source: R/rapfunc.R

Description

This function was adapted from its namesake in cupcake package (github.com/ollyburren/cupcake/) to no longer require allele frequencies.

Usage

1
wakefield_pp(beta, se, pi_i = 1e-04, sd.prior = 0.2)

Arguments

beta

a vector of effect sizes (β) from a quantitative trait GWAS

se

vector of standard errors of effect sizes (β)

pi_i

a scalar representing the prior probability (DEFAULT 1 \times 10^{-4})

sd.prior

a scalar representing our prior expectation of β (DEFAULT 0.2). The method assumes a normal prior on the population log relative risk centred at 0 and the DEFAULT value sets the variance of this distribution to 0.04, equivalent to a 95\ is in the range of 0.66-1.5 at any causal variant.

Value

a vector of posterior probabilities.

Author(s)

Olly Burren, Chris Wallace, Guillermo Reales


GRealesM/RapidoPGS documentation built on Feb. 2, 2021, 2:58 p.m.