[R] nls.control()
Dieter Menne
dieter.menne at menne-biomed.de
Wed Sep 3 08:00:04 CEST 2008
Ranney, Steven <steven.ranney <at> montana.edu> writes:
>
> and I'm trying to fit a simple Von Bertalanffy growth curve with program:
>
> VonB=nls(TL~Linf*(1-exp(-k*(age-t0))), data=box5.4,
> start=list(Linf=1000, k=0.1, t0=0.1), trace=TRUE)
...
> Everything works as it should right up until the confint(VonB) statement.
When I ask
> for the confidence intervals, R goes through the process then stops and gives
>
> "Error in prof$getProfile() :
> step factor 0.000488281 reduced below 'minFactor' of 0.000976562".
>
> However, when I use nls.control=(minFactor=0.000001) OR
nls.control(minFactor=0.000001),
In my experience, it never helps to reduce the control factors once this error
comes up; it is a problem that by profiling you hit a boundary where
convergence problems are serious.
Suggestions: first, try with some smaller p-value, e.g. 0.8 or 0.7 instead of
the default of 0.95; profiling will stay in safer region.
It may be possible to use algorithm "port" in combination with lower/upper, but
I am not sure if this is honoured by profile.
Dieter
More information about the R-help
mailing list