fitSystematicErrorModel function

Fit a systematic error model