Description Usage Arguments Value References See Also Examples

View source: R/locFDR_BF_theoretic.R

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`

.

1 | ```
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`

1 2 3 4 5 6 7 | ```
data(ExampleDataUncor)
BetaHat <- ExampleDataUncor$BetaHat
BetaHat
SE <- ExampleDataUncor$SE
SE
result <- analytic_locFDR_BF_uncor(BetaHat, SE)
str(result)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.