I need a nudge in the right direction to get started using NADA. I bought Helsel's second addition and am currently reading it; NADA is installed in R.
My data has been restructured with a couple of awk scripts. The data frame structure now has a flag if the quantity is censored (ceneq1 column) as well as a lower and upper limit for censored data. For present purposes, interval censoring can be ignored. The data frame structure is now: str(waterchem) 'data.frame': 46551 obs. of 7 variables: $ site : Factor w/ 126 levels "BC-0.5","BC-1",..: 22 22 22 13 3 13 ... $ sampdate: Date, format: "1996-05-22" "1996-07-19" ... $ param : Factor w/ 58 levels "-0.100","AGP",..: 47 58 10 16 16 26 ... $ quant : num 0.01 7.69 0.02 63.8 120 0.02 399 439 2 433 ... $ ceneq1 : int 1 0 0 0 0 0 0 0 0 0 ... $ low : num 0 7.69 0.02 63.8 120 0.02 399 439 2 433 ... $ high : num 0.01 7.69 0.02 63.8 120 0.02 399 439 2 433 ... What I want to first learn is how to specify a box plot (and whether I can use the lattice package) for specific chemicals. ?cenboxplot shows me the arguments, but I'm not entering them correctly, or there's a prerequisite step I need to take: cenboxplot(waterchem$quant, waterchem$ceneq1, group='SO4', log=T, range=1.5) Error in function (classes, fdef, mtable) : unable to find an inherited method for function "ros", for signature "numeric", "integer" Perhaps cenboxplot is looking for a separate data set and not a data frame? Or, perhaps I need to melt and re-cast the data frame to the wide format from the current narrow format? Pointers appreciated. Rich ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
