maxK=9 #最大的k值(最多可以将样品分成几个亚型)
data=t(data[,sigGenes])
results=ConsensusClusterPlus(data,
maxK=maxK,
reps=50,
pItem=0.8,
pFeature=1,
title=workDir,
clusterAlg="km",
distance="euclidean",
seed=123456,
plot="png")
运行代码出现