Mercurial > repos > stef > qdnaseq
diff QDNAseq-plot.R @ 88:d2ea2b842c21 draft default tip
Uploaded
author | stef |
---|---|
date | Mon, 06 Jul 2015 06:29:19 -0400 |
parents | 8391cdb7479c |
children |
line wrap: on
line diff
--- a/QDNAseq-plot.R Thu May 28 14:19:57 2015 -0400 +++ b/QDNAseq-plot.R Mon Jul 06 06:29:19 2015 -0400 @@ -51,13 +51,13 @@ if ( whatToPlot == 'everything' ){ catMsg( c( "Plotting all data in object" ) ) plot( qdnaseqObject[ ,sample ], ylab=ylab_text ) - abline( h=c(-2,-1,1,2,3,4), lty=1, lwd=0.5, col="grey" ) + abline( h=c(-4,-3,-2,-1,1,2,3,4), lty=1, lwd=0.5, col="grey" ) } else if( whatToPlot == 'chromosomes' ){ catMsg( c( "Plotting subset of chromosomes" ) ) fdata <- qdnaseqObject@featureData@data idx_region <- which( fdata$chromosome %in% chromosomesToPlot ) plot( qdnaseqObject[ idx_region, sample ], ylab=ylab_text ) - abline( h=c(-2,-1,1,2,3,4), lty=1, lwd=0.5, col="grey" ) + abline( h=c(-4,-3,-2,-1,1,2,3,4), lty=1, lwd=0.5, col="grey" ) } else if( whatToPlot == 'region' ){ regionC <- chrName regionS <- chrStart