Plot_cluster.R

Plotting_clustering_timeseries.R

par(mfrows=c(clustering$cluster,1))

for(j in 1:max(clustering$cluster))
{
  x <- 1:12
  d <- d[clustering$cluster == j,]
  plot.default(x = NA, xlim = c(1, 12), 
             ylim = c(min(y), max(y)), xlab = "Time", ylab = "Expression changes", 
             main = paste("Cluster", j), axes = FALSE)

  for(k in 1:nrows(d)) {
    lines(x, d[k, ])
  }
}