Re: [R] negative predicted values in poisson glm

From: Thomas Lumley <tlumley_at_u.washington.edu>
Date: Thu 19 Jan 2006 - 02:28:48 EST

On Wed, 18 Jan 2006, P. Olsson wrote:

> Dear R helpers,
> running the following code of a glm model of the family poisson, gives
> predicted values < 0. Why?

Look at the help page for predict.glm, particularly the "type" argument. Then exponentiate the results.

         -thomas

>
> library(MASS)
> library(stats)
> library(mvtnorm)
> library(pscl)
> data(bioChemists)
> poisson_glm <- glm(art ~ fem + mar + kid5 + phd + ment, data = bioChemists,
> family = poisson)
> predicted.values = predict(poisson_glm)
> range(predicted.values)
>
>
> Thank you in advance for any hints.
> Best regards,
> P. Olsson
>
> [[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
>

Thomas Lumley			Assoc. Professor, Biostatistics
tlumley@u.washington.edu	University of Washington, Seattle

______________________________________________
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 Jan 19 03:31:04 2006

This archive was generated by hypermail 2.1.8 : Thu 19 Jan 2006 - 06:12:21 EST