{

# read the files that contain the statistics and boxplot
#them

#par(mfrow=c(2,2))
#months = c("Oct","Nov","Dec","Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Ann")
months = c("Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec","Ann")
	
#stype = "means"
#stype = "sds"
#stype = "skews"
#stype = "corrs1"
#stype = "max"	
#stype = "min"

yrange = matrix(scan(paste("lfknn",stype, sep = "")),ncol=13,byrow=T)
	
source("myboxplots-stats.r")
source("myboxplots.r")	
			
## boxplot the modified KNN
xmeans = matrix(scan(paste("lfknn",stype, sep = "")),ncol=13,byrow=T)
n = length(xmeans[,1])
#the first row is the means of the original data
xmeans1 = xmeans[2:n,]
xs = 1:13
zz = myboxplot(split(t(xmeans1),xs), plot=F)
zz$names = rep("",length(zz$names))
z1 = bxp(zz,ylim=range(yrange),xlab="",ylab="",cex=.75)
axis(1,at=z1,labels=months, cex=.75, las=2)
points(z1,xmeans[1,],pch=16)
lines(z1[1:12],xmeans[1,1:12],pch=16)
title(main="Modified KNN")                                 

}