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)Run the code above in your browser using DataLab