predict from glmer
sj <ssj1364 <at> gmail.com> writes:
predPN <- model.matrix(terms(mod.PN),newDat) %*% fixef(mod.PN)
this seems to work fine and the I print the predictions
print(PredPN)
[1] [2]
2.358722 2.312340
However I am not certain what I am looking at here, my best guess is that
each of these could the odds (p/1-p). But I am not at all sure in my guess,
any insights would be appreciated
These are most likely to be log-odds rather than odds. (gmane wants me to say more or quote less ... you can use plogis(PredPN) to get the probabilities.)