I revisited the MACAU result again to:

• re-do heat maps with only loci that has 10x or greater coverage
• Generate heat maps of % methylation (and coverage >= 10x)
• Generate PC plot using count data (>= 10x coverage)

Check out my RMarkdown notebook, 06-Analyzing MACAU results