library(ConsensusClusterPlus)res <- ConsensusClusterPlus(es_max, # es_max为矩阵,行名为feature,列为samplemaxK = 10,reps = 1000, # 抽样次数pItem = 0.8, # 抽样比例pFeature = 1,clusterAlg = "pam", # 聚类方法corUse = "complete.obs",seed=456,plot="png", # 结果保存形式writeTable=T,innerLinkage = "complete")
遇到过的一些报错:
es_max一定要是矩阵
es_max里不可以有标注差为0的列 es_max = es_max[,apply(es_max, 2, function(x) sd(x)!=0)]
