plottab<-function(filename="mcmc.tab",lm=400,rows=4)
{
        op<-par(mfrow=c(rows,3))
        tab<-read.table(filename,header=TRUE,check.names=FALSE)
        names<-attr(tab,"names")
        count<-0
        for (i in names[2:length(names)]) {
                plot(ts(tab[[i]]),
                        main=paste("Trace plot for",i),
                        ylab="Value",xlab="Iteration",col=4)
                v<-var(tab[[i]])
                if (v>1e-10) {
                        acf(tab[[i]],lag.max=lm,ci=0,
                                main=paste("ACF plot for",i),col=2)
                        d<-density(tab[[i]])
                        plot(d,
                                main=paste("Density for",i),
                                xlab="Value",col=4,lwd=3)
                } else {
                        plot(0,0)
                        plot(0,0)
                }
                print(paste("Mean for",i,"is",mean(tab[[i]]),"SD
is",sqrt(v)))
                count<-count+1
                if (count %% rows == 0) {
                        readline("Press return to continue... ")
                }
        }
        par(op)
        NULL
}