[BioC] small bug in overview function of made4

Jenny Drnevich drnevich at illinois.edu
Wed Nov 11 21:23:11 CET 2009


Hi Aedin,

I noticed a change in the behavior of your overview function in the 
made4 package from BioC 2.4 to BioC 2.5. Before, when I only wanted 
to plot the hc dendrogram, it completely filled the plotting window. 
However, now no matter what combinations of hc, boxplot and hist, the 
window is always divided up into 3 sections. I looked at the code, 
and there seems to be one change that's causing the bug:

 >overview
function (dataset, labels = NULL, title = "", classvec = NULL,
     hc = TRUE, boxplot = TRUE, hist = TRUE, returnTree = FALSE)
{
     cols = NULL
     if (sum(boxplot, hist, hc) == 3)
         layout(matrix(c(1, 1, 2, 3), 2, 2, byrow = TRUE))
     if (sum(boxplot, hist, hc) == 2)
         layout(matrix(c(1, 1, 2, 2), 2, 2, byrow = TRUE))
     classvec = as.factor(classvec)

This last line is new, and changes the default NULL value into a 
factor, such that a few lines farther down in the code:

if (!is.null(classvec)) {
         cols = getcol(as.numeric(classvec))
         layout(matrix(c(1, 1, 2, 2, 3, 4), 3, 2, byrow = TRUE),
             heights = c(1.5, 0.1, 2))
     }

!is.null(classvec) is now TRUE, even if classvec was originally NULL.


And since I've e-mailed you, I'll say that I always thought it was 
weird that the main title for the hierarchical clustering dendrogram 
was "Histogram", not "Dendrogram":

plot(hc, hang = -1, labels = labels, main = paste("Histogram",
             title, sep = " "), sub = "", xlab = "")

Thanks!
Jenny

 > sessionInfo()
R version 2.10.0 (2009-10-26)
i386-pc-mingw32

locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252

attached base packages:
[1] grid      stats     graphics  grDevices datasets  utils     methods
[8] base

other attached packages:
  [1] limma_3.2.1          affyQCReport_1.24.0  lattice_0.17-26
  [4] xtable_1.5-5         simpleaffy_2.22.0    genefilter_1.28.0
  [7] made4_1.20.0         scatterplot3d_0.3-29 gplots_2.7.3
[10] caTools_1.10         bitops_1.0-4.1       gdata_2.6.1
[13] gtools_2.6.1         RColorBrewer_1.0-2   ade4_1.4-11
[16] affyPLM_1.22.0       preprocessCore_1.8.0 gcrma_2.18.0
[19] affycoretools_1.18.0 KEGG.db_2.3.5        GO.db_2.3.5
[22] RSQLite_0.7-3        DBI_0.2-4            AnnotationDbi_1.8.0
[25] affy_1.24.0          Biobase_2.6.0        RWinEdt_1.8-1

loaded via a namespace (and not attached):
  [1] affyio_1.14.0     annaffy_1.18.0    annotate_1.24.0   biomaRt_2.2.0
  [5] Biostrings_2.14.0 Category_2.12.0   GOstats_2.12.0    graph_1.24.0
  [9] GSEABase_1.8.0    IRanges_1.4.0     RBGL_1.20.0       RCurl_1.2-1
[13] splines_2.10.0    survival_2.35-7   tools_2.10.0      XML_2.6-0






Jenny Drnevich, Ph.D.

Functional Genomics Bioinformatics Specialist
W.M. Keck Center for Comparative and Functional Genomics
Roy J. Carver Biotechnology Center
University of Illinois, Urbana-Champaign

330 ERML
1201 W. Gregory Dr.
Urbana, IL 61801
USA

ph: 217-244-7355
fax: 217-265-5066
e-mail: drnevich at illinois.edu



More information about the Bioconductor mailing list