Bom dia a todos

Estou trabalhando em um modelo não-linear e necessito de calcular, para efeitos de programação da matriz de informação de fisher, o vetor gradiente e a hessiana de uma função de ligação não-linear do tipo eta(x,beta).
Por exemplo, eta(x,beta)=beta[1]/(1+beta[2]*exp(-beta[3]*x)).

Já tentei utilizar os pacotes/funções nls, optim, deriv e numericDeriv. Tem o pacote deSolve, mas não se dá pra utilizar.

Não preciso da solução analítica, mas sempre é bom, né? rssss

Att,
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
x Prof. Clécio da Silva Ferreira x
x Depto Estatística - UFJF x
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx