library(ConsensusClusterPlus)
res <- ConsensusClusterPlus(es_max, # es_max为矩阵,行名为feature,列为sample
maxK = 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)]