Hello All,
Is there a way to extract the plots separately from the ConsensusClusterPlus package? For example, using the example data from the package, I can print the last three plots as below,
library(ALL)
data(ALL)
d=exprs(ALL)
mads=apply(d,1,mad)
d=d[rev(order(mads))[1:5000],]
d = sweep(d,1, apply(d,1,median,na.rm=T))
library(ConsensusClusterPlus)
par(mfrow=c(1,3))
title=tempdir()
results = ConsensusClusterPlus(d,maxK=6,reps=50,pItem=0.8,pFeature=1,
title=title,clusterAlg="hc",distance="pearson",seed=1262118388.71279)
But, I wish to extract the consensus heatmap (say plot #2 which represents the optimal number of clusters) from the results and display separately.
Any help in this regard would be very helpful.
Thanks in advance,
M
Read more here: Source link