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 : Sélectionner tout - Visualiser dans une fenêtre à part
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 : Sélectionner tout - Visualiser dans une fenêtre à part
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 ?