methcounts will soon be updated to mirror bsrate's mechanism of keeping running tabulations of an entire chromosome's methylation counts and then computing all methylation levels when all the reads for that chromosome have been processed.
allelicmeth could benefit from this structure as well: instead of keeping tabulations of a chromosome's methylation counts, it could keep track of the 2-CpG methylation state counts ([0,0] [0,1] [1,0] [1,1]), and at the end of a chromosome could quickly calculate p-values for the allelic-ness of each pair of CpGs.
methcounts will soon be updated to mirror bsrate's mechanism of keeping running tabulations of an entire chromosome's methylation counts and then computing all methylation levels when all the reads for that chromosome have been processed.
allelicmeth could benefit from this structure as well: instead of keeping tabulations of a chromosome's methylation counts, it could keep track of the 2-CpG methylation state counts ([0,0] [0,1] [1,0] [1,1]), and at the end of a chromosome could quickly calculate p-values for the allelic-ness of each pair of CpGs.