visualization methods in ChIPseeker

April 30, 2014
By

(This article was first published on YGC » R, and kindly contributed to R-bloggers)

After two weeks developed, I have added/updated some plot functions in ChIPseeker (version >=1.0.1).

ChIP peaks over Chromosomes

?View Code RSPLUS
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
> files=getSampleFiles()
> peak=readPeakFile(files[[1]])
> peak
GRanges with 1331 ranges and 2 metadata columns:
         seqnames                 ranges strand   |             V4        V5
            <Rle>              <IRanges>  <Rle>   |       <factor> <numeric>
     [1]     chr1     [ 815092,  817883]      *   |    MACS_peak_1    295.76
     [2]     chr1     [1243287, 1244338]      *   |    MACS_peak_2     63.19
     [3]     chr1     [2979976, 2981228]      *   |    MACS_peak_3    100.16
     [4]     chr1     [3566181, 3567876]      *   |    MACS_peak_4    558.89
     [5]     chr1     [3816545, 3818111]      *   |    MACS_peak_5     57.57
     ...      ...                    ...    ... ...            ...       ...
  [1327]     chrX [135244782, 135245821]      *   | MACS_peak_1327     55.54
  [1328]     chrX [139171963, 139173506]      *   | MACS_peak_1328    270.19
  [1329]     chrX [139583953, 139586126]      *   | MACS_peak_1329    918.73
  [1330]     chrX [139592001, 139593238]      *   | MACS_peak_1330    210.88
  [1331]     chrY [ 13845133,  13845777]      *   | MACS_peak_1331     58.39
  ---
  seqlengths:
    chr1 chr10 chr11 chr12 chr13 chr14 ...  chr6  chr7  chr8  chr9  chrX  chrY
      NA    NA    NA    NA    NA    NA ...    NA    NA    NA    NA    NA    NA
> plotChrCov(peak, weightCol="V5")

chrCoverage

Heatmap of ChIP binding to TSS regions

?View Code RSPLUS
1
2
3
require(TxDb.Hsapiens.UCSC.hg19.knownGene)
txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
peakHeatmap(files, weightCol="V5", TranscriptDb=txdb, upstream=3000, downstream=3000, color=rainbow(length(files)))

peakHeat2

Average Profile of ChIP peaks binding to TSS region

?View Code RSPLUS
1
plotAvgProf(files, TranscriptDb=txdb, weightCol="V5", upstream=3000, downstream=3000)

avgProf2

Genomic Annotation

?View Code RSPLUS
1
2
peakAnnoList=lapply(files, annotatePeak)
plotAnnoPie(peakAnnoList[[1]])

pieAnno

?View Code RSPLUS
1
plotAnnoBar(peakAnnoList)

barAnno

Distance to TSS

?View Code RSPLUS
1
plotDistToTSS(peakAnnoList)

distTSS

Overlap of peak sets

?View Code RSPLUS
1
vennplot(peakAnnoList)

venn

In the future version, ChIPseeker will support statistical comparison among ChIP peak sets, and incorporate open access database GEO for users to compare their own dataset to those deposited in database. Signifciant overlap among peak sets can be used to infer cooperative regulation. This feature will soon be available.

Related Posts

To leave a comment for the author, please follow the link and comment on his blog: YGC » R.

R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...



If you got this far, why not subscribe for updates from the site? Choose your flavor: e-mail, twitter, RSS, or facebook...

Comments are closed.