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 <[email protected]>
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: [email protected]
To: [email protected]
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: [email protected]
__________________________________________
*Tenha Um Dia Abençoado Por DEUS!!!***
* Sl -- 34, 8*
_______________________________________________ R-br mailing list
[email protected]
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
[email protected]
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
[email protected]
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
[email protected]
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.
