[R] predict for a model with a subset

[Spencer Graves]

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
>  
>