Question: ChIPSeeker Gene ID not matched error
0
gravatar for Chenzhen Zhang
7 weeks 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 • 114 views
ADD COMMENTlink modified 7 weeks ago by ATpoint16k • written 7 weeks 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: 1408 users visited in the last hour