# NOT RUN {
library(org.Hs.eg.db)
data(rna)
cluster_which<-cbind(colnames(rna),c(rep(1,times=24),rep(2,times=24)))
me_x=rna
## calculate best number of clusters and
res<-top_supervised(me_x,TOP = 100,cluster_which)
me_TOP=res[[1]]
number_of_k=2
groups_men=res[[2]]
me_x=me_TOP
colnames(me_x)
o_g<-Supervised_Cluster_Heatmap(groups_men = groups_men, gene_matrix=me_x,
method="PAMR",show_sil=TRUE,print_genes=TRUE,threshold = 0,
TOP = 100,GSE=FALSE,plot_mean_sil=FALSE,stats_clust=res[[2]],
samples_data = as.data.frame(groups_men[,1,drop=FALSE]))
# }
Run the code above in your browser using DataLab