search for: minfactor

Displaying 20 results from an estimated 25 matches for "minfactor".

Did you mean: mfactor
2003 Aug 28
3
(no subject)
Dear All, A couple of questions about the nls package. 1. I'm trying to run a nonlinear least squares regression but the routine gives me the following error message: step factor 0.000488281 reduced below `minFactor' of 0.000976563 even though I previously wrote the following command: nls.control(minFactor = 1/4096), which should set the minFactor to a lower level than the default one, 1/1024=0.000976563. Is there any way of setting the new minfactor to a lower level? 2. Is it possible to set s...
2011 Mar 28
1
error in nls, step factor reduced below minFactor
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 last tried...
2013 Oct 03
2
SSweibull() : problems with step factor and singular gradient
...r messages: 1) # Example days <- c(163,168,170,175,177,182,185,189,196,203,211,217,224) height <- c(153,161,171,173,176,173,185,192,195,187,195,203,201) dat <- as.data.frame(cbind(days,height)) fit <- nls(y ~ SSweibull(x, Asym, Drop, lrc, pwr), data = dat, trace=T, control=nls.control(minFactor=1/100000)) Error in nls(y ~cbind(1, -exp(-exp(lrc)* x^pwr)), data = xy, algorithm = “plinear”, :                           step factor 0.000488281 reduced below `minFactor` of 0.000976562 I tried to avoid this error by reducing the step factor below the standard minFactor of 1/1024 usin...
2012 Jan 03
1
nls and rbinom function: step factor 0.000488281 reduced below 'minFactor' of 0.000976562
...           71            76           102  z <- nls(y ~ rbinom(1000,1,a+b*x),data=d,start= list(a =0.1,b=0.2),trace=T); 374 :  0.1 0.2  361 :  0.1 0.2  350 :  0.1 0.2  Error in nls(y ~ rbinom(1000, 1, a + b * x), data = d, start = list(a = 0.1,  :    step factor 0.000488281 reduced below 'minFactor' of 0.000976562 I have tried plinear, got the same.... Additionaly- why does the parameters not change with the iteration ? [[alternative HTML version deleted]]
2001 Aug 03
0
step factor below minimum
I am trying to fit the 4-parameter logistic model with SSfpl. The estimation is terminated with a message "step factor 0.000488281 reduced below 'minFactor' of 0.000976562". I tried to change the minFactor with nls.control(minFactor=yy) with no apparent success. A call of nls.control() without arguments always shows the value of 0.000976562. Also the estimation terminates as previously although I tried to use control=nls.control(minFactor=yy...
2007 Apr 26
1
gnls warning message
...NULL in: control$nlmStepMax What is new for me is the warning message. It appears when I try to specify different settings using the control statement.For example: > hm02<-gnls(ht~a1*hd*(1-a2*exp(-a3*(dbh/dq2))),data=hdat, + start=c(a1=1.04,a2=1.38,a3=2.88),weights=varExp(), + control=c(minFactor=1/2048)) Any comment will be appreciate Thanks PM
2007 Feb 07
1
Singular Gradient
I tried to fit data with the following function: fit<-nls(y~ Is*(1-exp(-l*x))+Iph,start=list(Is=-2e-5,l=2.3,Iph=-0.3 ),control=list(maxiter=500,minFactor=1/10000,tol=10e-05),trace=TRUE) But I get only a singular Gradient warning... the data can by found attached(there are two sampels of data col 1/2 and 3/4). I tried to fix it by chanching the start parameters but that didn't solve the problem. Would it be a possibiliti to use the s...
2010 Apr 19
3
nls for piecewise linear regression not converging to least square
....70471269, 21.92956392, 36.39401717, 32.43620195, 44.77442277) x1<-c(24.6, 28.9, 33.2, 37.6, 42.0, 46.4, 50.9, 55.3, 59.8, 64.3) dat <- data.frame(x1,y1) nlmod <- nls(y1 ~ ifelse(x1 < xint+(yint/slp), yint, yint + (x1-(xint+(yint/slp)))*slp), data=dat, control=list(minFactor=1e-5,maxiter=500,warnOnly=T), start=list(xint=39.27464924, yint=0.09795576, slp=2.15061064), na.action=na.omit, trace=T) ##plotting the function plot(dat$x1,dat$y1) segments(x0=0, x1=coef(nlmod)[1]+coef(nlmod)[2]*coef(nlmod)[3], y0=coef(nlmod)[2], y1=coef(nlmod)...
2007 Jun 14
0
nlsList problems: control option does not effect output and strange environment search
...the start, one at the end and three randomly in between. Now I get the following error messages: testnlslist2<-nlsList(SSfpl,groupedData(Y~Time|group, datset)) Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xy, : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 Error in nls(formula = formula, data = data, control = control) : singular gradient For the first error it could be that this is due to the fact that this group is simulated as a kind of placebo and therefore it may be hard to fit. To solve this I tried to specify cont...
2011 Jun 12
1
Error in NLS example in the documentation
...(*), data, true function and fit, n=100") curve(100 + 10 * exp(x / 2), col=4, add = TRUE) lines(x, predict(nlmod), col=2) However, when I try to execute this - I get the following error- Error in nls(y ~ Const + A * exp(B * x), trace = TRUE) : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 In addition: Warning message: In nls(y ~ Const + A * exp(B * x), trace = TRUE) : No starting values specified for some parameters Can someone propose how I can fix this error? Thanks. Diviya [[alternative HTML version deleted]]
2013 Oct 27
1
nls function error
data(Boston, package='MASS') y <- Boston$nox x <- Boston$dis nls(y~ A + B * exp(C * x), start=list(A=1, B=1, C=1)) Error in nls(y ~ A + B * exp(C * x), start = list(A = 1, B = 1, C = 1), : step factor 0.000488281 reduced below 'minFactor' of 0.000976562 I don't know how to fix this error. I think my problem is that I set the wrong start. Could somebody help please? [[alternative HTML version deleted]]
2013 May 28
6
Ocultar componentes de una lista
...ist(A = 5, B = 1), algorithm = "default", control = structure(list(maxiter = 50, tol = 1e-05, ... $ dataClasses: Named chr "numeric" ..- attr(*, "names")= chr "Samples" $ control :List of 5 ..$ maxiter : num 50 ..$ tol : num 1e-05 ..$ minFactor: num 0.000977 ..$ printEval: logi FALSE ..$ warnOnly : logi FALSE - attr(*, "class")= chr "nls" Como se observa muchos de estos elmentos no son mostrados en el resultado de la funcion a menos que se use la funcion "str" , quisiera saber como hacer esto, ocult...
2007 Oct 01
1
[nls] singular gradient
...2. 8,3.3,3.4,3.6,5.0) # set the data frame data1<- data.frame(cbind(N,CSR,ev)) # the initial values of parameters para1.st <- c(a=75.4,b=165.9,c=-22.4,d=0,e=0.8,f=-0.28,g=0) para2.st <- c(a=31.3,b=-176.9,c=53.41,d=0,e=75.2,f=-0.18,g=0) # call nls funciton try.control <- c(maxiter=100, minFactor=1/4096) out <- nls(N~CSR/(1/(a+b*ev+c*(ev^2))+CSR/(d+e*exp(f*ev)))-g, data1, start=para1.st, control=try.control, trace=T) The data is from experiments and the pattern of data is scatter quite a bit. I want to find the best fit coefficients (a,b,c,d,e,f,g) for the data. I...
2006 Sep 15
1
Formula aruguments with NLS and model.frame()
...argument to modify garchFit() to handle other variables in the mean equation besides just an arma(u,v) specification. My nonlinear model is y<-nls(t~a*sin(w*2*pi/365*id+p)+b*id+int,data=t1, start=list(w=.5,a=.1,p=.5,b=init.y$coef[2],int=init.y$coef[1] ), control=list(maxiter=1000000,minFactor=1e-18)) where t is change in daily temperatures, id is just a time trend and the a*sin is a one year fourier series. I have tried to debug the nls code using the following code t1<-data.frame(t=as.vector(x),id=index(x)) data=t1; formula <- as.formula(t ~ a *sin(w *2* pi/365 * id + p) + b...
2013 May 28
3
Ocultar componentes de una lista
...ist(A = 5, B = 1), algorithm = "default", control = structure(list(maxiter = 50, tol = 1e-05, ... $ dataClasses: Named chr "numeric" ..- attr(*, "names")= chr "Samples" $ control :List of 5 ..$ maxiter : num 50 ..$ tol : num 1e-05 ..$ minFactor: num 0.000977 ..$ printEval: logi FALSE ..$ warnOnly : logi FALSE - attr(*, "class")= chr "nls" Como se observa muchos de estos elmentos no son mostrados en el resultado de la funcion a menos que se use la funcion "str" , quisiera saber como hacer esto, oculta...
2013 May 28
3
Ocultar componentes de una lista
...list(A = 5, B = 1), algorithm = "default", control = structure(list(maxiter = 50, tol = 1e-05, ... $ dataClasses: Named chr "numeric" ..- attr(*, "names")= chr "Samples" $ control :List of 5 ..$ maxiter : num 50 ..$ tol : num 1e-05 ..$ minFactor: num 0.000977 ..$ printEval: logi FALSE ..$ warnOnly : logi FALSE - attr(*, "class")= chr "nls" Como se observa muchos de estos elmentos no son mostrados en el resultado de la funcion a menos que se use la funcion "str" , quisiera saber como hacer esto, ocultar...
2002 Mar 19
0
MASS Book Exercise
...* p["th"] c(p[c("b0", "b1")], d0 = as.vector(d0)) } However when running the above codes in R, it stopped at wtloss.gr with the error: Error in nls(Weight ~ expn(b0, b1, th, Days), data = wtloss, start = wtloss.st, : step factor 0.000488281 reduced below `minFactor' of 0.000976562 Execution halted Any hints on where I did wrong would be greatly appreciated, Ko-Kang Wang ------------------------------------------------------------------------------ Ko-Kang Kevin Wang Postgraduate PGDipSci Student Department of Statistics University of Auckland New Zeala...
2004 May 27
0
specifying starting values in nlsList/conversion problems
...guess, resulting in typical error messages as known from nls: Error in nls(formula = formula, data = data, start = start, control = control) : singular gradient Error in nls(formula = formula, data = data, start = start, control = control) : step factor 0.000488281 reduced below `minFactor' of 0.000976563 Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an Infinity produced when evaluating the model is there a possibility to specifiy group specific starting or an algorithm which tests a matrix of starting with a certain stepwidth to overcome the conv...
2011 Oct 24
1
nonlinear model
...0,16.633333,15.666667 ,17.883333,NaN) When using the nls command I recieve an error like this: >mod = nls(A~ Amax * (1-exp(-b*Pot)),start=list(Amax=10,b=0.003)) / / Error in nls(A~ Amax * (1 - exp(-b * Pot), start = list(Amax = 10, : / step factor 0.000488281 reduced below 'minFactor' of 0.000976562 /I suppose that the used model is wrong, specialy Amax. I studied somewhere I should use "plinear" model but I did not know. If you can guide me, any advices would be realy appreciated. Thank you in advance,/ Best regards, Mitra RAHMATI/ [[alternative...
2012 Oct 04
0
extract fit values from geom_smooth
...y", position = "stack", fill = 2 ) + geom_smooth( method = "nls", formula = y ~ (a/b)*exp(-(x-c)^2/(2*b^2)) + (d/e)*exp(-(x-f)^2/(2*e^2)), se=F, start=list(a=100, b=1, c=695, d=100, e=1, f=710), control=nls.control(tol=1E-5, minFactor=1/1024), n = 1000 ) Now I would like to extract some values for further calculations. These are: 1st: mean of the two gauss fits ( c and f ) 2nd: ?3 sigma integral of the two gaussians 3rd: minimum between the two gaussians How can I extract these values from the geom_smooth functi...