GenomeGraphs plot colours (bug?)
0
0
Entering edit mode
Tim Rayner ▴ 270
@tim-rayner-2913
Last seen 10.2 years ago
Hi, I'm using the GenomeGraphs package to plot expression data for multiple sample groups, and I've run into something that may or may not be classed as a bug. I've been using DisplayPars(color=c('blue','red')) to visually differentiate between sample groups, and that works well under most circumstances. However, I've found that when using e.g. makeGenericArray() that the addition of a probeEnd argument renders the graph using only the first colour in the DisplayPars list. Example code and sessionInfo() attached below. I only noticed this because I happened to try using DisplayPars(type="line") with probeStart and probeEnd to indicate the length of the interrogated genome region on the plots. Aside from this minor niggle, nice package! Cheers, Tim Rayner ## Example code follows: testGenomePlot <- function () { require('GenomeGraphs') data("exampleData", package = "GenomeGraphs") minBase <- 180292097 maxBase <- 180492096 ## Replace the default intensity data for example's sake (we need two samples) intensity <- matrix(rnorm(184), ncol=2) expres <- makeGenericArray(intensity = intensity, probeStart = exonProbePos, dp = DisplayPars( color = c("blue","red"), type = "point", lwd=2)) expres2 <- makeGenericArray(intensity = intensity, probeStart = exonProbePos, probeEnd = exonProbePos+200, # different dp = DisplayPars( color = c("blue","red"), type = "point", lwd=2)) gdPlot(list(good=expres, bad=expres2), minBase = minBase, maxBase = maxBase) } testGenomePlot() sessionInfo() R version 2.8.1 Patched (2009-01-19 r47650) i386-apple-darwin9.6.0 locale: en_GB.UTF-8/en_GB.UTF-8/C/C/en_GB.UTF-8/en_GB.UTF-8 attached base packages: [1] grid stats graphics grDevices utils datasets methods [8] base other attached packages: [1] GenomeGraphs_1.2.3 biomaRt_1.16.0 loaded via a namespace (and not attached): [1] RCurl_0.94-1 XML_1.99-0 tools_2.8.1
graph GenomeGraphs graph GenomeGraphs • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 723 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6