static void
update_dmr_name(Locus &dmr, double lod_sum, double mindiff_sum) {
ostringstream name_stream;
name_stream << "dmr:"
<< lod_sum / dmr.score() << ":"
<< mindiff_sum / dmr.score();
dmr.set_name(name_stream.str());
}
void
merge(istream &cpg_stream, ostream &dmr_stream, double cutoff) {
bool skipped_last_cpg;
Locus dmr;
if (!read_next_significant_cpg(cpg_stream, dmr, cutoff, skipped_last_cpg))
return;
dmr.set_score(1);
double lod_sum, mindiff_sum, unadjusted_pval;
extract_scores(dmr.name(), lod_sum, mindiff_sum, unadjusted_pval);
if (unadjusted_pval >= cutoff) {
lod_sum = 0;
mindiff_sum = 0;
}
Locus cpg;
while(read_next_significant_cpg(cpg_stream, cpg, cutoff, skipped_last_cpg)) {
double lod, mindiff;
extract_scores(cpg.name(), lod, mindiff, unadjusted_pval);
if (skipped_last_cpg || cpg.chrom() != dmr.chrom()) {
update_dmr_name(dmr, lod_sum, mindiff_sum);
if (lod_sum != 0) {
dmr_stream << dmr << std::endl;
}
dmr = cpg;
dmr.set_score(1);
if (unadjusted_pval >= cutoff) {
lod_sum = 0;
mindiff_sum = 0;
} else {
lod_sum = lod;
mindiff_sum = mindiff;
}
} else {
dmr.set_end(cpg.end());
dmr.set_score(dmr.score() + 1);
if (unadjusted_pval < cutoff) {
lod_sum += lod;
mindiff_sum += mindiff;
}
}
}
update_dmr_name(dmr, lod_sum, mindiff_sum);
dmr_stream << dmr << std::endl;
}