calculateThreshold {iBMQ} | R Documentation

## Calculate PPA significance threshold leading to a desired false discovery rate

### Description

In the context of multiple testing and discoveries, a popular approach is to use a common threshold leading to a desired false discovery rate (FDR). In the Bayesian paradigm, derivation of the PPA threshold is trivial and can be calculated using a direct posterior probability calculation as described in Newton et al. (2004).

### Usage

    calculateThreshold(prob, threshold)

### Arguments

`prob` | matrix or data frame that contains Posterior Probability of Association (output of eqtlMcmc function).
`threshold` | The desired false discovery rate.

### Value

`cutoff` | The significance threshold value

### References

Newton, MA., Noueiry, A., Sarkar, D. and Ahlquist, P. (2004): "Detecting differential gene expression with a semiparametric hierarchical mixture method."Biometrics, 5(2), 155-176

### Examples

    data(PPA.liver)
    cutoff.liver <- calculateThreshold(PPA.liver, 0.2)

[Package iBMQ version 1.22.0 Index]