> picloram.mod1 <- glm(formula = death ~ dose1, family=binomial("logit")) > > summary(picloram.mod1) Call: glm(formula = death ~ dose1, family = binomial("logit")) Deviance Residuals: Min 1Q Median 3Q Max -2.20648 -0.23398 0.01978 0.42836 1.44988 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -3.5845 0.4656 -7.698 1.38e-14 *** dose1 2.6941 0.3278 8.220 < 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: 424.89 on 312 degrees of freedom Residual deviance: 165.54 on 311 degrees of freedom AIC: 169.54 Number of Fisher Scoring iterations: 7 > > anova(picloram.mod1, test="Chisq") Analysis of Deviance Table Model: binomial, link: logit Response: death Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 312 424.89 dose1 1 259.36 311 165.54 < 2.2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > picloram.mod2 <- glm(formula = death ~ dose1 + replic2 + replic3, family=binomial("logit")) > > summary(picloram.mod2) Call: glm(formula = death ~ dose1 + replic2 + replic3, family = binomial("logit")) Deviance Residuals: Min 1Q Median 3Q Max -2.38737 -0.21818 0.02105 0.34525 1.57076 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -3.1452 0.5317 -5.915 3.32e-09 *** dose1 2.6979 0.3302 8.170 3.08e-16 *** replic2 -0.5809 0.5080 -1.144 0.253 replic3 -0.7119 0.4745 -1.500 0.134 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 424.89 on 312 degrees of freedom Residual deviance: 163.06 on 309 degrees of freedom AIC: 171.06 Number of Fisher Scoring iterations: 7 > > anova(picloram.mod2, test="Chisq") Analysis of Deviance Table Model: binomial, link: logit Response: death Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 312 424.89 dose1 1 259.356 311 165.54 <2e-16 *** replic2 1 0.194 310 165.34 0.6594 replic3 1 2.281 309 163.06 0.1309 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 >