Walmes, aí está o exemplo citado. Preciso derivar, em relação a beta, uma e duas vezes. Uma daquelas funções daria para ser utilizada, se eu escrevesse cada componente de beta como argumento, mas eu gostaria de uma solução mais prática. não tenho como escrever o CMR mais exato!
nlf<-function(x,beta){
resp<- beta[1]/(1 +beta[2]*exp(-beta[3]*x))
return(resp)
}
Clécio
Você se refere a informação de fisher, gradiente e hessiana da função de verossimilhança? a rootSolve::gradient() pode ser útil. Alguns exemplos de como obter essas quantidades você pode ver no material recente que apresentamos no SINAPE
http://www.leg.ufpr.br/doku.php/cursos:mcie
Forneça um CMR para que possamos contribuir mais.
À 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
skype: walmeszeviani
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.