[R] nls, convergence and starting values
Christian Ritz
ritz at life.ku.dk
Mon Mar 30 15:48:35 CEST 2009
Hi Patrick,
there exist specialized functionality in R that offer both automated calculation of
starting values and relatively robust optimization, which can be used with success in many
common cases of nonlinear regression, also for your data:
library(drc) # on CRAN
## Fitting 3-parameter logistic model
## (slightly different parameterization from SSlogis())
bdd.m1 <- drm(pourcma~transat, weights=sqrt(nbfeces), data=bdd, fct=L.3())
plot(bdd.m1, broken=TRUE, conLevel=0.0001)
summary(bdd.m1)
Of course, standard errors are huge as the data do not really support this model (as
already pointed out by other replies to this post).
Christian
More information about the R-help
mailing list