[R] predict for a model with a subset

Spencer Graves
Fri, 24 Oct 2003 16:56:29 -0000

Have you tried the following:

notmar1 <- glm(yprisx~age+harddrug+sex, subset = marcom == 0,
 family = quasipoisson), data=DF)

predict(notmar1, newdata=DF[DF$marcom==0,])

      I haven't tested this specific code, but I've gotten sensible
results from constructs like this in the past.

hope this helps.  spencer graves

Peter Flom wrote:

>Hello
>
>running R 1.7.1 on Windows 2000
>
>I have a model
>
>notmar1 <- glm(yprisx~age+harddrug+sex, subset = marcom == 0,
> family = quasipoisson)
>
>and summary(notmar1) gives (as it should) 433 df for the null model
>
>but when I run
>predict(notmar1 <- glm(yprisx~age+harddrug+sex, subset = marcom == 0,
>family = quasipoisson))
>
>I get preditions for 528 people (the full data set, not the subset)
>
>How do I get predict to work on just the subset of people for whom the
>model is estimated?
>
>Thanks
>
>Peter
>
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>
>