Gosto de usar expressoes nesses casos (obtendo assim o array) d3 <- deriv3(expression(x1/((1+((x1/x2)**(x3*t)-1)))),c("x1","x2","x3")) t = rnorm(10,5,3) x1=8500 x2=0.0156 x3=0.26 eval(d3) Em 15-04-2013 11:58, Rodrigo Coster escreveu:
A hessiana nao deveria ter retornado um array 3x3x3?
2013/4/13 Tiago Souza Marçal <tiagosouzamarcal@hotmail.com>
Carlos tente os seguintes comandos:
d3 <- deriv3(~x1/((1+((x1/x2)**(x3*t)-1))),c("x1","x2","x3"), function(t,x1,x2,x3){ NULL })
t = rnorm(10,5,3)
x1=8500
x2=0.0156
x3=0.26
hessiana<-d3(t,x1,x2,x3)
Espero ter contribuido
Att.
Tiago. ------------------------------ From: ccpsilva81@hotmail.com To: r-br@listas.c3sl.ufpr.br Date: Sat, 13 Apr 2013 05:31:38 +0300 Subject: [R-br] Matrix Hessiana - duvida
Oi pessoal to com um problema para calcular a matriz hessiana do seguinte modelo
modelo = x1/((1+((x1/x2)**(x3*t) - 1)))
sendo que t = rnorm(10,5,3) x1=8500 x2=0.0156 x3=0.26
Carlos Pereira da Silva Bacharel em Estatística -- UEPB
Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81@hotmail.com __________________________________________
*Tenha Um Dia Abençoado Por DEUS!!!*** * Sl -- 34, 8*
_______________________________________________ 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.
_______________________________________________ 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.
_______________________________________________ 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.