Comments (4)
Hi @muzyuri ,
sorry for the late reply. Did you fixed it in the meantime?
If not what does the output of the makeBindingSites() function looks like? Does it produce the same error?
It seems like there is an issue with the names in of your GRanges object and/ or the seqnames (names for the chromosomes). Can you once make sure that the seqnames in the GRanges object match with the names of the signal object. You can do this by doing sth by looking into head(names(getSignal(bds)[[1]][[1]]))
and head(seqnames(getRanges(bds)))
.
Besides that, maybe the '_' in the seqnames is causing the trouble.
I hope this is still helpful for you, even tough the reply was so late.
from bindingsitefinder.
Hi @MirkoBr, thank you for your kindly support!
I am still stacking here.
After your reply, I checked my GRanges object and the seqnames and made sure that the seqnames in the GRanges object matched with the names of the signal object.
Besides, I removed all '_' and '-' in the seqnames and retried to proceed with the workflow like this:
supportRatioPlot(bds, bsWidths = seq(from = 3, to = 19, by = 2))
but I got the following another error:
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'flesh' in selecting a method for function 'relist': subscript contains out-of-bounds ranges
I also got this same error in the makeBindingSites() function.
Do you have any idea what is wrong? Thanks a lot in advance.
from bindingsitefinder.
Hi @muzyuri,
the error message that you get is a very general one. So without seeing some of your objects it's hard to tell where it goes wrong.
Could you provide a small subset of your data. Just like the first 10-20 entries in the GRanges plus the iCLIP signal for these ranges? To get both sets, see the code snippet below.
library(BindingSiteFinder)
# subset the GRanges object
rangesSubset = head(getRanges(bds), 20)
# subset the RLE signal list
totalSignal = signal$signalPlus$`1_WT`$chr22
startAt = max(as.data.frame(ranges(totalSignal))$start)
endAt = min(as.data.frame(ranges(totalSignal))$start)
signalSubset = totalSignal[startAt:endAt]
# save both objects
saveRDS(rangesSubset, file = "./rangesSubset.rds")
saveRDS(signalSubset, file = "./signalSubset.rds")
from bindingsitefinder.
fixed in 10ce962
from bindingsitefinder.
Related Issues (10)
- ! Can't recycle `quantiles` (size 101) to size 0. HOT 2
- Fix version number in devel HOT 1
- mergeCrosslinksDiagnosticsPlot() Description HOT 1
- Problems if not all samples have the same chromosomes HOT 2
- Target gene assignment in vignette: Problem with "_PAR_Y" HOT 1
- Error when extending plots with theme from ggplot
- 3.3 Transcript region assignment - typo? HOT 1
- help with an erro please HOT 11
- Errors related to calculateBsBackground HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bindingsitefinder.