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