Chemical concentration data are in the columns that start with a P and are
followed by a number. Qualification codes for the concentration data are in
the columns that start with an R followed by the same numbers as the
associated concentration data. For example, column P04035 indicates
simazine data, 04035, being the U.S. Geological Survey parameter code
for simazine. The qualification codes for the simazine concentrations
are found in the column R04035, indicating a U.S. Geological Survey
remark code. Remark codes include _ or nothing, indicating no
qualification of the value in the associated concentration field; <,
indicating a censored value that is less than the number reported in the
associated concentration field; and E, indicating that the value has
been estimated. See Oblinger Childress and others (1999) for information
on the remark codes used by the U.S. Geological Survey.
Source
Data provided by Patrick Phillips, U.S. Geological Survey,
New York Water Science Center.
References
Oblinger Childress, C.J., Foreman, W.T., Connor, B.F., and Maloney, T.J.,
1999, New reporting procedures based on long-term method detection levels
and some considerations for interpretations of water-quality data provided
by the U.S. Geological Survey Open-File Report 99–193, 19 p. (Also
available at http://water.usgs.gov/owq/OFR_99-193/index.html.)
Examples
data(swData)
# summary of water-quality concentrations
apply(qwMoRivOmaha[,grep("P[[:digit:]]",dimnames(qwMoRivOmaha)[[2]])],
2, summary)
# simple boxplot of water-quality concentrations
rosBoxPlot(qwMoRivOmaha, qwcols=c("R", "P"))
# same boxplot function with many additional plotting arguments
rosBoxPlot(qwMoRivOmaha, site="06610000 Missouri River at Omaha, Nebr.",
log="y", yaxt="n", ylim=c(0.000001, 10), qwcols=c("R", "P"),
ylab=c("Concentration, micrograms per liter"), col="skyblue1",
cex.axis=0.7, cex.sub=0.8, par(tcl=0.5, las=1,
yaxs="i",
mgp=c(3,0.5,0),
mar=c(5,5,2,2),
cex.main=0.9))
axis(2, at=c(0.000001, 0.00001, 0.0001, 0.001, 0.01, 0.1, 1, 10),
labels=c("0.000001", "0.00001", "0.0001", "0.001", "0.01",
"0.1", "1", "10"), cex.axis=0.7)
# scatter plot of Simazine concentrations
cenScatPlot(qwMoRivOmaha, pname="04035")
# scatter plot with many additional plotting arguments
par(las=1, tcl=0.5)
cenScatPlot(qwMoRivOmaha, pname="04035",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Simazine concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0,0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"), as.Date("2004-01-01")),
xaxs="i", xaxt="n")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Prometon scatter plot
cenScatPlot(qwMoRivOmaha, pname="04037",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Prometon concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0,0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Cyanazine scatter plot
cenScatPlot(qwMoRivOmaha, pname="04041",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Cyanazine concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0.001,5), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n", log="y")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Metolachlor scatter plot
cenScatPlot(qwMoRivOmaha, pname="39415",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Metolachlor concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0.001,5), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n", log="y", legpos="bottomleft")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Alachlor scatter plot
cenScatPlot(qwMoRivOmaha, pname="46342",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Alachlor concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0, 0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Metribuzin scatter plot
cenScatPlot(qwMoRivOmaha, pname="82630",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Metribuzin concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0, 0.1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# Trifluralin scatter plot
cenScatPlot(qwMoRivOmaha, pname="82661",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="Trifluralin concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0,0.03), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)
# EPTC scatter plot
cenScatPlot(qwMoRivOmaha, pname="82668",
site="06610000 Missouri River at Omaha, Nebr.",
ylabel="EPTC concentration, in micrograms per liter",
legcex=0.7, qwcols=c("R", "P"),
ylim=c(0.001,1), yaxs="i", cex.lab=0.9, cex.axis=0.9,
xlim=c(as.Date("1996-01-01"),
as.Date("2004-01-01")), xaxs="i",
xaxt="n", log="y")
axdates<-c("1996-01-01", "1998-01-01", "2000-01-01", "2002-01-01",
"2004-01-01")
axis(1, as.Date(axdates), labels=c("1996", "1998", "2000", "2002",
"2004"), cex.axis=0.9)