Question: ChIPSeeker Gene ID not matched error
0
gravatar for Chenzhen Zhang
5 months ago by
Rutgers University
Chenzhen Zhang0 wrote:

Hi I am using the ChIPSeeker package to do gene annotation for ChIP-seq data. broadpeak files were generated with MACS2. Here is my code.

library('dplyr')
library("ChIPseeker")
library("org.Ce.eg.db")
library("TxDb.Celegans.UCSC.ce6.ensGene")
library("clusterProfiler")
txdb <- TxDb.Celegans.UCSC.ce6.ensGene

N2_K9_0.05<-readPeakFile("N2K9CallPeak_450_0.05_maxgap1000_peaks.broadPeak.bed")
peakAnno <- annotatePeak(N2_K9_0.05, tssRegion=c(-3000, 3000), 
                         TxDb=txdb, annoDb="org.Ce.eg.db")

I got this warning message:

Warning message:
In getGeneAnno(annoDb, peak.gr$geneId, type) :
  ID type not matched, gene annotation will not be added...

HELP! How can I match the ID type? I checked the txdb and and annoDb files, both of which have Ensembl gene ID.

chip-seq chipseeker • 234 views
ADD COMMENTlink modified 5 months ago by ATpoint22k • written 5 months ago by Chenzhen Zhang0
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 2.3.0
Traffic: 1141 users visited in the last hour