pdf("strplot.pdf", width=10, height=7)
par(mfrow=c(3,2), mar=c(2,2,2,2))

mcmcdat <- read.csv("mcmcdata", header=T)
attach(mcmcdat)

plot(Rep, Alpha, type="line", col="blue", xlab="MCMC Chains", ylab="alpha")
plot(Rep, F1, type="line", col="blue", xlab="MCMC Chains", ylab="F1")
plot(Rep, F2, type="line", col="blue", xlab="MCMC Chains", ylab="F2")
plot(Rep, D1_2, type="line", col="blue", xlab="MCMC Chains", ylab="D1_2")
plot(Rep, Ln_Like, type="line", col="blue", xlab="MCMC Chains", ylab="log(Likelihood)")
plot(Rep, Est_lnpd, type="line", col="blue", xlim=c(100000,200000), xlab="MCMC Chains", ylab="ln(Pr(X|K))")

dev.off()