A classical statistical estimation method for fitting IRT models by maximizing likelihood across all possible ability values.