Window size for MEDIPS
1
0
Entering edit mode
@anupriyaverma1408-13011
Last seen 6.5 years ago

Hello all,
I've been trying a package called 'MEDIPS' for MeDIP-seq analysis for few days now. The manual says : The genome will be divided into adjacent windows of length 100nt and all further calculations (short read coverage,
differential coverage between conditions etc.) will be applied to these windows, if ws =100. I've seen a paper (PMC4213696) apart from MEDIPS's original paper that researchers are using ws = 100 only. 

My question is what is the criteria for selecting window size? What if I increase the window size ? I guess the values such as 'counts', 'logFC' , 'rpkm' & eventually 'p-value' will change, is I do so.

Please help me out.

Thank you :)


package : https://www.bioconductor.org/packages/release/bioc/html/MEDIPS.html

medips medip-seq peak finding luckas • 1.4k views
ADD COMMENT
1
Entering edit mode
Lukas Chavez ▴ 570
@lukas-chavez-5781
Last seen 6.7 years ago
USA/La Jolla/UCSD
Dear anupriyaverma1408, you can certainly try other window sizes such es 300 or 500nt. This will have a strong impact on the number of significant cases after correction for multiple testing as you are massively reducing the number of windows tested for differential enrichment. Due to the same reason you should also make use of the minRowSum parameter (>=10) or focus on regions of interest such as previously calculated peak regions. For an improved MeDIP normalisation, please also consider our QSEA package. Best, Lukas > On 11. May 2017, at 11:39, anupriyaverma1408 [bioc] <noreply@bioconductor.org> wrote: > > Activity on a post you are following on support.bioconductor.org > User anupriyaverma1408 wrote Question: Window size for MEDIPS: > > > Hello all, > I've been trying a package called 'MEDIPS' for MeDIP-seq analysis for few days now. The manual says : The genome will be divided into adjacent windows of length 100nt and all further calculations (short read coverage, > differential coverage between conditions etc.) will be applied to these windows, if ws =100. I've seen a paper (PMC4213696) apart from MEDIPS's original paper that researchers are using ws = 100 only. > > My question is what is the criteria for selecting window size? What if I increase the window size ? I guess the values such as 'counts', 'logFC' , 'rpkm' & eventually 'p-value' will change, is I do so. > > Please help me out. > > Thank you :) > > > package : https://www.bioconductor.org/packages/release/bioc/html/MEDIPS.html > > > Post tags: medips, medip-seq, peak finding, luckas > > You may reply via email or visit Window size for MEDIPS >
ADD COMMENT

Login before adding your answer.

Traffic: 658 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6