Comments (3)
Hello @kbarrr the tutorial on using hmmratac
command in the MACS3 (version 3.0) can be found here:
https://macs3-project.github.io/MACS/docs/hmmratac.html
The command line options are different from the callpeak
command in MACS3.
Let us know if you find any issues when running hmmratac
in MACS3. We are still optimizing it. @philippadoherty is currently working on reducing the memory usage for this command.
from macs.
Thank you for your response.
I followed the instructions:
First i run
macs3 hmmratac -i M1_noM_sorted.bam -f BAMPE --cutoff-analysis-only
to know the upper and lower cut off and got this TSV file:
(macs3-env) *$ head *tsv
score npeaks lpeaks avelpeak
990.00 5 7609 1521.80
980.00 5 7643 1528.60
970.00 5 7650 1530.00
960.00 5 7654 1530.80
950.00 5 7659 1531.80
940.00 5 7662 1532.40
930.00 5 7665 1533.00
920.00 5 7673 1534.60
910.00 5 7679 1535.80
(macs3-env) *$ tail *tsv
90.00 171 313505 1833.36
80.00 201 353408 1758.25
70.00 263 415032 1578.07
60.00 385 523102 1358.71
50.00 615 695166 1130.35
40.00 1107 995113 898.93
30.00 2509 1774048 707.07
20.00 5513 3736145 677.70
10.00 14149 11493688 812.33
0.00 169628 2841257445 16749.93
i am expecting ~15-30 K peaks because they are CD4 cells.
would you help me choose my upper and lower cut off please? i am lost
from macs.
@kbarrr Thanks for the information! However, the resolution of the cutoff analysis report is not good enough for me to judge where the cutoff should be -- each step is 10 according to the first column. I would guess the settings of '-u 20 -l 10 -c 2' (or -c 5
) should work. We will release a new version soon to address 1) the resolution of cutoff analysis -- to make it with a step of 1 by default #642; 2) to reduce memory usage by a magnitude for large datasets #640 ; 3) to provide a simpler HMM model with Poisson Emission #635. If you can wait, please try the next release.
from macs.
Related Issues (20)
- Q: Defining the bin size when generating bedgraphs.
- Bug: test suite failure in MACS 3.0.0 HOT 18
- Bug: pip install error: subprocess-exited-with-error HOT 2
- What is the difference between "absolute peak summit" and "summit position" in narrowPeak format? HOT 1
- Q: conversion of output files to UCSC track hub format? HOT 1
- Feat: Enhanced clarification of the specific definition of BEDPE as used by MACS HOT 1
- Feat: Reduction of memory consumption HOT 3
- Q: HMMRATAC output and options for _summit.bed HOT 1
- Feat: Built-in scoring for HMMRATAC HOT 1
- Q: HMMRATAC producing too many peaks HOT 10
- Bug: undefined symbol: __pow_finite HOT 3
- Q: HMMRATAC reproducibility HOT 2
- Bug: Callvar throws an error when run on broadPeak file HOT 15
- Q: impact of sequencing throughput on peak calling HOT 4
- Setting cut-off when using logFE bigWig files from bdgcmp
- Q: How to use HMMRATAC to give us a narrow peak file and have signalValue, p and q-value ? HOT 1
- Bug: OverflowError when running hmmratac
- Number of de-duplicated reads indicated by filterdup log does not match number of lines in output BED file.
- Feat: Add --call-summits option in MACS3 bdgpeakcall?
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 macs.