<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Pessoal, boa tarde!</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Tenho interesse no tópico e fiz uns testes. Não sei se ajuda em algo, mas deixo o script pra consulta.</div><div class="gmail_default" style="font-family:verdana,sans-serif"><br></div><div class="gmail_default" style=""><div class="gmail_default"><font face="monospace, monospace">### <code r></font></div><div class="gmail_default"><font face="monospace, monospace">head(Orange, 3)</font></div><div class="gmail_default"><font face="monospace, monospace">mForm <- as.formula("circumference ~ A/(1+exp((B-age)/C))"); mForm; class(mForm)</font></div><div class="gmail_default"><font face="monospace, monospace">mExpr <- mForm[[3]]; mExpr; class(mExpr) </font></div><div class="gmail_default"><font face="monospace, monospace"># A/(1 + exp((B - age)/C))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">fm1  <- nls(mForm, Orange, start=list(A=160, B=700, C=350))</font></div><div class="gmail_default"><font face="monospace, monospace">mPar <- summary(fm1)$parameters[,1]</font></div><div class="gmail_default"><font face="monospace, monospace">xLim <- with(Orange, range(pretty(age)))</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">varx <- all.vars(mExpr)[!all.vars(mExpr) %in% names(mPar)]</font></div><div class="gmail_default"><font face="monospace, monospace">sapply(1:length(mPar), function(i) {assign(names(mPar)[i], mPar[i], envir=.GlobalEnv)})</font></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><span style="font-family:monospace,monospace">eval(call("curve", mExpr, xLim[1], xLim[2], xname=varx))</span><br></div><div class="gmail_default"><span style="font-family:monospace,monospace"># curve(A/(1 + exp((B - age)/C)), xLim[1], xLim[2], xname=varx)</span><span style="font-family:monospace,monospace"><br></span></div><div class="gmail_default"><font face="monospace, monospace"><br></font></div><div class="gmail_default"><font face="monospace, monospace">plot(circumference~age, Orange)</font></div><div class="gmail_default"><font face="monospace, monospace">eval(call("curve", mExpr, xLim[1], xLim[2], xname=varx, add=T, col=2))</font></div><div class="gmail_default"><font face="monospace, monospace">### </code></font></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><font face="arial, helvetica, sans-serif" style="font-size:small"><div style="font-family:'trebuchet ms',sans-serif;display:inline">​</div></font><font face="arial, helvetica, sans-serif" style="font-size:12.8px"><span style="font-size:small">================================================</span><br>Éder Comunello</font><div style="font-size:12.8px"><span style="font-family:arial,helvetica,sans-serif">Agronomist (UEM), MSc in Environ. Sciences (UEM)</span><font face="arial, helvetica, sans-serif"><br></font></div><div style="font-size:12.8px"><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">DSc in Agricultural Systems Engineering (USP/Esalq)</span><br></div><div><span style="font-size:small">Brazilian Agricultural Research Corporation (</span><font face="arial, helvetica, sans-serif" style="font-size:small">Embrapa)</font></div><div><div><font face="arial, helvetica, sans-serif" style="font-size:small">Dourados, MS, Brazil </font><font face="arial, helvetica, sans-serif" size="2">|<O>|</font></div><div style="font-size:small"><span style="font-family:arial,helvetica,sans-serif">================================================</span><br></div></div></div><div style="font-size:small"><font face="arial, helvetica, sans-serif">GEO, </font>-22.2752, -54.8182, 408m<span style="font-family:arial,helvetica,sans-serif"><br></span></div><div style="font-size:12.8px"><span style="font-size:small;font-family:arial,helvetica,sans-serif">UTC-04:00 / DST: UTC-03:00</span></div><div><div><div><br></div><div><br></div></div></div></div></div></div></div></div></div></div></div></div></div></div>