2017-09-04 57 views
1

我在使用泊松模型運行對數迴歸模型時遇到了問題。我應該如何防止這個警告信息?這也是我第一次使用泊松,所以我真的不知道該怎麼做。非常感謝Poisson模型的故障 - 整數

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 
+0

它給你一個警告,因爲泊松分佈只支持非負整數。 – James

+0

那麼,如果我沒有非負整數,我不能使用泊松分佈? –

+0

你應該更多地解釋你在做什麼(爲什麼你想使用泊松迴歸?)。這個問題可能更適合於[CrossValidated](http://stats.stackexchange.com) –

回答

0

這是很難確定,因爲你沒有足夠明確的有關模型,但它看起來像你需要簡單地砸在了log1p公式的左手邊; Poisson glm默認具有日誌鏈接功能(但您不需要添加1,因爲它是平均值,而不是將數據轉換爲線性預測器的比例)。

雖然你仍然需要它。

+0

好吧,我做到了,現在看來工作!萬分感謝 :) –