<div>Boa noite, </div>
<div> </div>
<div>Venho solicitar a ajuda de voces, pois estou tentando calcular o ponto de inflexão de um modelo logístico e não estou conseguindo, tenho realizado os seguintes comandos.</div>
<div> </div>
<div>> ajnll.gg1<-nls(PESO ~ A/(1+B*exp(-k*IDADE)),data=gg1.dat,start=c(A=900,B=1,k=0.01))<br>> s.ajnll.gg1<-summary(ajnll.gg1)<br>> s.ajnll.gg1</div>
<div>Formula: PESO ~ A/(1 + B * exp(-k * IDADE))</div>
<div>Parameters:<br> Estimate Std. Error t value Pr(>|t|) <br>A 9.284e+02 1.322e+01 70.202 < 2e-16 ***<br>B 1.801e+01 2.210e+00 8.148 2.68e-14 ***<br>k 2.826e-02 1.337e-03 21.140 < 2e-16 ***<br>---<br>
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </div>
<div>Residual standard error: 81.91 on 222 degrees of freedom</div>
<div>Number of iterations to convergence: 8 <br>Achieved convergence tolerance: 2.787e-06 </div>
<div> </div>
<div>> logi <- expression(A/(1+B*exp(-K*x)))<br>> D(logi,'x')<br>A * (B * (exp(-K * x) * K))/(1 + B * exp(-K * x))^2<br></div>
<div> </div>
<div>Substituindo na equação</div>
<div> > GGtc= 928.4*(18.01*(exp(-0.028*x)*0.028))/(1+18.01*exp(-0.028*x))^2<br>> GGtc<br>[1] 1.328331 1.361869 1.396139 1.431151 1.466914</div>
<div> </div>
<div> </div>