{

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

par(mfrow=c(2,2))

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=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)
points(z1,xmeans[1,],pch=16)
lines(z1,xmeans[1,],pch=16)
title(main="Means")                                 


## 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)
points(z1,xmeans[1,],pch=16)
lines(z1,xmeans[1,],pch=16)
title(main="Standard Deviation") 



## 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)
points(z1,xmeans[1,],pch=16)
lines(z1,xmeans[1,],pch=16)
title(main="Skews")

## 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)
points(z1,xmeans[1,],pch=16)
lines(z1,xmeans[1,],pch=16)
title(main="Lag-1 correlation")

                                
}