Hi, my problem sounds simple but I've been trying and looking for a solution, but just can not find out how to solve it. I have a dataset and want to do a GAM with binomial link function, random effects (for individuals as there are severeal measurements for everyone) and a time-spline. For this purpose I used the lmer-function from lme4 package. My code looks like this: mod <- lmer(aim ~ bs(time) + var1 + var2 + factor(var3) +... +factor(var12) + (1|ID) ,family=binomial(link = "logit"), data=data) It works, but what I also want to do is plotting the time-spline, and I just can't find out how to do this. I also tried a gamm4 and a amer function, but I think they can't handle all the NA-values in the dataset as I always get error messages when I put all the variables into the model ("Matrices must have same number of columns in rbind2(..1, r)"). I'd be really thankful if someone could help me on this topic... Best regards, Caro -- View this message in context: http://r.789695.n4.nabble.com/lmer-How-to-plot-a-spline-from-function-tp3603083p3603083.html Sent from the R help mailing list archive at Nabble.com.