Comments (3)
It’s impossible to call ‘depleted’ regions from ‘callpeak’ command since it is designed to find ‘enriched’ regions only. If you have MeDIP data from two conditions, you can consider ‘bdgdiff’ command after you get four bdg files from two separate runs of ‘callpeak’ on both conditions.
To analyze MeDIP data, you may consider more specific tools such as MEDIPS or MnM.
Tao
On Jun 23, 2014, at 7:28 AM, biocyberman [email protected] wrote:
I noticed that after callpeak command with treat vs control MeDIP samples, I don't have peaks with fold-changes < 1, which, I suppose, are demethylated spots. Is it correct to expect fold-change < 1 for demethylated peaks? In general, a plain English explanation of how peaks and fold-changes are calculated will be very helpful.
I tried looking for that information in the publication http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3120977 but could not find anything. The last resort could be to look into the source code, but I prefer to have the authors' words first.
—
Reply to this email directly or view it on GitHub.
from macs.
Hello Tao,
Thank you for very helpful hints. I will take a closer look at MEDIPS and MnM. However, just to think about a way (around) to use MACS for 'depleted' peaks also, would it still be OK to switch 'treated' and 'control' samples to call these peaks? That means, I first run
macs2 callpeak -t treat_sample -c control_sample
and then run macs2 callpeak -t control_sample -c treat_sample
Afterward I would do some preprocessing of peak_names, fold_changes before JOINING the two peak_list together so I will have both 'enriched' and 'depleted' peaks.
My concern is that it may not be correct to join the two output together because may be normalized differently. What is your thoughts about this?
Thanks
Vang
from macs.
Sorry for late reply! I am cleaning issues nowadays and just saw your question...
Yes. You can do so. By default, MACS2 uses size of smaller dataset to normalize, so in your case, this step is fine. But MACS2 will use control data to calculate so-called local background in 1,000bps, 10,000 bps surrounding regions. This may not be a good solution for your methylation data since you may want to compare two libraries at the same position. You can try to disable this by --slocal 0 --llocal 0.
I close this issue for now, since it's more about using MACS2 not about coding. If you still have question, please send email to MACS google user group or email me directly.
-T
from macs.
Related Issues (20)
- 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
- Q: how to use MACS3 for ATAC seq with "himmratac" option? 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?
- Q: using macs3 to detect peaks on long-read readings
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.