note C01.R # 31/95

cell.sdata
> pca <- princomp(cell.sdata)
> pca.dim1 <- pca$scores[,1]
> pca.dim2 <- pca$scores[,2]
> plot(pca.dim1, pca.dim2,  main="PCA for Cell Cycle Data with K-means Clustering", xlab="PCA-1", ylab="PCA-2", col=cell.kmeans$cluster)
> 

C02.R
# 25/109
no pca.dim1
 


cell.sdata <- scale(cell.data)
GBRcol <- color.Palette(low="green", mid="black", high="red")

hv <- heatmap(cell.sdata, col = GBRcol)

rc <- rainbow(5)[as.integer(gene.phase)]
cc <- rainbow(ncol(cell.sdata))




hv <- heatmap(cell.sdata, col = GBRcol, scale = "column", Colv=NA, Rowv=NULL,            
margins = c(5,10),
              xlab = "Times", ylab =  "Genes",
              main = "Heatmap of Microarray Data")