{

# read the files that contain the pdf statistics and plot them

#sets up the device to plot to
device = png
ext = "png"

#read observed flows
WYmonflow = scan("leesferry0695cy.txt")
#converts ac-ft/month to cms
WYmonflow = WYmonflow*.0004690502
obsmonflow = WYmonflow

#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")
	
source("myboxplots-stats.r")
source("myboxplots.r")	
source("simboxpdf3.r")
source("ypdfRange.r")

par(mfrow=c(2,2))	

yrange=0
ypdfRange(5,"lfknnmaysim",WYmonflow)	


# boxplot the modified KNN
simboxpdf(5,"lfknnmaysim",WYmonflow,yrange)
title(main="Modified KNN")                                 

yrange=0
ypdfRange(12,"lfknndecsim",WYmonflow)	
		
# boxplot the modified KNN
simboxpdf(12,"lfknndecsim",WYmonflow,yrange)
title(main="Modified KNN")                                 


}