## ------------------------------------------------------------------------
## Example1: CAGE
## ------------------------------------------------------------------------
require("GMD") # load library
data(cage) # load data
## construct a distance matrix and visualize it
short.labels <- gsub("(.+) \\(.+","\\1",names(cage)) # get short labels
x <- gmdm(cage[1:6],labels=short.labels[1:6])
plot(x)
## Not run:
# ## ------------------------------------------------------------------------
# ## Example2: ChIP-seq
# ## ------------------------------------------------------------------------
# data(chipseq_mES) # load data
# data(chipseq_hCD4T) # load data
#
# ## pairwise distance and alignment based on GMD metric
# plot(gmdm(chipseq_mES,sliding=FALSE))
#
# ## clustering on spatial distributions of histone modifications
# x <- gmdm(chipseq_hCD4T,sliding=FALSE,resolution=10)
# heatmap.3(x,revC=TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab