similar to: R-help Digest, Vol 112, Issue 6

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

The sample data sets that come with the NADA package are limited to one or two variables and a censored measurement indicator column. I try to mimic examples using my data but keep missing the target. My water chemistry data is available in two formats: long (as seen in a database table) and wide (as seen in a spreadsheet). The two structures are: str(chem) 'data.frame': 65349 obs. of

I've successfully reformatted data frames from long to wide with reshape2, but this time I'm getting errors that I want to understand and resolve. Here's the data frame structure and the results of the melt() and dcast() functions: str(waterchem) 'data.frame': 128412 obs. of 8 variables: $ site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1

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

One set of data has censored (less-than detection limits) water chemistry concentrations for 80-100% of all observations. My initial trial-and-error attempts to apply the cenboxplot() method suggests that it has an upper limit to the percentage of censored observations. I do not see this limit in Dennis Helsel's second edition. Has anyone experience plotting censored data and can provide

I have a data frame with 1 factor, one date, and 37 numeric values: str(waterchem) 'data.frame': 3525 obs. of 39 variables: site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1 1 1 1 ... $ sampdate : Date, format: "2007-12-12" "2008-03-15" ... $ CO3 : num 1 1 6.7 1 1 1 1 1 1 1 ... $ HCO3 : num 231 228 118 246

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

I'm brand new to Samba and I am trying to set in up on our network (linux OS). We have two windows clients: a notebook running win95 (soon to be win98) and a virtual machine running win98. Right now, the host is running RH 5.2 with the 2.0.36 kernel. Samba-2.04b was installed from the rpm. I'm using Carter and Sharpe's "Teach Yourself Samba in 24 Hours" as my guide. On

I am not sure if this is an R-users question, but since most of you here are statisticians, I decided to give it a shot. I am using the lm() function in R to fit a dependent variable to a set of 3 to 5 independent variables. For this, I used the following commands: >model1<-lm(function=PBW~SO4+NO3+NH4) Coefficients: (Intercept) SO4 NO3 NH4 0.01323 0.01968

Thanks to David's help I subset my large data set and produced a smaller one for a single stream and 7 factors of interest. The structure of this data frame is: str(burns.tds.anal) 'data.frame': 718 obs. of 4 variables: $ site : Factor w/ 143 levels "BC-0.5","BC-1",..: 1 1 4 6 4 4 4 5 5 5 $ sampdate: Date, format: "1996-06-02"

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

I have measured values for 47 chemicals in a stream. After processing the original data frame through reshape2, the recast data frame has this structure: 'data.frame': 256 obs. of 47 variables: $ site : Factor w/ 143 levels "BC-0.5","BC-1",..: 1 1 1 2 2 2 2 2 2 2 ... $ sampdate : Date, format: "1996-04-19" "1996-05-21" ... $ Acid :

I have water chemistry data with censored values (i.e., those less than reporting levels) in a data frame with a narrow (i.e., database table) format. The structure is: $ site : Factor w/ 64 levels "D-1","D-2","D-3",..: 1 1 1 1 1 1 1 1 ... $ sampdate: Date, format: "2007-12-12" "2007-12-12" ... $ preeq0 : logi TRUE TRUE TRUE TRUE TRUE

This is my first excursion into using reshape2 and I want to ensure that the melt() function call is syntactically correct. The unmodifed data frame is organized this way: head(tds.anal) site sampdate param quant 1 UDS-O 2006-12-06 TDS 10800 4 STC-FS 1996-06-14 Cond 280 7 UDS-O 2007-10-04 Mg 1620 9 UDS-O 2007-10-04 SO4 7580 19 JCM-10B 2007-06-21 Ca 79 20

I would appreciate pointers on what I should read to understand this output: summary(lm(TDS ~ Cond + Ca + Cl + Mg + Na + SO4)) Call: lm(formula = TDS ~ Cond + Ca + Cl + Mg + Na + SO4) Residuals: ALL 1 residuals are 0: no residual degrees of freedom! Coefficients: (6 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 125 NA

hi there, trying to follow the steps of Yan & Rossini 2003, I have two problems : first when I 'make CrossCompileBuild', I get : ****************************************************************************** make[1]: Entering directory `/home/ros/RWORK/CROSSSOMP/WinR/R-1.8.0/src/gnuwin32' make -f Makefile.docfiles make[2]: Entering directory

All but one of the summaries of multiple linear regressions in this analysis set present the residuals by min, 1Q, median, 3Q, and max. Example: lm(formula = TDS ~ Cond + Ca + Cl + Mg + Na + SO4, data = snow.cast) Residuals: Min 1Q Median 3Q Max -277.351 -32.551 -2.621 40.812 245.272 The one that doesn't has only a small number of rows (23) and presents the

I have read ?zoo but am not sure how to relate the parameters (x, order.by, frequency, and style) to my data.frame. The structure of the data.frame is 'data.frame': 11169 obs. of 4 variables: $ stream : Factor w/ 37 levels "Burns","CIL",..: 1 1 1 1 1 1 1 1 1 1 ... $ sampdate: Date, format: "1987-07-23" "1987-09-17" ... $ param : Factor w/

Hello all, I'm quite a newbie to wine, and have some basic things running under it like WT Library, & Seamonkey. I need to install the local Revenue software because they sent me all tax documentation by email now :-(. My PC: Athlon 2.6Ghz with 1 gig, 40G hd, running Fedora Core 5 ROS ((Irish) Revenue Online Service) System Requirements

I'm stuck getting Samba to permit my copy of win98 (running in a VMware virtual machine) to connect to the linux host. Specs: Samba-2.04b installed via rpm on a Red Hat 5.2 system (kernel 2.0.36). For whatever reason, win98 is asking for a network password when I click on the icon for the linux box in Network Neighborhood. When I installed samba via the rpm, I have no recollection of