Matt:LabNotes/2015-12-11

From ZhangLabWiki
Revision as of 20:26, 11 December 2015 by >Mzcai (Created page with " load.packages("gplots") library(gplots) tileA <- read.table("C:\\Users\\Matt\\Dropbox\\GradZhangLab\\SpatialSearching\\20151029_20XDecoding\\DARTFISH_BA8_suppv2_20151029_Ti...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search


load.packages("gplots") library(gplots)

tileA <- read.table("C:\\Users\\Matt\\Dropbox\\GradZhangLab\\SpatialSearching\\20151029_20XDecoding\\DARTFISH_BA8_suppv2_20151029_TileA_240genes21bins.txt",header=TRUE,sep="\t",row.names=1) tileB <- read.table("C:\\Users\\Matt\\Dropbox\\GradZhangLab\\SpatialSearching\\20151029_20XDecoding\\DARTFISH_BA8_suppv2_20151029_TileB_240genes21bins.txt",header=TRUE,sep="\t",row.names=1) y <- as.matrix(tileA) z <- as.matrix(tileB)


my_palette <- colorRampPalette(c("turquoise","yellow","red"))(n=299) col_breaks=c(seq(-4,-1,length=100),seq(-1,1,length=100),seq(1,4,length=100))

heatmap.2(y,density.info="none",trace="none",col=my_palette,breaks=col_breaks)