Asymptotic Posterior Normality of Multivariate Latent Traits in an IRT Model

Abstract

The asymptotic posterior normality (APN) of the latent variable vector in an item response theory (IRT) model is a crucial argument in IRT modeling approaches. In case of a single latent trait and under general assumptions, Chang and Stout (Psychometrika, 58(1):37–52, 1993) proved the APN for a broad class of latent trait models for binary items. Under the same setup, they also showed the consistency of the latent trait’s maximum likelihood estimator (MLE). Since then, several modeling approaches have been developed that consider multivariate latent traits and assume their APN, a conjecture which has not been proved so far. We fill this theoretical gap by extending the results of Chang and Stout for multivariate latent traits. Further, we discuss the existence and consistency of MLEs, maximum a-posteriori and expected a-posteriori estimators for the latent traits under the same broad class of latent trait models.

1. Introduction

In the context of item response theory (IRT) methodology, statistical inference for the examinee’s ability relies often on the assumption that its posterior distribution given the test response is a normal distribution. As this is usually hard to justify and in contradiction to common models of the examinees abilities distribution in the population, it can assumed to be, for a long test, well approximated by a normal distribution. This assumption of asymptotic posterior normality (APN) is part of the famous Dutch identity conjecture of Holland (1990), who mentioned then that he was not aware of a thorough discussion of APN of latent variables and this would be an interesting area for future research. Shortly after, Chang and Stout (1993) proved the APN for univariate latent traits (LTs), mentioning that APN for multivariate LTs can be proved, but without providing further details or discussing the associated regularity conditions required.

As far as we know, APN of multivariate LTs has not been proved so far for IRT models of a general context, although posterior normality or APN is assumed quite often under various IRT setups (e.g., Anderson & Vermunt, 2000; Anderson & Yu, 2017; Anderson et al., 2007; Hessen, 2012; Li, 2010; Paek, 2016). For example, Pelle et al. (2016) assume posterior multivariate normality for the latent variable vector of a log-linear multidimensional Rasch model for capture–recapture analysis of registration data.

Sometimes the APN-assumption is justified by the APN in a Bayesian framework (pointing to the “Bernstein–von Mises Theorem”) without however proceeding to further details (e.g., the computationally efficient adaptive quadrature methods for high-dimensional item factor analysis (Schilling & Bock, 2005) and for generalized linear mixed models (Rabe-Hesketh et al., 2002) are based on the APN assumption).

In this work, we study the APN for multivariate latent trait models, focusing on models for dichotomous items and targeting at conditions that are tailored to IRT models and thus simpler to verify. APN of LTs, univariate or multivariate, is related to Bayesian asymptotics. In the light of this connection, we deepen in the approach of Ghosal et al. (1995), who discussed asymptotic posterior distributions in a very general setup that includes the regular cases and some non-regular cases as well. They also proved a general result on the asymptotic equivalence of the Bayes and maximum likelihood estimators, a well-known result for the regular cases. In particular, we generalize the approach and results of Chang and Stout (1993), CS hereafter, linking them to the semiproper centering concept of Ghosal et al. (1995), GGS hereafter, and embedding them in their approach. We provide conditions for multivariate APN that correspond one to one to the conditions of CS for univariate LTs, which is the standard approach for IRT models, as alternatives to the conditions imposed in Ghosal et al. (1995). Even for the case of univariate LTs, the proposed approach could be an interesting alternative to that of CS, since it has the advantage of applying also to models with non-monotone item response functions, which is not the case in the CS setup. Furthermore, we discuss conditions under which the existence of the maximum likelihood estimators (MLEs) for latent variable vectors is ensured. The consistency of MLEs under mild conditions, which was indicated as an open issue by Sinharay (2015), follows as a natural consequence of the proof of the APN. Finally, we prove the consistency of maximum a-posteriori (MAP) and expected a-posteriori (EAP) estimators for multivariate LTs.

The paper is organized as follows. Basic notation and the adopted IRT framework is set in Sect. 2, while the CS-theory for a univariate LT is briefly reviewed in Sect. 3. The approach of Ghosal et al. (1995) is discussed and linked to the APN of LTs and the CS-results in Sect. 4. The CS-conditions are generalized for the multivariate case and commented in Sect. 5 while they are verified for characteristic examples in Sect. 6. The main result on APN for multivariate LTs and properties of the MLEs, MAPs and EAPs of LTs are provided in Sect. 7 and supported by a simulation study in Sect. 8. Finally, the results are summarized in Sect. 9. A brief version of the proofs of the results of Sect. 7 is given in “Appendix” while their extended version can be found in the web-appendix. For a preliminary version of these results, see also Chapter 3 in Kornely (2021).

2. Preliminaries

Consider a test consisting of d binary response variables $Y_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i$$\end{document}, $i\in \left[d\right]:=\{1,\dots ,d\}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in [d]:=\{1,\ldots ,d\}$$\end{document}, with $Y_i\in \{0,1\}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i\in \{0,1\}$$\end{document} for the i-th item, where 1 (0) denotes a correct (incorrect) response, and defined over a probability space $(\Omega ,A,P)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\Omega ,\mathcal {A},\mathsf {P})$$\end{document}. Consider further the response vector $Y^{\left(d\right)}={(Y_1,\dots ,Y_d)}^{{}^{⊺}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {Y}}}^{(d)}=(Y_1, \ldots , Y_d)^{^{\intercal }}$$\end{document}, with superscript ${}^{{}^{⊺}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{^{\intercal }}$$\end{document} denoting the transpose of a vector. Thus, the manifest probability for a specific response pattern $y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d)}$$\end{document} is given by $P\left(y^{\left(d\right)}\right)=P(Y^{\left(d\right)}=y^{\left(d\right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P({\varvec{\mathsf {y}}}^{(d)})=\mathsf {P}({\varvec{\mathsf {Y}}}^{(d)}= {\varvec{\mathsf {y}}}^{(d)})$$\end{document}. In an multidimensional IRT (MIRT) modeling framework, manifest probabilities are derived via conditioning on an absolutely continuous latent variable vector $\eta ={({\eta }_1,\dots ,{\eta }_q)}^{{}^{⊺}}\in \Theta \subseteq R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\eta }}=(\eta _1,\ldots ,\eta _q)^{^{\intercal }}\in \Theta \subseteq \mathbb {R}^q$$\end{document}, defined over the same probability space as the binary items with probability density function (pdf) and cumulative distribution function (cdf) $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document} and $H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}$$\end{document}, respectively. In particular, the conditional probability mass function (pmf) of $Y_i|\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i\,|\,\varvec{\eta }$$\end{document} is thus given by

(1) $\begin{array}{c}P(Y_i=y|\eta )=P_i{\left(\eta \right)}^y{(1-P_i\left(\eta \right))}^{1-y},y\in \{0,1\},i\in \left[d\right],\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}(Y_i=y\,|\,\varvec{\eta })=P_i(\varvec{\eta })^y(1-P_i(\varvec{\eta }))^{1-y},y\in \{0,1\}, \ \ i\in [d], \end{aligned}$$\end{document}

with $P_i\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i(\varvec{\eta })$$\end{document} being known as the i-th item response function. In MIRT modeling, specific assumptions are usually imposed on the conditional distribution $P(Y^{\left(d\right)}=y^{\left(d\right)}|\eta )$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {P}({\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)}\,|\,\varvec{\eta })$$\end{document}; namely the assumption of local independence

(2) $\begin{array}{cc}{P}^{\left(d\right)}\left(y^{\left(d\right)}\right|\eta ):=P(Y^{\left(d\right)}=y^{\left(d\right)}|\eta )& =\prod _{i=1}^{d}P(Y_i=y_i|\eta ),y^{\left(d\right)}\in {\{0,1\}}^d,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P^{(d)}({\varvec{\mathsf {y}}}^{(d)}\,|\,\varvec{\eta }):=\mathsf {P}({\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)}\,|\,\varvec{\eta })&=\prod _{i=1}^d\mathsf {P}(Y_i=y_i\,|\,\varvec{\eta }),\quad {\varvec{\mathsf {y}}}^{(d)}\in \{0,1\}^d, \end{aligned}$$\end{document}

and that of monotonicity for the item response functions $P_i\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i(\varvec{\eta })$$\end{document}, i.e., for $i\in \left[d\right]$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in [d]$$\end{document}

(3) $\begin{array}{c}{P}_i\left(\eta \right)\text{is strictly monotonic in every dimension of}\eta \text{being measured.}\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P_i(\varvec{\eta }) \ \text {is strictly monotonic in every dimension of} \ \varvec{\eta }\ \text {being measured.} \end{aligned}$$\end{document}

Note that assumption (3), which is required in the CS-approach, is relaxed in our setup. In the sequel, we denote by ${\left\{Y_i\right\}}_{i\in N}\sim P\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}\sim {\mathcal {P}}(\varvec{\eta })$$\end{document} a sequence of Bernoulli random variables that fulfill (1) and (2) for all $d\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\in \mathbb {N}$$\end{document}.

Due to assumption (2) and using (1), the manifest probabilities are derived through the following integral

(4) $\begin{array}{c}P\left(y^{\left(d\right)}\right)=\int \dots \int (\prod _{i=1}^d{P}_i{\left(\eta \right)}^{y_i}{(1-P_i\left(\eta \right))}^{1-y_i})h\left(\eta \right)d\eta .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P({\varvec{\mathsf {y}}}^{(d)})=\int \ldots \int \Bigg (\prod _{i=1}^d P_i(\varvec{\eta })^{y_i}(1-P_i(\varvec{\eta }))^{1-y_i}\Bigg ) \mathfrak {h}(\varvec{\eta })\mathrm {d}\varvec{\eta }. \end{aligned}$$\end{document}

Remark 1.

For simplicity of notation, we use $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} to denote the random latent variable vector as well as a realization of it. If not clear from the context, we write explicitly $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document} for a realization or $\eta \sim H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\sim \mathcal {H}$$\end{document} for the random vector with values in $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document}. In the sequel, we abbreviate the term latent variable vector to latent vector.

The posterior density of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document}, given an observed response $y^{\left(d\right)}\in {\{0,1\}}^d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d)}\in \{0,1\}^d$$\end{document}, is then given by

(5) $\begin{array}{c}h\left(\eta \right|y^{\left(d\right)}):=h\left(\eta \right|Y^{\left(d\right)}=y^{\left(d\right)})=\frac{P(Y^{\left(d\right)}=y^{\left(d\right)}|\eta )h\left(\eta \right)}{P(Y^{\left(d\right)}=y^{\left(d\right)})}=\frac{exp\left({\ell }^{\left(d\right)}\left(\eta \right|y^{\left(d\right)})\right)h\left(\eta \right)}{P\left(y^{\left(d\right)}\right)},\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} h(\varvec{\eta }\,|\,{\varvec{\mathsf {y}}}^{(d)}):=\mathfrak {h}(\varvec{\eta }\,|\,{\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)})=\frac{\mathsf {P}({\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)}\,|\,\varvec{\eta })\mathfrak {h}(\varvec{\eta })}{\mathsf {P}({\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)})} =\frac{\exp (\ell ^{(d)}(\varvec{\eta }\,|\,{\varvec{\mathsf {y}}}^{(d)}))\mathfrak {h}(\varvec{\eta })}{P({\varvec{\mathsf {y}}}^{(d)})}, \end{aligned}$$\end{document}

where ${\ell }^{\left(d\right)}(·|y^{\left(d\right)}))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell ^{(d)}(\cdot \,|\,{\varvec{\mathsf {y}}}^{(d)}))$$\end{document} is the log-likelihood corresponding to (1), given by

(6) $\begin{array}{c}{\ell }^{\left(d\right)}\left(\eta \right|y^{\left(d\right)}))=\sum _{i=1}^d\left(y_i,{\lambda }_i,\left(\eta \right),-,\psi ,\left({\lambda }_i\left(\eta \right)\right)\right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \ell ^{(d)}(\varvec{\eta }\,|\,{\varvec{\mathsf {y}}}^{(d)}))= \sum _{i=1}^d\left( y_i\lambda _i(\varvec{\eta })-\psi (\lambda _i(\varvec{\eta }))\right) , \end{aligned}$$\end{document}

with ${\lambda }_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _i$$\end{document} denoting the item logit, i.e.,

