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 quet = rnorm(10,5,3)x1=8500x2=0.0156x3=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

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 quet = rnorm(10,5,3)x1=8500x2=0.0156x3=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.

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.

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.

Não por que você usou 10 valores sortiados pela rnorm() a matriz hessiana mostrada é a substituição destes valores na função derivada. Att. Tiago. Date: Mon, 15 Apr 2013 11:58:12 -0300 From: rcoster@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Matrix Hessiana - duvida 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 quet = rnorm(10,5,3)x1=8500x2=0.0156x3=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.
participantes (4)
-
Carlos Pereira
-
Elias Krainski
-
Rodrigo Coster
-
Tiago Souza Marçal