fit_LME_landmark function

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