[R] drop1: warning message
Martin Bulla
bulla at centrum.cz
Tue Jan 12 08:17:56 CET 2010
Dear R collegues,
I am new to R and do not understad the warning messages, which appear when I am trying to apply type III procedure to my glm:
> inkm$Gamie glm.inkmZ drop1(glm.inkmZ,test="Ch")
Single term deletions
Model:
InkMb ~ IncStart + Vol + Gamie + PC1 + PC2 + PC3 + SpotsN + PC1:SpotsN
Df Deviance AIC LRT Pr(Chi)
9.3879 27.388
IncStart 1 14.6274 30.627 5.2395 0.0220791 *
Vol 1 15.9723 31.972 6.5844 0.0102876 *
Gamie 1 13.8659 29.866 4.4780 0.0343330 *
PC2 1 9.6899 25.690 0.3020 0.5826585
PC3 1 10.7326 26.733 1.3447 0.2462124
PC1:SpotsN 1 21.6517 37.652 12.2638 0.0004618 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Warning messages:
1: In glm.fit(x[, jj, drop = FALSE], y, wt, offset = object$offset, :
fitted probabilities numerically 0 or 1 occurred
2: In glm.fit(x[, jj, drop = FALSE], y, wt, offset = object$offset, :
fitted probabilities numerically 0 or 1 occurred
Where did PC1 disappiered? What is the meaning of these warning messages?
Similarly, I do not understand where did PCs dissapeared in the following calculations:
> inkm$Gamie glm.inkm0 drop1(glm.inkm0,test="Ch")
Single term deletions
Model:
InkMb ~ IncStart + Vol + Gamie + PC1 + PC2 + PC3 + SpotsN + Gamie:PC1 +
Gamie:PC2 + Gamie:PC3
Df Deviance AIC LRT Pr(Chi)
18.368 40.368
IncStart 1 26.078 46.078 7.7106 0.005490 **
Vol 1 26.515 46.515 8.1475 0.004312 **
SpotsN 1 18.663 38.664 0.2958 0.586537
Gamie:PC1 1 20.418 40.418 2.0505 0.152155
Gamie:PC2 1 18.460 38.460 0.0919 0.761818
Gamie:PC3 1 19.832 39.832 1.4647 0.226180
---
Any suggestions would help me greatly.
Best regards,
Martin
More information about the R-help
mailing list