powered by
data(venice)
data(venice) attach(venice) y <- y1[year>1930] # for analysis in Section 5 of Davison (2003) x <- year[year>1930]-1956 plot(x+1956,y,ylab="Sea level (cm)",xlab="Year") lm(y~x)
Run the code above in your browser using DataLab