particle <- read.table("particle.csv",header=T,sep=",")

ps.options(horizontal=FALSE,onefile=FALSE)
ps.options(pagecentre=TRUE)
ps.options(paper="letter")
#ps.options(width=7.0,height=10.0)
ps.options(width=7.0,height=5.0)

mean1 <- particle$mean1
mean2 <- particle$mean2
ulim1 <- mean1 + 2.0*particle$sdev1
llim1 <- mean1 - 2.0*particle$sdev1
ulim2 <- mean2 + 2.0*particle$sdev2
llim2 <- mean2 - 2.0*particle$sdev2
true1 <- particle$x1
true2 <- particle$x2
data <- particle$y
time <- seq(1,length(mean1))

postscript(file="particle1.eps")

par(mfrow=c(2,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(b,l,t,r)

xbot <- min(time)
xtop <- max(time)

ybot <- min(c(mean1,ulim1,llim1,true1))
ytop <- max(c(mean1,ulim1,llim1,true1))

plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
  lines(time,mean1,lty="solid",col="red")
  lines(time,ulim1,lty="dashed")
  lines(time,llim1,lty="dashed")

ybot <- min(c(mean2,ulim2,llim2,true2))
ytop <- max(c(mean2,ulim2,llim2,true2))

plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
  lines(time,mean2,lty="solid",col="red")
  lines(time,ulim2,lty="dashed")
  lines(time,llim2,lty="dashed")

dev.off()

postscript(file="particle2.eps")

par(mfrow=c(2,1),mar=c(2.5,4,1.5,2)+0.1) # mar=c(b,l,t,r)

xbot <- min(time)
xtop <- max(time)

ybot <- min(c(mean1,true1))
ytop <- max(c(mean1,true1))

plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
  lines(time,mean1,lty="solid",col="red")
  lines(time,true1,lty="solid",col="blue")

ybot <- min(c(mean2,true2))
ytop <- max(c(mean2,true2))

plot(c(xbot,xtop),c(ybot,ytop),type="n",ylab="",xlab="")
  lines(time,mean2,lty="solid",col="red")
  lines(time,true2,lty="solid",col="blue")

dev.off()