t = rnorm(10,5,3)
x1=8500
x2=0.0156
x3=0.26

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.