Function libprosic::estimation::fdr_bh::control_fdr

pub fn control_fdr<E: Event, W: Write>(calls: &mut Reader, null_calls: &mut Reader, writer: &mut W, event: &E, vartype: &VariantType) -> Result<(), Box<Error>>

Print thresholds to control FDR of given calls at multiple levels.

Arguments

• calls - BCF reader with prosic calls
• null_calls - calls under the null model, e.g. obtained by swapping tumor and normal sample
• writer - writer for resulting thresholds
• event - the event to control
• vartype - the variant type to consider