similar to: Accessing Results from cenmle function in NADA package

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

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

Hi, Has anyone ever used the censtats from the nada package and carried it out per group on a dataframe? I have it working on the entire dataframe but I need to do it by group. Thanks -- Le gach dea ghui, *Shane Carey* *GIS and Data Solutions Consultant* [[alternative HTML version deleted]]

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

Dear Users Regarding the NADA package, would anyone be able to help me understand what values are actually plotted on the Y axis of the plot obtained by using the *ros* function on the data and plotting the result with the plot() function? The Y axis is labeled "Values". According to the NADA user manual, ros performs a log transformation of the data by default, but the user can specify

I want to get access to the code of an S4 method in order to possibly modify a function to accomplish my particular needs: in my case the function in is cenfit() from the package NADA Now, given my reproducible example: my.ex<-structure(list(TEC = c(0.21, 0.077, 0.06, 0.033, 0.014, 0.0072), LR = c(0L, 0L, 1L, 0L, 1L, 0L)), .Names = c("TEC", "LR"), class =

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

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

Hello Instead of reading in group1.txt I want to read in groups1 for the first iteration of i, then groups2 for the second and so on. Obviously I can't use groups(i) but assume there is a way to do this. group<-read.table("C:/Data/April 2005/group1.txt",header=T) thanks in advance Meredith

Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide

I have some questions regarding Zero Inflation Poisson models. I am using count data to analyze abundance trends of salamanders. However, I have surveys which differ in the amount of effort (i.e. the number of people searching and amount of time - I am using a museum database so not all surveys were conducted by me). Therefore I need to account for the effort. If change the count (response

Hi Ansely, As Jim points out we really need some sample data to go with the code. Have a look at ?dput which is the best way to supply sample data here or have a look at http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example and/or http://adv-r.had.co.nz/Reproducibility.html for some general suggestions on asking questions here---including discussions of using

Hello, *Error in tx %*% comb : non-conformable arguments* Suggestions greatly appreciated. I am a beginner and this is my first time posting. I would like to get the summary for indicator species analysis, using package indicspecies with multipatt. I am getting errors, I believe, do to my data organization. After reorganizing and reorganizing, nothing has helped. > data<-

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................

Hi all, I try to represent a multiple curve graphic where the x-axis is the temperature and the different y-axes are the different X (X22,X43,X44...) some X corresponds to the same molecule (22 and 44 are for CO2 for instance) so I use the same colour for them. I wanna mix the linetype with the colour to be able to visually see the difference between X43 and X45 The best I have done up to now

All, I have some code that I have been working on that crashes R (1.8.1 on MacOS X installed from source). Here are excerpts of the code and data structures. When the lists that I am using to construct the data.frame are truncated (shortened), R catches the differing rows and warns. However, when the entire lists are used (I haven't looked for an exact breakpoint or for a particular list

Hi, I'm attempting to add a "Julian Day" column to a data frame. Here is my code and the resulting data frame: vic.data <- read.table("C:/VIC/data/vic.data.csv", header=F) names(vic.data) <- c("year", "month", "day", "precip", "evap", "runoff", "baseflow", "Tsup",

Hi, I wonder if anyone can answer the following or point me in the direction of how to obtain answers to the questions. Below is Output from R and further down are the questions raised and explanation of the study. Output from R: glm(formula = CB95TO00 ~ URB + INC, family = poisson) Deviance Residuals: Min 1Q Median 3Q Max -1.2272 -1.1290 0.2709 0.4272 2.1376

Dear R-helpers, When I try > anova(txtE2.lme, txtE2.lme1) Model df AIC BIC logLik Test L.Ratio p-value txtE2.lme 1 10 8590 8638 -4285 txtE2.lme1 2 7 8591 8624 -4288 1 vs 2 6.79 0.0789 > latex(anova(txtE2.lme, txtE2.lme1)) Error: object "n.group" not found I don't even see n.group as one of the arguments of latex() I checked to see >

I would be very grateful to anyone who could point to the error of my ways in the following. I have a dataframe called net1, as such: > str(net1) `data.frame': 114192 obs. of 9 variables: $ server : Factor w/ 122 levels "AB93-99","AMP93-1",..: 1 1 1 1 1 1 1 1 1 1 ... $ ts :'POSIXct', format: chr "2006-06-30 12:31:44"