Skip to content
Prev 392346 / 398502 Next
Original Raw

Predicted values from glm() when linear predictor is NA.

Rolf Turner 
I have a data frame with a numeric ("TrtTime") and a categorical
("Lifestage") predictor.

Level "L1" of Lifestage occurs only with a single value of TrtTime,
explicitly 12, whence it is not possible to estimate a TrtTime "slope"
when Lifestage is "L1".

Indeed, when I fitted the model

    fit <- glm(cbind(Dead,Alive) ~ TrtTime*Lifestage, family=binomial,
               data=demoDat)

I got:
That is, TrtTime:LifestageL1 is NA, as expected.

I would have thought that fitted or predicted values corresponding to
Lifestage = "L1" would thereby be NA, but this is not the case:
That is, the predicted values on the scale of the linear predictor are
large and positive, rather than being NA.

What this amounts to, it seems to me, is saying that if the linear
predictor in a Binomial glm is NA, then "success" is a certainty.
This strikes me as being a dubious proposition.  My gut feeling is that
misleading results could be produced.

Can anyone explain to me a rationale for this behaviour pattern?
Is there some justification for it that I am not currently seeing?
Any other comments?  (Please omit comments to the effect of "You are as
thick as two short planks!". :-) )

I have attached the example data set in a file "demoDat.txt", should
anyone want to experiment with it.  The file was created using dput() so
you should access it (if you wish to do so) via something like

    demoDat <- dget("demoDat.txt")

Thanks for any enlightenment.

cheers,

Rolf Turner

Thread (8 messages)

Rolf Turner Predicted values from glm() when linear predictor is NA. Jul 27 Andrew Robinson Predicted values from glm() when linear predictor is NA. Jul 27 Ebert,Timothy Aaron Predicted values from glm() when linear predictor is NA. Jul 27 David Winsemius Predicted values from glm() when linear predictor is NA. Jul 27 Rolf Turner Predicted values from glm() when linear predictor is NA. Jul 27 Rolf Turner Predicted values from glm() when linear predictor is NA. Jul 27 Rolf Turner Predicted values from glm() when linear predictor is NA. Jul 27 Ebert,Timothy Aaron Predicted values from glm() when linear predictor is NA. Jul 28