今天做展示PPT的時候發(fā)現(xiàn)可以用一種更簡單的方式計算之前計算過的內(nèi)容泡徙,因此更新一下之前的結(jié)果。
1.相關(guān)性分析結(jié)果
> 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.組織特異性甲基化的啟動子區(qū)域
> nrow(methpromspec)
[1] 2006
> nrow(exprsetspec)
[1] 12059
> nrow(merge_2)
[1] 24
得到2006個組織特異性甲基化的啟動子區(qū)域。(由17702個啟動子區(qū)域得到)
overlap數(shù)目為24個羹蚣。