pirtmodel function

Fits a polytomous IRT model.