similar to: labRow/labCol options in heatmap()

Dear All, There might be a bug in heatmap(): x <- matrix(rnorm(10*5), 10, 5) heatmap(x, labCol=1:5) Error in axis(1, 1L:nc, labels = labCol, las = 2, line = -0.5, tick = 0, : 'at' and 'labels' lengths differ, 5 != 10 Works fine under 4.0.1. Looking at the code from 4.0.1 vs. Rdevel, esp. this part sticks out: labRow <- labRow[rowInd] %||% rownames(x) %||%

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

I'm trying to display the probability space of a function and wanted to see specifically where the maximum posterior probability is. data=read.table("PosteriorData",header=F) mat=as.matrix(data) heatmap.2(mat/max(mat,na.rm=T),dendrogram=c("none"),trace=c("none"), Rowv=F,Colv=F,labRow=3*c(10:-10),labCol=3*c(-10:10),symm=T,col=rainbow(100)) As the

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,

Full_Name: Anja von Heydebreck Version: 1.3.0 OS: Alpha Unix Submission from: (NULL) (141.14.19.61) Hi, I repeatedly obtained meaningless results from the function 'cutree' in the 'mva' package, when the argument 'h' was greater or equal to the maximum height occuring: > library('mva') > y [,1] [,2] [,3] [,4] [1,] 0 1 -1 1 [2,] 0 -1

Hi, I am using heatmap.2 of gplots to make heatmaps of my the attached file. I am giving my code for the same .. library(gplots) x=read.table("1.txt", header=TRUE) mat=data.matrix(x) heatmap.2(mat, col=greenred(75), Rowv=TRUE, Colv=TRUE, distfun = dist, hclustfun = hclust, dendrogram = c("both"), scale = c("row"), na.rm=TRUE, trace="none",

I'm suppressing the labeling of my rows and columns in heatmap.2 using the commands: labRow = " ", labCol = " " But I'd like to annotate them again using the axis command: axis(1, at=seq(500, 1000, 500)) mtext("Group 1", "Group 2") For some reason however it appears that the axis command is having no effect. My complete

Full_Name: Anja von Heydebreck Version: 1.9.1 OS: Linux Submission from: (NULL) (155.250.128.25) I obtained a negative p-value from the fisher.test() function: > fisher.test(matrix(c(14576,3023,89,68),2), alternative="g")$p.value [1] -8.426593e-13 With R 1.8.1, I got a p-value of -6.239231e-12 for this example. Anja von Heydebreck

Dear R users, I'm trying to create images of 2 symmetrical matrices using heatmap.2 of the 'gplots' library. Both matrices have the same row and column names For the first matrix I'm using the default clustering and column denrogram options: heatmap.2(Matrix1, symm=T, breaks=20, col=rich.colors(19), main="matrix1", density.info="histogram",

Dear Experts, I am using the below script to generate the heat map of gene expression data. I am using Hierarchical Clustering (hclust) for clustering. Now I want to compare different clustering parameters such as *K-means* clustering, Model Based Clustering, I have two queries: 1. How to incorporate different clustering method in the same code? 2. Is this possible to implement pvclust in the

Dear R users, need help with my heatmap. I will really approciate some help. Given the matrix: > head(x) A B C D time [1,] 0 8 0 0 1 [2,] 0 160 0 0 2 [3,] 0 175 0 0 3 [4,] 0 253 0 0 4 [5,] 79 212 0 0 5 [6,] 6 105 0 0 6 and call: ##------------ Heatmap -------------- library(RColorBrewer) rc <- grey((nrow(x):1)/nrow(x)) heatmap(data.matrix(x[,1:4]),

Hi, I am drawing a heat map using gplots of R package. The file (new4) is having values 0 and 1 only. library(gplots) library(marray) x=read.table("new4", header=TRUE) mat=data.matrix(x) heatmap.2(mat, col=colorRampPalette(c("black","red"))(256), #col=greenred(75), Rowv=TRUE, Colv=TRUE, distfun = dist, hclustfun = hclust, dendrogram = c("row"), scale =

Hi, Can anyone tell me how to set Z-score according to my own requirement as the below code is taking as per the file entries. Any help would be appreciable. library(gplots) x=read.table("final.txt", header=TRUE) mat=data.matrix(x) heatmap.2(mat, col=colorRampPalette(c("green","white","red"))(256), #col=greenred(75), Rowv=TRUE, Colv=FALSE, distfun = dist,

I would like to use gridBase to place four separate heatmaps (actually, a stripped-down heatmap.2 from ght gregmisc package that contains only the "image" part) into four different viewports. I can get the placement correct, but I keep 'losing' the previous plot. Any suggestions? Here is some quick example code trying to put a heatmap into the left viewport and then

Dear R Community, I am trying to create an heatmap for the following set of data: ##example of data matrix o4 V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 green 27 28 29 29 28 28 26 25 25 23 23 22 22 21 21 22 22 22 yellow 6 8 8 7 7 7 6 6 6 7 7 7 6 6 6 6 6 5 red 15 15 15 15 15 15 14 13 12 11 12 10 9 8 7 6 8 9 pink

Hi, I have searched for a solution but I failed to find an answer. I am hoping you may be able to help me. I have a data set where I have observations for a number of units (n =~40) over a period of time (t =~100) and I have a variable (Z) that codes a categorical variable for each observation. I want to produce a 2D plot where time is on the x-axis and units are on the y-axis. Then each block

Hello R Gurus, I have a matrix for which I am doing a heatmap using heatmap.2. I want to put the rownames on the lefthand side instead of the right side of the heatmap. how can i put the rownames on left hand side: I have already tried axis but could not make it work ccc<-structure(c(1, 0.283300333755851, 0.237863231117007, 0.0148696794159253, -0.0780756406815149, -0.106161465097659,

Hi, Is there a way I can draw an abline on a heatmap? I try the abline function, but don't get the line. My sample code is: mat <- matrix(sample(1:100,40),nrow=5) heatmap(mat,col=greenred(75),trace="none", dendrogram = "column",labCol = NULL) abline(h=5,v=4) thanks! [[alternative HTML version deleted]]

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 ==

Dear All: I tried to use heatmap.2 to generate hierarchical clustering using the following command: heatmap.2(datamatrix, scale="row", trace="none", col=greenred(256), labRow=genelist[,1], margins=c(10,10), Rowv=TRUE, Colv=TRUE) datamatrix is subset of a RMA normalized data subset by a genelist. The problem is a lot of times, the z-score in key are from, like -5 to 15 or