Window size for MEDIPS
0
0
Entering edit mode
3.9 years ago
anu014 ▴ 180

Hello Biostars!

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

UPDATE :

The author has answered on Bioconductor:

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

medips medip-seq peakcalling sequencing • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1259 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