Struct libprosic::model::PairCaller

pub struct PairCaller<A: AlleleFreqs, B: AlleleFreqs, P: PairModel<A, B>> { /* fields omitted */ }

Joint variant calling model, combining two latent variable models.

Methods

impl<A: AlleleFreqs, B: AlleleFreqs, P: PairModel<A, B>> PairCaller<A, B, P>

fn new(case_sample: Sample, control_sample: Sample, prior_model: P) -> Self

Create new PairCaller.

Arguments

• case_sample - case sample
• control_sample - control sample
• prior_model - prior model

fn pileup(
    &self,
    chrom: &[u8],
    start: u32,
    variant: Variant,
    chrom_seq: &[u8]
) -> Result<PairPileup<A, B, P>, Box<Error>>

Calculate pileup and marginal probability for given variant.

Arguments

• chrom - the chromosome of the variant
• start - the starting position of the variant
• variant - the variant

Returns

The PairPileup, or an error message.