今天做展示PPT的时候发现可以用一种更简单的方式计算之前计算过的内容，因此更新一下之前的结果。
1.相关性分析结果

> cor.test(merge_1$value, merge_1$rate,method = c("pearson"))

    Pearson's product-moment correlation

data:  merge_1$value and merge_1$rate
t = -6.7244, df = 198313, p-value = 1.769e-11
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 -0.01949810 -0.01069772
sample estimates:
       cor 
-0.0150982 

2.组织特异性甲基化的启动子区域

> nrow(methpromspec)
[1] 2006
> nrow(exprsetspec)
[1] 12059
> nrow(merge_2)
[1] 24

得到2006个组织特异性甲基化的启动子区域。（由17702个启动子区域得到）
overlap数目为24个。