关于 ChIPseeker 注释转录因子 Peak 的一点思考

个人观点：将 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)