Senhores, bom dia! 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> ================================================ Éder Comunello Agronomist (UEM), MSc in Environ. Sciences (UEM) DSc in Agricultural Systems Engineering (USP/Esalq) Brazilian Agricultural Research Corporation (Embrapa) Dourados, MS, Brazil |<O>| ================================================ GEO, -22.2752, -54.8182, 408m UTC-04:00 / DST: UTC-03:00