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

[Walmes Zeviani]

Adriano,

Aplicação sequencial das regras do cálculo. Derive f(x) em relação à x
obtendo f'(x), depois derive f'(x) em relação à x obtendo f''(x), Iguale
f''(x) à zero e resolva para x. Isso é um excelente exercício, acredite.
Você pode usar funções que encontram raízes para ajudar na tarefa de
encontrar a raíz de f''(x), veja as funções do pacote rootSolve. Veja as
funções do pacote MASS, lá há uma função para cálculo dos "tempos de vida"
da logística (no caso glm), o tempo de 50% de vida é a inflexão. Seu código
não é reproduzível.

Outra opção é usar essa parametrização da logística

f(x) = Asym/(1+exp(-(x-xmid)/Scal))

que a paremetrização do R na função SSlogis(), que possui até selfStart, o
parâmetro xmid já é o tempo de meia vida. Estudos do qual participei
apontam que essa é a melhor parametrização do modelo logístico para
inferência estatística (contornos elípticos ortogonais para regiões de
confiança, intervalos simétricos, baixa não linearidade e viés).

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120120/98a962e7/attachment.html>