ChIP-seq processing pipeline
This page last changed on Jun 11, 2009 by peterk.

The page provides download links to an R package for processing ChIP-seq data. The rational behind most of the methods is described in the following paper: Kharchenko PK, Tolstorukov MY, Park PJ "Design and analysis of ChIP-seq experiments for DNA-binding proteins" Nat. Biotech. doi:10.1038/nbt.1508

Releases and updates
version 1.7

Download. Changes: fixed a step filtering problem with get.mser.interpolation (thanks to Joe Luquette), fixed remove.tag.anomalies flag in get.binding.characteristics which wasn't being checked.

version 1.6

Download. Changes: Changed header files to accommodate newer GCC versions; Fixed a big in write.broadpeak.info.

version 1.5

Download. Changes: Added support for the extended Eland format (extended=T option for read.eland.tags).

version 1.4

Download. Changes: Added support for bowtie output, basic routines for broader regions of enrichment.

version 1.3

Download. Changes: fixed topN threshold is used together with background anomaly subtraction.

version 1.2

Download. Changes: added ability to read in MAQ alignments, skip binding.characteristics calculations for select.informative.tag call, improved effective dataset size estimation.

version 1.1

Download. Change: capped reported FDR estimates for peaks taller than any control peaks at the lowest FDR reached by the dataset

version 1.0

Initial package release. Download Unix version.

Features

  • Assess overall DNA-binding signals in the data and select appropriate quality of tag alignment.
  • Discard or restrict positions with abnormally high number of tags.
  • Calculate genome-wide profiles of smoothed tag density and save them in WIG files for viewing in other browsers.
  • Calculate genome-wide profiles providing conservative statistical estimates of fold enrichment ratios along the genome. These can be exported for browser viewing, or thresholded to determine regions of significant enrichment/depletion.
  • Determine statistically significant point binding positions
  • Assess whether the set of point binding positions detected at a current sequencing depth meets saturation criteria, and if does not, estimate what sequencing depth would be required to do so.

Installation instructions

Under unix use "R CMD INSTALL spp_1.7.tar.gz" to install the package.

Brief Tutorial

A brief tutorial showing use of the main methods can be found here.
Document generated by Confluence on Jun 11, 2009 20:14