[R] Extracting mean level of factor for all levels from a glm

[Ruben Roa Ureta]

ComRades:

I have a glm where interest lies in the estimated mean levels of the
response Y for every level of factor Year (regardeless of statistical
significance of the factor), when the effect of other factors have been
accounted for. The model is of the form:

fit<-glm(Y~Year+Mo+Area+Ship+Eff,family=gaussian(link="log"),data=data)

where Year is categorical and has T levels (labeled 1 to T in the data).
Help much appreciated on how to extract these 'standardised' annual
estimates of Y and their standard errors.

Rubén