dat = read.table("sspot.dat",head=T) plot(ts(dat[,1],freq=12,start=c(1749,1)),ty="l",ylab="sunspots (monthly means)")