(7) $\begin{array}{c}{\lambda }_i\left(\eta \right):=log\left(\frac{P_i\left(\eta \right)}{1-P_i\left(\eta \right)}\right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lambda _i(\varvec{\eta }):=\log \left( \frac{P_i(\varvec{\eta })}{1-P_i(\varvec{\eta })}\right) , \end{aligned}$$\end{document}

and the function $\psi (·)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi (\cdot )$$\end{document} being defined as $\psi \left(x\right)=log\left(1,+,exp,\left(,x,\right)\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\psi (x)=\log \left( 1+\exp (x)\right) $$\end{document}, $x\in R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x\in \mathbb {R}$$\end{document}.

Let ${\hat{\eta }}_d=\hat{\eta }\left(y^{\left(d\right)}\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d = {\varvec{\hat{\eta }}}({\varvec{\mathsf {y}}}^{(d)})$$\end{document} denote the MLE of the true value of the latent vector ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}, based on a test realization $y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d)}$$\end{document}. Furthermore, the Fisher information matrix of the test at point $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} is given by

(8) $\begin{array}{c}{I}^{\left(d\right)}\left(\eta \right):=E_{\eta }\left(\nabla ,{\ell }^{\left(d\right)},(\eta \mid Y^{\left(d\right)}),,,{\nabla }^{⊺},{\ell }^{\left(d\right)},(\eta \mid Y^{\left(d\right)})\right)=\sum _{i=1}^d{I}_i\left(\eta \right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {I}^{(d)}(\varvec{\eta }):=\mathsf {E}_{\varvec{\eta }}\left( \nabla \ell ^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)})\,\,\nabla ^{\intercal }\ell ^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)})\right) =\sum _{i=1}^d\mathcal {I}_i(\varvec{\eta }), \end{aligned}$$\end{document}

where $I_i(·)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}_i(\cdot )$$\end{document} is the i-th item information matrix

$\begin{array}{cc}{I}_i\left(\eta \right):=& E_{\eta }\left(\nabla ,log,(,P_i,{\left(\eta \right)}^{Y_i},{(1-P_i\left(\eta \right))}^{1-Y_i},),,,{\nabla }^{⊺},log,(,P_i,{\left(\eta \right)}^{Y_i},{(1-P_i\left(\eta \right))}^{1-Y_i},)\right)\\ =& P_i\left(\eta \right)(1-P_i\left(\eta \right))\nabla {\lambda }_i\left(\eta \right)\nabla {\lambda }_i{\left(\eta \right)}^{⊺},\eta \in \Theta .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathcal {I}_i(\varvec{\eta }):=&\,\mathsf {E}_{\varvec{\eta }}\left( \nabla \log \Big (P_i(\varvec{\eta })^{Y_i}(1-P_i(\varvec{\eta }))^{1-Y_i}\Big )\,\,\nabla ^{\intercal }\log \Big (P_i(\varvec{\eta })^{Y_i}(1-P_i(\varvec{\eta }))^{1-Y_i}\Big ) \right) \\ =&P_i(\varvec{\eta })(1-P_i(\varvec{\eta }))\nabla \lambda _i(\varvec{\eta })\nabla \lambda _i(\varvec{\eta })^{\intercal },\quad \varvec{\eta }\in \Theta . \end{aligned}$$\end{document}

This work studies the APN of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} for $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document}, based on a sequence of random variables ${\left\{Y_i\right\}}_{i\in N}\sim P\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}\sim {\mathcal {P}}(\varvec{\eta })$$\end{document}, as defined above. Particularly, we shall prove that, under certain conditions, (8) is invertible at ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d$$\end{document} and $\eta |Y^{\left(d\right)}=y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\,|\,{\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}}^{(d)}$$\end{document} is approximately normal distributed, $N({\hat{\eta }}_d,{\left[I^{\left(d\right)}\left(\widehat{{\eta }_d}\right)\right]}^{-1})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {N}}({\varvec{\hat{\eta }}}_d, [\mathcal {I}^{(d)}({\hat{\varvec{\eta }_d}})]^{-1})$$\end{document}, for a realization $y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d)}$$\end{document} of $Y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {Y}}}^{(d)}$$\end{document}. This enables the approximation of probabilities of the type

(9) $\begin{array}{c}P\left(\left(I^{\left(d\right)},{\left({\hat{\eta }}_d\right)}^{1/2},\left(\eta ,-,{\hat{\eta }}_d\right),\in ,B,\right),,Y^{\left(d\right)}\right),B\in B^q,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}\left( \left. \mathcal {I}^{(d)}({\varvec{\hat{\eta }}}_d)^{1/2}\left( \varvec{\eta }-{\varvec{\hat{\eta }}}_d\right) \in B \,\right| \, {\varvec{\mathsf {Y}}}^{(d)}\right) , \ \ B\in \mathcal {B}^q, \end{aligned}$$\end{document}

where $B^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {B}^q$$\end{document} denotes the Borel- $\sigma$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}-algebra of $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document}. Practically speaking, a set B can be any countable union or intersection of q-dimensional real cubes.

Next, we define some functions that are useful for the sequel derivations. For all $d\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\in \mathbb {N}$$\end{document}, set $Z^{\left(d\right)}:\Theta \times \Theta \to R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^{(d)}:\Theta \times \Theta \rightarrow \mathbb {R}$$\end{document} with

(10) $\begin{array}{c}{Z}^{\left(d\right)}(\eta ,{\eta }^{\prime }):=\prod _{i=1}^d{Z}_i(\eta ,{\eta }^{\prime }),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Z^{(d)}(\varvec{\eta },\varvec{\eta }'):= \prod _{i=1}^d Z_i(\varvec{\eta },\varvec{\eta }'), \end{aligned}$$\end{document}

where $Z_i:\Theta \times \Theta \to R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z_i:\Theta \times \Theta \rightarrow \mathbb {R}$$\end{document}, $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, are defined as

(11) $\begin{array}{c}{Z}_i(\eta ,{\eta }^{\prime }):=\frac{P_i{\left(\eta \right)}^{Y_i}{(1-P_i\left(\eta \right))}^{1-Y_i}}{P_i{\left({\eta }^{\prime }\right)}^{Y_i}{(1-P_i\left({\eta }^{\prime }\right))}^{1-Y_i}}.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Z_i(\varvec{\eta },\varvec{\eta }'):=\frac{P_i(\varvec{\eta })^{Y_i}(1-P_i(\varvec{\eta }))^{1-Y_i}}{P_i(\varvec{\eta }')^{Y_i}(1-P_i(\varvec{\eta }'))^{1-Y_i}} \ . \end{aligned}$$\end{document}

Note that for given d and $\eta ,{\eta }^{\prime }\in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta },\,\varvec{\eta }'\in \Theta $$\end{document}, (10) is the likelihood ratio of the likelihoods for $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} and ${\eta }^{\prime }$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }'$$\end{document}. Furthermore,

(12) $\begin{array}{c}log\left(Z^{\left(d\right)},(\eta ,{\eta }^{\prime })\right)=\sum _{i=1}^{d}log{Z}_i(\eta ,{\eta }^{\prime })={\ell }^{\left(d\right)}(\eta \mid Y^{\left(d\right)})-{\ell }^{\left(d\right)}({\eta }^{\prime }\mid Y^{\left(d\right)}),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \log \left( Z^{(d)}(\varvec{\eta },\varvec{\eta }')\right) =\sum _{i=1}^d\log Z_i(\varvec{\eta },\varvec{\eta }')=\ell ^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)})-\ell ^{(d)}(\varvec{\eta }'\mid {\varvec{\mathsf {Y}}}^{(d)}), \end{aligned}$$\end{document}

while $-E_{{\eta }_0}(log{Z}_i(\eta ,{\eta }_0))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\mathsf {E}_{\varvec{\eta }_0}(\log Z_i(\varvec{\eta },\varvec{\eta }_0))$$\end{document} is the Kullback–Leibler divergence between the conditional distributions of $Y_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i$$\end{document} given $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} and ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}, respectively. A basic approach for deriving APN results relies on a quadratic approximation of (12).

3. Review of APN for Univariate Latent Traits

In case of a single latent variable ( $q=1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=1$$\end{document}, $\eta =\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }=\eta $$\end{document}), Chang and Stout (1993) proved the APN of the univariate latent trait, adopting the approach of Walker (1969) for binary $Y_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i$$\end{document}, $i\in \left[d\right]$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in [d]$$\end{document}, that are independent but not identically distributed (inid). We briefly review their results, so that we can extend in the sequel their approach to the multivariate case ( $q>1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q>1$$\end{document}).

Additional to the general assumptions (2) and (3), they also introduced the following regularity conditions.

1. (CS1)

   1. [i] Let $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \in \Theta $$\end{document}, where $\Theta \subseteq (-\infty ,\infty )$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta \subseteq (-\infty ,\infty )$$\end{document} is a bounded or unbounded interval.

   2. [ii] Let the prior density $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document} be continuous and positive at the true value ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document}.

2. (CS2) $P_i\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i(\eta )$$\end{document} is twice continuously differentiable with the first two derivatives being uniformly bounded in absolute value with respect to both $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} and i in some closed interval ${\Theta }_0\subset \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta _0\subset \Theta $$\end{document} around ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document}.

3. (CS3) For every fixed $\eta \ne {\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \ne \eta _0$$\end{document}, $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \in \Theta $$\end{document}, there is a $c\left(\eta \right)>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c(\eta )>0$$\end{document} such that

   (13) $\begin{array}{cc}\underset{d⟶\infty }{lim\; sup}\frac{1}{d}\sum _{i=1}^d{E}_{{\eta }_0}log{Z}_i(\eta ,{\eta }_0)\le & -c\left(\eta \right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \limsup _{d\longrightarrow \infty }\frac{1}{d}\sum _{i=1}^d\mathsf {E}_{\eta _0}\log Z_i(\eta ,\eta _0)\le&-c(\eta ), \end{aligned}$$\end{document}
   and ${sup}_{i\in N}\left|{\lambda }_i\left(\eta \right)\right|<\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sup _{i\in \mathbb {N}}|\lambda _i(\eta )|<\infty $$\end{document}.

4. (CS4) If restricted to ${\Theta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta _0$$\end{document}, the following sets of functions are uniformly bounded:

   $\begin{array}{c}\left(\left(\frac{d{I}_i}{d\eta }\right),\mid ,i,\in ,N\right),\left(\left(\frac{d^2{\lambda }_i}{d{\eta }^2}\right),\mid ,i,\in ,N\right),\left(\left(\frac{d^3{\lambda }_i}{d{\eta }^3}\right),\mid ,i,\in ,N\right).\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \left\{ \left| \frac{\mathrm {d}\mathcal {I}_i}{\mathrm {d}\eta }\right| \mid i\in \mathbb {N}\right\} ,\quad \left\{ \left| \frac{\mathrm {d}^2 \lambda _i}{\mathrm {d}\eta ^2}\right| \mid i\in \mathbb {N}\right\} ,\quad \left\{ \left| \frac{\mathrm {d}^3 \lambda _i}{\mathrm {d}\eta ^3}\right| \mid i\in \mathbb {N}\right\} . \end{aligned}$$\end{document}

5. (CS5) Asymptotically, the average information at ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} is bounded away from 0, i.e.,

   $\begin{array}{c}\underset{d⟶\infty }{lim\; inf}\frac{I^{\left(d\right)}\left({\eta }_0\right)}{d}>0.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \liminf _{d\longrightarrow \infty }\frac{\mathcal {I}^{(d)}(\eta _0)}{d}>0. \end{aligned}$$\end{document}

Remark 2.

With respect to the prior of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}, additional to (CS1[ii]), Chang and Stout (1993) implicitly assumed its properness, which was stated explicitly in the earlier associated technical report (Chang & Stout, 1991, p. 15).

Remark 3.

Reasonable models for applications do not depend on a specific compact interval in $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document} since usually ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} is unknown. For this, also the conditions depending on ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} should be satisfied for almost all ${\eta }_0\in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0\in \Theta $$\end{document} and for almost each ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} there should be some (arbitrary small) interval ${\Theta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta _0$$\end{document}. In the usual models these conditions are satisfied.

Chang and Stout (1993) argued convincingly that conditions (CS1)–(CS5) are realistic and non-restrictive in practice for commonly used IRT models of well-designed tests. They particularly commented condition (CS3) and (13), which plays an important role in the proof of their main theorem. (CS3) is required when the item responses ${\left\{Y_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}$$\end{document} are independent but not identically distributed. If they are iid, (CS3) is automatically satisfied, which however is not necessarily the case in IRT models. Their main results are expressed in the three theorems given below.

Theorem 1.

(Chang & Stout, 1993, Theorem 1) Suppose that conditions (CS1) through (CS5) hold for a fixed ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document}. Let ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\eta }}_d$$\end{document} be the MLE of ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} and ${\hat{\sigma }}_d={\left(I^{\left(d\right)}\left({\hat{\eta }}_d\right)\right)}^{-1/2}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\sigma }}_d=(\mathcal {I}^{(d)}({\hat{\eta }}_d))^{-1/2}$$\end{document}. Then, for $-\infty \le a<b\le \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\infty \le a<b\le \infty $$\end{document}, the posterior probability of ${\hat{\eta }}_d+a{\hat{\sigma }}_d<\eta <{\hat{\eta }}_d+b{\hat{\sigma }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\eta }}_d+a{\hat{\sigma }}_d<\eta <{\hat{\eta }}_d+b{\hat{\sigma }}_d$$\end{document} approaches the probability of $Z\in (a,b)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z\in (a,b)$$\end{document} in $P_{{\eta }_0}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {P}_{\eta _0}$$\end{document} for $Z\sim N(0,1)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z\sim \mathcal {N}(0,1)$$\end{document}, that means

