# plot cholestrole values over time for the paper

postscript(file="seize.ps", horizontal = F, 
width=8.5, height=6.5)
par(mfrow=c(1,2), pty="s")

out <- read.table(file="seize.dat", 
col.names=c("id", "seize", "visit", "trt", "age"))

plot(out$visit[out$trt==1], out$seize[out$trt==1], pch=" ",
xlab="Order of visit", ylab="Seizure counts", cex=0.8)
title("Seizure counts for progabide arm", cex=0.7)
points(out$visit[out$trt==1], out$seize[out$trt==1])
lines(out$visit[out$id==103], out$seize[out$id==103])
lines(out$visit[out$id==128], out$seize[out$id==128])
lines(out$visit[out$id==214], out$seize[out$id==214])

plot(out$visit[out$trt==0], out$seize[out$trt==0], pch=" ",
xlab="Order of visit", ylab="Seizure counts", cex=0.8)
title("Seizure counts for control arm", cex=0.7)
points(out$visit[out$trt==0], out$seize[out$trt==0])
lines(out$visit[out$id==107], out$seize[out$id==107])
lines(out$visit[out$id==130], out$seize[out$id==130])
lines(out$visit[out$id==205], out$seize[out$id==205])

dev.off()