Struct libprosic::model::SingleCaller

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

Joint variant calling model, combining two latent variable models.

Implementations

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

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

Create new JointModel.

Arguments

• sample - sample
• prior_model - prior model

pub fn pileup(
    &self,
    chrom: &[u8],
    start: u32,
    variant: Variant,
    chrom_seq: &[u8]
) -> Result<SinglePileup<'_, A, P>, Box<dyn 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.