$\begin{array}{c}{A}_d\equiv {\int }_{{\hat{\eta }}_d+a{\hat{\sigma }}_d}^{{\hat{\eta }}_d+b{\hat{\sigma }}_d}h(\eta \mid Y^{\left(d\right)})d\eta \stackrel{P_{{\eta }_0}}{⟶}\frac{1}{\sqrt{2\pi }}{\int }_a^{b}exp\left(-,\frac{1}{2},{\eta }^2\right)d\eta \equiv A,d\to \infty .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} A_d \equiv \int _{{\hat{\eta }}_d+a{\hat{\sigma }}_d}^{{\hat{\eta }}_d+b{\hat{\sigma }}_d}h(\eta \mid {\varvec{\mathsf {Y}}}^{(d)})\,\mathrm {d}\eta \overset{\mathsf {P}_{\eta _0}}{\longrightarrow } \frac{1}{\sqrt{2\pi }}\int _a^b\exp \left( -\frac{1}{2}\eta ^2\right) \,\mathrm {d}\eta \equiv A,\quad \quad d\rightarrow \infty . \end{aligned}$$\end{document}

Theorem 2.

(Chang & Stout, 1993, Theorem 2) Suppose that conditions (CS1) through (CS5) hold for fixed ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} and let ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\eta }}_d$$\end{document} and ${\hat{\sigma }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\sigma }}_d$$\end{document} be defined as in Theorem 1. Then, for $-\infty \le a<b\le \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\infty \le a<b\le \infty $$\end{document}, the posterior probability $A_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_d$$\end{document} approaches A $P_{{\eta }_0}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {P}_{\eta _0}$$\end{document}-almost surely, as $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document}.

Theorem 3.

(Chang & Stout, 1993, Theorem 3) Assume $\Theta ={\Theta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta =\Theta _0$$\end{document}, a finite interval. Suppose that conditions (CS1) through (CS5) hold for all ${\eta }_0\in {\Theta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0\in \Theta _0$$\end{document} and let ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\eta }}_d$$\end{document} and ${\hat{\sigma }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\sigma }}_d$$\end{document} be defined as in Theorem 1. Then, for $-\infty \le a<b\le \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-\infty \le a<b\le \infty $$\end{document}, the posterior probability $A_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_d$$\end{document} approaches A in manifest probability $P$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {P}$$\end{document}, as $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document}.

The result of Theorem 3 does not depend on the true value ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} and is thus of special practical interest for estimation and prediction purposes. As Chang and Stout (1993) comment, Theorems 1 and 3 treat sampling from a fixed ability sub-population and from the whole population, respectively. An important by-product of the proof of the APN of latent variables distributions was the establishment of the weak and strong consistency of the MLE of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} under milder conditions than Lord (1983).

Due to the theorems above, the following approximation for a large d and any observed response pattern $y\in {\{0,1\}}^d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}\in \{0,1\}^d$$\end{document}, i.e., the construction of asymptotic credible intervals, is justified

(14) $\begin{array}{c}P(a\le \eta \le b\mid Y^{\left(d\right)}=y)\approx {\Phi }_1(b;{\hat{\eta }}_d,{\hat{\sigma }}^2)-{\Phi }_1(a;{\hat{\eta }}_d,{\hat{\sigma }}^2),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}(a\le \eta \le b\mid {\varvec{\mathsf {Y}}}^{(d)}={\varvec{\mathsf {y}}})\approx \Phi _1(b\,;\,{\hat{\eta }}_d,{\hat{\sigma }}^2)-\Phi _1(a\,;\,{\hat{\eta }}_d,{\hat{\sigma }}^2), \end{aligned}$$\end{document}

for $-\infty \le a\le b\le \infty \in R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ -\infty \le a\le b \le \infty \in \mathbb {R}$$\end{document}, where ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\eta }}_d$$\end{document} is the MLE of ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} based on the sample $y$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}$$\end{document}, ${\hat{\sigma }}^2={\left(I^{\left(d\right)}\left({\hat{\eta }}_d\right)\right)}^{-1}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{\sigma }}^2=(\mathcal {I}^{(d)}({\hat{\eta }}_d))^{-1}$$\end{document} and ${\Phi }_1(·;{\hat{\eta }}_d,{\hat{\sigma }}^2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Phi _1(\cdot \,;\,{\hat{\eta }}_d,{\hat{\sigma }}^2)$$\end{document} denotes the cdf of $N({\hat{\eta }}_d,{\hat{\sigma }}^2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathcal {N}}({\hat{\eta }}_d,{\hat{\sigma }}^2)$$\end{document}. Approximation (14) is of special practical importance in the context of long tests where the exact computation of posterior probabilities for latent variables is commonly intractable. Furthermore, (14) allows the approximation of the posterior if the exact distribution $H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}$$\end{document} of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} is unavailable or uncertain.

Finally, Chang and Stout (1993) noted that their theory, under suitable regularity conditions, can be extended to prove the APN for latent vectors of general multidimensional IRT models, without however commenting further the proving procedure or the regularity conditions required. Next, we discuss the asymptotic posterior distribution of multivariate latent traits in the context of MIRT.

4. APN for Multivariate Latent Traits

The theory of APN of the latent variables is naturally linked to Bayesian procedures and results on the convergence of posterior distributions. In particular, interesting and inspiring is the fundamental contribution by GGS (Ghosal et al., 1995), who consider asymptotic multivariate posterior distributions (not necessarily normal) in a very general and flexible framework discussing different types of convergence, relying on earlier works by Ghosh et al. (1994) and Ibragimov and Has’minskii (1981), denoted as IH hereafter. In particular, they studied posterior convergence of suitably centered and normalized posteriors. Their results provide a very general framework, which can be adopted for the APN in the IRT setup. Next, we adjust the GGS approach for MIRT models and discuss their conditions, embedding the CS approach in the GGS framework.

Following Ghosal et al. (1995, Definition 2), we distinguish two types of APN and link them to the statistic used for the centering of the posterior distribution of the latent vector.

Definition 1.

Let $Z\sim N_q(0,I_q)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {Z}}}\sim \mathcal {N}_q({\varvec{\mathsf {0}}}, {\varvec{\mathsf {I}}}_q)$$\end{document} be a q-variate standard normal distributed random vector. A $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document}-valued statistic ${\tilde{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\tilde{\varvec{\eta }}}_d$$\end{document} is called a proper centering (with limiting normal distribution) if

(15) $\begin{array}{c}\underset{A\in B^q}{sup}\left(P,\left(I^{\left(d\right)},{\left({\tilde{\eta }}_d\right)}^{1/2},(\eta -{\tilde{\eta }}_d),\in ,A,\left(,Y^{\left(d\right)}\right)\right),-,P,(Z\in A)\right)\stackrel{P_{{\eta }_0}}{⟶}0,\text{as}d\to \infty .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sup _{A\in \mathcal {B}^q}\left| \mathsf {P}\left( \mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2}(\varvec{\eta }-{\tilde{\varvec{\eta }}}_d)\in A\left| \ {\varvec{\mathsf {Y}}}^{(d)}\right. \right) -\mathsf {P}({\varvec{\mathsf {Z}}}\in A)\right| \overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }0, \ \text {as} \ d\rightarrow \infty . \end{aligned}$$\end{document}

A statistic ${\tilde{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\tilde{\varvec{\eta }}}_d$$\end{document} is called semiproper centering (with limiting normal distribution) if, for all $A\in B^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A\in \mathcal {B}^q$$\end{document},

(16) $\begin{array}{c}P\left(I^{\left(d\right)},{\left({\tilde{\eta }}_d\right)}^{1/2},(\eta -{\tilde{\eta }}_d),\in ,A,\left(,Y^{\left(d\right)}\right)\right)\stackrel{P_{{\eta }_0}}{⟶}P(Z\in A),\text{as}d\to \infty .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}\left( \mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2}(\varvec{\eta }-{\tilde{\varvec{\eta }}}_d)\in A\left| \ {\varvec{\mathsf {Y}}}^{(d)}\right. \right) \overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }\mathsf {P}({\varvec{\mathsf {Z}}}\in A), \ \text {as} \ d\rightarrow \infty . \end{aligned}$$\end{document}

A statistic ${\tilde{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\tilde{\varvec{\eta }}}_d$$\end{document} is called compatible (with the posterior), if

$\begin{array}{c}\left(I^{\left(d\right)},{\left({\tilde{\eta }}_d\right)}^{1/2},({\tilde{\eta }}_d-{\eta }_0),,,h^{\ast },\left(,·,,\left(,Y^{\left(d\right)}\right)\right)\right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \left( \mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2}({\tilde{\varvec{\eta }}}_d-\varvec{\eta }_0),h^{*}\left( \,\cdot \,\left| \,{\varvec{\mathsf {Y}}}^{(d)}\right) \right. \right) , \end{aligned}$$\end{document}

as a random element in $R^q\times L^1\left(R^q\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q\times L^1(\mathbb {R}^q)$$\end{document}, converges in distribution for $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document}, where $h^{\ast }(·\mid Y^{\left(d\right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{*}(\,\cdot \mid {\varvec{\mathsf {Y}}}^{(d)})$$\end{document} denotes the density of the posterior distribution of $I^{\left(d\right)}{\left({\eta }_0\right)}^{1/2}(\eta -{\eta }_0)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}^{(d)}(\varvec{\eta }_0)^{1/2}(\varvec{\eta }-\varvec{\eta }_0)$$\end{document} and $L^1\left(R^q\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^1(\mathbb {R}^q)$$\end{document} stands for the space of all q-variate Lebesgue-integrable real functions on $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document}.

Proper and semiproper centering correspond to uniform and pointwise convergence of the posterior of the standardized latent vector $I^{\left(d\right)}{\left({\tilde{\eta }}_d\right)}^{1/2}(\eta -{\tilde{\eta }}_d)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2}(\varvec{\eta }-{\tilde{\varvec{\eta }}}_d)$$\end{document}, respectively. Hence, proper centering is a stronger property than semiproper centering, and is consequently expected to require stronger assumptions.

Under this view, one can easily recognize that Theorem 1 of Chang and Stout (1993) is the semiproper centering of the MLE, since it can be formulated as

$\begin{array}{c}P\left(I^{\left(d\right)},{\left({\hat{\eta }}_d\right)}^{1/2},\left(\eta ,-,{\hat{\eta }}_d\right),\in ,[a,b],\left(,Y^{\left(d\right)}\right)\right)\stackrel{P_{{\eta }_0}}{⟶}P(Z\in [a,b]),d⟶\infty ,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}\left( \mathcal {I}^{(d)}({\hat{\eta }}_d)^{1/2}\left( \eta -{\hat{\eta }}_d\right) \in [a,b]\left| \ {\varvec{\mathsf {Y}}}^{(d)}\right. \right) \overset{\mathsf {P}_{\eta _0}}{\longrightarrow }\mathsf {P}\big (Z\in [a,b]\,\big ),\quad \quad d\longrightarrow \infty , \end{aligned}$$\end{document}

for $Z\sim N(0,1)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z\sim \mathcal {N}(0,1)$$\end{document}, $A=[a,b]\subset R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A=[a,b]\subset \mathbb {R}$$\end{document} and ${\tilde{\eta }}_d={\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{{\tilde{\eta }}}}_d={\hat{\eta }}_d$$\end{document} being the MLE of ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta _0$$\end{document} based on $Y^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {Y}}}^{(d)}$$\end{document}. Thus, for the extension of the CS-theory for multivariate LTs, we focus on semiproper centering.

