# NOT RUN {
data(redfieldNP, package="ocedata")
plot(redfieldNP$PO4, redfieldNP$NO3, xlab=expression(PO[4]), ylab=expression(NO[3]))
abline(0, 20, lwd=3) # line N:P = 20:1, as shown by Redfield
# }
Run the code above in your browser using DataLab