2017-09-04 32 views
1

포아송 로그 로그 회귀 모델을 실행하는 데 문제가 있습니다. 이 경고 메시지를 어떻게 방지해야합니까? 내가 처음부터 포아송을 사용하기 때문에 어떻게하는지 모릅니다. 덕분에 많은포아송 모델 문제 - 정수 번호

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 
summary(sardegnalog.lm) 
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

좋아, 이제 해냈어! 정말 고마워 :) –