The asymptotic results of GGS adjusted in our setup, primarily focus on the convergence of the posterior distribution of the standardized latent vector

(17) $\begin{array}{c}{\eta }^{\ast }=I^{\left(d\right)}{\left({\tilde{\eta }}_d\right)}^{1/2}(\eta -{\eta }_0),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \varvec{\eta }^* = \mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2} (\varvec{\eta }-\varvec{\eta }_0), \end{aligned}$$\end{document}

with ${\eta }^{\ast }\in {\Theta }_d:=I^{\left(d\right)}{\left({\tilde{\eta }}_d\right)}^{1/2}(\Theta -{\eta }_0)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }^*\in \Theta _d:= \mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2} (\Theta -\varvec{\eta }_0)$$\end{document}. We need the likelihood ratio (10) expressed in terms of ${\eta }^{\ast }$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }^*$$\end{document}, which is denoted by

(18) $\begin{array}{c}{Z}^{\ast \left(d\right)}\left({\eta }^{\ast }\right):=Z^{\left(d\right)}({\eta }_0+I^{\left(d\right)}{\left({\tilde{\eta }}_d\right)}^{-1/2}{\eta }^{\ast },{\eta }_0).\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Z^{*(d)}(\varvec{\eta }^*):=Z^{(d)}(\varvec{\eta }_0+\mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{-1/2}\varvec{\eta }^*,\varvec{\eta }_0). \end{aligned}$$\end{document}

In our setup, for binary response variables $Y_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i$$\end{document}, $i\in \left[d\right]$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in [d]$$\end{document}, and log-likelihoods given by (6), the likelihood ratio $Z^{\ast \left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^{*(d)}$$\end{document} takes the form

$\begin{array}{cc}{Z}^{\ast \left(d\right)}\left({\eta }^{\ast }\right)=& exp(\sum _{i=1}^d(Y_i\left({\lambda }_i,\left({\eta }_0,+,I^{\left(d\right)},{\left({\tilde{\eta }}_d\right)}^{-1/2},{\eta }^{\ast }\right),-,{\lambda }_i,(,{\eta }_0,)\right)\\ & -\left(\psi ,\left({\lambda }_i,\left({\eta }_0,+,I^{\left(d\right)},{\left({\tilde{\eta }}_d\right)}^{-1/2},{\eta }^{\ast }\right)\right),-,\psi ,\left({\lambda }_i,\left({\eta }_0\right)\right)\right)),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} Z^{*(d)}(\varvec{\eta }^*)=&\exp \Bigg (\sum _{i=1}^d\Bigg (Y_i\left( \lambda _i\left( \varvec{\eta }_0+\mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{-1/2}\varvec{\eta }^*\right) -\lambda _i\Big (\varvec{\eta }_0\Big )\right) \\&-\left( \psi \left( \lambda _i\left( \varvec{\eta }_0+\mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{-1/2}\varvec{\eta }^*\right) \right) -\psi \left( \lambda _i\left( \varvec{\eta }_0\right) \right) \right) \Bigg ), \end{aligned}$$\end{document}

with the item logits ${\lambda }_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _i$$\end{document} provided in (7).

The primary conditions of GGS for APN are given as follows:

1. (GGS1) For some $M>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M>0$$\end{document}, $m_1\ge 0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_1\ge 0$$\end{document} and $\alpha >0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha >0$$\end{document} holds

   $\begin{array}{c}{E}_{{\eta }_0}\left({\left(Z^{\ast \left(d\right)},{\left({\eta }_1^{\ast }\right)}^{1/2},-,Z^{\ast \left(d\right)},{\left({\eta }_2^{\ast }\right)}^{1/2}\right)}^2\right)\le M(1+R^{m_1}){\Vert {\eta }_1^{\ast }-{\eta }_2^{\ast }\Vert }^{\alpha },\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {E}_{\varvec{\eta }_0}\left( \left| Z^{*(d)}(\varvec{\eta }^*_1)^{1/2}-Z^{*(d)}(\varvec{\eta }^*_2)^{1/2}\right| ^2\right) \le M(1+R^{m_1})\Vert \varvec{\eta }_1^*-\varvec{\eta }_2^*\Vert ^\alpha , \end{aligned}$$\end{document}
   for all ${\eta }_j^{\ast }\in {\Theta }_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_j^*\in \Theta _d$$\end{document}, satisfying $\Vert {\eta }_j^{\ast }\Vert \le R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Vert \varvec{\eta }_j^*\Vert \le R$$\end{document}, $j=1,2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j=1,2$$\end{document}, where $\Vert ·\Vert$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Vert \cdot \Vert $$\end{document} is the Euclidean norm.

2. (GGS2) For all ${\eta }^{\ast }\in {\Theta }_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }^*\in \Theta _d$$\end{document} holds

   $\begin{array}{c}{E}_{{\eta }_0}\left(Z^{\ast \left(d\right)},{\left({\eta }^{\ast }\right)}^{1/2}\right)\le exp(-g_d(\Vert {\eta }^{\ast }\Vert )),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {E}_{\varvec{\eta }_0}\left( Z^{*(d)}(\varvec{\eta }^*)^{1/2}\right) \le \exp \Big (-g_d\big (\Vert \varvec{\eta }^*\Vert \big )\Big ), \end{aligned}$$\end{document}
   where ${\left\{g_d\right\}}_{d\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{g_d\}_{d\in \mathbb {N}}$$\end{document} is a sequence of real-valued functions on $[0,\infty )$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[0,\infty )$$\end{document} satisfying the following: (a) for a fixed $d\ge 1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\ge 1$$\end{document}, ${lim}_{x\to \infty }{g}_d\left(x\right)=\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lim _{x\rightarrow \infty }g_d(x)=\infty $$\end{document}; (b) for any $N>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N>0$$\end{document},
   $\begin{array}{c}\underset{x\to \infty }{lim}\underset{d\to \infty }{lim}{x}^{N}exp(-g_d\left(x\right))=\underset{d\to \infty }{lim}\underset{x\to \infty }{lim}{x}^{N}exp(-g_d\left(x\right))=0.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lim _{x\rightarrow \infty }\lim _{d\rightarrow \infty } x^N\exp (-g_d(x))=\lim _{d\rightarrow \infty }\lim _{x\rightarrow \infty } x^N\exp (-g_d(x))=0. \end{aligned}$$\end{document}

3. (GGS3) For all $n\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} and ${\eta }_1^{\ast },\dots ,{\eta }_n^{\ast }\in R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }^*_1,\ldots ,\varvec{\eta }^*_n\in \mathbb {R}^q$$\end{document}, the vector of the likelihood-ratios, defined in (18), satisfies

   $\begin{array}{c}(Z^{\ast \left(d\right)}\left({\eta }_1^{\ast }\right),\dots ,Z^{\ast \left(d\right)}\left({\eta }_n^{\ast }\right))\stackrel{D}{⟶}(Z\left({\eta }_1^{\ast }\right),\dots ,Z\left({\eta }_n^{\ast }\right)),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} (Z^{*(d)}(\varvec{\eta }^*_1),\ldots ,Z^{*(d)}(\varvec{\eta }^*_n))\overset{\mathcal {D}}{\longrightarrow }(Z(\varvec{\eta }^*_1),\ldots ,Z(\varvec{\eta }^*_n)), \end{aligned}$$\end{document}
   for $d⟶\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\longrightarrow \infty $$\end{document}, where $\stackrel{D}{⟶}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overset{\mathcal {D}}{\longrightarrow }$$\end{document} denotes convergence in distribution and $Z\left({\eta }^{\ast }\right)=exp({\xi }^T{\eta }^{\ast }-\frac{1}{2}\Vert {\eta }^{\ast }{\Vert }^2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z(\varvec{\eta }^*)=\exp (\varvec{\xi }^T\varvec{\eta }^*-\frac{1}{2}\Vert \varvec{\eta }^*\Vert ^2)$$\end{document}, ${\eta }^{\ast }\in R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }^*\in \mathbb {R}^q$$\end{document}, where $\xi \sim N_q(0,I_q)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\xi }\sim \mathcal {N}_q(\pmb 0,\mathrm {I}_q)$$\end{document}.

Under these conditions, Ghosal et al. (1995) provided the following general result. Notice that they discussed a far more general framework, allowing further distributions for the response variable and considering cases for which the posterior may converge to another distribution than a normal. We refer to GGS for further details regarding these cases.

Theorem 4.

(Ghosal et al., 1995, Theorem 1) Assume that conditions (GGS1) through (GGS3) hold. If either a proper centering or a semiproper compatible centering sequence ${\left\{{\tilde{\eta }}_d\right\}}_{d\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{{\tilde{\varvec{\eta }}}_d\}_{d\in \mathbb {N}}$$\end{document} exists, then it exists a random vector $W$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb W$$\end{document}, such that (a) $I^{\left(d\right)}{\left({\tilde{\eta }}_d\right)}^{1/2}({\tilde{\eta }}_d-{\eta }_0)\stackrel{D}{⟶}W$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}^{(d)}({\tilde{\varvec{\eta }}}_d)^{1/2}({\tilde{\varvec{\eta }}}_d-\varvec{\eta }_0)\overset{\mathcal {D}}{\longrightarrow }\pmb W$$\end{document} for $d⟶\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\longrightarrow \infty $$\end{document} and (b) for almost all $x\in R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb x\in \mathbb {R}^q$$\end{document}, $\frac{Z(x-W)}{{\int }_{R^q}Z(x^{\ast }-W)d{x}^{\ast }}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{Z(\pmb x-\pmb W)}{\int _{\mathbb {R}^q}Z(\pmb x^*-\pmb W)\,\mathrm {d}\pmb x^*}$$\end{document} is nonrandom, where Z is as defined in condition (GGS3). Conversely, if (b) holds for a random vector $W$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pmb W$$\end{document}, then any Bayes estimator (with respect to a prior and loss considered by Ghosal et al. (1995)) is a compatible proper centering.

Applying Theorem 4 for an appropriate Bayes estimator for ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}, the APN of an MIRT model under conditions (GGS1) to (GGS3) is derived. The extension of Theorem 4 for an MLE, i.e., for ${\tilde{\eta }}_d={\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\tilde{\varvec{\eta }}}_d = {\varvec{\hat{\eta }}}_d$$\end{document}, is based on its asymptotic equivalence to an arbitrary Bayes estimator, which has been proved by Ghosal et al. (1995, cf. Corollary 1) under (GGS2)–(GGS3) and the following strengthened form of (GGS1):

1. (GGS1') For some $M>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M>0$$\end{document}, $m_1\ge 0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_1\ge 0$$\end{document} and $m\ge \alpha >q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m\ge \alpha >q$$\end{document} holds

   $\begin{array}{c}{E}_{{\eta }_0}\left({\left(Z^{\ast \left(d\right)},{\left({\eta }_1^{\ast }\right)}^{1/m},-,Z^{\ast \left(d\right)},{\left({\eta }_2^{\ast }\right)}^{1/m}\right)}^m\right)\le M(1+R^{m_1}){\Vert {\eta }_1^{\ast }-{\eta }_2^{\ast }\Vert }^{\alpha },\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {E}_{\varvec{\eta }_0}\left( \left| Z^{*(d)}(\varvec{\eta }^*_1)^{1/m}-Z^{*(d)}(\varvec{\eta }^*_2)^{1/m}\right| ^m\right) \le M(1+R^{m_1})\Vert \varvec{\eta }_1^*-\varvec{\eta }_2^*\Vert ^\alpha , \end{aligned}$$\end{document}
   for all ${\eta }_j^{\ast }\in {\Theta }_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_j^*\in \Theta _d$$\end{document}, satisfying $\Vert {\eta }_j^{\ast }\Vert \le R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Vert \varvec{\eta }_j^*\Vert \le R$$\end{document}, $j=1,2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j=1,2$$\end{document}.

Remark 4.

Alternatively to the GGS conditions discussed above, one could consider the conditions of Ibragimov & Has’minskii (1981, Section III.4) for general regular models for independent non-necessarily identical distributed (inid) random variables. They proved that these conditions are sufficient for the set of conditions N1–N4 of IH, Section III.1, where N1 is the uniform asymptotic normality and corresponds to (GGS3), while N3 and N4 correspond to (GGS1) and (GGS2), respectively.

5. Regularity Conditions for Asymptotic Properties of Latent Vectors

Aiming to generalize the CS approach, we provide conditions for APN of (multivariate) LTs that correspond one to one to the conditions of CS for univariate LTs, which is the standard approach for IRT models, as alternatives to the conditions imposed in Ghosal et al. (1995). Throughout, we assume that ${\left\{Y_i\right\}}_{i\in N}\sim P\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}\sim {\mathcal {P}}(\varvec{\eta })$$\end{document}, i.e., ${\left\{Y_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}$$\end{document} are Bernoulli random variables fulfilling (1) and (2), and that the true latent vector ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} lies in the interior of the parameter space, i.e., ${\eta }_0\in \Theta \backslash \partial \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0\in \Theta \setminus \partial \Theta $$\end{document}, where $\partial \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\partial \Theta $$\end{document} denotes the boundary of $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document}. The asymptotic results of Sect. 7 rely on the following regularity conditions.

