[R-br] Matrix Hessiana - duvida

[Elias Krainski]

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 em 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 em hotmail.com
>> To: r-br em 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 em hotmail.com
>> __________________________________________
>>
>> *Tenha Um Dia Abençoado Por DEUS!!!***
>> * Sl -- 34, 8*
>>
>>
>> _______________________________________________ R-br mailing list
>> R-br em 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 em 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 em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130415/074d492a/attachment.html>