teethplot: produces teeth plot useful for demonstating a grouping on...

In bclust: Bayesian Hierarchical Clustering Using Spike and Slab Models

Description

The function alone is useless but can be attached to a horizontal dendrogram, a profile plot, or an image plot to show a specified partitioning.

Usage

teethplot(x, teeth.space = 0.25, teeth.lwd = 1)

Arguments

x	A bclustvs object.
teeth.space	The space between two teeth, a value between 0 and 0.25.
teeth.lwd	The thickness of the lines used to draw the teeth.

Details

The teeth plot for the moment shows the grouping vertically.

See Also

profileplot, ditplot, dptplot.

Examples

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),var.select=TRUE)
#start plotting
layout(matrix(c(1,2),1,2,byrow=TRUE), c(9,1),10, respect=TRUE) 
# divide plot space into two unequal parts
par(mar=c(0,0,0,2)) 
# preserve some space for labels in dendrogram plot
plot(as.dendrogram(gaelle.bclust),
horiz=TRUE,yaxs="i") #plot the dendrogram
abline(v=gaelle.bclust$cut) 
#show the optimal allocation by a line on the dendrogram
par(mar=c(0,0,0,0)) # we need no space for teeth plot
teethplot(gaelle.bclust) #show the optimal allocation using teeth plot

Example output

2/14
3/14
4/14
5/14
6/14
7/14
8/14
9/14
10/14
11/14
12/14
13/14
14/14

bclust documentation built on May 2, 2019, 3:01 a.m.