use scorefunctions with regularization

E.g. use BIC when fitting with lm_fit