includes CpG density and H3K79me2 in region of relevance (in gene body)
get_covar is now called within rate_analysis