Question: Control data for chip-seq analysis of SRX681547, SRX681548 data
0
gravatar for gskbioinfo143
2.8 years ago by
India
gskbioinfo14350 wrote:

Dear all,

im working on chip-seq analysis of SRX681547, SRX681548 data using galaxy suite for peak calling, motif finding and differential expression analysis. i have query regarding control data set while doing MACS analysis. for this dataset can we do MACS analysis without control data?. if not then what has to be used as control data.

Need Help.

Hope to hear from you

Thank you

chip-seq • 756 views
ADD COMMENTlink modified 2.8 years ago by Devon Ryan90k • written 2.8 years ago by gskbioinfo14350
2
gravatar for Satyajeet Khare
2.8 years ago by
Satyajeet Khare1.3k
Pune, India
Satyajeet Khare1.3k wrote:

You can perform MACS peak calling without control. But since you are referring to a published dataset, you can refer to the publication or contact the authors for their analysis pipeline. I am not sure how you will perform differential expression analysis with MedIP data.

ADD COMMENTlink written 2.8 years ago by Satyajeet Khare1.3k
1
gravatar for Devon Ryan
2.8 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

MACs won't be very useful, you're not looking for peaks, you're looking in differences in methylation, which generally looks more like a dip in the signal instead. You're better off using something like batman than either MACs or homer (what they did in the paper that created those samples).

ADD COMMENTlink written 2.8 years ago by Devon Ryan90k

Hi im actually looking for unique enriched peaks with p values from treatment (SRX681547) and control (SRX681548) samples using peak calling. I have used MACS in galaxy to get peaks (.bed) from both samples then i inputted into SICER to get enriched regions but analysis fails.

Is there any tool which gives unique enriched peaks taking treatment and control samples .bed peak files

need help thank you

ADD REPLYlink written 2.7 years ago by gskbioinfo14350

Bed intersect should give you unique peaks in treatment unique peaks in control, overlapping peaks in treatment and control, or all peaks in both treatment and control.

ADD REPLYlink written 2.7 years ago by Satyajeet Khare1.3k

Hi,

in MACS output files, the peak interval file has -10*log10(pvalue) values for each peaks.

know how to get the original pvalues for each peak data. In result files it gave calculated value using ablove formula. shall i reverse the formula to get that pvalue?

thank you

ADD REPLYlink written 2.7 years ago by gskbioinfo14350

Hi Sir,

in MACS output files, the peak interval file has -10*log10(pvalue) values for each peaks.

know how to get the original pvalues for each peak data. In result files it gave calculated value using ablove formula. shall i reverse the formula to get that pvalue?

thank you

ADD REPLYlink written 2.7 years ago by gskbioinfo14350
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1486 users visited in the last hour