Re: [R-br] Obtenção do ponto de inflexão em modelo logístico [OFF TOPIC]
Walmes, esta função SSlogis() está em qual pacote? Não achei no MASS . Gostei desta parametrização, pode ser útil em um trabalho que estou desenvolvendo. Abs e obrigado Em 20 de janeiro de 2012 10:57, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Fernando A.B. Colugnati
Fernando, SSlogis é do pacote stats (já instalado e carregado em toda sessão R) ?SSlogis Exitem mais funções do tipo SS (SelfStart) no pacote stats, veja
apropos("^SS") [1] "SSasymp" "SSasympOff" "SSasympOrig" "SSbiexp" "SSD" [6] "SSfol" "SSfpl" "SSgompertz" "SSlogis" "SSmicmen" [11] "SSweibull"
Esse número aumenta se você carregar o pacote drc (dose response curve), nlstools, enfim, vários pacotes contem mais funções SS. À 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Obrigado! Em 20 de janeiro de 2012 11:27, Walmes Zeviani <walmeszeviani@gmail.com>escreveu:
Fernando,
SSlogis é do pacote stats (já instalado e carregado em toda sessão R)
?SSlogis
Exitem mais funções do tipo SS (SelfStart) no pacote stats, veja
apropos("^SS") [1] "SSasymp" "SSasympOff" "SSasympOrig" "SSbiexp" "SSD" [6] "SSfol" "SSfpl" "SSgompertz" "SSlogis" "SSmicmen" [11] "SSweibull"
Esse número aumenta se você carregar o pacote drc (dose response curve), nlstools, enfim, vários pacotes contem mais funções SS.
À 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Fernando A.B. Colugnati
find("SSlogis") [1] "package:stats"
On Fri, 20 Jan 2012, Fernando Colugnati wrote:
Walmes, esta função SSlogis() está em qual pacote? Não achei no MASS . Gostei desta parametrização, pode ser útil em um trabalho que estou desenvolvendo.Abs e obrigado
Em 20 de janeiro de 2012 10:57, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: 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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Fernando A.B. Colugnati
participantes (3)
-
Fernando Colugnati -
Paulo Justiniano -
Walmes Zeviani