YuLab-SMU / ChIPseeker

:dart: ChIP peak Annotation, Comparison and Visualization
https://onlinelibrary.wiley.com/share/author/GYJGUBYCTRMYJFN2JFZZ?target=10.1002/cpz1.585
219 stars 74 forks source link

error with Ensembl-based annotation #233

Open yfan18 opened 5 months ago

yfan18 commented 5 months ago

Describe you issue

Ensembl-based annotation,

edb <- EnsDb.Hsapiens.v97 seqlevelsStyle(edb) <- "UCSC" peakAnno.edb <- annotatePeak(myFile, tssRegion=c(-3000, 3000), TxDb=edb, annoDb="org.Hs.eg.db")

seqnames did not match geneChr. e.g. chr20 vs 2

<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns="http://www.w3.org/TR/REC-html40">

seqnames | start | end | annotation | geneChr | geneStart | geneEnd | geneLength | geneStrand -- | -- | -- | -- | -- | -- | -- | -- | -- chr20 | 58888051 | 58890615 | Promoter (<=1kb) | 2 | 58888413 | 58900410 | 11998 | 1

The same error for the example in you tutorial where ensDb.Hsapiens.v75 was used.

Ask in right place