Realmente o script disponibilizado pelo Walmes é muito bom. É pra aprender o caminho das pedras!
Mas me surgiu uma ideia... Supondo que se queira apenas visualizar o ajuste do modelo, talvez bastasse utilizar predict() e depois adicionar os resultados ao gráfico.
### <code r>
head(Orange, 3)
mForm <- as.formula("circumference ~ A/(1+exp((B-age)/C))"); mForm; class(mForm)
mExpr <- mForm[[3]]; mExpr; class(mExpr) # A/(1 + exp((B - age)/C))
fm1 <- nls(mForm, Orange, start=list(A=160, B=700, C=350))
new <- seq(xLim[1], xLim[2], len=101)
pre1 <- predict(fm1, newdata=list(age=new))
plot(circumference~age, Orange)
lines(new, pre1, col=2)
# title(mForm)
title(mExpr)
### </code>