A well-known person fit statistic in the item response theory (IRT) literature is the \(l_{z}\) statistic (Drasgow et al. in Br J Math Stat Psychol 38(1):67-86, 1985). Snijders (Psychometrika 66(3):331-342, 2001) derived \(l_{z}^{*}\) , which is the asymptotically correct version of \(l_{z}\) when the ability parameter is estimated. However, both statistics and other extensions later developed concern either only the unidimensional IRT models or multidimensional models that require a joint estimate of latent traits across all the dimensions. Considering a marginalized maximum likelihood ability estimator, this paper proposes \(l_{zt}\) and \(l_{zt}^{*}\) , which are extensions of \(l_{z}\) and \(l_{z}^{*}\) , respectively, for the Rasch testlet model. The computation of \(l_{zt}^{*}\) relies on several extensions of the Lord-Wingersky algorithm (1984) that are additional contributions of this paper. Simulation results show that \(l_{zt}^{*}\) has close-to-nominal Type I error rates and satisfactory power for detecting aberrant responses. For unidimensional models, \(l_{zt}\) and \(l_{zt}^{*}\) reduce to \(l_{z}\) and \(l_{z}^{*}\) , respectively, and therefore allows for the evaluation of person fit with a wider range of IRT models. A real data application is presented to show the utility of the proposed statistics for a test with an underlying structure that consists of both the traditional unidimensional component and the Rasch testlet component.