list <- c("01","02","03","04","05","06","07","08","09","10","11","12","13")

list <- c("01")

ps.options(horizontal=FALSE,onefile=FALSE)
ps.options(pagecentre=TRUE)
ps.options(paper="letter")
ps.options(width=7.0,height=10.0)
  
for (i in list) {

  run <- paste("0",i,"001",sep="") 
  
  filename <- paste("theta",run,".dat",sep="") 
  tmp <- scan(filename)
  rows <- tmp[1]
  cols <- tmp[2]
  
  post <- matrix(tmp[3:(2+rows*cols)],nrow=rows,ncol=cols)
  
  idx <- 1:cols
  
  theta <- c("g","R11","R12","R22","phi","delta","gamma")
  prior.lo <- c(-0.005,  0.00,  0.1,  0.00,  0.85, 0.97,  1.0 )
  prior.hi <- c( 0.005,  0.04, -0.1,  0.25,  1.00, 1.00,  3.0 )
  
  filename <- paste("thser",run,".eps",sep="")   
  postscript(file=filename)
  par(mfrow=c(rows,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(bottom,left,top,right)
  for (j in 1:rows)
  {
    plot(idx,post[j,],type="n",ylab="")
      lines(idx,post[j,],lty="solid")
      title(theta[j])
  }
  dev.off()
  
  filename <- paste("thden",run,".eps",sep="")
  postscript(file=filename)
  par(mfrow=c(rows,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(bottom,left,top,right)
  for (j in 1:rows)
  {
    den <- density(post[j,], from=prior.lo[j], to=prior.hi[j],adjust=3)
    plot(den,main=theta[j])
  }
  dev.off()
  
  filename <- paste("thcor",run,".eps",sep="")
  postscript(file=filename) 
  par(mfrow=c(rows,rows),mar=c(4.5,4.0,1.5,1)+0.1,pty="s")
  for (j in 1:rows)
  {
    for (i in 1:rows)
    {
       plot(x=post[i,],y=post[j,],xlab=theta[i],ylab=theta[j],pch=".") 
    }
  }
  dev.off()

  tpost <- t(post);
  dimnames(tpost) <- list(NULL,theta);
  print(cor(tpost)) 
}

for (i in list) {

  run <- paste("0",i,"001",sep="") 
  
  filename <- paste("stats",run,".dat",sep="") 
  tmp <- scan(filename)
  rows <- tmp[1]
  cols <- tmp[2]
  
  post <- matrix(tmp[3:(2+rows*cols)],nrow=rows,ncol=cols)
  
  idx <- 1:cols
  
  stats <- c("bond mean","bond variance","stock mean","stock variance") 
  prior.lo <- c( 0.99, 0,     0.99, 0)  
  prior.hi <- c( 1.05, 0.015, 1.05, 0.025)
  
  filename <- paste("stser",run,".eps",sep="")   
  postscript(file=filename)
  par(mfrow=c(rows,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(bottom,left,top,right)
  for (j in 1:rows)
  {
    plot(idx,post[j,],type="n",ylab="")
      lines(idx,post[j,],lty="solid")
      title(stats[j])
  }
  dev.off()
  
  filename <- paste("stden",run,".eps",sep="")
  postscript(file=filename)
  par(mfrow=c(rows,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(bottom,left,top,right)
  for (j in 1:rows)
  {
    den <- density(post[j,], from=prior.lo[j], to=prior.hi[j],adjust=3)
    plot(den,main=stats[j])
  }
  dev.off()
  
  filename <- paste("stcor",run,".eps",sep="")
  postscript(file=filename) 
  par(mfrow=c(rows,rows),mar=c(4.5,4.0,1.5,1)+0.1,pty="s")
  for (j in 1:rows)
  {
    for (i in 1:rows)
    {
       plot(x=post[i,],y=post[j,],xlab=stats[i],ylab=stats[j],pch=".") 
    }
  }
  dev.off()
  
  spost <- t(post)
  dimnames(spost) <- list(NULL,stats);
  print(cor(spost)) 
}