Syntaxe gradient dans maxLik

Je cherche à préciser la fonction de gradient d'une régression poisson avec offset.
J'ai déclaré ma fonction objectif (log vraisemblance) :

Code:

---

1 2 3 4 5

loglik <- function(param) { b<-param m=as.vector(of+z%*%b) ll <- sum(v*m-exp(m)) }

---

En csq, à partir d'une lm puis maxlik, je retrouve les estimateurs d'une glm avec family poisson donc c'est ok.

Ensuite, je veux préciser le gradient en ajoutant :

Code:

---

1 2 3 4

gradlik <- function(param) { b<-param m=as.vector(of+z%*%b) gg <- (v-exp(m))*z) }

---

Et là ca se corse car maxLik(loglik,grad=gradlik,start=sc) me donne rien.

Quelqu'un aurait une idée pour mieux écrire ce gradient ?