Alternative model for multiallelic and rare-variant calling. Proof-of-principle imple...

[samtools.git] / bcftools / prob1.h

diff --git a/bcftools/prob1.h b/bcftools/prob1.h
index 0a51a0ae2650f377f4fc48665109487687eb69cb..88ed5bf0861ad709364b507fdfa39a24b825eca0 100644 (file)
--- a/bcftools/prob1.h
+++ b/bcftools/prob1.h
@@ -26,7 +26,9 @@ extern "C" {
        void bcf_p1_init_prior(bcf_p1aux_t *ma, int type, double theta);
        void bcf_p1_init_subprior(bcf_p1aux_t *ma, int type, double theta);
        void bcf_p1_destroy(bcf_p1aux_t *ma);
+    void bcf_p1_set_ploidy(bcf1_t *b, bcf_p1aux_t *ma);
        int bcf_p1_cal(const bcf1_t *b, int do_contrast, bcf_p1aux_t *ma, bcf_p1rst_t *rst);
+    int call_multiallelic_gt(bcf1_t *b, bcf_p1aux_t *ma, double threshold);
        int bcf_p1_call_gt(const bcf_p1aux_t *ma, double f0, int k);
        void bcf_p1_dump_afs(bcf_p1aux_t *ma);
        int bcf_p1_read_prior(bcf_p1aux_t *ma, const char *fn);