1. (CS1')

   1. [i] The set $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document} is closed, convex and has non-empty interior.

   2. [ii] The prior density $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document} of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} is proper and continuous at ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} with $h\left({\eta }_0\right)>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}(\varvec{\eta }_0)>0$$\end{document}.

2. (CS2') $P_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i$$\end{document} is thrice continuously differentiable, $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}. If restricted to a compact subset $K\subseteq \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K\subseteq \Theta $$\end{document}, all $\left(\frac{\partial P_i}{\partial {\eta }_k}\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| \frac{\partial P_i}{\partial \eta _k} \right| $$\end{document} and $\left(\frac{{\partial }^2{P}_i}{\partial {\eta }_k\partial {\eta }_j}\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left| \frac{\partial ^2 P_i}{\partial \eta _k\partial \eta _j} \right| $$\end{document} are uniformly bounded for all $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, $1\le j,k\le q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\le j,k\le q$$\end{document}. Moreover, there exist constants $0<{\zeta }_0\left(K\right)<{\zeta }_1\left(K\right)<1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0<\zeta _0(K)<\zeta _1(K)<1$$\end{document}, which are independent of $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, such that

   (19) $\begin{array}{c}{\zeta }_0\left(K\right)\le \underset{(i,\eta )\in N\times K}{inf}{P}_i\left(\eta \right)\le \underset{(i,\eta )\in N\times K}{sup}{P}_i\left(\eta \right)\le {\zeta }_1\left(K\right).\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _0(K)\le \inf _{(i,\varvec{\eta })\in \mathbb {N}\times K} P_i(\varvec{\eta })\le \sup _{(i,\varvec{\eta })\in \mathbb {N}\times K} P_i(\varvec{\eta })\le \zeta _1(K). \end{aligned}$$\end{document}

3. (CS3') For each $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document}, $\eta \ne {\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\ne \varvec{\eta }_0$$\end{document}, there is a $c\left(\eta \right)<0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c(\varvec{\eta })<0$$\end{document} such that

   $\begin{array}{c}\underset{d\to \infty }{lim\; sup}\frac{1}{d}\sum _{i=1}^d{E}_{{\eta }_0}(log{Z}_i(\eta ,{\eta }_0))=\underset{d\to \infty }{lim\; sup}\frac{1}{d}{E}_{{\eta }_0}({\ell }^{\left(d\right)}\left(\eta \right|Y^{\left(d\right)})-{\ell }^{\left(d\right)}\left({\eta }_0\right|Y^{\left(d\right)}))\le c\left(\eta \right),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \limsup _{d\rightarrow \infty }\frac{1}{d}\sum _{i=1}^d\mathsf {E}_{\varvec{\eta }_0}(\log Z_i(\varvec{\eta },\varvec{\eta }_0))=\limsup _{d\rightarrow \infty }\frac{1}{d}\mathsf {E}_{\varvec{\eta }_0}(\ell ^{(d)}(\varvec{\eta }\,|\,{\varvec{\mathsf {Y}}}^{(d)})-\ell ^{(d)}(\varvec{\eta }_0\,|\,{\varvec{\mathsf {Y}}}^{(d)}))\le c(\varvec{\eta }), \end{aligned}$$\end{document}
   and if $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document} is unbounded holds additionally
   (20) $\begin{array}{c}\underset{\eta \in \Theta \backslash B_{\delta }\left({\eta }_0\right)}{sup}c\left(\eta \right)<0,\text{for all}\delta >0,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sup _{\varvec{\eta }\in \Theta \setminus B_\delta (\varvec{\eta }_0)}c(\varvec{\eta })<0 , \quad \quad \text {for all }\delta >0 , \end{aligned}$$\end{document}
   where $B_{\delta }\left({\eta }_0\right):=\{\eta \in R^q:\Vert \eta -{\eta }_0\Vert <\delta \}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_\delta (\varvec{\eta }_0):=\{\varvec{\eta }\in \mathbb {R}^q:\Vert \varvec{\eta }-\varvec{\eta }_0\Vert <\delta \}$$\end{document} is the open ball of radius $\delta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document} and center ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}.

4. (CS4') If restricted to any compact set $K\subseteq \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$K\subseteq \Theta $$\end{document}, the following set of functions is uniformly bounded

   $\begin{array}{cc}\{\left(\frac{{\partial }^3{P}_i}{\partial {\eta }_k\partial {\eta }_g\partial {\eta }_u}\right)& :i\in N,1\le k,g,u\le q\}.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Bigg \{\left| \frac{\partial ^3 P_i}{\partial \eta _k\partial \eta _g\partial \eta _u}\right| \,&:\,i\in \mathbb {N},1\le k,g,u\le q \Bigg \}. \end{aligned}$$\end{document}

5. (CS5') For all $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document} holds

   (21) $\begin{array}{c}\underset{d\to \infty }{lim\; inf}{\nu }_{\text{min}}\left(\frac{1}{d},\sum _{i=1}^d,\nabla ,{\lambda }_i,\left(\eta \right),\nabla ,{\lambda }_i,{\left(\eta \right)}^{⊺}\right)>0,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \liminf _{d\rightarrow \infty }\,\nu _\text {min}\left( \frac{1}{d}\sum _{i=1}^d\nabla \lambda _i(\varvec{\eta })\nabla \lambda _i(\varvec{\eta })^{\intercal }\right) >0, \end{aligned}$$\end{document}
   where ${\nu }_{\text{min}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nu _\text {min}$$\end{document} denotes the smallest eigenvalue.

These regularity conditions correspond one to one to conditions (CS1)–(CS5), given in Sect. 3. For the comparison of these conditions, have in mind that convexity and connectivity are equivalent properties in $R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}$$\end{document}. The convexity condition in (CS1’) is at first place stronger but it does not impose a real practical restriction, since non-convex $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document} are only rarely needed in MIRT. Analogue to the CS-theory (s. Remark 3), conditions involving ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}, like $h\left({\eta }_0\right)>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}(\varvec{\eta }_0)>0$$\end{document}, should be interpreted as $h>0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}>0$$\end{document} almost surely. Condition (CS1’[ii]) on $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document} seems more strict than (CS1[ii]). However, Chang and Stout (1993) require additional a proper prior (s. Remark 2). Thus, under the consideration that ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} is still unknown and we consider $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document} instead of $R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}$$\end{document}, the requirements on proper priors in (CS1’[ii]) are analogue to (CS1[ii]). Finally note that in the generalization of conditions (CS3) and (CS4), some requirements have been removed as the remaining requirements on ${\lambda }_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda _i$$\end{document}, $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, and its derivatives are implied by conditions (CS2’) and (CS4’).

A common assumption in one-dimensional IRT models ( $q=1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=1$$\end{document}) is the strict monotonicity assumption (3) of $P_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i$$\end{document} in $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document}, for all $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}. Conceptually, this represents the notion that a more able subject has a higher probability of responding correct in any item of an educational test. Thus, models fulfilling this strict monotonicity assumption are easier to interpret. However, models with non-generalized-linear latent variable effects can be more adequate in practice. For example, Rizopoulos and Moustaki (2008) considered IRT models with possibly non-monotonic latent variable dependencies (like polynomial effects). Due to this reason, in order to allow for more flexible modeling options, in our semiproper centering theory, we abandon the requirement on strict monotonicity of $\eta \mapsto P_i\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\mapsto P_i(\varvec{\eta })$$\end{document}, for all $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, in each component. Since the results of Chang and Stout (1993) rely on this monotonicity assumption, the merit of the current contribution is not only the extension of the CS-results for latent vectors ( $q>1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q>1$$\end{document}) but also for univariate latent variables in case of a non-monotonic latent variable effect.

If all $P_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i$$\end{document}, $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}, are strictly monotonic in each component, then requirement (19) of condition (CS2’) is satisfied as in the univariate case. Otherwise, the requirement in (19) is generally not really restrictive; it is the technical formulation of the notion that the response probabilities can (but not necessarily have to) approach zero or one only if $\Vert \eta \Vert ⟶\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Vert \varvec{\eta }\Vert \longrightarrow \infty $$\end{document}. Assumption (20) in (CS3’) serves for ensuring the identifiability of the latent vector in case of $\Vert \eta \Vert ⟶\infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Vert \varvec{\eta }\Vert \longrightarrow \infty $$\end{document}. Hence, this condition is quite natural for a statistical model. In the univariate case, (20) of (CS3’) is implied by the strict monotonicity, too. But in contrast to (19), (20) cannot be concluded directly from the strict monotonicity of all $P_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i$$\end{document} in each component if $q>1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q>1$$\end{document}. Moreover, while a single item can suffice in the univariate case for identifiability, there are always at least q needed in the q-dimensional one. Similarly, the average test information $\frac{1}{d}{I}^{\left(d\right)}\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{d}\mathcal {I}^{(d)}(\varvec{\eta })$$\end{document} is always singular for $d<q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d<q$$\end{document}, since it is a sum of d rank-one matrices. Condition (CS5’) ensures that $\frac{1}{d}{I}^{\left(d\right)}\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{d}\mathcal {I}^{(d)}(\varvec{\eta })$$\end{document} becomes regular for $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document} and can be interpreted as a condition to ensure that the asymptotic posterior of $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} is regular q-dimensional distributed and does not have a lower dimensional support (cf. Lemma W.6 in the web-appendix).

To get a better impression of the conditions, we exemplary discuss them next for model (22), the multidimensional version of a model of Lee and Bolt (2018) and a logit model with an interaction of the latent variables (Rizopoulos, 2006).

6. Verification of the CS Regularity Conditions for Multidimensional IRT Models

We shall verify the proposed conditions (CS1’) to (CS5’) for a multidimensional version of a model by Lee and Bolt (2018) and discuss them also for the models of Pelle et al. (2016) and a logit model with interaction of the latent variables (Rizopoulos, 2006).

Consider first the IRT model by Lee and Bolt (2018) or its multidimensional version

(22) $\begin{array}{c}{P}_i\left(\eta \right)={\Phi }_1\left(\frac{{\alpha }_i^{{}^{⊺}}\eta +{\beta }_i}{\sqrt{2}{(1+exp(-{\delta }_i^{{}^{⊺}}\eta ))}^{-1/2}}\right),\eta \in \Theta ,i\in N.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} P_i(\varvec{\eta })=\Phi _1\left( \frac{\varvec{\alpha }_i^{^{\intercal }}\varvec{\eta }+\beta _i}{\sqrt{2}(1+ \exp (-{\varvec{\delta }}_i^{^{\intercal }}\varvec{\eta }))^{-1/2}}\right) ,\quad \quad \varvec{\eta }\in \Theta ,\,i\in \mathbb {N}. \end{aligned}$$\end{document}

