Description Usage Arguments Value References See Also Examples
After fiting the CENTIPEDE model this function is useful for visualizing the footprint cleavage pattern.
1 2 | imageCutSites(data, main="", xlab = "Dist. to motif (bp)")
imageCutSitesCombined(data, main="",xlab= "Dist. to motif (bp)",breaks)
|
data |
A matrix with the data supplied to CENTIPEDE |
xlab |
Label of the x axis |
main |
Title for the plot |
breaks |
Breaks for the color assigments |
No return value, just the plot
Roger Pique-Regi, Jacob F. Degner, Athma A. Pai, Daniel J. Gaffney, Yoav Gilad, Jonathan K. Pritchard. "Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data", Genome Research, Submitted Aug 2010
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ## Not run:
#GETS EXAMPLE DATA FOR NRSF
data(NRSFcuts, package='CENTIPEDE')
data(NRSF_Anno, package='CENTIPEDE')
#FITS THE CENTIPEDE MODEL
centFit <- fitCentipede(Xlist = list(as.matrix(NRSFcuts)), Y=cbind(rep(1, dim(NRSF_Anno)[1]), NRSF_Anno[,5], NRSF_Anno[,6]))
#PLOTS IMAGE OF CUTSITES RANKED BY CENTIPEDE POSTERIORS
imageCutSites(NRSFcuts[order(centFit$PostPr),][c(1:100, (dim(NRSFcuts)[1]-100):(dim(NRSFcuts)[1])),])
#PLOT ESTIMATED FOOTPRINT
plotProfile(centFit$LambdaParList[[1]],Mlen=21)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.