[R] gnlsControl

From: Daniel Coleman <coleman.daniel_at_gene.com>
Date: Fri 04 Aug 2006 - 11:47:20 EST


When I run gnls I get the error:  

Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xy, :

        step factor 0.000488281 reduced below 'minFactor' of 0.000976563  

My first thought was to decrease minFactor but gnlsControl does not contain minFactor nor nlsMinFactor (see below). It does however contain nlsMaxIter and nlsTol which I assume are the analogs of maxiter and tol in nls.control. I would be happy to hear from anyone who has an idea on what parameters in gnlsControl to change to get convergence.  

Cheers  

Dan Coleman

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, msScale = msScale, returnObject = returnObject,

        msVerbose = msVerbose, apVar = apVar, nlmStepMax = nlmStepMax,

        opt = match.arg(opt), optimMethod = optimMethod, .relStep = .relStep,

        minAbsParApVar = minAbsParApVar)

}  

> nls.control

$maxiter

[1] 50  

$tol

[1] 1e-05  

$minFactor

[1] 0.0009765625

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch 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. Received on Fri Aug 04 11:51:08 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 05 Aug 2006 - 00:17:49 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.