Comments (4)
Hi,
I was wondering if you could provide more information on how to use the different parameters in callStripes to optimize the calls. I run it on my data using default parameters and while I notice there are some calls that are accurate, many of them are not visually obvious. The algorithm is also missing stripes that tend to be longer and wider.
I am attaching some screenshots.
Dear User,
The parameters as default were determined by HiChIP H3K27ac data, of which the stripes could be more kind of "sharp" compared to Hi-C data I guess.
Could you also please show me some examples the stripes were well called?
For how to tune key parameters , hope following explanations could help.
Increase the key parameter -eps from 20000 to 50000 may leads to longer and wider stripes.
Additionally, I would suggest ignore the significant flag as the last column of .strip file (example ), convert all called potential stripes to juicebox visible file and determine your cutoffs, and then modify the code of callStripes from line 226-255 to set your cutoffs.
Hope this will be helpful.
Could you please email me some example data if you think it's better for me to determine the parameters? BEDPE format for small chromosome such as chr21 would be great.
Best,
Yaqiang
from cloops.
Thanks a lot Yaqiang for your suggestions and help! I am attaching the jd file for chr19 (mm10) from calling loops that I used for calling stripes.
I would also try increasing the eps parameter and see if there is any improvement. Could you also explain what the -ext mean and whether increasing could also help ?
from cloops.
Hi Yaqiang,
I tried running callStripes with the -eps parameter, but it is crashing no matter what value I give to this parameter (even providing the default value of 20000 fails) . I get a long list of errors, but I am taking a screenshot of what I think are the key errors. This does not happen if I don't the eps parameter, and I can successfully get stripe calls with the default parameters.
I am using:
$ conda list -n cLoops | grep joblib
joblib 0.11 py27_0 conda-forge
$ cLoops -v
cLoops v0.92
( I think the cLoops version is actually 0.93 according to the setup.py file, but your -v parameter calls a different version)
from cloops.
-ext option means to enlarge the another axis, to make the stripe region as a rectangle (a loop), it will help.
For the bug, I have made some modifications to the code (add the type=int for input eps and minPts) and push it to the repo. Could you please clone the new code and try? For your data, other eps can run successfully with modifications.
from cloops.
Related Issues (20)
- Input format HOT 1
- Some details about the parameter 'eps' HOT 4
- adjust parameters to increase/decrease stringency HOT 5
- calling loops on a .hic file HOT 5
- HiC input from TADbit HOT 3
- Distance cut-off HOT 2
- joblib error HOT 5
- Conda env create issue HOT 8
- OpenSSL HOT 1
- cLopps hang for more than 24 hours HOT 1
- How to plot the loops in WashU HOT 1
- hicpro2bedpe bug HOT 9
- ERROR: no inter-ligation PETs detected for eps 5000 minPts 50,can't model the distance cutoff,continue anyway HOT 17
- Extracting enrichment score from a matrix given a list of loops HOT 1
- call stripe example file has error HOT 1
- cLoops differential loop calling HOT 1
- using cloops from Juicer HOT 2
- Loops Visualization HOT 4
- Issue with file preprocessing. HOT 8
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 cloops.