Systems and methods are provided for implementing a parallel Expectation Minimization algorithm for generalized latent variable models. Item response data that is based on responses to items from multiple respondents is accessed. The item response data includes data for multiple response variables. The item response data is analyzed using a generalized latent variable model, and the analysis includes an application of a Parallel-E Parallel-M (PEPM) algorithm. In a parallel Expectation step of the PEPM algorithm, the respondents are subdivided into N groups of respondents, and computations for the N groups are performed in parallel using the N processor cores. In a parallel Maximization step of the PEPM algorithm, the response variables are subdivided into N groups of response variables, and computations for the N groups of response variables are performed in parallel using the N processor cores.