# NOT RUN {
data(redfieldPlankton, package="ocedata")
par(mfrow=c(1,2))
boxplot(redfieldPlankton$Nitrogen)
abline(h=16.7, col='gray')
boxplot(redfieldPlankton$Phosphorus)
abline(h=1.85, col='gray')
# }
Run the code above in your browser using DataLab