pub fn control_fdr<R, W>(
    inbcf: R,
    outbcf: Option<W>,
    events: &[SimpleEvent],
    vartype: Option<&VariantType>,
    alpha: LogProb,
    local: bool,
    smart: bool
) -> Result<()>
where R: AsRef<Path>, W: AsRef<Path>,
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 sample
  • writer - writer for resulting thresholds
  • events - the set of events to control (sum of the probabilities of the individual events at a site)
  • vartype - the variant type to consider
  • alpha - the FDR threshold to control for