plmm_fit function

PLMM fit: a function that fits a PLMM using the values returned by plmm_prep()