[R-br] transferir modelo da função nls para a função curve() automaticamente

Éder Comunello comunello.eder em gmail.com
Sábado Fevereiro 27 09:16:52 BRT 2016


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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160227/ce7ddfcf/attachment.html>


Mais detalhes sobre a lista de discussão R-br