Function varlociraptor::filtration::fdr::control_fdr
source · pub fn control_fdr<R, W>(
inbcf: R,
outbcf: Option<W>,
events: &[SimpleEvent],
vartype: Option<&VariantType>,
alpha: LogProb,
local: bool,
smart: bool
) -> Result<()>
Expand description
Print thresholds to control FDR of given calls at multiple levels.
§Arguments
inbcf
- path to BCF with varlociraptor calls (None for stdin)outbcf
- path to BCF with filtered varlociraptor calls (None for stdout)null_calls
- calls under the null model, e.g. obtained by swapping tumor and normal samplewriter
- writer for resulting thresholdsevents
- the set of events to control (sum of the probabilities of the individual events at a site)vartype
- the variant type to consideralpha
- the FDR threshold to control for