I used generalized linear model with the commands as follows:
 
For S-Plus:
>attach(data)
>options(contrasts=c("contr.treatment","contr.poly"))
>model1<-glm(logcatchrate~Year+Month+Vessclass+Livebaitday+....+Vessclass:Livebaitday,family=gaussian)
>anova(model1,test="Chisq")
 
For R :
>data<-read.table("data.dat", header=T)
>fix(data)
>attach(data)
>model1<-glm(logcatchrate~Year+Month+Vessclass+Livebaitday+...+Vessclass:Livebaitday,family=gaussian)>anova(model1,test="Chisq")
 
 
 

Prof Brian Ripley <ripley@stats.ox.ac.uk> wrote:
Depends on the `GLM' model and what you mean by that (general linear 
model, generalized linear model?). And by `Anova' (capitalized), for that 
matter.

S-PLUS (sic) does some inconsistent things with dispersion parameters in 
glm() fits, and it is likely that you are using the wrong test. But P 
values that near one should arouse your suspicions.

Please read the R posting guide, and then tell us enough of what you did
so we can help you. 




	[[alternative HTML version deleted]]


