When dealing with longitudinal data, linear mixed-effects models (LMMs) are often used by researchers. However, LMMs are not always the most adequate models, especially if we expect a nonlinear relationship between the outcome and a continuous covariate. To allow for more flexibility, we propose the use of a semiparametric mixed-effects model to evaluate the overall treatment effect on the hemodynamic responses during bone graft healing and build a prediction model for the healing process. The model relies on a closed-form expectation-maximization algorithm, where the unknown nonlinear function is estimated using a Lasso-type procedure. Using this model, we were able to estimate the effect of time for individual mice in each group in a nonparametric fashion and the effect of the treatment while accounting for correlation between observations due to the repeated measurements. The treatment effect was found to be statistically significant, with the autograft group having higher total hemoglobin concentration than the allograft group.