[R-br] Obtenção do ponto de inflexão em modelo logístico

Adriano Carvalho Costa acarvalhocosta em gmail.com
Sexta Janeiro 20 02:13:18 BRST 2012


Boa noite,

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.

> ajnll.gg1<-nls(PESO ~
A/(1+B*exp(-k*IDADE)),data=gg1.dat,start=c(A=900,B=1,k=0.01))
> s.ajnll.gg1<-summary(ajnll.gg1)
> s.ajnll.gg1
Formula: PESO ~ A/(1 + B * exp(-k * IDADE))
Parameters:
   Estimate Std. Error t value Pr(>|t|)
A 9.284e+02  1.322e+01  70.202  < 2e-16 ***
B 1.801e+01  2.210e+00   8.148 2.68e-14 ***
k 2.826e-02  1.337e-03  21.140  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 81.91 on 222 degrees of freedom
Number of iterations to convergence: 8
Achieved convergence tolerance: 2.787e-06

> logi <- expression(A/(1+B*exp(-K*x)))
>  D(logi,'x')
A * (B * (exp(-K * x) * K))/(1 + B * exp(-K * x))^2

Substituindo na equação
 > GGtc= 928.4*(18.01*(exp(-0.028*x)*0.028))/(1+18.01*exp(-0.028*x))^2
> GGtc
[1] 1.328331 1.361869 1.396139 1.431151 1.466914
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120120/5c15c0af/attachment.html>


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