similar to: heatmap.2() problems with re-ordering of rows and columns

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

Hi STATISTICA has a function called "Two-way joining" (see http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the reference material states that this is based on the method as published by Hartigan (found this paper: http://www.jstor.org/pss/2284710 through wikipedia). What is the relationship (if any) between the "heatmap" function in R and this technique? Is there an

Hi I am using the heatmap function in package mva to look at large correlation matrices visually. Is there any way to "flip" the output of heatmap plot left-right so that, if presented with a correlation matrix, it plots the unity elements in the correlation matrix along a diagonal from top left to bottom right? For example: library(mva) x = matrix(rnorm(1000), ncol=10) z = cor(x)

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

Dear R community! I am trying to create a heatmap based on the following data. As you can see the diagonal (0,0 to 10,10) is always 0). If I run the heatmap command like i posted it I get a graph with horizontal lines. What is my mistake? > heatmap(activity.matrix, Rowv = NA, Colv = NA, scale="column", xlab="x Compounds", ylab="y Compounds",

I have a bunch of lines I want to plot using plotCI() What Id like to know is, how can I connect the points with a line and how can I print multiple lines on the same graph? -- View this message in context: http://www.nabble.com/plotCI----multiple-plots-on-same-graph-tp19435198p19435198.html Sent from the R help mailing list archive at Nabble.com.

I want to make sure this piece of code I wrote is doing what I want it to do. ll<-function(string) { grep(string,dir(),value=T) } subdir = ll("Coverage_[1-9][0-9]$") I basically wrote a little function that would grab all the files of form Coverage_[0-99] The way I wrote it, will it grab Coverage_5 or does it have to have 2 numbers (10-99)? -- View this message in context:

I have a huge dataset(15k X 18) and tried to use the heatmap in R to examine the patterns. However, I found that heatmap and heatmap.2 gave me completely different outputs. Here are the codes: ------------ > dim(as.matrix(data.dcpm)) [1] 15462??? 18 > > heatmap(as.matrix(data.dcpm), col=topo.colors(100)) > heatmap.2(as.matrix(data.dcpm), col=topo.colors(100), key=TRUE,

I am having trouble displaying a dendrogram of evolutionary relationships (a phylogram imported from the ape package) as the vertical component of a heatmap, but keeping the hierarchical clustering of the horizontal component. The relationships of the vertical component in the generated heatmap are not that of the dendrogram, although the ordering is. In more detail, I am attempting to generate

Hi, I'm kinda new of the R world. I need to use the heatmap or heatmap.2 function to plot correlation values. 2 questions: - how can I specify a color palette? (would like to have a transition from blue(-1) to white(0) to red(1)) - how can I use heatmap.2 imposing a specific order (like hclust.order)? I mean, I don't want heatmap.2 to re-cluster/re-order my data but just to plot them

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 Sorry to revisit an old problem, I seemed to solve this in 2004, only for it to resurface :-S I am trying to plot a heatmap, and I don't want the columns of my matrix re-ordered. The function doesn't seem to behave as the help would have you believe: a <- matrix(rnorm(100),nr=20) a.d <- dist(a) a.hc <- hclust(a.d) a.de <- as.dendrogram(a.hc) # columns are re-ordered

Hi, I was trying to use par(mfrow) to put 4 heatmaps on a single page. However, I get one plot per page and not one page with 4 plots. What should I modify? Test code is given below: test = matrix(rnorm(60), 20, 3) pdf(file='test.pdf',width=10,height=8) par(mfrow=c(2,2)) heatmap(test) heatmap(test) heatmap(test) heatmap(test) dev.off() thanks! [[alternative HTML version deleted]]

Hi, I have what would seem to be a very simple problem but have been unable to find any way to go about solving it: I am using heatmap.2 to create a heatmap and by default the x axis label is at the bottom of the heatmap, I simply want to move the label to the top of the heatmap. Any help with this would be greatly appreciated. Thanks, Mark

Hi,   I need to generate a heatmap on a square matrix and wouldn't want to reorder the columns and the rows on the heatmap display.    I have used the options Rowv=NULL and Colv=NULL but doesn't seem to work. Following is a snippet of the heatmap function i am using. args <- commandArgs(); inputfile <- args[2] imgfile   <- args[3] bitmap(imgfile, height=15, width=15, res=100,

Howdy, I am trying to make a simple monochrome heatmap from 1 row of data. Essentially, I just want a long bar where black represents the max value in the data, white is the minimum, and all values in between are interpolated appropriately. I have tried using heatmap and heatmap.2, but both have issues. heatmap seems to be restricted to only drawing a square plot, and this data requires a

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

Hi I am having some problems getting my heatmap to be the right size! Let me explain. I am experienced at getting an hclust or a dendrogram object to be the right size. For example, I have a dataset which has 4000 rows, which I clustered using hclust and I wanted to plot it as a horizontal dendrogram. So I used jpeg(), set the image height to be 4000 and plotted the dendrogram and got the

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

Good afternoon, I ran heatmap and hclust on the same matrix x (strictly, I ran heatmap(x), and hclust(dist(t(x))), and realized that the two dendrograms were slightly different, in that the left-right arrangement of one pair of subclusters (columns) was reversed in the two functions (but all individual columns were grouped correctly). Looking through the code for heatmap as a most definite