fit_lm_metamodel function

Fit linear metamodel