个人观点：将 ignoreDownstream 和 ignoreUpstream 的结果取并集更好的能代表 TF 结合的靶基因启动子区？？希望明天服务器早点恢复，看另外一个结果怎么样。
upstream = 1000, downstream = 200
peak2Gene_GO_igD(peak, ignoreDownstream = TRUE, GO = FALSE, Promoter = TRUE)
peak2Gene_GO_igU(peak, ignoreUpstream = TRUE, GO = FALSE, Promoter = TRUE)
# 合并两个结果
Promoter_Gene_y <- read.csv("Promoter_u1000_d200_igD.csv")
Promoter_Gene_x <- read.csv("Promoter_u1000_d200_igU.csv")
end_binding_Promoter <- rbind(Promoter_Gene_y[, c(1:4, 10)],
Promoter_Gene_x[, c(1:4, 10)]
) %>%
distinct(geneId, .keep_all = TRUE)
网友评论