# critical values for cloning data with all 13 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$all$means ))
tmpy <- range( unlist( c(cloning.crit$all,cloning.mst$all[,-1]) ) )
tmp <- round( mean( tmpx ) )
zaxis <- seq(1.5,5,.5)
zmult <- sqrt(2/5)*239

plot( tmpx, tmpy, type="n", xlab="", ylab="" )
axis(4,zmult*zaxis,lab=as.character(zaxis))
#axis(4,zmult*zaxis,lab=F,las=0)
mtext( "(a) based on F test", 1, 2 )
for ( i in c("LSD","BON","SCHEFFE","WALLER") )
{
	abline( cloning.crit$all[i], 0, lty=2 )
	text( tmp, cloning.crit$all[i], i )
}
for ( i in c("REGWF") )
{
	lines( cloning.mst$all$means, cloning.mst$all[[i]] )
	text( tmp, cloning.mst$all[[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$all[i]
for ( i in c("REGWQ","SNK","DUNCAN") )
	tmpp[i] <- cloning.mst$all[[i]][tmp]
tmpp["SNK"] <- tmpp["SNK"] - 25
tmpp["HSD"] <- tmpp["HSD"] + 25

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