Struct libprosic::model::SingleCaller

pub struct SingleCaller<A: AlleleFreqs, P: Model<A>> { /* fields omitted */ }

Joint variant calling model, combining two latent variable models.

Methods

impl<A: AlleleFreqs, P: Model<A>> SingleCaller<A, P>

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

Create new JointModel.

Arguments

• sample - sample
• prior_model - prior model

fn pileup(
    &self,
    chrom: &[u8],
    start: u32,
    variant: Variant,
    chrom_seq: &[u8]
) -> Result<SinglePileup<A, 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 SinglePileup, or an error message.