Thank you I got it!
features = attributes(peakAnno)
head(features['detailGenomicAnnotation'][[1]])
the link is: https://github.com/GuangchuangYu/ChIPseeker/issues/43
How to get the all the overlap information of a peak in ChIPseeker
Hi guys, I think there must be all the annotation information of one peak in ChIPseeker analysis, for example, a peak could overlap with an exon or an intron. But when I look into the annotation file using as.GRanges method, I could only see one feature corresponding to one peak. I use this command line:
anno<-as.data.frame(as.GRanges(peakAnno))
write.table(anno,"anno.txt",sep="\t",quote=F)
Maybe there is other files that storing all the information? How to get this? Thanks!
1 answer
peakAnno@detailGenomicAnnotation
• 1 views
•
link
written
by
Guangchuang Yu
Log in to answer this question.
More posts like this
-
line 1: syntax error near unexpected token `Rsubread
written by sehriban.buyukkilic •
-
Why not use ONLY promoter-bound peaks when testing for enrichment in differentially-bound regions?
written by e.r.zakiev •
-
Question about ChIPseeker parameters
written by Jen •
-
Statistical test for genomic annotation by ChIPseeker
written by arisarkar88 •
-
What is the codification in genestrand 1 and 2?
written by Lila M
-
Chip-seq peak annotation [Chipseeker]
written by anu014
-
s pombe non-coding regions annotation
written by Lila M
-
warning trying to create a TxDB for S.pombe
written by Lila M
-
Process multiples files
written by Lila M
-
Question About Chippeakanno And Iranges
written by e.karasmani •
overlap-of-peaks-and-annotated-genes section given in ChIPseeker manual may help you.
Thank you my friend, and I found another answer below.