# critical values for cloning data with 9 groups
tmpar <- par( mar=c(3.1,3.1,0,2.1), mfrow=c(1,2) )
on.exit( { par( tmpar ); rm( tmp, tmpx, tmpy, tmpp, zaxis, zmult, tmpar ) } )

tmpx <- c( 2, max( cloning.mst$bal$means) )
tmpy <- range( unlist( c( cloning.crit$bal, cloning.mst$bal[,-1] ) ) )
tmp <- round( mean( tmpx ) )
zaxis <- seq( 1.5, 5, .5 )
zmult <- sqrt( 2 / 4.91) * 198

plot( tmpx, tmpy, type="n", xlab="", ylab="" )
axis(4, zmult * zaxis, lab = as.character( zaxis ) )
mtext( "(a) based on F test", 1, 2 )
for ( i in c("LSD","BON","SCHEFFE","WALLER") )
{
	abline( cloning.crit$bal[i], 0, lty=2 )
	text( tmp, cloning.crit$bal[i], i )
}
for ( i in c("REGWF") )
{
	lines( cloning.mst$bal$means, cloning.mst$bal[[i]] )
	text( tmp, cloning.mst$bal[[i]][tmp], i )
}

plot( tmpx, tmpy, type="n", xlab="", ylab="" )
axis(4, zmult * zaxis,lab=as.character(zaxis))
mtext( "(b) based on range", 1, 2 )
mtext( "critical value", 2, 2.5 )

tmpp <- NULL
for ( i in c("LSD","HSD") )
	tmpp[i] <- cloning.crit$bal[i]
for ( i in c("REGWQ","SNK","DUNCAN") )
	tmpp[i] <- cloning.mst$bal[[i]][tmp]
tmpp["SNK"] <- tmpp["SNK"] - 25
tmpp["HSD"] <- tmpp["HSD"] + 25

for ( i in c("LSD","HSD") )
{
	abline( cloning.crit$bal[i], 0, lty=2 )
	text( tmp, tmpp[i], i )
}
for ( i in c("REGWQ","SNK","DUNCAN") )
{
	lines( cloning.mst$bal$means, cloning.mst$bal[[i]] )
	text( tmp, tmpp[i], i )
}