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