今天做展示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个。
网友评论