Problème avec le numéro de modèle de poisson-entier

Question

J'ai des problèmes avec l'exécution d'un modèle de régression log-log avec Poisson. Comment devrais-je empêcher ce message d'avertissement? C'est aussi la première fois que j'utilise Poisson donc je ne sais vraiment pas comment faire. Merci beaucoup

sardegnalog.lm <-glm(log1p(fulldata[381:400,1])~log1p(fulldata[381:400,2])+log1p(fulldata[381:400,3])+log1p(fulldata[381:400,4])+log1p(fulldata[381:400,8]), family="poisson") 
Warning messages: 
1: In dpois(y, mu, log = TRUE) : non-integer x = 8.868132 
2: In dpois(y, mu, log = TRUE) : non-integer x = 9.885069 
3: In dpois(y, mu, log = TRUE) : non-integer x = 9.410911 
4: In dpois(y, mu, log = TRUE) : non-integer x = 7.876259 
5: In dpois(y, mu, log = TRUE) : non-integer x = 11.826326 
6: In dpois(y, mu, log = TRUE) : non-integer x = 9.632728 
7: In dpois(y, mu, log = TRUE) : non-integer x = 9.872616 
8: In dpois(y, mu, log = TRUE) : non-integer x = 6.899723 
9: In dpois(y, mu, log = TRUE) : non-integer x = 9.027379 
10: In dpois(y, mu, log = TRUE) : non-integer x = 16.733528 

Call: 
glm(formula = log1p(fulldata[381:400, 1]) ~ log1p(fulldata[381:400, 
2]) + log1p(fulldata[381:400, 3]) + log1p(fulldata[381:400, 
4]) + log1p(fulldata[381:400, 8]), family = "poisson") 

Deviance Residuals: 
Min  1Q Median  3Q  Max 
-3.267 -2.082 -1.093 1.085 3.123 

Coefficients: 
         Estimate Std. Error z value Pr(>|z|)  
(Intercept)     -17.5129  5.2594 -3.330 0.000869 *** 
log1p(fulldata[381:400, 2]) 1.3144  0.4632 2.838 0.004544 ** 
log1p(fulldata[381:400, 3]) 0.7884  0.2384 3.307 0.000944 *** 
log1p(fulldata[381:400, 4]) -0.1477  0.2613 -0.565 0.571836  
log1p(fulldata[381:400, 8]) -0.7765  0.2960 -2.623 0.008715 ** 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for poisson family taken to be 1) 

Null deviance: 144.602 on 19 degrees of freedom 
Residual deviance: 80.231 on 15 degrees of freedom 
AIC: Inf 

Number of Fisher Scoring iterations: 6 

Answer 1

Il est difficile d'être certain que vous n'avez pas été suffisamment explicite sur le modèle, mais il vous ressemble besoin de déposer simplement le log1p sur la côté gauche de la formule; le glm de Poisson a déjà une fonction de lien de journal par défaut (mais vous n'avez pas besoin d'ajouter 1, car c'est la moyenne, pas les données qui sont transformées à l'échelle du prédicteur linéaire).

Vous en auriez toujours besoin sur la droite.