fit_domains function

Estimate the Rasch and the Interaction model per domain