# NOT RUN {
data(QWstacked)
# Plot the data for whole-water phosphorus
with(subset(QWstacked, parm_cd == "00665"), plot(sample_dt, result_va))
# Overlay the data for dissolved phosphorus (should be smaller values)
with(subset(QWstacked, parm_cd == "00666"), points(sample_dt, result_va, col='red'))
# }
Run the code above in your browser using DataLab