ChIPseeker icon indicating copy to clipboard operation
ChIPseeker copied to clipboard

Published 20 hours ago • verified publisher icon YuLab-SMU

bug fix #201

Open MingLi-929 opened this issue 2 years ago • 0 comments

this commit fix bug #201, after updating,

library(TxDb.Hsapiens.UCSC.hg19.knownGene)

files <- getSampleFiles()
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene

peakAnno <- annotatePeak(files[[4]], 
                         tssRegion=c(-3000, 3000),
                         TxDb=txdb, annoDb="org.Hs.eg.db")
plotAnnoPie(peakAnno)

1

nrow <- nrow(peakAnno@annoStat)

col <- c("#8dd3c7", "#ffffb3", "#bebada",
         "#fb8072", "#80b1d3", "#fdb462",
         "#b3de69", "#fccde5", "#d9d9d9",
         "#bc80bd", "#ccebc5", "#ffed6f")

cols <- col[1:nrow]

plotAnnoPie(peakAnno,col = col[1:nrow])

2

MingLi-929 avatar Oct 15 '22 12:10 MingLi-929