Muito bom, Walmes. Não havia tentado este pacote, dá para generalizar.<br><br>Abs,<br>Clécio<br><br><div class="gmail_quote">Em 20 de setembro de 2012 16:49, Walmes Zeviani <span dir="ltr"><<a href="mailto:walmeszeviani@gmail.com" target="_blank">walmeszeviani@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><font face="trebuchet ms,sans-serif">Veja o exemplo de como obter o gradiente e hessianos numericamente<br><br><span style="font-family:courier new,monospace">f <- function(beta, t){<br>
with(as.list(beta),<br> b1/(1+b2*exp(-b3*t)))<br>
}<br><br>f(c(b1=1, b2=0.5, b3=-0.1), t=1:3)<br><br>require(rootSolve)<br><br>help(gradient, , help_type="html")<br>gradient(f, x=c(b1=1, b2=0.5, b3=-0.1), t=1)<br>gradient(f, x=c(b1=1, b2=0.5, b3=-0.1), t=1:3)<br>
<br>help(hessian, help_type="html")<br>hessian(f, x=c(b1=1, b2=0.5, b3=-0.1), t=1)<br>hessian(f, x=c(b1=1, b2=0.5, b3=-0.1), t=1:2)</span><br><br>Dependendo das operações que vai fazer, é melhor organizar o hessiano em um array.<br>
É possível obter analiticamente usando a deriv3(). Para isso veja esse post<br><br><a href="http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-de-regressao-nao-linear/" target="_blank">http://ridiculas.wordpress.com/2011/05/19/bandas-de-confianca-para-modelo-de-regressao-nao-linear/</a><br>
<br>À disposição.<span class="HOEnZb"><font color="#888888"><br>Walmes.<br><br clear="all"></font></span></font><div class="HOEnZb"><div class="h5"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: <a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
</div></div><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><br>-- <br><br>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx<br>
x Prof. Clécio da Silva Ferreira x<br>x Depto Estatística - UFJF x<br>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx<br><br>