Question: How to extract plots from ConsensusClusterPlus package
0
gravatar for manali.rupji
3 days 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

ADD COMMENTlink written 3 days 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 2 days ago by arta40
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: 957 users visited in the last hour