## Not run:
# library(cluster)
# data(subset.loring.eset)
# clustObj <- agnes(as.dist(1-t(cor(exprs(subset.loring.eset)))))
# cfmod.vals <- NULL
# for( i in 1:10 ){
# cfmod.vals <- c(cfmod.vals, calcModfstat(exprs(subset.loring.eset), cutree(clustObj,i), pData(subset.loring.eset)$celltype))
# }
# k <- (1:10)[cfmod.vals==max(cfmod.vals)]
# ## End(Not run)
Run the code above in your browser using DataLab