# NOT RUN {
## using TROM to map developmental stages of D. melanogaster and C. elegans
## The .rda files used in this example can be downloaded and unzipped from
## http://www.stat.ucla.edu/~jingyi.li/packages/TROM/TROM_Rdata.zip.
# }
# NOT RUN {
load("dm_gene_expr.rda")
load("ce_gene_expr.rda")
load("dm_ce_orthologs.rda")
dm_ce_trom <- bs.trom(sp1_gene_expr = dm_gene_expr,
sp2_gene_expr = ce_gene_expr,
sp1_sp2_orthologs = dm_ce_orthologs, z_thre=1.5,
provide=FALSE)
heatmap.3( dm_ce_trom,
max_score = 6,
Rowv = NULL,
Colv = NULL,
dendrogram = c("none"),
distfun = dist,
hclustfun = hclust,
xlab = "worm stages",
ylab = "fly stages",
main = "D. melanogaster vs. C. elegans Stage Mapping",
key = TRUE,
keysize = 1,
trace = "none",
density.info = c("none"),
col = terrain.colors(120),
)
# }
Run the code above in your browser using DataLab