data(wilson, package="ocedata")
plot(wilson$Age, wilson$Distance, xlab="Age [My]", ylab="Distance [km]",
xlim=c(125,0), ylim=c(0,4000))
abline(0,4000/125)
m <- lm(Distance~Age-1, data=wilson)
# Wilson reports in cm/year; 1cm/y = 1km/1My * 1e5 / 1e6
spreadingRate <- as.numeric(coef(m)) / 10
print(spreadingRate)
abline(m, col=2)
Run the code above in your browser using DataLab