View Single Post
Old 05-09-2011, 11:26 AM   #3
NearyJL78
Member
 
Location: San Antonio, Texas

Join Date: Feb 2011
Posts: 14
Default

EDIT - the paragraph below is referring to the next steps in MEDIPS analysis, sorry! I'm typing the answer to your question (AS I UNDERSTAND IT) after the ***


I am just starting with MEDIPS also, also looking at MeDIP-seq data. I am still somewhat puzzled by the analysis, but I know the higher the correlation values the better. Roughly, they represent how much your sample has been enriched for CpGs if you imagine the genome to have a value of 1. The higher above one you are, the more likely the beads are pulling out what you want (methylated CpGs). That being said, I am not sure which correlation value is the more important one. The author of the program is quick to respond and was helpful for some of my questions.
Best of luck!



***
You are looking at genomic coverage essentially - whether the entire genome is covered by your reads or not. Again, the higher the value the better, but you'd want at least one. The way the math is done in the program, your reads are divided into two groups, A and B. The true correlation is the coverage using half your reads. The estimated correlation is half your reads but doubled... so essentially you would expect your true genome coverage to be somewhere in between. Most alignment programs will also give you a coverage value. The main thing is you want to make sure you are covering the reference ENOUGH, without overkill because that wastes money on sequencing. What value range you need really depends on your experiment. In our case we like to see at least 10x but less than 50x coverage, and if we are seeing too much coverage we will multiplex several samples into one lane of the sequencing plate.

I hope this helps! Sorry for the mix-up.

Last edited by NearyJL78; 05-09-2011 at 12:21 PM.
NearyJL78 is offline   Reply With Quote