Making Heatmaps: Difference between revisions

From Powers Wiki
No edit summary
No edit summary
Line 52: Line 52:
  **the script file will be opened in the R Editor window
  **the script file will be opened in the R Editor window
  ***If necessary changes can be made on the R Editor window  
  ***If necessary changes can be made on the R Editor window  
  ===To run the script===
  ===To run the script===
* click Run all or Run line or selection form the Edit menu
* click Run all or Run line or selection form the Edit menu

Revision as of 05:05, 4 December 2012

Load the data

The data can be prepared in csv or text file to load the data in to R. For example, to load list.csv:

heatmapname <- read.csv (“list.csv”)

Name the rows

You can then name the rows of the heatmap like so:

row.names(heatmapname) <- data$Name

To exclude the first column from the heat map, use a command of similar form to the following:

heatmapname <- heatmapname [,2: n]

In the above command, n is the number of columns to be included in the heat map.

Build a data matrix

Use the following command to build a data matrix for making the heat map:

heatmapname_matrix <- data.matrix (heatmapname)

Plot the heat map

Note: Gnuplot package must be installed in R before heat maps may be displayed! Run the following command to load in the gplots library.

library("gplots")

To make heat map, run the following command:

heatmap.2 (heatmapname_matrix, dendrogram="row", col= redgreen (75), scale="none",
key = TRUE, keysize = 1.0, margins = c(4,30),
density.info="none", trace="none")

Note: To rescale the color key, add break function to the heatmap.2.

breaks=c(seq(-1,0.8,length=10),seq(0.8,1.2,length=10),seq(1.2,3,length=10),

This will define the color range:

"red=[-1,0.8]

black=[0.8,1.2]

green=[1.2,3]"

It is also possible to write all the steps in txt file as a single script.

open the script

*In the file menu click open script
**choose the proper directory and script file name
**the script file will be opened in the R Editor window
***If necessary changes can be made on the R Editor window 
===To run the script===
  • click Run all or Run line or selection form the Edit menu