search for: nlmstepmax

...ends on how the method (ML/REM) and which (and how much) parameters will depend randomly on the cluster-variable. How get the bist fit without convergence? I set the parameters msVerbose and returnObject to TRUE: lmeControl(maxIter=50000, msMaxIter=200, tolerance=1e-4, niter=50, msTol=1e-5, nlmStepMax=500, ,msVerbose=TRUE ,returnObject=TRUE ) However, the lme-functions does not produce verbose output, nor does it return the best fit if lme is not converging. It returns only an error: Error in lme.formula(y ~ lndbh + I(lndbh^2) + lnh + I(lnh^2), random = ~lndbh + : iteration limi...

...eman Genentech Inc. > gnlsControl function (maxIter = 50, nlsMaxIter = 7, msMaxIter = 50, minScale = 0.001, tolerance = 1e-06, nlsTol = 0.001, msTol = 1e-07, msScale = lmeScale, returnObject = FALSE, msVerbose = FALSE, apVar = TRUE, .relStep = (.Machine$double.eps)^(1/3), nlmStepMax = 100, opt = c("nlminb", "optim"), optimMethod = "BFGS", minAbsParApVar = 0.05) { list(maxIter = maxIter, nlsMaxIter = nlsMaxIter, msMaxIter = msMaxIter, minScale = minScale, tolerance = tolerance, nlsTol = nlsTol, msTol = msTol, msScal...

...inXP) and I got the following error and warning message: Error in gnls(ht ~ a1 * hd * (1 - a2 * exp(-a3 * (dbh/dq2))), data = hdat, : Step halving factor reduced below minimum in NLS step In addition: Warning message: $ operator is deprecated for atomic vectors, returning 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...

...`R CMD check nlme' fails on my machine. The final output in nlme-Ex.Rout is: > library(nlme) > data(Soybean) > fm1 <- nlme(weight ~ SSlogis(Time, Asym, xmid, scal), data = Soybean, + fixed = Asym + xmid + scal ~ 1, start = c(18, 52, 7.5), + control = list(nlmStepMax = 1.0)) Error in nlme.formula(weight ~ SSlogis(Time, Asym, xmid, scal), data = Soybean, : Maximum number of iterations reached without convergence Execution halted I presume that `R CMD check MASS' fails for some similar reason. The final output in MASS-Ex.Rout is: > library(nlme) >...

...Below are the defaults.. lmeControl function (maxIter = 50, msMaxIter = 50, tolerance = 1e-06, niterEM = 25, msTol = 1e-07, msScale = lmeScale, msVerbose = FALSE, returnObject = FALSE, gradHess = TRUE, apVar = TRUE, .relStep = (.Machine$double.eps)^(1/3), minAbsParApVar = 0.05, nlmStepMax = 100, optimMethod = "BFGS", natural = TRUE) I was reading on the R listserve that lmer from the lme4 package may be preferable to lme (for convergence problems) but lmer seems to need you to put in starting values and I'm not sure how to go about chosing them. I was wonder...

...)) Maximum number of PQL iterations: 200 iteration 1 iteration 2 Error in MEestimate(lmeSt, grps) : NA/NaN/Inf in foreign function call (arg 1) > gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~1),family="poisson", niterPQL=200,control=lmeControl(nlmStepMax=0.00000000000000000001)) Maximum number of PQL iterations: 200 iteration 1 iteration 2 Error in MEestimate(lmeSt, grps) : NA/NaN/Inf in foreign function call (arg 1) > gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~1),family="poisson", niterPQL=...

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems