# plot cholestrole values over time for the paper

postscript(file="cholst.ps", horizontal = F, font=3, pointsize=14)
par(mfrow=c(1,1), pty="s")

out <- read.table(file="cholst.dat", col.names=c("newid", "id", "cholst",
 "sex", "age",  "year"))

plot(out$year, out$cholst, xlab="Time in years", ylab="Cholesterol level",
cex=0.8)

title("Cholesterol levels over time", cex=0.8)
lines(out$year[out$newid==2], out$cholst[out$id==2])
lines(out$year[out$newid==47], out$cholst[out$id==47])
lines(out$year[out$newid==74], out$cholst[out$id==74])
lines(out$year[out$newid==107], out$cholst[out$id==107])
lines(out$year[out$newid==171], out$cholst[out$id==171])	
	
dev.off()