search for: ceneq1

...Date, format: "2007-12-12" "2007-12-12" ... $ preeq0 : logi TRUE TRUE TRUE TRUE TRUE TRUE ... $ param : Factor w/ 37 levels "Ag","Al","Alk_tot",..: 1 2 8 17 3 4 ... $ quant : num 0.005 0.106 1 231 231 0.011 0.001 0.002 0.001 100 ... $ ceneq1 : logi TRUE FALSE TRUE FALSE FALSE FALSE ... $ floor : num 0 0.106 0 231 231 0.011 0 0 0 100 ... $ ceiling : num 0.005 0.106 1 231 231 0.011 0.001 0.002 0.001 100 ... > chem.melt <- melt(waterchem, idvars = c('site', 'sampdate', 'preeq0', 'param', ...

...07-12-12" "2007-12-12" ... $ era : Factor w/ 2 levels "Post","Pre": 1 1 1 1 1 1 1 1 1 1 ... $ param : Factor w/ 64 levels "AgDis","AgTot",..: 2 4 5 7 11 15 25 ... $ quant : num 1.30e-04 1.06e-01 2.31e+02 1.13e-02 5.00e-03 ... $ ceneq1 : logi TRUE FALSE FALSE FALSE TRUE FALSE ... $ floor : num 0 0.106 231 0.0113 0 100 0 1.43 0 0.0239 ... $ ceiling : num 1.30e-04 1.06e-01 2.31e+02 1.13e-02 5.00e-03 2.39e-02 ... and str(chem.cast) 'data.frame': 56938 obs. of 70 variables: $ site : Factor w/ 64 levels &quo...

Before reading water chemistry into a data frame I removed all missing data. Yet when I try to run cenros() to summarize a specific chemical I get an error that I do not understand: with( subset(chem, param=='Ag'), cenros(quant,ceneq1) ) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : NA/NaN/Inf in 'y' I would like to learn what I did incorrectly so I can avoid these errors in the future. The data frame structure is str(chem) 'data.frame': 120309 obs. of 8 variables: $ site...

...te, format: "2007-12-12" "2007-12-12" ... $ preeq0 : logi TRUE TRUE TRUE TRUE TRUE TRUE ... $ param : Factor w/ 37 levels "Ag","Al","Alk_tot",..: 1 2 8 17 3 4 9 ... $ quant : num 0.005 0.106 1 231 231 0.011 0.001 0.002 0.001 100 ... $ ceneq1 : logi TRUE FALSE TRUE FALSE FALSE FALSE ... $ floor : num 0 0.106 0 231 231 0.011 0 0 0 100 ... $ ceiling : num 0.005 0.106 1 231 231 0.011 0.001 0.002 0.001 100 ... The logical 'preeq0' separates sampdate into two groups; 'ceneq1' indicates censored/uncensored values...

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...

...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...

...A values, and neither does the data frame resulting from applying the reshape2 melt() function to it. However, the data frame produced by the dcast() function does contain NAs for all chemicals. I assume this is because of the syntax I used: chem.cast <- dcast(chem.melt, site + sampdate + era + ceneq1 + floor + ceiling ~ param) How should I reshape the data frame from long to wide without adding these spurious NAs? Rich

Hi, I am using the cenmle function from the NADA package and some of my data are negative values. As a result the cenmle function will not work and NaN's are produced. I try to change the distribution to Gaussian, but it still will not run. Could somebody please help me with this? Thanks -- Shane [[alternative HTML version deleted]]

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 -------------- stru...

I've searched the Web with Google and do not find what might cause this particular error from an invocation of cenboxplot: cenboxplot(cu.t$quant, cu.t$ceneq1, cu.t$era, range=1.5, main='Total Recoverable Copper', ylab='Concentration (mg/L)', xlab='Time Period') Error in if (n > 0) (1L:n - a)/(n + 1 - 2 * a) else numeric() : argument is of length zero I do not an observation with a concentration of zero; the only zeros...