Description Usage Arguments Details References See Also Examples
This is a plot function suitable for visualisation of a bclustvs
object when the bclust
function is used on replicated data, but can be applied for unreplicated data too.
The function attaches a coloured horizontal dendrogram to the left side of a profile plot with the optimal grouping highlighted by a teethplot on the right.
1 2 3 4 5 6 7 8 9 | dptplot(x, scale = 1, xlab = colnames(x$data), ylab = x$labels,
xlab.cex = 1, ylab.cex = 1, dendrogram.lwd = 1,
dendrogram.size = 2, xlab.mar = 3, ylab.mar = 3,
horizbar.plot = FALSE,
horizbar.col = rev(c(heat.colors(5)[-4], "white")),
horizbar.distance = 4, varimp = rep(0, ncol(x$data)),
horizbar.size = 0.5, vertbar = NULL,
vertbar.col = rainbow(max(vertbar)),
teeth.size = 0.25, plot.width = 10)
|
x |
A |
scale |
A positive value. If the profile plots are too noisy or too flat, adjust the scaling factor. |
xlab |
A string vector, the labels for the clustering types. |
ylab |
A string vector, the variable labels. |
xlab.cex |
A positive value, the magnitude of the type labels. |
ylab.cex |
A positive value, the magnitude of the variable labels. |
dendrogram.lwd |
A positive value, the thickness of lines used to plot the dendrogram. |
dendrogram.size |
A positive value, the size of the dendrogram plot. |
xlab.mar |
A positive value, the margin reserved to write variable labels. |
ylab.mar |
A positive value, the margin reserved to write type labels. |
horizbar.plot |
A logical value. If |
horizbar.col |
Colours of the horizontal bar. |
horizbar.distance |
A positive value, the distance between the horizontal bar and the profile plot. |
varimp |
A numerical vector denoting the importance of variables. We propose to use |
horizbar.size |
A positive value, the size of the horizontal bar. |
vertbar |
A positive integer vector that may be used to draw an additional vertical bar on the right of the teethplot. This may be helpful to represent another optional grouping on the data. |
vertbar.col |
The colours used to plot the additional vertical bar. |
teeth.size |
A positive integer, the size of the teeth plot. |
plot.width |
A positive integer, the width of the whole plot. If the plot region is unbalanced in width and height, adjust this value. |
The varimp
is assumed to be the log Bayes factor and therefore categorised according to Kass and Raftery (1995) for a better visualisation. If varimp
is not
result of the imp
function, keep horizbar.plot=FALSE
.
Kass and Raftery (1995) Bayes Factors, Journal of the American Statistical Association, Vol. 90, pp. 773–795.
ditplot, profileplot, teethplot, viplot.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 | data(gaelle)
gaelle.id<-rep(1:14,c(3,rep(4,13)))
# first 3 rows replication of ColWT, 4 for the rest
gaelle.lab<-c("ColWT","d172","d263","isa2",
"sex4","dpe2","mex1","sex3","pgm","sex1","WsWT","tpt","RLDWT","ke103")
gaelle.bclust<-bclust(gaelle,rep.id=gaelle.id,labels=gaelle.lab,
transformed.par=c(-1.84,-0.99,1.63,0.08,-0.16,-1.68))
dptplot(gaelle.bclust,scale=5,varimp=imp(gaelle.bclust)$var,
horizbar.plot=TRUE)
#replicated clustering
gaelle.bclust<-bclust(gaelle,
transformed.par=c(-1.84,-0.99,1.63,0.08,-0.16,-1.68))
dptplot(gaelle.bclust,scale=10,varimp=imp(gaelle.bclust)$var,
horizbar.plot=TRUE,plot.width=5,horizbar.size=0.2,ylab.mar=4)
#unreplicated clustering
wildtype<-rep(1,55) #initiate a vector
wildtype[c(1:3,48:51,40:43)]<-2 #associate 2 to wildtypes
dptplot(gaelle.bclust,scale=10,varimp=imp(gaelle.bclust)$var,
horizbar.plot=TRUE,plot.width=5,horizbar.size=0.2,vertbar=wildtype,
vertbar.col=c("white","violet"),ylab.mar=4)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.