pscript(horizo=F,file='fig1.ps')
par(mfrow=c(2,2),bty='l',las=1,oma=c(0,0,0,0))
x<-matrix(scan(file='fig1.rsd'),byrow=T,ncol=3)
plot(x[,3],x[,2],xlab='Predicted Value',ylab='Residual',main='(a)')
y<-x[,2]
n<-length(y)
sy<-sqrt(sum(y*y)/(n-2))
y<-y/sy
qqnorm(y,xlab='Expected Value',ylab='Observed Value',main='(b)')
lines(c(-1000,1000),c(-1000,1000))
dev.off()