similar to: suggestions for nls error: false convergence

Hi, I'm trying to fit some data using a logistic function defined as y ~ a * (1+m*exp(-x/tau)) / (1+n*exp(-x/tau) My data is below: x <- 1:100 y <- c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 0,0,0,0,0,1,1,1,2,2,2,2,2,3,4,4,4,5, 5,5,5,6,6,6,6,6,8,8,9,9,10,13,14,16,19,21, 24,28,33,40,42,44,50,54,69,70,93,96,110,127,127,141,157,169,

Dear list, I do have a problem with nls. I use the following data: >test time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1 1.75 18.00 1 2.00 17.00 1 2.50 13.90 1 3.00 11.20 1 3.50 9.90 1 4.00 4.70 1 5.00 5.00 1 6.00 1.90 1 7.00 1.90 1 9.00 1.10 1 12.00 0.95 1 14.00

A couple of weeks ago there was a question regarding apparent convergence in nls when using the SSfol selfStart model for fitting a first-order pharmacokinetic model. I can't manage to find the original message either in my archive or in the list archives but the data were time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1

Dear Colleagues, Our group is also working on implementing the use of R for pharmacokinetic compartmental analysis. Perhaps I have missed something, but > fit <- nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=0.5, k2=0.5)), + data=C1.lsoda, + start=list(K1=0.3, k2=0.7), + trace=T + ) Error in eval(as.name(varName), data) : Object

Dear Friends. I tried to use nls.control() to change the 'minFactor' in nls( ), but it does not seem to work. I used nls( ) function and encountered error message "step factor 0.000488281 reduced below 'minFactor' of 0.000976563". I then tried the following: 1) Put "nls.control(minFactor = 1/(4096*128))" inside the brackets of nls, but the same error message

Hi, Parameters assessment in R with nls doesn't work, though it works fine with MS Excel with the internal solver :( I use nls in R to determine two parameters (a,b) from experimental data. m V C0 Ce Qe 1 0.0911 0.0021740 3987.581 27.11637 94.51206 2 0.0911 0.0021740 3987.581 27.41915 94.50484 3 0.0911 0.0021740 3987.581 27.89362

I am relatively new to R and am looking for advice, ideas or both... I have a data set that consists of pathogen population sizes on individual plant units in an experimental field plot. However, in order to estimate the pathogen population sizes I had to destroy the plant unit and could not determine if that plant unit became diseased or to what extent it would have become diseased. I

Hi I want to change a control parameter for an nls () as I am getting an error message "step factor 0.000488281 reduced below 'minFactor' of 0.000976562". Despite all tries, it seems that the control parameter of the nls, does not seem to get handed down to the function itself, or the error message is using a different one. Below system info and an example highlighting the

Hi, I have some data I want to fit with a non-linear function using nls, but it won't solve. > regresjon<-nls(lcfu~lN0+log10(1-(1-10^(k*t))^m), data=cfu_data, > start=(list(lN0 = 7.6, k = -0.08, m = 2))) Error in nls(lcfu ~ lN0 + log10(1 - (1 - 10^(k * t))^m), data = cfu_data, : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 Tried to increase minFactor

All - I have data: TL age 388 4 418 4 438 4 428 5 539 10 432 4 444 7 421 4 438 4 419 4 463 6 423 4 ... [truncated] and I'm trying to fit a simple Von Bertalanffy growth curve with program: #Creates a Von Bertalanffy growth model VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4, start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE) #Scatterplot of the data plot(TL~age, data=box5.4,

Hello, I've seen various threads on people reporting: step factor 0.000488281 reduced below `minFactor' of 0.000976563 While I know how to set the minFactor, what I'd like to have happen is for nls to return to me, the last or closest fitted parameters before it errors out. In other words, so I don't get convergence, I'd still like to acquire the values of the parameters

Hello, I am trying to estimate 4 parameters of a non-linear model using nls. My model function is a Fourier integral and is very expensive to calculate. I get the following error: > theta0 <- c(0.045, 1.02*10^(-4), 0.00169, 5.67*10^(-4)) > res <- nls(log(y) ~ log(model(theta,r,t)), data=dataModel, + start=list(theta=theta0), trace=TRUE, + control=nls.control(tol=1e-2))

Hello, I have been successfully using nls to fit a non-linear, self-limiting function to several sets of data collected in 2010 (example found below). To generate confidence intervals for parameter estimates, I've been attempting to bootstrap my sample. Unfortunately, I have meet with singularity gradients and failures to converge. I have altered maxiter and minFactor statements, used

Hi I?m trying to fit a nonlinear model to a derivative of the logistic function y = a/(1+exp((b-x)/c)) (this is the parametrization for the SSlogis function with nls) The derivative calculated with D function is: > logis<- expression(a/(1+exp((b-x)/c))) > D(logis, "x") a * (exp((b - x)/c) * (1/c))/(1 + exp((b - x)/c))^2 So I enter this expression in the nls function:

Hello. I am trying to fit an asymptotic relationship (nonlinear) to some ecological data, and am having problems. I am interested in the upper bound on the data (i.e. if there is an upper limit to 'y' across a range of 'x'). As such, I am using the nonlinear quantile regression package (nlrq) to fit a michaelis mention type model. The errors I get (which are dependant on

Hello colleagues, I am attempting to determine the nonlinear least-squares estimates of the nonlinear model parameters using nls. I have come across a common problem that R users have reported when I attempt to fit a particular 3-parameter nonlinear function to my dataset: Error in nls(r ~ tlm(a, N.fix, k, theta), data = tlm.data, start = list(a = a.st, : step factor 0.000488281

Hi R experts, I'm trying to use nls() for a piecewise linear regression with the first slope constrained to 0. There are 10 data points and when it does converge the second slope is almost always over estimated for some reason. I have many sets of these 10-point datasets that I need to do. The following segment of code is an example, and sorry for the overly precise numbers, they are just

I'm trying to use nls() to fit an exponential decay with an unknown offset in the time (independent variable). (Perhaps this is inherently very difficult?). > decay.pl <- nls (amp ~ expn(b0,b1,tau,t0,t), data = decay, + start = c(b0=1, b1=7.5, tau=3.5, t0=0.1), trace=T) Error in nlsModel(formula, mf, start) : singular gradient matrix at initial parameter estimates

This data are kilojoules of energy that are consumed in starving fish over a time period (Days). The KJ reach a lower asymptote and level off and I would like to use a non-linear plot to show this leveling off. The data are noisy and the sample sizes not the largest. I have tried selfstarting weibull curves and tried the following, both end with errors. Days<-c(12, 12, 12, 12, 22, 22, 22,

I wish to fit a dynamical model in R and I am running in a problem that requires some of your wisdom to solve. For SAS users I am searching for the equivalent of the retain statement. For people that want to read complicated explanations to help me: I have a system of two equations written as difference equations here. To boil it down. I have a dataframe with three variables y, X1, X2 which