thr3ads.net - R help - [R] cenboxplot(): Reporting Limit Twice Correct Concentration [Sep 2012]

If this information is useful, please help other people find it:
Share via:

Rich Shepard

2012-Sep-04 23:14 UTC

[R] cenboxplot(): Reporting Limit Twice Correct Concentration

I've gone over the data and do not see my error; the dput() output of the
data frame and the pdf output of cenboxplot() are attached.

   The command used:

cenboxplot(sb.t$quant, sb.t$ceneq1, range=1.5, main='Total Recoverable
Antimony', xlab='Pre-Mining Era', ylab='Concentration (log
mg/L)')

(on a single line in emacs).

   The RL on the plot is drawn at 0.01 rather than at 0.005 and I'd like to
learn why.

Rich
-------------- next part --------------
structure(list(site = structure(c(12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 
12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 16L, 16L, 17L, 
17L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 
22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 22L, 
22L, 22L, 22L, 22L, 22L, 22L, 22L, 31L, 31L, 31L, 31L, 31L, 31L, 
31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 
31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 31L, 
32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 
32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 32L, 
32L, 32L, 32L, 32L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 
33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 
33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 33L, 35L, 35L, 35L, 
35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 
35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 35L, 
35L, 35L, 35L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 
36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 
36L, 36L, 36L, 36L, 36L, 36L, 36L, 36L, 37L, 37L, 37L, 37L, 37L, 
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 
37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 37L, 
37L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 
47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 47L, 
47L, 47L, 47L, 47L, 47L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 
48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 48L, 49L, 49L, 49L, 49L, 
49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 49L, 
49L, 49L, 49L, 49L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
50L, 50L, 50L, 50L, 50L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 
51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 51L, 
51L, 51L, 51L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 
52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 52L, 57L, 57L, 
57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 
57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 57L, 
59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 
59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 59L, 
59L, 59L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 
60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 60L, 
60L, 60L, 60L, 60L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 
61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 61L, 62L, 62L, 
62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 
62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 62L, 
62L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 
63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 
63L, 63L, 63L, 63L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 64L, 
64L, 64L, 64L, 64L, 64L), .Label = c("D-1", "D-2",
"D-3", "D-4",
"D-5", "D-6", "D-7", "D-8",
"D-9", "Effluent", "GB-ES1", "GW-2",
"Gold Bowl Sump", "Influent", "JJ-16",
"JJ-18", "JJ-20", "JJ-22",
"MC-1", "MC-2", "MC-3", "MW-1",
"MW-10", "MW-11", "MW-12", "MW-13",
"MW-14", "MW-15", "MW-16", "MW-2R",
"MW-3", "MW-4", "MW-5", "MW-6R",
"MW-7", "MW-8", "MW-9", "Mine Domestic
Well", "P-10", "P-5",
"P-6", "P-8", "P-9", "RAS-001",
"SN-22", "SW Zone Sump", "SW-1",
"SW-10", "SW-11", "SW-12", "SW-13",
"SW-14", "SW-15", "SW-16",
"SW-17", "SW-18", "SW-2", "SW-3",
"SW-4", "SW-5", "SW-6", "SW-7",
"SW-8", "SW-9"), class = "factor"), sampdate =
structure(c(8421,
8442, 8477, 8503, 8537, 8567, 8595, 8596, 8636, 8666, 8686, 8713, 
8741, 8782, 8810, 8833, 8861, 8889, 8924, 8957, 8988, 9022, 9049, 
9114, 9142, 9246, 9274, 9295, 9407, 8573, 8693, 8573, 8693, 8363, 
8390, 8421, 8442, 8476, 8505, 8532, 8567, 8604, 8638, 8665, 8686, 
8717, 8741, 8778, 8805, 8833, 8860, 8887, 8924, 8959, 8987, 9022, 
9050, 9113, 9141, 9169, 9205, 9244, 9274, 9296, 9408, 8364, 8391, 
8421, 8441, 8476, 8505, 8531, 8566, 8594, 8636, 8665, 8686, 8714, 
8741, 8777, 8805, 8832, 8861, 8888, 8923, 8959, 8987, 9022, 9050, 
9112, 9141, 9168, 9205, 9244, 9274, 9293, 9408, 8367, 8392, 8420, 
8441, 8475, 8504, 8531, 8566, 8594, 8637, 8664, 8687, 8715, 8742, 
8777, 8806, 8832, 8861, 8888, 8923, 8958, 8987, 9021, 9140, 9168, 
9204, 9246, 9274, 9295, 9406, 8364, 8392, 8421, 8441, 8476, 8504, 
8531, 8566, 8594, 8637, 8664, 8687, 8715, 8747, 8777, 8806, 8832, 
8861, 8888, 8923, 8958, 8987, 9021, 9049, 9112, 9140, 9168, 9204, 
9246, 9274, 9296, 9406, 8364, 8397, 8419, 8441, 8474, 8505, 8531, 
8565, 8594, 8638, 8664, 8686, 8714, 8741, 8777, 8805, 8832, 8861, 
8888, 8923, 8959, 8987, 9022, 9050, 9112, 9141, 9168, 9204, 9246, 
9274, 9295, 9406, 8363, 8391, 8419, 8440, 8475, 8503, 8530, 8565, 
8593, 8637, 8663, 8684, 8713, 8740, 8776, 8803, 8831, 8860, 8887, 
8922, 8957, 8986, 9020, 9048, 9113, 9140, 9204, 9244, 9267, 9296, 
9412, 8363, 8390, 8421, 8442, 8476, 8505, 8531, 8567, 8594, 8638, 
8665, 8686, 8714, 8741, 8778, 8805, 8833, 8861, 8888, 8924, 8959, 
8987, 9022, 9050, 9112, 9141, 9168, 9205, 9244, 9274, 9295, 9408, 
8420, 8441, 8475, 8504, 8532, 8568, 8595, 8638, 8664, 8687, 8717, 
8747, 8785, 8810, 8833, 8859, 8862, 8888, 8923, 8958, 8988, 9021, 
9049, 9112, 9147, 9178, 9247, 9268, 9300, 9406, 8506, 8537, 8568, 
8596, 8638, 8665, 8687, 8714, 8889, 8924, 8959, 8986, 9022, 9051, 
9261, 9307, 9413, 8422, 8443, 8477, 8506, 8537, 8569, 8596, 8638, 
8666, 8691, 8713, 8743, 8783, 8889, 8924, 8959, 9023, 9051, 9265, 
9309, 9414, 8422, 8443, 8477, 8506, 8537, 8569, 8596, 8638, 8666, 
8691, 8713, 8743, 8783, 8811, 8834, 8862, 8889, 8924, 8959, 9023, 
9051, 9147, 9178, 9211, 9265, 9309, 9414, 8422, 8443, 8477, 8506, 
8537, 8569, 8596, 8638, 8666, 8691, 8713, 8743, 8783, 8811, 8834, 
8862, 8889, 8924, 8959, 9023, 9051, 9265, 9309, 9414, 8422, 8443, 
8477, 8537, 8569, 8596, 8638, 8666, 8691, 8713, 8743, 8783, 8862, 
8889, 8924, 8959, 9023, 9051, 9265, 9309, 9414, 8418, 8440, 8475, 
8502, 8530, 8565, 8593, 8637, 8684, 8713, 8740, 8783, 8810, 8833, 
8860, 8887, 8922, 8957, 8986, 9020, 9048, 9114, 9147, 9178, 9247, 
9267, 9300, 9412, 8422, 8443, 8477, 8506, 8537, 8569, 8596, 8640, 
8665, 8687, 8714, 8782, 8810, 8834, 8862, 8889, 8924, 8959, 8985, 
9023, 9051, 9114, 9147, 9178, 9246, 9261, 9307, 9413, 8422, 8443, 
8477, 8507, 8537, 8568, 8596, 8640, 8666, 8691, 8715, 8747, 8783, 
8811, 8834, 8862, 8889, 8924, 8959, 8988, 9023, 9051, 9147, 9181, 
9211, 9261, 9309, 9415, 8421, 8441, 8476, 8504, 8537, 8568, 8594, 
8638, 8664, 8687, 8715, 8782, 8861, 8888, 8923, 8958, 9021, 9274, 
9296, 9406, 8420, 8441, 8475, 8504, 8532, 8568, 8594, 8637, 8664, 
8687, 8715, 8742, 8777, 8806, 8832, 8861, 8888, 8923, 8958, 8987, 
9021, 9049, 9112, 9140, 9168, 9246, 9274, 9295, 9406, 8418, 8440, 
8475, 8502, 8530, 8565, 8593, 8637, 8663, 8687, 8713, 8740, 8783, 
8810, 8834, 8860, 8887, 8922, 8957, 8986, 9020, 9048, 9114, 9147, 
9178, 9247, 9267, 9300, 9412, 8505, 8537, 8568, 8595, 8636, 8665, 
8687, 8714, 8889, 8924, 8958, 8988, 9268, 9296), class = "Date"), 
    era = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L), .Label = c("Mining", "Pre-mining"),
class = "factor"),
    param = structure(c(50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 
    50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L, 50L), .Label =
c("AgDis",
    "AgTot", "AlDis", "AlTot", "Alk",
"AsDis", "AsTot", "BaDis",
    "BaTot", "BeDis", "BeTot", "BiDis",
"BiTot", "CaDis", "CaTot",
    "CdDis", "CdTot", "ClDis", "ClTot",
"CoDis", "CoTot", "CrDis",
    "CrTot", "CuDis", "CuTot", "DO",
"FeDis", "FeTot", "HgDis",
    "HgTot", "MgDis", "MgTot", "MnDis",
"MnTot", "MoDis", "MoTot",
    "N", "NH4", "NaDis", "NaTot",
"NiDis", "NiTot", "PDis", "PTot",
    "PbDis", "PbTot", "SC", "SO4",
"SbDis", "SbTot", "SeDis",
    "SeTot", "SrDis", "SrTot", "TDS",
"TSS", "TlDis", "TlTot",
    "Vdis", "Vtot", "ZnDis", "ZnTot",
"oilgrease", "pH"), class = "factor"),
    quant = c(0.005, 0.005, 0.005, 0.005, 0.005, 0.007, 0.005, 
    0.009, 0.006, 0.005, 0.007, 0.005, 0.007, 0.006, 0.005, 0.009, 
    0.01, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 
    0.007, 0.009, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 
    0.009, 0.008, 0.005, 0.008, 0.008, 0.005, 0.006, 0.012, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.006, 0.008, 0.008, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.008, 0.005, 0.006, 
    0.006, 0.005, 0.008, 0.011, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.007, 0.006, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.009, 
    0.005, 0.005, 0.008, 0.005, 0.005, 0.009, 0.011, 0.008, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.007, 0.008, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.008, 0.007, 
    0.01, 0.008, 0.011, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.01, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.007, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.006, 0.008, 0.006, 0.006, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.007, 0.01, 0.005, 0.007, 0.005, 0.005, 0.005, 
    0.009, 0.009, 0.008, 0.007, 0.005, 0.005, 0.005, 0.005, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.007, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.01, 0.009, 0.007, 0.006, 0.005, 0.006, 0.008, 0.008, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.008, 0.005, 0.007, 0.005, 0.005, 
    0.005, 0.005, 0.007, 0.005, 0.005, 0.005, 0.005, 0.006, 0.008, 
    0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.006, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 0.009, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.009, 0.005, 0.005, 0.006, 0.005, 0.006, 
    0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.008, 0.008, 0.005, 
    0.005, 0.005, 0.008, 0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.015, 0.007, 0.005, 0.005, 
    0.005, 0.005, 0.007, 0.006, 0.006, 0.008, 0.005, 0.005, 0.007, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.007, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.006, 0.007, 
    0.005, 0.005, 0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 0.009, 0.005, 
    0.005, 0.005, 0.005, 0.011, 0.007, 0.007, 0.005, 0.005, 0.005, 
    0.005, 0.007, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.01, 0.005, 0.005, 0.005, 0.005, 0.005, 0.008, 0.005, 
    0.005, 0.007, 0.005, 0.007, 0.008, 0.008, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.008, 0.005, 0.005, 0.005, 0.006, 0.005, 
    0.006, 0.005, 0.005, 0.005, 0.01, 0.006, 0.007, 0.005, 0.005, 
    0.005, 0.01, 0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.012, 0.006, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.009, 0.007, 0.005, 0.006, 0.007, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.009, 0.006, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.006, 0.005, 0.011, 0.005, 0.005, 0.005, 
    0.006, 0.01, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.012, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.014, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.01, 0.009, 0.006, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.008, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005), ceneq1 = c(TRUE, TRUE, 
    TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, FALSE, TRUE, FALSE, 
    TRUE, FALSE, FALSE, TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, FALSE, TRUE, TRUE, FALSE, FALSE, TRUE, FALSE, FALSE, 
    TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, FALSE, FALSE, TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, FALSE, TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, FALSE, 
    FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, 
    FALSE, TRUE, TRUE, TRUE, FALSE, FALSE, FALSE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, FALSE, FALSE, FALSE, TRUE, FALSE, FALSE, 
    FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    FALSE, TRUE, TRUE, FALSE, TRUE, FALSE, FALSE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, FALSE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, FALSE, FALSE, 
    TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, 
    FALSE, TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, FALSE, TRUE, 
    FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, TRUE, TRUE, TRUE, FALSE, 
    FALSE, FALSE, TRUE, TRUE, TRUE, FALSE, TRUE, TRUE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, FALSE, 
    FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, FALSE, TRUE, TRUE, 
    TRUE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, FALSE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, FALSE, FALSE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE, FALSE, TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE), floor = c(0, 0, 0, 0, 0, 0.007, 0, 0.009, 0.006, 
    0, 0.007, 0, 0.007, 0.006, 0, 0.009, 0.01, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.006, 0.007, 0.009, 0, 
    0, 0, 0, 0.006, 0, 0, 0.009, 0.008, 0, 0.008, 0.008, 0, 0.006, 
    0.012, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.006, 0, 0, 0, 0, 0.006, 
    0.008, 0.008, 0, 0, 0, 0, 0, 0, 0, 0, 0.008, 0, 0.006, 0.006, 
    0, 0.008, 0.011, 0, 0, 0, 0, 0, 0, 0, 0, 0.007, 0.006, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.009, 0, 0, 0.008, 0, 0, 
    0.009, 0.011, 0.008, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0.007, 0.008, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.008, 
    0.007, 0.01, 0.008, 0.011, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0.007, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0.006, 0.008, 0.006, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0.006, 0, 0, 0, 0, 0, 0.007, 0.01, 0, 0.007, 
    0, 0, 0, 0.009, 0.009, 0.008, 0.007, 0, 0, 0, 0, 0.006, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.007, 0, 0, 0, 0, 0, 0.01, 
    0.009, 0.007, 0.006, 0, 0.006, 0.008, 0.008, 0.006, 0, 0, 
    0, 0, 0, 0, 0, 0, 0.006, 0, 0, 0, 0, 0.008, 0, 0.007, 0, 
    0, 0, 0, 0.007, 0, 0, 0, 0, 0.006, 0.008, 0.006, 0, 0, 0, 
    0, 0, 0.006, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.006, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0.006, 0, 0, 0, 0, 0.009, 0, 0, 
    0, 0, 0, 0.009, 0, 0, 0.006, 0, 0.006, 0.006, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.008, 0.008, 0, 0, 0, 0.008, 
    0.008, 0, 0, 0, 0, 0, 0, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0.015, 0.007, 0, 0, 0, 0, 0.007, 0.006, 0.006, 0.008, 
    0, 0, 0.007, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.007, 0, 0, 0, 0, 
    0.006, 0, 0.006, 0.007, 0, 0, 0.008, 0, 0, 0, 0, 0, 0, 0, 
    0, 0.006, 0, 0, 0, 0, 0.009, 0, 0, 0, 0, 0.011, 0.007, 0.007, 
    0, 0, 0, 0, 0.007, 0, 0, 0, 0, 0, 0, 0, 0, 0.01, 0, 0, 0, 
    0, 0, 0.008, 0, 0, 0.007, 0, 0.007, 0.008, 0.008, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.008, 0, 0, 0, 0.006, 
    0, 0.006, 0, 0, 0, 0.01, 0.006, 0.007, 0, 0, 0, 0.01, 0, 
    0, 0.006, 0, 0, 0, 0, 0, 0, 0, 0.012, 0.006, 0, 0, 0, 0, 
    0, 0, 0.009, 0.007, 0, 0.006, 0.007, 0, 0, 0, 0, 0, 0, 0, 
    0.009, 0.006, 0, 0, 0, 0, 0, 0.006, 0, 0.011, 0, 0, 0, 0.006, 
    0.01, 0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.012, 
    0, 0, 0, 0, 0, 0.014, 0.006, 0, 0, 0, 0, 0, 0.01, 0.009, 
    0.006, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0.006, 0, 0.008, 0, 0, 0, 0, 0, 0, 0), ceiling = c(0.005, 
    0.005, 0.005, 0.005, 0.005, 0.007, 0.005, 0.009, 0.006, 0.005, 
    0.007, 0.005, 0.007, 0.006, 0.005, 0.009, 0.01, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.007, 0.009, 0.005, 
    0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 0.009, 0.008, 0.005, 
    0.008, 0.008, 0.005, 0.006, 0.012, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.006, 0.008, 0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.008, 0.005, 0.006, 0.006, 0.005, 0.008, 
    0.011, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.007, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.009, 0.005, 0.005, 0.008, 
    0.005, 0.005, 0.009, 0.011, 0.008, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.007, 0.008, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.008, 0.007, 0.01, 0.008, 0.011, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.01, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.007, 0.006, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.006, 0.008, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.007, 
    0.01, 0.005, 0.007, 0.005, 0.005, 0.005, 0.009, 0.009, 0.008, 
    0.007, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.007, 0.005, 0.005, 0.005, 0.005, 0.005, 0.01, 0.009, 0.007, 
    0.006, 0.005, 0.006, 0.008, 0.008, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 0.005, 
    0.005, 0.008, 0.005, 0.007, 0.005, 0.005, 0.005, 0.005, 0.007, 
    0.005, 0.005, 0.005, 0.005, 0.006, 0.008, 0.006, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.006, 0.006, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 
    0.005, 0.005, 0.005, 0.009, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.009, 0.005, 0.005, 0.006, 0.005, 0.006, 0.006, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.008, 0.008, 0.005, 0.005, 0.005, 0.008, 
    0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.015, 0.007, 0.005, 0.005, 0.005, 0.005, 0.007, 
    0.006, 0.006, 0.008, 0.005, 0.005, 0.007, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.007, 0.005, 0.005, 
    0.005, 0.005, 0.006, 0.005, 0.006, 0.007, 0.005, 0.005, 0.008, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.009, 0.005, 0.005, 0.005, 0.005, 
    0.011, 0.007, 0.007, 0.005, 0.005, 0.005, 0.005, 0.007, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.01, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.008, 0.005, 0.005, 0.007, 0.005, 
    0.007, 0.008, 0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.008, 0.005, 0.005, 0.005, 0.006, 0.005, 0.006, 0.005, 0.005, 
    0.005, 0.01, 0.006, 0.007, 0.005, 0.005, 0.005, 0.01, 0.005, 
    0.005, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.012, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.009, 
    0.007, 0.005, 0.006, 0.007, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.009, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.006, 0.005, 0.011, 0.005, 0.005, 0.005, 0.006, 0.01, 0.006, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.012, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.014, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.01, 
    0.009, 0.006, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.006, 0.005, 0.008, 0.005, 0.005, 0.005, 0.005, 0.005, 
    0.005, 0.005)), .Names = c("site", "sampdate",
"era", "param",
"quant", "ceneq1", "floor", "ceiling"),
row.names = c(10064L,
10100L, 10136L, 10172L, 10208L, 10244L, 10280L, 10313L, 10347L, 
10383L, 10419L, 10455L, 10491L, 10527L, 10563L, 10616L, 10657L, 
10693L, 10729L, 10765L, 10801L, 10837L, 10873L, 10909L, 10945L, 
10984L, 11020L, 11056L, 11092L, 16420L, 16456L, 17759L, 17795L, 
20141L, 20176L, 20212L, 20248L, 20284L, 20319L, 20355L, 20391L, 
20427L, 20463L, 20499L, 20535L, 20571L, 20607L, 20643L, 20679L, 
20715L, 20751L, 20787L, 20823L, 20859L, 20895L, 20931L, 20967L, 
21003L, 21039L, 21075L, 21111L, 21147L, 21183L, 21219L, 21255L, 
28175L, 28211L, 28247L, 28282L, 28318L, 28353L, 28389L, 28425L, 
28461L, 28497L, 28533L, 28569L, 28605L, 28641L, 28677L, 28713L, 
28749L, 28785L, 28821L, 28857L, 28893L, 28929L, 28965L, 29001L, 
29036L, 29072L, 29108L, 29144L, 29180L, 29216L, 29252L, 29288L, 
30743L, 30779L, 30815L, 30851L, 30887L, 30922L, 30958L, 30994L, 
31030L, 31066L, 31102L, 31138L, 31174L, 31210L, 31246L, 31282L, 
31318L, 31354L, 31390L, 31426L, 31462L, 31498L, 31534L, 31573L, 
31608L, 31644L, 31680L, 31716L, 31752L, 31788L, 33270L, 33306L, 
33342L, 33378L, 33414L, 33449L, 33485L, 33521L, 33557L, 33593L, 
33629L, 33665L, 33701L, 33737L, 33773L, 33809L, 33845L, 33881L, 
33917L, 33953L, 33989L, 34025L, 34061L, 34097L, 34133L, 34169L, 
34205L, 34241L, 34277L, 34313L, 34349L, 34385L, 34849L, 34885L, 
34921L, 34956L, 34992L, 35027L, 35063L, 35099L, 35135L, 35171L, 
35207L, 35243L, 35279L, 35315L, 35351L, 35387L, 35423L, 35459L, 
35495L, 35531L, 35567L, 35603L, 35639L, 35675L, 35711L, 35747L, 
35783L, 35819L, 35855L, 35891L, 35927L, 35963L, 37438L, 37474L, 
37510L, 37546L, 37582L, 37618L, 37654L, 37690L, 37726L, 37762L, 
37798L, 37834L, 37870L, 37906L, 37942L, 37978L, 38014L, 38050L, 
38086L, 38122L, 38158L, 38194L, 38230L, 38266L, 38302L, 38338L, 
38377L, 38413L, 38449L, 38485L, 38521L, 38642L, 38677L, 38713L, 
38749L, 38785L, 38820L, 38856L, 38892L, 38928L, 38964L, 39000L, 
39036L, 39072L, 39108L, 39144L, 39180L, 39216L, 39252L, 39288L, 
39324L, 39360L, 39396L, 39432L, 39468L, 39503L, 39539L, 39575L, 
39611L, 39647L, 39683L, 39719L, 39755L, 43504L, 43540L, 43576L, 
43611L, 43647L, 43683L, 43719L, 43755L, 43791L, 43827L, 43863L, 
43899L, 43935L, 43971L, 44024L, 44065L, 44099L, 44134L, 44170L, 
44206L, 44242L, 44278L, 44314L, 44350L, 44386L, 44421L, 44457L, 
44493L, 44529L, 44565L, 45060L, 45096L, 45132L, 45168L, 45204L, 
45240L, 45276L, 45312L, 45348L, 45384L, 45420L, 45456L, 45492L, 
45528L, 45564L, 45600L, 45636L, 46273L, 46309L, 46345L, 46380L, 
46416L, 46452L, 46488L, 46524L, 46560L, 46596L, 46632L, 46668L, 
46704L, 46740L, 46776L, 46812L, 46848L, 46884L, 46920L, 46956L, 
46992L, 47744L, 47780L, 47816L, 47851L, 47887L, 47923L, 47959L, 
47995L, 48031L, 48067L, 48103L, 48139L, 48175L, 48211L, 48264L, 
48305L, 48341L, 48377L, 48413L, 48449L, 48485L, 48521L, 48557L, 
48593L, 48629L, 48665L, 48701L, 49604L, 49640L, 49676L, 49711L, 
49747L, 49783L, 49819L, 49855L, 49891L, 49927L, 49963L, 49999L, 
50035L, 50071L, 50124L, 50165L, 50201L, 50237L, 50273L, 50309L, 
50345L, 50381L, 50417L, 50453L, 51352L, 51388L, 51424L, 51460L, 
51496L, 51532L, 51568L, 51604L, 51640L, 51676L, 51712L, 51748L, 
51784L, 51820L, 51856L, 51892L, 51928L, 51964L, 52000L, 52036L, 
52072L, 53379L, 53415L, 53451L, 53487L, 53523L, 53559L, 53595L, 
53631L, 53670L, 53706L, 53742L, 53778L, 53814L, 53867L, 53908L, 
53944L, 53980L, 54016L, 54052L, 54088L, 54124L, 54160L, 54196L, 
54232L, 54268L, 54304L, 54340L, 54376L, 54905L, 54941L, 54977L, 
55012L, 55048L, 55084L, 55120L, 55156L, 55192L, 55228L, 55264L, 
55300L, 55336L, 55389L, 55430L, 55466L, 55502L, 55538L, 55574L, 
55610L, 55646L, 55682L, 55718L, 55754L, 55790L, 55826L, 55862L, 
55898L, 56835L, 56871L, 56907L, 56942L, 56978L, 57014L, 57050L, 
57086L, 57122L, 57158L, 57194L, 57230L, 57266L, 57302L, 57355L, 
57396L, 57432L, 57468L, 57504L, 57540L, 57576L, 57612L, 57648L, 
57684L, 57720L, 57756L, 57792L, 57828L, 58629L, 58665L, 58701L, 
58736L, 58772L, 58808L, 58844L, 58880L, 58916L, 58952L, 58988L, 
59024L, 59060L, 59096L, 59132L, 59168L, 59204L, 59240L, 59276L, 
59312L, 59576L, 59612L, 59648L, 59683L, 59719L, 59755L, 59791L, 
59827L, 59863L, 59899L, 59935L, 59971L, 60007L, 60043L, 60095L, 
60136L, 60172L, 60208L, 60244L, 60280L, 60316L, 60352L, 60388L, 
60424L, 60460L, 60510L, 60546L, 60582L, 60618L, 61997L, 62033L, 
62069L, 62105L, 62141L, 62177L, 62213L, 62249L, 62285L, 62321L, 
62357L, 62393L, 62429L, 62465L, 62518L, 62559L, 62595L, 62631L, 
62667L, 62703L, 62739L, 62775L, 62811L, 62847L, 62883L, 62919L, 
62955L, 62990L, 63025L, 64108L, 64144L, 64180L, 64216L, 64252L, 
64288L, 64324L, 64360L, 64396L, 64432L, 64468L, 64504L, 64540L, 
64576L), class = "data.frame")
-------------- next part --------------
A non-text attachment was scrubbed...
Name: total.sb.pdf
Type: application/pdf
Size: 4758 bytes
Desc: 
URL:
<https://stat.ethz.ch/pipermail/r-help/attachments/20120904/0e56b5d2/attachment-0002.pdf>

