Question: How to extract plots from ConsensusClusterPlus package
0
gravatar for manali.rupji
11 months ago by
United States
manali.rupji0 wrote:

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

plots consensusclusterplus • 669 views
ADD COMMENTlink written 11 months ago by manali.rupji0

Hey,

Why do not use plot="pdf" in ConsensusClusterPlus function, so you will get 6 individual plots for each number of cluster, later you can take interest of yours.

results = ConsensusClusterPlus(d,maxK=6,reps=50,pItem=0.8,pFeature=1,
                           title=title,clusterAlg="hc",distance="pearson",seed=1262118388.71279, plot="pdf")
ADD REPLYlink written 11 months ago by arta540
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1733 users visited in the last hour