pim.fit function

Fitter function for a probabilistic index model