Annotating peak using Chipseeker for non-model organism
I have a chip-seq data analyzed using macs3 to get the peak files. I would like to annotate those files using the ChIPSeeker package. Unfortunately, I'm working aon an organism without a TxDB or org*.db packages.
I know I can create a TxDb object from a GTF file using the GenomicFeatures:: makeTxDbFromGFF function.
The function annotatePeaks needs though an annoDB object (like org.Hs.eg.db for the human genome).
Is there a way to create this object from scratch? If so, what files does one needs to have for that?
thanks
Assa
• 1,333 views
•
link
updated
by
rpolicastro
•
written
by
Assa Yeroslaviz
•
written
by
Assa Yeroslaviz
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Question about ChIPseeker parameters
written by Jen •
-
How to include the metadata when we build the txDb object?
written by qwesxadzc9 •
-
Creating taxonomic database using GFF3 file ?
written by lokraj2003
-
Chip-seq peak annotation [Chipseeker]
written by anu014
-
chip-seq analysis and peak annotation
written by tarek.mohamed
-
how to exclude set of genes from a txdb
written by tarek.mohamed
-
Chip-Seq peaks annotation
written by tarek.mohamed
-
ChIPSeeker - problem with GO enrichment analysis
written by jk •
-
Process multiples files
written by Lila M
-
Peak annotation using ChIPseeker (bioconductor package)
written by s.singh •
You don't actually need to provide the annoDB object, it's optional. Most of the extra information it would have added is present in the GTF file to, so you can just merge the information into the results after you run the function.