> cas.mod1 <- glm(probcaslty ~ 1 , binomial("logit"), weights=totsoldier, data=g1) > > summary(cas.mod1) Call: glm(formula = probcaslty ~ 1, family = binomial("logit"), data = g1, weights = totsoldier) Deviance Residuals: Min 1Q Median 3Q Max -78.916 -18.848 -2.277 15.314 75.238 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.1437072 0.0075644 -151.20 < 2.2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 18665.9 on 17 degrees of freedom Residual deviance: 18665.9 on 17 degrees of freedom AIC: 18802.5 Number of Fisher Scoring iterations: 5 > > cas.mod2 <- glm(probcaslty ~ 1 + officer, binomial("logit"), weights=totsoldier) > > summary(cas.mod2) Call: glm(formula = probcaslty ~ 1 + officer, family = binomial("logit"), weights = totsoldier) Deviance Residuals: Min 1Q Median 3Q Max -78.644 -18.682 -3.027 14.134 75.534 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.1502230 0.0078302 -146.8950 < 2.2e-16 *** officer 0.1000510 0.0303317 3.2986 0.0009718 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 18665.9 on 17 degrees of freedom Residual deviance: 18655.2 on 16 degrees of freedom AIC: 18793.8 Number of Fisher Scoring iterations: 5 > > cas.mod3 <- glm(probcaslty ~ 1+ officer + division, binomial("logit"), weights=totsoldier) > > summary(cas.mod3) Call: glm(formula = probcaslty ~ 1 + officer + division, family = binomial("logit"), weights = totsoldier) Deviance Residuals: Min 1Q Median 3Q Max -3.73404 -0.38461 -0.05373 1.02485 2.36496 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.421062 0.020565 20.4751 <2e-16 *** officer 0.049521 0.033560 1.4756 0.1400 division2 -1.092001 0.027636 -39.5136 <2e-16 *** division3 -1.029677 0.028009 -36.7617 <2e-16 *** division4 -1.976044 0.031170 -63.3964 <2e-16 *** division5 -4.571900 0.067933 -67.2998 <2e-16 *** division6 -0.886714 0.029091 -30.4808 <2e-16 *** division7 -2.361214 0.038404 -61.4835 <2e-16 *** division8 -2.949276 0.041053 -71.8409 <2e-16 *** division9 -2.676410 0.070603 -37.9079 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 18665.9088 on 17 degrees of freedom Residual deviance: 28.5736 on 8 degrees of freedom AIC: 183.187 Number of Fisher Scoring iterations: 4 > > cas.mod4 <- glm(probcaslty ~ 1+officer + division + officer:division, binomial("logit"), weights=totsoldier) > > summary(cas.mod4) Call: glm(formula = probcaslty ~ 1 + officer + division + officer:division, family = binomial("logit"), weights = totsoldier) Deviance Residuals: [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.443009 0.021234 20.8633 < 2.2e-16 *** officer -0.256025 0.078423 -3.2647 0.0010959 ** division2 -1.121497 0.028738 -39.0254 < 2.2e-16 *** division3 -1.058943 0.029089 -36.4031 < 2.2e-16 *** division4 -2.005011 0.032359 -61.9618 < 2.2e-16 *** division5 -4.591188 0.070325 -65.2857 < 2.2e-16 *** division6 -0.904857 0.030091 -30.0709 < 2.2e-16 *** division7 -2.375369 0.039678 -59.8662 < 2.2e-16 *** division8 -2.980073 0.042423 -70.2458 < 2.2e-16 *** division9 -2.734870 0.073106 -37.4094 < 2.2e-16 *** officer:division2 0.406079 0.105260 3.8579 0.0001144 *** officer:division3 0.409924 0.108259 3.7865 0.0001528 *** officer:division4 0.408276 0.120918 3.3765 0.0007342 *** officer:division5 0.266042 0.272290 0.9771 0.3285422 officer:division6 0.239187 0.118902 2.0116 0.0442589 * officer:division7 0.178787 0.158979 1.1246 0.2607588 officer:division8 0.462066 0.169112 2.7323 0.0062892 ** officer:division9 1.024391 0.287555 3.5624 0.0003675 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 1.86659e+04 on 17 degrees of freedom Residual deviance: -1.79323e-12 on 0 degrees of freedom AIC: 170.613 Number of Fisher Scoring iterations: 3 >