## 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),
# )
# ## End(Not run)
Run the code above in your browser using DataLab