fdrControl: Classification of genomic windows given an FDR control...

Description Usage Arguments Value Author(s) References Examples

View source: R/fdrControl.R

Description

For a given FDR control threshold and a set of posterior probabilities, 'fdrControl()' classifies genomic windows pertaining to a particular HMM state.

Usage

1
fdrControl(prob, fdr = 0.05)

Arguments

prob

numeric vector of posterior probabilities of genomic windows pertaining to a particular HMM state

fdr

FDR threshold (default is 0.05)

Value

A logical vector indicating whether genomic windows pertain to a particular HMM state (TRUE) or not (FALSE), under a given FDR control threshold.

Author(s)

Pedro L. Baldoni, pedrobaldoni@gmail.com

References

https://biostats.bepress.com/jhubiostat/paper115/

https://doi.org/10.1093/biostatistics/5.2.155

Examples

1
2
3
4
data(ENCODE)
ENCODE <- createOffset(object = ENCODE,type = 'loess',span = 1)
## Not run: ENCODE <- mixNBHMM(object = ENCODE,control = controlEM())
## Not run: fdrControl(metadata(ENCODE)$prob$Differential,fdr = 0.05)

plbaldoni/mixNBHMM documentation built on Dec. 24, 2019, 1:31 p.m.