Question: ChIPseeker genomicAnnotationPriority
0
gravatar for dolata
4 months ago by
dolata0
dolata0 wrote:

Hi, I'm trying to annotate peaks using annotatepeak (ChIPseeker) There is no problem using default genomicAnnotationPriority

peakAnno <- annotatePeak(files[[1]], tssRegion = c(-1000, 0), TxDb = txdb, level = "transcript", assignGenomicAnnotation = TRUE, genomicAnnotationPriority = c("Promoter", "5UTR", "3UTR", "Exon", "Intron", "Downstream", "Intergenic"), annoDb = "org.At.tair.db", addFlankGeneInfo = FALSE, flankDistance = 5000, sameStrand = FALSE, ignoreOverlap = FALSE, ignoreUpstream = FALSE, ignoreDownstream = FALSE, overlap = "TSS", verbose = TRUE)

However, when I move "Promoter" behind "Intron", promoters are gone from "Genomic Annotation Summary" table. I mean that promoters are not used for annotation and not that "Frequency" is 0.

Have anybody idea how to fix it?

Jakub

annotation chipseeker • 156 views
ADD COMMENTlink modified 4 months ago • written 4 months ago by dolata0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 165 users visited in the last hour