代码库12-ChIPseeker注释基因组区域
作者:
江湾青年 | 来源:发表于
2024-01-06 14:02 被阅读0次library(ChIPseeker)
library(TxDb.Hsapiens.UCSC.hg38.knownGene)
library(org.Hs.eg.db)
library(clusterProfiler)
options(ChIPseeker.ignore_1st_exon = T)
options(ChIPseeker.ignore_1st_intron = T)
options(ChIPseeker.ignore_downstream = T)
options(ChIPseeker.ignore_promoter_subcategory = T)
peak_regions <- annotatePeak(SeuratObj_ATAC@assays$ATAC@ranges, level = 'gene',tssRegion = c(-1000, 1000),
TxDb=TxDb.Hsapiens.UCSC.hg38.knownGene, annoDb = "org.Hs.eg.db")
head(peak_regions@anno)
plotAnnoPie(peak_regions,legend.position = 'top')
image.png
本文标题：代码库12-ChIPseeker注释基因组区域
本文链接：https://www.haomeiwen.com/subject/gmbindtx.html
网友评论