fit_LME_longitudinal function

Fit a landmarking model using a linear mixed effects (LME) model for the longitudinal submodel