Nordlund, Dan (DSHS/RDA)

2012-Sep-05 01:07 UTC

head link

[R] cenboxplot(): Reporting Limit Twice Correct Concentration

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Rich Shepard
> Sent: Tuesday, September 04, 2012 4:15 PM
> To: r-help at r-project.org
> Subject: [R] cenboxplot(): Reporting Limit Twice Correct Concentration
> 
>    I've gone over the data and do not see my error; the dput() output
> of the
> data frame and the pdf output of cenboxplot() are attached.
> 
>    The command used:
> 
> cenboxplot(sb.t$quant, sb.t$ceneq1, range=1.5, main='Total Recoverable
> Antimony', xlab='Pre-Mining Era', ylab='Concentration (log
mg/L)')
> 
> (on a single line in emacs).
> 
>    The RL on the plot is drawn at 0.01 rather than at 0.005 and I'd
> like to
> learn why.
> 
> Rich
Rich,

The line is drawn at the maximum censored value.  
>max(sb.t$quant[sb.t$ceneq1==TRUE])[1] 0.01

The line is drawn at 0.01 because there is a censored record where quant=0.01,
i.e.
>sb.t[which(sb.t$quant==0.01 & sb.t$ceneq1==TRUE),]      site   sampdate        era param quant ceneq1 floor ceiling
34169 MW-5 1995-01-10 Pre-mining SbTot  0.01   TRUE     0    0.01



Hope this is helpful,

Dan

Daniel J. Nordlund
Washington State Department of Social and Health Services
Planning, Performance, and Accountability
Research and Data Analysis Division
Olympia, WA 98504-5204

Seemingly Similar Threads

Search for more possibly parallel threads

R help - Sep 2012 - cenboxplot(): Reporting Limit Twice Correct Concentration

[R] cenboxplot(): Reporting Limit Twice Correct Concentration

[R] cenboxplot(): Reporting Limit Twice Correct Concentration

Seemingly Similar Threads