minor tweaks

[samtools.git] / bam_mcns.c

diff --git a/bam_mcns.c b/bam_mcns.c
index d98d48030560a42fd65c00bda4f1633bc90ac169..d337e5a290aa76a1f1faa2c8a687db4639962a16 100644 (file)
--- a/bam_mcns.c
+++ b/bam_mcns.c
@@ -264,8 +264,8 @@ static void mc_add_afs(mc_aux_t *ma, double PD, double *f_map, double *p_map)
                        if (ma->afs1[k] > max) max = ma->afs1[k], max_k = k;
                        e += k * ma->afs1[k];
                }
-               *f_map = .5 * max_k / ma->n; *p_map = max;
-               printf(" * %.3lg:%.3lg:%.3lg:%.3lg * ", sum, 1.-.5*max_k/ma->n, max, 1.-.5*e/ma->n);
+               *f_map = .5 * max_k / ma->n; *p_map = max; // e should equal mc_rst_t::f_exp
+//             printf(" * %.3lg:%.3lg:%.3lg:%.3lg * ", sum, 1.-.5*max_k/ma->n, max, 1.-.5*e/ma->n);
        }
 }
 
@@ -297,3 +297,13 @@ int mc_cal(int ref, int *n, const bam_pileup1_t **plp, mc_aux_t *ma, mc_rst_t *r
                mc_add_afs(ma, rst->PD, &rst->f_map, &rst->p_map);
        return tot;
 }
+
+void mc_dump_afs(mc_aux_t *ma)
+{
+       int k;
+       fprintf(stderr, "[afs]");
+       for (k = 0; k <= ma->M; ++k)
+               fprintf(stderr, " %d:%.3lf", k, ma->afs[ma->M - k]);
+       fprintf(stderr, "\n");
+       memset(ma->afs, 0, sizeof(double) * (ma->M + 1));
+}