fit.FirmaModel function

Estimates the model parameters