If $H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}$$\end{document} is one of the usual structural models or any other regular distribution, for example $N_q(0,I_q)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {N}_q({\varvec{\mathsf {0}}},{\varvec{\mathsf {I}}}_q)$$\end{document}, a mixture of normals or a uniform distribution on some compact set, then (CS1’) is directly satisfied. Further, considering the model parameters as random variables, we assume that the parameter sequences ${\left\{{\alpha }_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\varvec{\alpha }_i\}_{i\in \mathbb {N}}$$\end{document} and ${\left\{{\delta }_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\varvec{\delta }_i\}_{i\in \mathbb {N}}$$\end{document} behave as they were two independent iid sequences drawn from absolutely continuous regular distributions in some bounded region in $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document} and the sequence ${\left\{{\beta }_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\beta _i\}_{i\in \mathbb {N}}$$\end{document} is in an arbitrary bounded subset of $R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}$$\end{document}, then conditions (CS2’) and (CS4’) are directly satisfied. The assumption of regular distributions with a bounded support for the model parameters is reasonable, since in IRT practice items with arbitrarily large discrimination are not realistic and items of arbitrarily high or low difficulty are avoided. Furthermore, in almost all cases, the latent vector is identifiable if q arbitrary items are given. Hence, (CS3’) is satisfied, too. The gradient of the response probabilities is given by

(23) $\begin{array}{c}\begin{array}{cc}\nabla P_i\left(\eta \right)={\varphi }_1(& \frac{{\alpha }_i^{⊺}\eta +{\beta }_i}{\sqrt{2}{(1+exp(-{\delta }_i^{⊺}\eta ))}^{-1/2}})\\ & \times (\sqrt{\frac{1+exp(-{\delta }_i^{⊺}\eta )}{2}}{\alpha }_i-\frac{exp(-{\delta }_i^{⊺}\eta )}{2\sqrt{2}}\frac{{\alpha }_i^{⊺}\eta +{\beta }_i}{\sqrt{1+exp(-{\delta }_i^{⊺}\eta )}}{\delta }_i),\end{array}\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{aligned} \nabla P_i(\varvec{\eta })=\phi _1\bigg (&\frac{\pmb \alpha _i^{\intercal }\varvec{\eta }+\beta _i}{\sqrt{2}(1+\exp (-\pmb \delta _i^{\intercal }\varvec{\eta }))^{-1/2}}\bigg )\\&\times \Bigg (\sqrt{\frac{1+\exp (-\pmb \delta _i^{\intercal }\varvec{\eta })}{2}}\pmb \alpha _i-\frac{\exp (-\pmb \delta _i^{\intercal }\varvec{\eta })}{2\sqrt{2}}\frac{\pmb \alpha _i^{\intercal }\varvec{\eta }+\beta _i}{\sqrt{1+\exp (-\pmb \delta _i^{\intercal }\varvec{\eta })}}\pmb \delta _i\Bigg ), \end{aligned} \end{aligned}$$\end{document}

for all $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document} and $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document}, where ${\varphi }_q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _q$$\end{document} denotes the pdf of $N_q(0,I_q)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {N}_q(\pmb 0,\mathrm {I}_q)$$\end{document}. In particular, we see from (23), that ${\left\{\nabla P_i\left(\eta \right)\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{\nabla P_i(\varvec{\eta })\}_{i\in \mathbb {N}}$$\end{document} behaves in almost all cases for all $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document} as an iid sequence drawn from a regular distribution with bounded support in $R^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^q$$\end{document}, since the parameters are iid distributed for all items and every $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document} is considered separately, i.e., $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} is held fixed. Exceptions are pathological cases like the one in which zero belongs to the support of the distributions of all model parameters and all model parameters equal zero, i.e., $P_i\left(\eta \right)=0.5$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_i(\varvec{\eta })=0.5$$\end{document} for all $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document} and $i\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \mathbb {N}$$\end{document}. However, the subset of such cases is of zero probability for regular continuous distributions, i.e., is a null-set. Thus,

$\begin{array}{c}\frac{1}{d}\sum _{i=1}^d\nabla P_i\left(\eta \right)\nabla P_i{\left(\eta \right)}^{⊺}\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{1}{d}\sum _{i=1}^d\nabla P_i(\varvec{\eta })\nabla P_i(\varvec{\eta })^{\intercal }\end{aligned}$$\end{document}

converges to the second moment of the distribution of $\nabla P_1\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\nabla P_1(\varvec{\eta })$$\end{document}, as a random vector formed by the multivariate transformation of the randomly selected parameter values described directly after equation (23), and is thus positive definite.

With respect to (CS5’), note that (19) in (CS2’) implies that (21) is equivalent to

(24) $\begin{array}{c}\underset{d\to \infty }{lim\; inf}{\nu }_{\text{min}}\left(\frac{1}{d},\sum _{i=1}^d,I_i,\left(\eta \right)\right)>0\text{and to}\underset{d\to \infty }{lim\; inf}{\nu }_{\text{min}}\left(\frac{1}{d},\sum _{i=1}^d,\nabla ,P_i,\left(\eta \right),\nabla ,P_i,{\left(\eta \right)}^{⊺}\right)>0,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \liminf _{d\rightarrow \infty }\,\nu _\text {min}\left( \frac{1}{d}\sum _{i=1}^d\mathcal {I}_i(\varvec{\eta })\right)>0\quad \text {and to}\quad \liminf _{d\rightarrow \infty }\,\nu _\text {min}\left( \frac{1}{d}\sum _{i=1}^d\nabla P_i(\varvec{\eta })\nabla P_i(\varvec{\eta })^{\intercal }\right) >0, \end{aligned}$$\end{document}

which in our case ensures that (CS5’) is satisfied (cf. Lemma W.6 in the web-appendix).

For illustrative purposes, consider an example with $d=30$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=30$$\end{document} and $q=2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=2$$\end{document} and model parameter values, as given in Table 1, which are independently drawn from a uniform distribution on $(-2,2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(-2,2)$$\end{document} for ${\beta }_i$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _i$$\end{document}, and on $(-0.5,1)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(-0.5,1)$$\end{document} for all other parameters.

Table 1 Hypothetical parameter values for model (22) and the first 30 items.

In Fig. 1 (top) visualizations of $E_{{\eta }_0}(log\left(Z_5(\eta ,{\eta }_0)\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}_{\varvec{\eta }_0}(\log (Z_5(\varvec{\eta },\varvec{\eta }_0)))$$\end{document} are provided for two exemplary values of ${\eta }_0\in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0\in \Theta $$\end{document} and the parameters in Table 1. In particular, we can recognize lines in $\Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta $$\end{document}, for which $E_{{\eta }_0}(log\left(Z_5(\eta ,{\eta }_0)\right))=0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}_{\varvec{\eta }_0}(\log (Z_5(\varvec{\eta },\varvec{\eta }_0)))=0$$\end{document} holds. In Fig. 1 (bottom), surfaces of $\frac{1}{30}{\sum }_{i=1}^{30}{E}_{{\eta }_0}(log\left(Z_i(\eta ,{\eta }_0)\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{30}\sum _{i=1}^{30}\mathsf {E}_{\varvec{\eta }_0}(\log (Z_i(\varvec{\eta },\varvec{\eta }_0)))$$\end{document} are illustrated for further two exemplary ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} values. The surfaces are drawn over ${[-1,1]}^2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[-1,1]^2$$\end{document}. There is a nearly parabolic surface, which illustrates that there is no reason to doubt for (20) in (CS3’).

Figure. 1 Top left: $E_{(-0.5,0.1)}(log\left(Z_5(\eta ,(-0.5,0.1))\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}_{(-0.5,0.1)}(\log (Z_5(\varvec{\eta },(-0.5,0.1))))$$\end{document}, top right: $E_{(-0.1,0.3)}(log\left(Z_9(\eta ,(-0.1,0.3))\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}_{(-0.1,0.3)}(\log (Z_9(\varvec{\eta },(-0.1,0.3))))$$\end{document}, bottom left: $\frac{1}{30}{\sum }_{i=1}^{30}{E}_{(0,0)}(log\left(Z_i(\eta ,(0,0))\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{30}\sum _{i=1}^{30}\mathsf {E}_{(0,0)}(\log (Z_i(\varvec{\eta },(0,0))))$$\end{document}, bottom right: $\frac{1}{30}{\sum }_{i=1}^{30}{E}_{(0.5,0.5)}(log\left(Z_i(\eta ,(0.5,0.5))\right))$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{30}\sum _{i=1}^{30}\mathsf {E}_{(0.5,0.5)}(\log (Z_i(\varvec{\eta },(0.5,0.5))))$$\end{document}.

Figure. 2 Minimal eigenvalue of $\frac{1}{d}{\sum }_{i=1}^d\nabla {\lambda }_i\left(\eta \right)\nabla {\lambda }_i{\left(\eta \right)}^{⊺}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{d}\sum _{i=1}^d\nabla \lambda _i(\varvec{\eta })\nabla \lambda _i(\varvec{\eta })^{\intercal }$$\end{document} for $\eta \in {[-1,1]}^2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in [-1,1]^2$$\end{document}.

Figure 2 provides the minimal smallest eigenvalue of $\frac{1}{d}{\sum }_{i=1}^d\nabla {\lambda }_i\left(\eta \right)\nabla {\lambda }_i{\left(\eta \right)}^{⊺}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{1}{d}\sum _{i=1}^d\nabla \lambda _i(\varvec{\eta })\nabla \lambda _i(\varvec{\eta })^{\intercal }$$\end{document} on ${[-1,1]}^2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[-1,1]^2$$\end{document} for $d\in \{1,\dots ,30\}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\in \{1,\ldots ,30\}$$\end{document}, cf. (CS5’). Overall, in this case, the regularity conditions can be considered as justified to apply the APN for arbitrary response patterns on the illustrated 30 items.

Conditions (CS1’) to (CS5’) for other models can be verified similarly. For example, the multidimensional Rasch model implemented in Pelle et al. (2016) is the logit model

(25) $\begin{array}{c}\mathrm{logit}\left(P(Y_i=1\mid \eta )\right):=log\left(\frac{P(Y_i=1\mid \eta )}{1-P(Y_i=1\mid \eta )}\right)={\alpha }_{i0}+{\alpha }_i^{⊺}\eta ,\eta \in R^q,i=1,\dots ,d,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {logit}(\mathsf {P}(Y_i=1\mid \varvec{\eta })):=\log \left( \frac{\mathsf {P}(Y_i=1\mid \varvec{\eta })}{1-\mathsf {P}(Y_i=1\mid \varvec{\eta })}\right) =\alpha _{i0}+\pmb \alpha _i^{\intercal }\varvec{\eta },\quad \quad \varvec{\eta }\in \mathbb {R}^q,\quad i=1,\ldots ,d, \end{aligned}$$\end{document}

where $P(Y_i=1\mid \eta )$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {P}(Y_i=1\mid \varvec{\eta })$$\end{document} is the probability of inclusion in registration i, given the vector of latent variables. In this case,

$\begin{array}{c}\nabla {\lambda }_i\left(\eta \right)={\alpha }_i,i=1,\dots ,d,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \nabla \lambda _i(\varvec{\eta })=\pmb \alpha _i,\quad \quad i=1,\ldots ,d, \end{aligned}$$\end{document}

would replace (23) for $\eta \in \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }\in \Theta $$\end{document}, while the subsequent arguments are the same as above.

Another example is the two-dimensional model

$\begin{array}{c}{\lambda }_i\left(\eta \right)={\alpha }_{i0}+{\alpha }_{i1}{\eta }_1+{\alpha }_{i2}{\eta }_2+{\alpha }_{i3}{\eta }_1{\eta }_2,i=1,\dots ,d,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lambda _i(\varvec{\eta })=\alpha _{i0}+\alpha _{i1}\eta _1+\alpha _{i2}\eta _2+\alpha _{i3}\eta _1\eta _2,\quad \quad i=1,\ldots ,d, \end{aligned}$$\end{document}

which is a logit model that contains an interaction term between the two latent variables and was considered by Rizopoulos (2006) (see Section 4). While it is still logit-linear in the model parameters ${\alpha }_{\mathrm{ij}}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _{ij}$$\end{document}, $i\in \left[d\right]$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in [d]$$\end{document}, $j\in \left[4\right]$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j\in [4]$$\end{document}, it is no longer linear in the latent variables. However, with

$\begin{array}{c}\nabla {\lambda }_i\left(\eta \right)=\left(\begin{array}{c}{\alpha }_{i1}+{\alpha }_{i3}{\eta }_2\\ {\alpha }_{i2}+{\alpha }_{i3}{\eta }_1\end{array}\right),\eta \in R^2,i=1,\dots ,d,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \nabla \lambda _i(\varvec{\eta })=\begin{pmatrix} \alpha _{i1}+\alpha _{i3}\eta _2 \\ \alpha _{i2}+\alpha _{i3}\eta _1\end{pmatrix},\quad \quad \varvec{\eta }\in \mathbb {R}^2,\quad i=1,\ldots ,d, \end{aligned}$$\end{document}

the same arguments still apply (compare also to Rizopoulos and Moustaki (2008), who discuss MIRT models within a more general form of the generalized latent variable model, allowing nonlinear effects of latent variables).

7. Main Results

Our main contribution is the generalization of Theorems 1 and 3 of Chang and Stout (1993) for $q>1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q>1$$\end{document}, under the assumptions (CS1’) to (CS5’). Furthermore, we embed the CS-approach in the GGS framework (see Theorem 5 (iii)). Similarly to Chang and Stout (1993), the consistency of the MLE is received as a by-product, along with an assertion on its existence. Additionally, the consistency of a penalized MLE is derived. The results are provided in the next theorem, while their proofs along with some preliminary required lemmas are given in appendix.

Theorem 5.

Let $Z\sim N_q(0,I_q)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {Z}}}\sim \mathcal {N}_q({\varvec{\mathsf {0}}}, {\varvec{\mathsf {I}}}_q)$$\end{document} be a q-variate standard normal distributed random vector and ${\left\{Y_i\right\}}_{i\in N}\sim P\left({\eta }_0\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}\sim {\mathcal {P}}(\varvec{\eta }_0)$$\end{document} is a sequence of binary response variables for a sequence of item response functions ${\left\{Y_i\right\}}_{i\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{Y_i\}_{i\in \mathbb {N}}$$\end{document} satisfying (CS1’[i]), (CS2’) and (CS3’) for ${\eta }_0\in \Theta \backslash \partial \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0\in \Theta \setminus \partial \Theta $$\end{document}. Then, the following statements holds:

1. (i) There is a sequence ${\left\{{\hat{\eta }}_d\right\}}_{d\in N}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\{{\varvec{\hat{\eta }}}_d\}_{d\in \mathbb {N}}$$\end{document} of measurable mappings so that

   $\begin{array}{cc}& \underset{d\to \infty }{lim}{P}_{{\eta }_0}\left(\nabla ,{\ell }^{\left(d\right)},\left({\hat{\eta }}_d\right|Y^{\left(d\right)}),=,0\right)=1,\\ & \underset{d\to \infty }{lim}{P}_{{\eta }_0}\left({\ell }^{\left(d\right)},({\hat{\eta }}_d\mid Y^{\left(d\right)}),=,\underset{\eta \in \Theta }{max},{\ell }^{\left(d\right)},(\eta \mid Y^{\left(d\right)})\right)=1\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\lim _{d\rightarrow \infty }\mathsf {P}_{\varvec{\eta }_0}\left( \nabla \ell ^{(d)}({\varvec{\hat{\eta }}}_d\,|\,{\varvec{\mathsf {Y}}}^{(d)})=\pmb 0\right) =1,\\&\lim _{d\rightarrow \infty }\mathsf {P}_{\varvec{\eta }_0}\left( \ell ^{(d)}({\varvec{\hat{\eta }}}_d\mid {\varvec{\mathsf {Y}}}^{(d)}) = \max _{\varvec{\eta }\in \Theta }\ell ^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)})\right) =1 \end{aligned}$$\end{document}
   and ${\hat{\eta }}_d\stackrel{P_{{\eta }_0}}{⟶}{\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d\overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }\varvec{\eta }_0$$\end{document} for $d\to \infty$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\rightarrow \infty $$\end{document}.

2. (ii) Statement (i) remains valid if ${\ell }^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell ^{(d)}$$\end{document} is replaced by the penalized log-likelihood

   $\begin{array}{c}{\tilde{\ell }}^{\left(d\right)}(\eta \mid Y^{\left(d\right)}))={\ell }^{\left(d\right)}(\eta \mid Y^{\left(d\right)})+log\left(W\left(\eta \right)\right),\eta \in \Theta ,d\in N,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} {\tilde{\ell }}^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)}))=\ell ^{(d)}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)}) + \log (\mathcal {W}(\varvec{\eta })),\quad \quad \varvec{\eta }\in \Theta ,\, d\in \mathbb {N}, \end{aligned}$$\end{document}
   for some continuously differentiable, positive and bounded function $W$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {W}$$\end{document}.

