[R] Understanding cenros Error

MacQueen, Don macqueen1 at llnl.gov
Wed Jul 11 21:50:30 CEST 2012


An my "easy" but not very useful answer is that this particular subset
probably violates some assumption of the cenros() model. I myself would
start with simple inspections of the data, such as

  with( subset(chem, param=='Ag'),  table(ceneq1) )
  with( subset(chem, param=='Ag'), qqnorm(quant) )
  with( subset(chem, param=='Ag'), range(quant) )


in the hopes that something pops out.

Do you have any zeros in quant?  (see ?cenros)

-Don

-- 
Don MacQueen

Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062





On 7/10/12 12:54 PM, "Rich Shepard" <rshepard at appl-ecosys.com> wrote:

>   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    : Factor w/ 65 levels ";Influent","D-1",..: 2 2 2 2 2 2 ...
>  $ sampdate: 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 9 ...
>  $ quant   : num  0 0.106 1 231 231 0.011 0 0.002 0 100 ...
>  $ ceneq1  : logi  FALSE FALSE TRUE FALSE FALSE FALSE ...
>  $ floor   : num  0 0.106 0 231 231 0.011 0 0 0 100 ...
>  $ ceiling : Factor w/ 3909 levels "0.000","0.000)",..: 1 116 841 1771
>...
>
>   I ran dput() on the data frame but cannot make sense of the output (a
>5.5M
>ASCII text file).
>
>   Pointers appreciated.
>
>Rich
>
>______________________________________________
>R-help at 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.



More information about the R-help mailing list