pub fn control_fdr<E: Event, R, W>(
inbcf: R,
outbcf: Option<W>,
events: &[E],
vartype: Option<&VariantType>,
alpha: LogProb,
local: 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