Structural equation modeling (SEM) with ordinal indicators rely on an assumption of categorized normality. This assumption may be tested for pairs of variables using the likelihood ratio G2 or Pearson’s X2 statistics. For increased computational efficiency, SEM programs usually estimate polychoric correlations in two stages. However, two-stage polychoric estimates are not asymptotically efficient and G2 and X2 need not be asymptotically chi-square when the estimator is not efficient. Recently, Maydeu-Olivares and Joe (2005) have introduced a new statistic, Mn , that is asymptotically chi-square even for estimators that are not efficient. We investigate the behavior of G2, X2, and Mn when testing underlying bivariate normality with polychoric correlations estimated in two stages.