[R] How to shift the heat map

Muhammad Hafiz muhdhafizali at gmail.com
Thu Jun 9 10:19:34 CEST 2011


Hi, what i need is as the example i attached in the box plot. As you can
see, the boxplot have margins and it kind
    of follow the coding I've added, but, it does not happen for the heat
map. Could you help me with this.. Thanks in advance.
the following is the coding i've used:

# Draw the heatmap
#    heatmap.r
#
#    Purpose:    Create a heatmap
#
#    Input:    Data matrix as specified in Data-matrix-format.pdf
#    Output:    Heatmap with dendrograms showing clustering of both
#               samples and variables.


#    Load necessary libraries, and install them if they are missing
#
tryCatch(library(gplots), error=function(err)
    # if this produces an error:
    install.packages("gplots",repos="http://cran.ms.unimelb.edu.au/"))


#    Prepare the data matrix

par(omi=c(1,1,1,1))
# Read in the .csv file
data<-read.csv("MEAN PERAK RT100 300-399.csv", header = TRUE, row.names=1)
boxplot(data)
heatmap_data<- data[,-1]

colnames(heatmap_data) <- if
    (length(grep("^X[\\d]",colnames(heatmap_data),perl=TRUE)) != 0) # then
    {gsub("^X([\\d].*)","\\1",colnames(heatmap_data),perl=TRUE)} else
    {colnames(heatmap_data)}

# Define palette colours
rch<-rich.colors(256)
rwb<-colorRampPalette(c("red","white","blue"))(256)
rbb<-colorRampPalette(c("red","black","blue"))(256)
pwg<-colorRampPalette(c("purple","white","green"))(256)
pbg<-colorRampPalette(c("purple","black","green"))(256)

pic_onscr<-function(matrix, title="", cex_val=1, pal=rch)
    {x11()
    heatmap.2(matrix,                              # matrix to use
        #Colv=FALSE,                             # reorder dendrogram
        dendrogram="both",                         # dendrograms to draw
        col=pal,                                   # palette (see above)
        trace="none",                              # trace line (show
distance)
        cexCol=cex_val,                            # size of column labels
        cexRow=cex_val,                            # size of row labels
        key=TRUE,                                  # show the colour key
        keysize=cex_val,                           # key size
        density.info="none",                       # superimpose "histogram"
or a
                                                   # "density" plot on
colour key
        margins=c(10,10),                          # add space for labels
        main=title,                                # title
        xlab="",                                   # x-axis title
        ylab=""                                    # y-axis title
        )
    }
# Draw the heatmap
box("plot", col="red")
box("figure", col="green")
par(mar=c(3,2,2,1) + 0.1)
par(omi=c(1,1,1,1))
pic_onscr(t(heatmap_data))
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Heat map.png
Type: image/png
Size: 44268 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110609/e3379c3f/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: BOXPLOT.png
Type: image/png
Size: 13418 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110609/e3379c3f/attachment-0001.png>


More information about the R-help mailing list