Comments (5)
Hi,
I tested the command on some sample data, and it works fine in my setup. Is it possible that your file is, in fact, empty, because of the filtering you applied? To check, could you send me the output of
python -c "import fanc; loops = fanc.load('chrX.filtered.noSinglets.merged.loops'); print(len(loops))"
please?
from fanc.
Thanks for the quick reply, I have attached my sample loops file below.
Thanks[
chrX.filtered.noSinglets.merged.loops.gz
](url)
from fanc.
Yes, there are no valid loops in that file. I would advise you to try adjusting your filter settings to be less strict. Getting the settings right can be really difficult, in our experience.
from fanc.
I'm working on Drosophila genome and this is the command (write below) how I get the loops. Could you please suggest me if I can try some other parameters.
3.7.0/bin/fanc loops chrX.hic chrX.loops -t 16
3.7.0/bin/fanc loops chrX.loops chrX.filtered.loops --rh-filter -d 5 -o 5
3.7.0/bin/fanc loops chrX.filtered.loops chrX.filtered.noSinglets.merged.loops -j --remove-singlets
3.7.0/bin/fanc loops chrX.filtered.noSinglets.merged.loops -b chrX.filtered.noSinglets.merged.loops.bedpe
Note: I have also tried using default -d and -o values but still my loops.bedpe file is empty.
Thanks
from fanc.
First of all, I'd like to note that the default loops
filtering parameters have been optimised for human and mouse genomes. Drosophila loops can have a very different appearance in Hi-C maps, and will likely require some experimentation to get the settings right.
If you haven't done so, read the original HICCUPS article to really understand the different neighbourhoods and procedures.
My recommendation is
- Run
fanc loops chrX.hic chrX.loops -t 16
as previously. This just annotates each pixel with information we need for filtering later on - Choose a set of loops from the original Hi-C matrices that you can clearly see by eye. You are going to use these to check if your settings work (you can use
fancplot
for this) - Run
fanc loops chrX.loops chrX.filtered.loops
with very lenient parameters. For example--enrichment 1.25 --fdr 0.2 --observed 5 --min-distance 5
- Run
fancplot <region> -p square chrX.filtered.loops
for each of the loops regions you identified in 2. This allows you to see which pixels have passed the filtering step - Adjust the parameters in 3. If you have removed the "loop" pixels, choose less stringent cutoffs. If there are too many pixels outside of your identified loops, especially clusters of pixels, choose more stringent cutoffs.
- Repeat 3-5 until you found a suitable set of parameters. When you are getting close, you can also fine-tune individual cutoffs. The most important ones will likely be
--enrichment-donut
and--fdr-donut
, but also have a look at--mappability
to reduce noise in your results. - Once you have a set of parameters that will filter out non-loop pixels, but leave in the loops you are interested in, you can run the merge and singlet removal steps (which will further reduce noise)
I hope that helps!
from fanc.
Related Issues (20)
- How to retrieve and save the compartment domain and values in python API HOT 1
- Issue with fancplot HOT 4
- 'Figure' object has no attribute 'draw_without_rendering' HOT 2
- Per chromosome compartment strength with different eigenvectors HOT 2
- is it SAM/BAM format?
- Issue with installation HOT 2
- fanc boundaries: AttributeError: 'Bed' object has no attribute 'close' HOT 2
- hic files format after the fanc auto pipeline HOT 1
- fancplot layer: AttributeError: 'Bed' object has no attribute '_parameters' HOT 1
- cooler version difference in parameters HOT 2
- Is it possible A/B compartment analysis entire genome
- Problem with fanc auto
- Load .insulation file (File type not recognised)
- Statistical methods employed comparing insulation scores HOT 1
- Fanc insulation hangs?
- Specific package versions for matplotlib and pandas? HOT 1
- Wrong name for parameter ylim
- How to do plotting.TriangularMatrixPlot without ice balancing HOT 1
- The command fanc compartments -r 5 encountered an error
- illegal hardware instruction fanc –version
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 fanc.