soi = ts(scan("/mydata/soi.dat"), start=1950, frequency=12)
 rec = ts(scan("/mydata/recruit.dat"), start=1950, frequency=12)
 par(mfrow=c(2,1), mar=c(3, 4, 3, 2))   # sets up the graphics - ?par for info
 plot(soi, ylab="", xlab="", main="Southern Oscillation Index")
 plot(rec, ylab="", xlab="", main="Recruitment")