search

aertslab / cisTopic

cisTopic: Probabilistic modelling of cis-regulatory topics from single cell epigenomics data
135 stars 29 forks source link

getBigwigFiles functions error #51

Closed ytang0831 closed 2 years ago

ytang0831 commented 4 years ago

Hi! I tried to generate BW file after get region scores.

library(TxDb.Mmusculus.UCSC.mm9.knownGene) txdb<-TxDb.Mmusculus.UCSC.mm9.knownGene getBigwigFiles(cisTopicObject, path='output/cisTopics_asBW', seqlengths=seqlengths(txdb))

However, I get the error message below: Error in FUN(extractROWS(unlisted_X, IRanges(X_elt_start[i], X_elt_end[i])), : BigWig ranges cannot overlap

Looking forward to your Suggestions！

cbravo93 commented 4 years ago

Hi @ytang0831 !

Do your regions overlap?

C

luzgaral commented 4 years ago

Hi @cbravo93

I'm facing the same issue. Because of the way I call peaks, they may overlap in 1bp. Any suggestion?

Thanks!