# NOT RUN {
library(oce)
data(redfieldNC, package="ocedata")
plot(redfieldNC$CO2, redfieldNC$NO3, xlab=expression(CO[2]), ylab=expression(NO[3]))
slope <- 1/7
abline(-275.5, slope, lwd=3)
abline(-288.5, slope, lwd=3)
# }
Run the code above in your browser using DataLab