get_log10bf: Return the Bayes Factor for each effect

Description Usage Arguments Value Examples

View source: R/get_functions.R

Description

Return the Bayes Factor for each effect

Usage

1

Arguments

m

the mash result (from joint or 1by1 analysis); must have been computed using usepointmass=TRUE

Value

if m was fitted using usepointmass=TRUE then returns a vector of the log10(bf) values for each effect. That is, the jth element lbf[j] is log10(Pr(Bj | g=ghat-nonnull)/Pr(Bj | g = 0)) where ghat-nonnull is the non-null part of ghat. Otherwise returns NULL.

Examples

1
2
3
4
simdata = simple_sims(50,5,1)
data = mash_set_data(simdata$Bhat, simdata$Shat)
m = mash(data, cov_canonical(data))
get_log10bf(m)

mashr documentation built on May 24, 2021, 1:06 a.m.