1. (iii) If additional (CS1’[ii]), (CS4’) and (CS5’) are satisfied, then the following statement holds: If ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} is held fix, then, for all $B\in B^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\in \mathcal {B}^q$$\end{document},

   (26) $\begin{array}{c}P\left(\left(I^{\left(d\right)},{\left({\hat{\eta }}_d\right)}^{1/2},\left(\eta ,-,{\hat{\eta }}_d\right),\in ,B,\right),,Y^{\left(d\right)}\right)\stackrel{P_{{\eta }_0}}{⟶}P(Z\in B).\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}\left( \left. \mathcal {I}^{(d)}({\varvec{\hat{\eta }}}_d)^{1/2}\left( \varvec{\eta }-{\varvec{\hat{\eta }}}_d\right) \in B \,\right| \, {\varvec{\mathsf {Y}}}^{(d)}\right) \overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }\mathsf {P}({\varvec{\mathsf {Z}}}\in B). \end{aligned}$$\end{document}
   That is, the MLE ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d$$\end{document} is a semiproper centering (cf. Definition 1). If ${\eta }_0\sim G$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0\sim \mathcal {G}$$\end{document}, where $G$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {G}$$\end{document} is an absolutely continuous proper distribution with $\mathrm{supp}\left(G\right)\subseteq \Theta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {supp}(\mathcal {G})\subseteq \Theta $$\end{document}, then furthermore
   (27) $\begin{array}{c}P\left(\left(I^{\left(d\right)},{\left({\hat{\eta }}_d\right)}^{1/2},\left(\eta ,-,{\hat{\eta }}_d\right),\in ,B,\right),,Y^{\left(d\right)}\right)\stackrel{P}{⟶}P(Z\in B),\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}\left( \left. \mathcal {I}^{(d)}({\varvec{\hat{\eta }}}_d)^{1/2}\left( \varvec{\eta }-{\varvec{\hat{\eta }}}_d\right) \in B \,\right| \, {\varvec{\mathsf {Y}}}^{(d)}\right) \overset{\mathsf {P}}{\longrightarrow }\mathsf {P}({\varvec{\mathsf {Z}}}\in B), \end{aligned}$$\end{document}
   for all $B\in B^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\in \mathcal {B}^q$$\end{document}.

Remark 5.

For $W=h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {W}=\mathfrak {h}$$\end{document}, the penalized MLE in part Theorem 5 (ii) becomes the maximum a-posteriori estimator (MAP), which is an important estimator for $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document} in IRT, also because it ensures the existence of estimates in cases the MLE becomes infinite (for example when ${\sum }_{i=1}^d{y}_i=0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sum _{i=1}^d y_i=0$$\end{document} or d). The restriction on $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document} in part (ii) is stronger than (CS1’[ii]), but still mild.

As already noted in Sect. 3, Theorem 5 (iii) can be used for the construction of credible regions for $\eta$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }$$\end{document}. Additionally, it allows the interpretation of the MLE as a Bayesian estimator of ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} and thus enables the use of ${\hat{\eta }}_d$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d$$\end{document} to derive some kind of objective posterior, in the sense that it is prior-free constructed.

An important concept in the asymptotic analysis of Bayesian procedures is the consistency of the posterior distribution, which forms a basis for the asymptotic validity of inferential methods, and is proved in Theorem 6 (i). The consistency of the EAP is stated in Theorem 6 (ii).

Theorem 6.

Consider the setup and the assumptions of Theorem 5(iii), the following statements hold:

