Correlation Coefficient Between Chip-Seq Replicates
5
2
Entering edit mode
12.4 years ago
Shumon ▴ 110

Any one knows any command line tools/methods/scripts to calculate correlation coefficient from two BED files of two ChIP-seq replicates?

thanks

chip-seq correlation • 12k views
ADD COMMENT
5
Entering edit mode
12.4 years ago

A very similar question was asked previously. I suggested a recent publication by Zhao and Sandelin (2012). They provide an R package called 'GMD' for calculating the "similarity between spatial distributions of read-based sequencing data such as ChIP-seq and RNA-seq". They also provide a detailed set of case studies (vignettes) in CRAN using some ChIP-seq data to illustrate. The authors seem responsive and showed up here on Biostars with some help for the previous poster.

ADD COMMENT
1
Entering edit mode
12.4 years ago
Ian 6.1k

Not quite what you've asked for, but htSeqTools (Bioconductor) can produce a PCA plot of sequence coverage similarity between samples.

ADD COMMENT
1
Entering edit mode
12.4 years ago
Ying W ★ 4.3k

A quick and dirty (though space consuming) way is to create per-base counts for the whole genome and then find correlation when you exclude all the places in the genome with read count of 0. This method is outline in in this paper "A computational pipeline for comparative ChIP-seq analyses (2011)" and the shell script to do the correlation can be found here under 'correlation.awk'

ADD COMMENT
0
Entering edit mode
10.9 years ago
Matt • 0

BiSA (http://bisa.sourceforge.net) statistical module does this, PloS One manuscript is in press. It is called overlap correlation value (OCV) in BiSA.

ADD COMMENT
0
Entering edit mode
10.7 years ago
dnaseiseq ▴ 220

abotu PCC

Bailey T, Krajewski P, Ladunga I, Lefebvre C, Li Q, Liu T, Madrigal P, Taslim C, Zhang J (2013) Practical Guidelines for the Comprehensive Analysis of ChIP-seq Data. PLoS Comput Biol 9(11): e1003326. doi:10.1371/journal.pcbi.1003326

http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.1003326

ADD COMMENT

Login before adding your answer.

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