> ships1 shiptype constprd operprd aggsrvc dmgincd 1 1 1 1 127 0 2 1 1 2 63 0 3 1 2 1 1095 3 4 1 2 2 1095 4 5 1 3 1 1512 6 6 1 3 2 3353 18 7 1 4 2 2244 11 8 2 1 1 44882 39 9 2 1 2 17176 29 10 2 2 1 28609 58 11 2 2 2 20370 53 12 2 3 1 7064 12 13 2 3 2 13099 44 14 2 4 2 7117 18 15 3 1 1 1179 1 16 3 1 2 552 1 17 3 2 1 781 0 18 3 2 2 676 1 19 3 3 1 783 6 20 3 3 2 1948 2 21 3 4 2 274 1 22 4 1 1 251 0 23 4 1 2 105 0 24 4 2 1 288 0 25 4 2 2 192 0 26 4 3 1 349 2 27 4 3 2 1208 11 28 4 4 2 2051 4 29 5 1 1 45 0 30 5 2 1 789 7 31 5 2 2 437 7 32 5 3 1 1157 5 33 5 3 2 2161 12 34 5 4 2 542 1 > > ships1a.poireg <- glm(formula = dmgincd ~ offset(log(aggsrvc)) + shiptype + + constprd, family=poisson("log")) > > summary(ships1a.poireg) Call: glm(formula = dmgincd ~ offset(log(aggsrvc)) + shiptype + constprd, family = poisson("log")) Deviance Residuals: Min 1Q Median 3Q Max -1.9538 -1.0525 -0.4846 0.5779 2.4183 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -6.26892 0.21241 -29.513 < 2e-16 *** shiptype2 -0.55861 0.17767 -3.144 0.00167 ** shiptype3 -0.68388 0.32908 -2.078 0.03769 * shiptype4 -0.07032 0.29071 -0.242 0.80887 shiptype5 0.30425 0.23581 1.290 0.19698 constprd2 0.75347 0.14855 5.072 3.94e-07 *** constprd3 0.96530 0.16369 5.897 3.70e-09 *** constprd4 0.70828 0.22153 3.197 0.00139 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 146.328 on 33 degrees of freedom Residual deviance: 49.355 on 26 degrees of freedom AIC: 163.22 Number of Fisher Scoring iterations: 5 > > ships2a.poireg <- glm(formula = dmgincd ~ offset(log(aggsrvc)) + shiptype + + operprd, family=poisson("log")) > > summary(ships2a.poireg) Call: glm(formula = dmgincd ~ offset(log(aggsrvc)) + shiptype + operprd, family = poisson("log")) Deviance Residuals: Min 1Q Median 3Q Max -3.4278 -1.0990 -0.3862 0.6206 3.2326 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -5.8023 0.1784 -32.519 < 2e-16 *** shiptype2 -0.7430 0.1691 -4.393 1.12e-05 *** shiptype3 -0.7545 0.3276 -2.303 0.0213 * shiptype4 -0.1845 0.2876 -0.642 0.5211 shiptype5 0.3885 0.2348 1.654 0.0980 . operprd2 0.5030 0.1115 4.511 6.46e-06 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 146.328 on 33 degrees of freedom Residual deviance: 70.103 on 28 degrees of freedom AIC: 179.97 Number of Fisher Scoring iterations: 5 > > > ships3a.poireg <- glm(formula = dmgincd ~ offset(log(aggsrvc)) + + constprd + operprd, family=poisson("log")) > > summary(ships3a.poireg) Call: glm(formula = dmgincd ~ offset(log(aggsrvc)) + constprd + operprd, family = poisson("log")) Deviance Residuals: Min 1Q Median 3Q Max -2.5072 -0.7776 -0.2107 0.9251 3.4695 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -6.9477 0.1269 -54.733 < 2e-16 *** constprd2 0.7542 0.1488 5.070 3.99e-07 *** constprd3 1.0509 0.1576 6.669 2.57e-11 *** constprd4 0.7041 0.2203 3.196 0.00139 ** operprd2 0.3875 0.1181 3.281 0.00104 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 146.328 on 33 degrees of freedom Residual deviance: 62.365 on 29 degrees of freedom AIC: 170.23 Number of Fisher Scoring iterations: 4 > > ships3.poireg <- glm(formula = dmgincd ~ offset(log(aggsrvc)) + + constprd + operprd + shiptype, family=poisson("log")) > > summary(ships3.poireg) Call: glm(formula = dmgincd ~ offset(log(aggsrvc)) + constprd + operprd + shiptype, family = poisson("log")) Deviance Residuals: Min 1Q Median 3Q Max -1.6768 -0.8293 -0.4370 0.5058 2.7912 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -6.40590 0.21744 -29.460 < 2e-16 *** constprd2 0.69714 0.14964 4.659 3.18e-06 *** constprd3 0.81843 0.16977 4.821 1.43e-06 *** constprd4 0.45343 0.23317 1.945 0.05182 . operprd2 0.38447 0.11827 3.251 0.00115 ** shiptype2 -0.54334 0.17759 -3.060 0.00222 ** shiptype3 -0.68740 0.32904 -2.089 0.03670 * shiptype4 -0.07596 0.29058 -0.261 0.79377 shiptype5 0.32558 0.23588 1.380 0.16750 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for poisson family taken to be 1) Null deviance: 146.328 on 33 degrees of freedom Residual deviance: 38.695 on 25 degrees of freedom AIC: 154.56 Number of Fisher Scoring iterations: 5 >