{

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

par(mfrow=c(2,2))

 source("http://civil.colorado.edu/~balajir/CVEN5454/R-sessions/sess1/myboxplot.r")
 source("http://civil.colorado.edu/~balajir/CVEN5454/R-sessions/sess1/myboxplot-stats.r")


months=c("Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec")


xmeans=matrix(scan("armeans"),ncol=12,byrow=T)
n=length(xmeans[,1])
#the first row is the means of the original data
xmeans1=xmeans[2:n,]
xs=1:12
zz=myboxplot(split(xmeans1,col(xmeans1)), plot=F, cex=1.0)
zz$names=rep("",length(zz$names))
z1=bxp(zz,ylim=range(xmeans),xlab="",ylab="",cex=1.00,axes=F)
axis(1,at=z1,labels=months)
axis(2)
points(z1,xmeans[1,],pch=16,col="red")
lines(z1,xmeans[1,],pch=16)
title(main="Means")                                 
box()


## boxplot the standard deviations.
xmeans=matrix(scan("arstdevs"),ncol=12,byrow=T)
n=length(xmeans[,1])
#the first row is the means of the original data
xmeans1=xmeans[2:n,]
xs=1:12
zz=boxplot(split(xmeans1,col(xmeans1)), plot=F, cex=1.0)
zz$names=rep("",length(zz$names))
z1=bxp(zz,ylim=range(xmeans),xlab="",ylab="",cex=1.00,axes=F)
axis(1,at=z1,labels=months)
axis(2)
points(z1,xmeans[1,],pch=16,col="red")
lines(z1,xmeans[1,],pch=16)
title(main="Standard Deviation") 
box()


## boxplot the skews
xmeans=matrix(scan("arskews"),ncol=12,byrow=T)
n=length(xmeans[,1])
#the first row is the means of the original data
xmeans1=xmeans[2:n,]
xs=1:12
zz=boxplot(split(xmeans1,col(xmeans1)), plot=F, cex=1.0)
zz$names=rep("",length(zz$names))
z1=bxp(zz,ylim=range(xmeans),xlab="",ylab="",cex=1.00,axes=F)
axis(1,at=z1,labels=months)
axis(2)
points(z1,xmeans[1,],pch=16,col="red")
lines(z1,xmeans[1,],pch=16)
title(main="Skews")
box()

## boxplot the lag-1 correlations
xmeans=matrix(scan("arcorrs1"),ncol=12,byrow=T)
n=length(xmeans[,1])
#the first row is the means of the original data
xmeans1=xmeans[2:n,]
xs=1:12
zz=boxplot(split(xmeans1,col(xmeans1)), plot=F, cex=1.0)
zz$names=rep("",length(zz$names))
z1=bxp(zz,ylim=range(xmeans),xlab="",ylab="",cex=1.00)
axis(1,at=z1,labels=months)
axis(2)
points(z1,xmeans[1,],pch=16,col="red")
lines(z1,xmeans[1,],pch=16)
title(main="Lag-1 correlation")
box()
                                
}