similar to: Questions about 'heatmap.2' function

I have written a wrapper for heatmap.2 called heatmap.w.row.and.col.clust which auto-generates breaks using breaks<-round((c(seq(from=(-20 * stddev), to=(20 * stddev))))/20, digits = 2) #(stddev in this case = 2.5) This has always worked well in the past but now I am getting an error that non-finite breaks are being generated. Drilling down, it seems that my wrapper is generating finite

Hello again, lately I have been playing a lot with heatmap functions, in particular I complemented heatmap.2 with a nice feature from heatmap.plus, i.e. the possibility of assigning a matrix to ColSideColors so to take into account more annotations. Below is the code I introduced in heatmap.2, everything else is maintained equal to heatmap.2: if (!missing(ColSideColors)) { #if

Hi Allen and list, See the code below.? I've tried it on R2.13 and R2.8.0 using either heatmap.plus 1.3 or the latest.? All gave the same results.? The problem is in the last line: when I tried to plot two different color bars, the one corresponding to "cm.colors(10)" is not correct (it starts with one black and one red.? Not sure where they're from?) Any ideas? Thanks!

Dear R-Helpers, I have a short question related to heatmap.2 function. I wanted to change the dimension of each cell in the heatmap drawn by this function. here is a reproducible example from the help page library(gplots) data(mtcars) x <- as.matrix(mtcars) rc <- rainbow(nrow(x), start=0, end=.3) cc <- rainbow(ncol(x), start=0, end=.3) hv <- heatmap.2(x, col=cm.colors(255),

Hello, I have come across some unexpected behaviour of the function heatmap in the stats package. This looks like a bug to me, but I might have misunderstood something. When calling the function in symmetric mode, the ColSideColors are plotted correctly, but RowSideColors appear in reverse order. This code (modified from the example on the help page) demonstrates the problem: cU <-

Can someone recommend a more sophisticated way to annotate heatmaps than the ColSideColors argument of heatmap and heatmap.2? In particular, I would like to be able to annotate columns with more than one piece of information, like in Figure 1 of the article at http://www3.interscience.wiley.com/cgi-bin/fulltext/117905619/HTMLSTART / doi:10.1002/hep.22256. Some example data and a heatmap:

I am using heatmap with the arguments below. The title size stays the same no matter what I set cex.main to. Is this expected? Can I adjust the title size in heatmap? Also, the position of the main title is at the very upper edge of the output and if I use a "\n" to stack the title the upper line is out of bounds and doesn't show up. I am outputting to pdf. Any help? Thanks,

Hi, I wanted to make more than one side color bar. For example, I can make one side color bar (col1) with the following code: ----------------------------------- library(gplots) mat <- matrix(sample(1:100,40),nrow=5) class1 <- c(rep(0,4),rep(1,4)) col1 <- ifelse(class1 == 0,"blue","red") class2 <- c(rep(1,3),rep(2,5)) col2 <- ifelse(class2 ==

I am using attributes colsep and sepwidth from heatmap.2 function, however when using larger values for sepwidth the following columns are not moved to the right but the seperator is covering the column(s) with a white bar and the column labels and ColSideColors remain at their positions, is this intended? I expected the following columns to be shifted to the right by sepwidth. Is there a way

Dear list members, I have been using R to create a heatmap where my data has continous variables from 0 to 100. When I create the heatmap, although the branches are correct, they do not order themselves so that the row with the most zeros is at one end and the row with the most 100s is at the other, which is what I would like them to do, so as to create a colour gradient down the graphic. I have

Hi, I want to generate a heatmap for my data (in a matrix). However, the data has some missing values (represented as blank). I get the following errors (with the blanks and with blanks replaced by NA and including the option rm.na = TURE): > filename = "input_heatmap.txt" > g <- as.matrix(filedata) > fg <- rainbow(nrow(g), start=0, end=.3) > gg <-

Dear list, I have a question on moving color keys when side color bars are added to a heatmap. The R code below produces the color key in the upper left corner. Notice I have added side bars to the heatmap, but how could I move the color key below the image? ------ library(gplots) data(mtcars) x <- as.matrix(mtcars) rc <- rainbow(nrow(x), start=0, end=.3) cc <- rainbow(ncol(x),

I am in a situation needing 3 ColSideColors for heatmap.2. How can I do that? TIA Daren

Hello, I have a large data set from RNA sequencing and I am trying to make a heat map of my data. I have am having issues formatting my heat map figure. My data set is large with the log2 fold change for over 6oo genes across 4 treatments. My csv file is formatted as such: Gene Drought Ozone Temp1 Temp2 Glyma0041s00260

Hi, I have a question about the heatmap dendrogram in R. I loaded my data matrix in command heatmap(), and it gave me a heatmap accordingly, and the rows (genes) also were clustered accordingly. But now I don't want the genes clustered in that way, I have a new order of these genes, and want the rows of the heatmap are drawed as the order I give. Could anybody help me on that? Thank you!

Hi all, Does anyone know of a fairly easy way to "stretch" a heatmap vertically? I've got 42 arrays and would like to be able to see as many significant genes as possible (right now I can only get 50 genes with it still being readable). In some comparisons there are several hundred significant genes. I've fiddled with the "asp" argument, but that doesn't

Hello! I looked for help through google and the help-files and spend several hours with trial and error, but didn't find a correct way. It's all about lines in a heatmap to separate different data block from each other to underline the significance of the found clusters! The heatmap is build like that: heatmap(X, Rowv=NA, Colv=NA, symm=TRUE, cexRow=0.3, cexCol=0.3) I switched the

Hi All Not sure if this a bioconductor question or general R mailing list so apologies if this has gone to the wrong one................. When plotting dendrograms created by hclust you can "identify" clusters by clicking on the graphics and returning a list of what is contained in each cluster. However I'd like to be able to "zoom in" on specific clusters and plot

Hello, I've been clustering my data using hclust and cutting the resulting tree with cutree. Separately, I visualize the clusterings with heatmap. Is it possible to have the dendrogram on the heatmap reflect the cutree results? That is, instead of having one large dendrogram, it would have 4 or 25 in the example below. Any guidance on if that's possible or not, and what kinds of

Hi, I am using the function heatmap(stats) to draw a microarray heatmap, columns are samples and rows are gene features. I did a 2D clustering during the heatmap drawing. The features and samples indeed cluster into several blocks both vertically and horizontally. I can get the index of re-ordered rows and columns after the heatmap drawing by typing the the return variable of the heatmap