copepod.lst<-scan(what=list(stage="",dose=0,y=0,nsurv=0))
 A    0    23    25    C    0    25    25    N    0    19    20
I other data omitted for space considerations I
 A   129    0    25                          N    84    1    20

 copepod.df<-data.frame(copepod.lst)

 copepod.df$iA<-as.numeric(copepod.df$stage=="A")
 copepod.df$iC<-as.numeric(copepod.df$stage=="C")

 copepod.fit1<-glm(cbind(y,nsurv-y)~ I(dose-mean(dose)) + iA + iC,
               family=binomial(link=logit),data=copepod.df)

 copepod.fit2<-glm(cbind(y,nsurv-y)~ I(dose-mean(dose))*(iA + iC),
               family=binomial(link=logit),data=copepod.df)

 copepod.fit3<-glm(cbind(y,nsurv-y)~ (I(dose-mean(dose)) +
               I((dose-mean(dose))^2)) *(iA + iC),
               family=binomial(link=logit),data=copepod.df)

 anova(copepod.fit1,copepod.fit2,copepod.fit3,test="Chisq")

.i.Analysis of deviance table:in S-plus anova;
 Analysis of Deviance Table
 Response: cbind(y, nsurv - y)
 model  Resid. Df Resid. Dev Df Deviance      Pr(Chi)
      1        83  183.9291
      2        81  154.1944  2 29.73473 3.492900e-07
      3        78  129.3039  3 24.89052 1.627607e-05


 summary(copepod.fit3)


 Coefficients:
                                     Value   Std. Error   t value
                 (Intercept) -0.5071041791 0.1356760682 -3.737610
        I(dose - mean(dose)) -0.0938132223 0.0096979248 -9.673536
    I((dose - mean(dose))^2) -0.0008770721 0.0002366171 -3.706715
                          iA  1.4914504520 0.1937024537  7.699698
                          iC  1.5284498745 0.1904798914  8.024206
     I(dose - mean(dose)):iA  0.0481300815 0.0102795837  4.682104
     I(dose - mean(dose)):iC  0.0587773040 0.0101854209  5.770729
 I((dose - mean(dose))^2):iA  0.0006605230 0.0002521107  2.619972
 I((dose - mean(dose))^2):iC  0.0007538493 0.0002485815  3.032604

 (Dispersion Parameter for Binomial family taken to be 1 )

     Null Deviance: 1156.764 on 86 degrees of freedom

 Residual Deviance: 129.3039 on 78 degrees of freedom