[R] prediction from glm

From: Hugues Santin-Janin <hsantin-janin_at_gct.org.uk>
Date: Thu 04 Aug 2005 - 21:32:47 EST


Hello r-help,  

I try to fit birds counts over years using glm. I have done (with Estate and year as factors):

Model1 <- glm(Females~Estate+Year+offset = log(area)), family = quasipoisson(link = log), na.action = "na.exclude")

After I have calculated the prediction using:

Pred1 <- predict(Model1, type = "response", na.action = "na.exclude")  

My question is: How can I obtain predictions for Females in each year that are standardized by averaging over the levels of Estate?

I have been advised to use log(area) as an offset because we are actually interested in density of birds but would appreciate any good comparisons between offset and weight.  

Thanks in advance for any help you can give.  

Hugues SANTIN-JANIN

hugues_sj@yahoo.fr

"Game Conservancy Trust"

**Win some fantastic prizes in our 2005 prize draw. Visit www.gct.org.uk, click on the membership menu and choose 2005 Prize Draw and purchase your draw tickets.**

E-mail Disclaimer
This communication from Hugues Santin-Janin contains information which is CONFIDENTIAL and may also be privileged. It is for the exclusive use of the addressee. If you are not the addressee please note that any distribution, copying or use of this communication or the information in it is prohibited. If you have received this communication in error, please telephone us immediately to arrange for its return.

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Aug 04 21:46:58 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:39:40 EST