averagepeaksignal-methods: Retrieve average peak profiles
In ChIPQC: Quality metrics for ChIPseq data

Description Methods Author(s) See Also Examples

Retrieve the average peak profile for a sample or set of samples.

signature(object = "ChIPQCexperiment"): Retrieve a matrix of the average peak profiles for all of the samples in an ChIP-seq experiment. Each column represents a sample, and each row a base pair position, centered on peak summits.
signature(object = "list"): Retrieve a matrix of the average peak profiles for all of the samples in a list of ChIPQC samples. Each column represents a sample, and each row a base pair position, centered on peak summits.
signature(object = "ChIPQCsample"): Retrieve a vector representing the average peak profile for a sample. Each column represents a basepair position, centered on the peak summits.

Thomas Carroll and Rory Stark

ChIPQC-package, ChIPQCexperiment, ChIPQCsample

data(example_QCexperiment)
CTCFprofile = averagepeaksignal(QCsample(exampleExp,1))
length(CTCFprofile)
plot(CTCFprofile,type='l',ylab="mean pileup")

allprofiles = averagepeaksignal(exampleExp)
dim(allprofiles)
for(i in 1:ncol(allprofiles)) lines(allprofiles[,i],col=i)