I have fitted the following model:
fitmixed2 <-lmer(Type~rcs(Pt.Age,3)+rcs(Time.min,3)+Pt.Gender+ASA+School+(1|provider.id))
and thought I could use the predict function with it, but it does not work.
I really want to get the predicted probability plots for some fixed values for age (say, mean(age)), gender=F, ASA=1, and three probability plots for the three levels of school on same axis and with confidence intervals.
lme4you are using? AFAIK, only the developement version on R-forge has apredictfunction. – chl♦ Aug 8 '12 at 1:19