analytic_locFDR_BF_uncor {CPBayes} | R Documentation |
Run the analytic_locFDR_BF_uncor
function to analytically compute the local FDR & Bayes factor (BF)
that quantifies the evidence of aggregate-level pleiotropic association for uncorrelated summary statistics.
Here a fixed value of slab variance is considred instead of a range of it in cpbayes_uncor
.
analytic_locFDR_BF_uncor(BetaHat, SE, SpikeVar = 1e-04, SlabVar = 0.8)
BetaHat |
A numeric vector of length K where K is the number of phenotypes. It contains the beta-hat values across studies/traits. No default. |
SE |
A numeric vector with the same dimension as BetaHat providing the standard errors corresponding to BetaHat. Every element of SE must be positive. No default. |
SpikeVar |
Variance of spike (normal distribution with small variance) representing the null effect distribution. Default is 10^(-4). |
SlabVar |
Variance of slab normal distribution representing the non-null effect distribution. Default is 0.8. |
The output produced by the function is a list which consists of the local FDR and log10(Bayes factor).
locFDR |
It provides the analytically computed local false discovery rate (posterior probability of null association) under CPBayes model (a Bayesian analog of the p-value) which is a measure of the evidence of the aggregate-level pleiotropic association. Bayes factor is adjusted for prior odds, but locFDR is solely a function of the posterior odds. |
log10_BF |
It provides the analytically computed log10(Bayes factor) produced by CPBayes that measures the evidence of the overall pleiotropic association. |
Majumdar A, Haldar T, Bhattacharya S, Witte JS (2018) An efficient Bayesian meta analysis approach for studying cross-phenotype genetic associations. PLoS Genet 14(2): e1007139.
cpbayes_uncor
, analytic_locFDR_BF_cor
, cpbayes_cor
, estimate_corln
, post_summaries
, forest_cpbayes
data(ExampleDataUncor) BetaHat <- ExampleDataUncor$BetaHat BetaHat SE <- ExampleDataUncor$SE SE result <- analytic_locFDR_BF_uncor(BetaHat, SE) str(result)