Function libprosic::utils::filter_by_threshold
source · pub fn filter_by_threshold<E: Event>(
calls: &mut Reader,
threshold: Option<LogProb>,
out: &mut Writer,
events: &[E],
vartype: &VariantType
) -> Result<(), Box<dyn Error>>
Expand description
Filter a VCF record stream by a minimum threshold on the sum of posterior probabilities of a given set of Events. The threshold should be an informative false discovery rate (FDR) threshold, e.g. determined with the libprosic FDR control functionality.
Arguments
calls
- BCF reader with libprosic callsthreshold
- minimum threshold for the sum of posterior probabilities of the set of Events consideredcalls
- BCF writer for the filtered libprosic callsevents
- the set of Events to filter onvartype
- the variant type to consider