1. (i) If ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} is held fix, then

   $\begin{array}{c}P(\eta \in B\mid Y^{\left(d\right)})\stackrel{P_{{\eta }_0}}{⟶}{\delta }_{{\eta }_0}\left(B\right):=\left(\begin{array}{cc}1,& {\eta }_0\in B\\ 0,& {\eta }_0\notin B\end{array}\right),d\to \infty ,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {P}(\varvec{\eta }\in B\mid {\varvec{\mathsf {Y}}}^{(d)})\overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }\delta _{\varvec{\eta }_0}(B):=\left\{ \begin{matrix} 1,&{}\varvec{\eta }_0\in B\\ 0,&{}\varvec{\eta }_0\not \in B \end{matrix}\right. ,\quad \quad d\rightarrow \infty , \end{aligned}$$\end{document}
   for all Borel-sets $B\in B^q$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\in \mathcal {B}^q$$\end{document} with ${\eta }_0\notin \partial B$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0\not \in \partial B$$\end{document}.

2. (ii) Suppose that ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} is held fix and that there is a continuous mapping $f:\Theta \to R$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f:\Theta \rightarrow \mathbb {R}$$\end{document} so that ${\int }_{\Theta }f\left(\eta \right)h\left(\eta \right)d\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\int _\Theta f(\varvec{\eta })\mathfrak {h}(\varvec{\eta })\,\mathrm {d}(\varvec{\eta })$$\end{document} exists. Then, the posterior expected value $E(f\left(\eta \right)\mid Y^{\left(d\right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}(f(\varvec{\eta })\mid {\varvec{\mathsf {Y}}}^{(d)})$$\end{document} exists for all $d\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\in \mathbb {N}$$\end{document} and is weakly consistent for $f\left({\eta }_0\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\varvec{\eta }_0)$$\end{document}, i.e.,

   $\begin{array}{c}E(f\left(\eta \right)\mid Y^{\left(d\right)})\stackrel{P_{{\eta }_0}}{⟶}f\left({\eta }_0\right),\text{for}d\to \infty .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathsf {E}(f(\varvec{\eta })\mid {\varvec{\mathsf {Y}}}^{(d)})\overset{\mathsf {P}_{\varvec{\eta }_0}}{\longrightarrow }f(\varvec{\eta }_0),\quad \quad \text { for }d\rightarrow \infty . \end{aligned}$$\end{document}
   If in particular $E\left(\eta \right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}(\varvec{\eta })$$\end{document} exists, then the posterior expected value $E(\eta \mid Y^{\left(d\right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathsf {E}(\varvec{\eta }\mid {\varvec{\mathsf {Y}}}^{(d)})$$\end{document} exists for all $d\in N$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\in \mathbb {N}$$\end{document} and is weakly consistent for ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document}.

8. Simulation Study

The simulation study that follows examines the convergence to zero of the error for the approximation of the MLE-centered normalized posterior by a standard normal distribution and its relation to the convergence of the MLE, for the case of a bivariate latent variable vector ( $q=2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=2$$\end{document}). Convergences are evaluated based on the following measures. For the MLE, we use the root-mean-square error

$\begin{array}{c}\mathrm{RMSE}({\hat{\eta }}_d,{\eta }_0)=\sqrt{\frac{1}{2}\left({({\hat{\eta }}_{d1}-{\eta }_{01})}^2,+,{({\hat{\eta }}_{d2}-{\eta }_{02})}^2\right)}.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {RMSE}({\varvec{\hat{\eta }}}_d, \varvec{\eta }_0) = \sqrt{\frac{1}{2}\left( ({\hat{\eta }}_{d1}-\eta _{01})^2+({\hat{\eta }}_{d2}-\eta _{02})^2\right) }. \end{aligned}$$\end{document}

For the approximation of the normalized posterior density $h^{\ast }$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^*$$\end{document} by a bivariate normal pdf ${\varphi }_2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _2$$\end{document}, we compute the density approximation error (also known as $L^1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^1$$\end{document}-distance)

$\begin{array}{c}\mathrm{DAE}(h^{\ast },{\varphi }_2)={\int }_{R^2}\left(h^{\ast },\left(\eta \right),-,{\varphi }_2,\left(\eta \right)\right)d\eta ,\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {DAE}(h^*, \phi _2)=\int _{\mathbb {R}^2}\left| h^*(\varvec{\eta })-\phi _2(\varvec{\eta })\right| \,\mathrm {d}\varvec{\eta }, \end{aligned}$$\end{document}

the Hellinger-distance

$\begin{array}{c}\mathrm{HD}(h^{\ast },{\varphi }_2)=\sqrt{\frac{1}{2}{\int }_{R^2}{\left(\sqrt{h^{\ast }\left(\eta \right)},-,\sqrt{{\varphi }_2\left(\eta \right)}\right)}^{2}d\eta },\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {HD}(h^*, \phi _2)=\sqrt{\frac{1}{2}\int _{\mathbb {R}^2}\left( \sqrt{h^*(\varvec{\eta })}-\sqrt{\phi _2(\varvec{\eta })}\right) ^2\,\mathrm {d}\varvec{\eta }} , \end{aligned}$$\end{document}

and the Kullback–Leibler divergence

$\begin{array}{c}\mathrm{KLD}(h^{\ast },{\varphi }_2)={\int }_{R^2}log\left(\frac{{\varphi }_2\left(\eta \right)}{h^{\ast }\left(\eta \right)}\right){\varphi }_2\left(\eta \right)d\eta .\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathrm {KLD}(h^*, \phi _2)=\int _{\mathbb {R}^2}\log \left( \frac{\phi _2(\varvec{\eta })}{h^*(\varvec{\eta })}\right) \phi _2(\varvec{\eta })\,\mathrm {d}\varvec{\eta }. \end{aligned}$$\end{document}

The simulation study is based on model (22) with the same item parameters across all replications, to mimic the situation that different persons respond on the same test. These are generated as in Sect. 6. For the structural model we assume $H=N_2(0,I_2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}=\mathcal {N}_2(\pmb 0,\mathrm {I}_2)$$\end{document}, resulting in $\Theta =R^2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Theta =\mathbb {R}^2$$\end{document}. The number of items d varies from 10 to 70 in steps of ten items, to mimic the asymptotic behavior with test lengthening. All involved integrals are approximated using an importance sampling Monte Carlo (MC) approximation with $N_2(0,I_2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {N}_2(\pmb 0,\mathrm {I}_2)$$\end{document} being the importance distribution.

We replicate 1000 times ( $\ell =1,\dots ,1000$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell =1,\ldots , 1000$$\end{document}) the following procedure.

1. 1. Draw ${\eta }_0^{\left(\ell \right)}\sim H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0^{(\ell )}\sim \mathcal {H}$$\end{document}.

2. 2. Draw $y^{(70,\ell )}=(y_1^{\left(\ell \right)},\dots ,y_{70}^{\left(\ell \right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(70, \ell )}=(y_1^{(\ell )},\ldots ,y_{70}^{(\ell )})$$\end{document} from model (22) with underlying true latent variable vector ${\eta }_0^{\left(\ell \right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0^{(\ell )}$$\end{document} and item parameter values as described above (setting $y_i^{\left(\ell \right)}=1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_i^{(\ell )}=1$$\end{document} if $y_i^{\ast \left(\ell \right)}<P_i\left({\eta }_0^{\left(\ell \right)}\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_i^{*(\ell )}<P_i(\varvec{\eta }_0^{(\ell )})$$\end{document} and $y_i=0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_i=0$$\end{document} otherwise, where $y_i^{\ast \left(\ell \right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$y_i^{*(\ell )}$$\end{document} is drawn from iid $U(0,1)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {U}(0,1)$$\end{document}, $i=1,\dots ,70$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=1,\ldots , 70$$\end{document}). Then set $y^{(d,\ell )}=(y_1^{\left(\ell \right)},\dots ,y_d^{\left(\ell \right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d,\ell )}=(y_1^{(\ell )},\ldots ,y_{d}^{(\ell )})$$\end{document}, for $d=10,20,\dots ,70$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=10, 20, \ldots , 70$$\end{document}.

3. 3. Compute the MLE ${\hat{\eta }}_d^{\left(\ell \right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\hat{\eta }}}_d^{(\ell )}$$\end{document} and the test information matrix $I^{\left(d\right)}\left({\hat{\eta }}_d^{\left(\ell \right)}\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}^{(d)}({\varvec{\hat{\eta }}}_d^{(\ell )})$$\end{document}, based on $y^{(d,\ell )}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varvec{\mathsf {y}}}^{(d,\ell )}$$\end{document}, for $d=10,20,\dots ,70$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=10, 20, \ldots , 70$$\end{document}.

4. 4. Derive the posterior pdf $h^{\ast \left(\ell \right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$h^{*(\ell )}$$\end{document} of the normalized latent vector $I^{\left(d\right)}\left({\hat{\eta }}_d\right)(\eta -{\hat{\eta }}_d)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {I}^{(d)}({\varvec{\hat{\eta }}}_d)\big (\varvec{\eta }-{\varvec{\hat{\eta }}}_d\big )$$\end{document}, estimating its normalization constant by a MC quadrature.

5. 5. Compute ${\mathrm{RMSE}}_{\ell }^{\left(d\right)}=\mathrm{RMSE}({\hat{\eta }}_d^{\left(\ell \right)},{\eta }_0^{\left(\ell \right)})$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {RMSE}_{\ell }^{(d)} = \mathrm {RMSE}({\varvec{\hat{\eta }}}_d^{(\ell )}, \varvec{\eta }_0^{(\ell )})$$\end{document}, ${\mathrm{DAE}}_{\ell }^{\left(d\right)}=\mathrm{DAE}(h^{\ast \left(\ell \right)},{\varphi }_2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {DAE}_{\ell }^{(d)} = \mathrm {DAE}(h^{*(\ell )}, \phi _2)$$\end{document}, ${\mathrm{HD}}_{\ell }^{\left(d\right)}=\mathrm{HD}(h^{\ast \left(\ell \right)},{\varphi }_2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {HD}_{\ell }^{(d)} = \mathrm {HD}(h^{*(\ell )}, \phi _2)$$\end{document} and ${\mathrm{KLD}}_{\ell }^{\left(d\right)}=\mathrm{KLD}(h^{\ast \left(\ell \right)},{\varphi }_2)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {KLD}_{\ell }^{(d)} = \mathrm {KLD}(h^{*(\ell )}, \phi _2)$$\end{document}, for $d=10,20,\dots ,70$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=10, 20, \ldots , 70$$\end{document}.

Our results are visualized in Fig. 3, where the box-plots of the RMS, DAE, HD and KLD values computed above are pictured, for all d values considered. As expected, all evaluation measures and their range are decreasing in d.

Figure. 3 Box-Plots of the RMSE, DAE, HD and KLD simulated values for $d=10,20,\cdots ,70$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=10, 20, \dots , 70$$\end{document}. Under every d value on the horizontal axis, the percentage of points located outside the corresponding whiskers is given.

Table 2 provides the average values of the evaluation measures, i.e., ${\overline{\mathrm{RMSE}}}^{\left(d\right)}=\frac{{\sum }_{\ell }{\mathrm{RMSE}}_{\ell }^{\left(d\right)}}{1000}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {RMSE}}^{(d)} = \frac{\sum _{\ell } {\mathrm {RMSE}_{\ell }^{(d)}}}{1000}$$\end{document}, and ${\overline{\mathrm{DAE}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {DAE}}^{(d)}$$\end{document}, ${\overline{\mathrm{HD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {HD}}^{(d)}$$\end{document}, and ${\overline{\mathrm{KLD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {KLD}}^{(d)}$$\end{document}, defined analogously. Notice that in our simulation study in case of relatively small number of items ( $d\le 30$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\le 30$$\end{document}), we observed simulation cycles for which the Kullback–Leibler divergence was numerically infinite (due to floating point arithmetics), indicating that the divergence between the two compared distributions for these cases was extremely large. In particular, this occurred in 112 cases for $d=10$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=10$$\end{document}, 20 cases for $d=20$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=20$$\end{document} and one case for $d=30$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=30$$\end{document} (out of 1000). These cases were excluded from the calculation of the corresponding average KLD-values reported in Table 2.

Table 2 Average values of the root-mean-square error (RMSE) for the MLE along with the density approximation error (DAE), the Hellinger distance (HD) and the Kullback–Leibler divergence (KLD) between the density of the MLE-centered normalized posterior distribution and a bivariate standard normal density, based on 1000 simulations of model (19) with $q=2$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q=2$$\end{document}, for different numbers of items d.

Figure 4 visualizes the relation of the divergence measures of the normalized posterior from the standardized normal distribution to the RMSE of the MLE for the values of Table 2, pictured for $d\ge 30$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\ge 30$$\end{document}. Observe that as d increases, ${\overline{\mathrm{DAE}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {DAE}}^{(d)}$$\end{document} and ${\overline{\mathrm{HD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {HD}}^{(d)}$$\end{document} are linear in ${\overline{\mathrm{RMSE}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {RMSE}}^{(d)}$$\end{document}, while ${\overline{\mathrm{KLD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {KLD}}^{(d)}$$\end{document} is linear in ${\overline{\mathrm{RMSE}}}^{\left(d\right)}/\sqrt{d}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {RMSE}}^{(d)}/\sqrt{d}$$\end{document}. This is an indication that DAE and HD have the same rate of convergence to zero as RMSE while that of KLD is scaled by $d^{-1/2}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d^{-1/2}$$\end{document}.

Figure. 4 Linear regression of ${\overline{\mathrm{DAE}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {DAE}}^{(d)}$$\end{document} and ${\overline{\mathrm{HD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {HD}}^{(d)}$$\end{document} on ${\overline{\mathrm{RMSE}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {RMSE}}^{(d)}$$\end{document} and of ${\overline{\mathrm{KLD}}}^{\left(d\right)}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {KLD}}^{(d)}$$\end{document} on ${\overline{\mathrm{RMSE}}}^{\left(d\right)}/\sqrt{d}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{\mathrm {RMSE}}^{(d)}/\sqrt{d}$$\end{document} for $d\ge 30$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\ge 30$$\end{document} (s. Table 2).

Notice that the convergence of the DAE to zero in probability is equivalent to proper centering. Thus, our simulation results suggest that the MLE is a proper centering for the multivariate version of the model of Lee and Bolt (2018) with the parameters we consider.

Simulation studies for other IRT models can be conducted similarly, expecting analogous results.

9. Discussion

In this work, we proved the APN of LTs under mild conditions that are fulfilled by a broad class of MIRT models for binary items. Furthermore, we obtained as by-products the existence and consistency of the MLE and the MAP estimator. Note that though the MLE is commonly known as consistent in IRT and MIRT settings, Sinharay (2015) indicated the lack of asymptotic results under milder conditions than some of the usual ones (such as test lengthening by strictly parallel forms). Thus, Theorem 5 (i) is a contribution toward this direction.

The distribution $G$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {G}$$\end{document} in Theorem 5 (iii) can be different from $H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}$$\end{document} used in the model. Hence, the asymptotic result above is robust to misspecifications of $H$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {H}$$\end{document} as long as the support is sufficiently large. An interesting task for further investigation, pointed out by one of the reviewers, is the study of the effect of misspecified item response functions.

Under similar mild conditions we provided results on the weak consistency of posterior distributions. In Theorem 6 (ii), we get the existence and consistency of the expected a-posteriori estimator (EAP) for estimating ${\eta }_0$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\eta }_0$$\end{document} as well as for estimating $f\left({\eta }_0\right)$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f(\varvec{\eta }_0)$$\end{document}. To the best of our knowledge, a proof of these properties in such a general setup and under comparably mild or milder conditions on the MIRT model does not exist in the related literature.

Our results are under the assumption of a proper prior $h$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathfrak {h}$$\end{document}. This is appropriate in IRT settings, where the prior is a model of the population distribution of the latent traits. However, in a Bayesian framework, improper priors can also be considered. If this is the case, the proper prior assumption in (CS1’[ii]) can be replaced by the following condition if the posterior is still proper

(28) $\begin{array}{c}\frac{1}{P^{\left(d\right)}(Y^{\left(d\right)}\mid {\eta }_0)}{\int }_{\Theta \backslash B_{\delta }\left({\eta }_0\right)}{P}^{\left(d\right)}(Y^{\left(d\right)}\mid \eta )h\left(\eta \right)d\eta =o_{P_{{\eta }_0}}(d^{-q/2}),\text{for all}\delta >0.\end{array}$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \frac{1}{P^{(d)}({\varvec{\mathsf {Y}}}^{(d)}\mid \varvec{\eta }_0)}\int _{\Theta \setminus B_\delta (\varvec{\eta }_0)}P^{(d)}({\varvec{\mathsf {Y}}}^{(d)}\mid \varvec{\eta })\mathfrak {h}(\varvec{\eta })\,\mathrm {d}\varvec{\eta }=o_{\mathsf {P}_{\varvec{\eta }_0}}\big (d^{-q/2}\big ), \ \ \ \text {for all} \ \ \ \delta >0 . \end{aligned}$$\end{document}

Condition (28) is sufficient for the derivation of the results stated here and is satisfied by a proper prior.

The APN for a univariate LT for polytomous items was discussed by Chang (1996). The extension of the results for MIRT models with polytomous items is the subject of our current research.

Here, we derived conditions for APN of LTs by generalizing the contribution of Chang and Stout (1993) to $q>1$ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q>1$$\end{document}. The methodology of GGS/IH, discussed in Sect. 4, provides a general framework for APN in various contexts, including MIRT. The results of Ghosal (1997, 1999) are helpful in deriving alternative conditions for APN tailored for MIRT models.