Question: ChIPseeker genomicAnnotationPriority
0
gravatar for dolata
11 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 • 298 views
ADD COMMENTlink modified 11 months ago • written 11 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: 177 users visited in the last hour