Disclosed herein is a method of analyzing examinee item response data comprising constructing a diagnosis model for reporting skill profiles of examinees, wherein the diagnosis model comprises at least a variable representing unobserved subpopulations, creating an item design matrix, distributing examinees across the unobserved subpopulations, iteratively estimating values for a plurality of variables within the diagnosis model, and reporting the estimated values to a user.