ChIPQC (version 1.8.2)

averagepeaksignal-methods: Retrieve average peak profiles

Description

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

Methods

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 = "ChIPQCsample")

Retrieve a vector representing the average peak profile for a sample. Each column represents a basepair position, centered on the peak summits.

See Also

ChIPQC-package, ChIPQCexperiment, ChIPQCsample

Examples

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)