A Generalized Two-Dimensional Index to Measure the Degree of Deviation from Double Symmetry in Square Contingency Tables

Abstract

The double symmetry model satisfies both the symmetry and point symmetry models simultaneously. To measure the degree of deviation from the double symmetry model, a two-dimensional index that can concurrently measure the degree of deviation from symmetry and point symmetry is considered. This two-dimensional index is constructed by combining two existing indexes. Although the existing indexes are constructed using power divergence, the existing two-dimensional index that can concurrently measure both symmetries is constructed using only Kullback-Leibler information, which is a special case of power divergence. Previous studies note the importance of using several indexes of divergence to compare the degrees of deviation from a model for several square contingency tables. This study, therefore, proposes a two-dimensional index based on power divergence in order to measure deviation from double symmetry for square contingency tables. Numerical examples show the utility of the proposed two-dimensional index using two datasets.

1. Introduction

Consider an $r\times r$ square contingency table that has the same row and column classifications with nominal categories. Let ${\pi }_{ij}$ denote the probability that an observation will fall in the ith row and jth column of the table ($i=1,\dots ,r;j=1,\dots ,r$).

The symmetry (S) model proposed by Bowker [1] is defined by

$${\pi }_{ij}={\pi }_{ji}\mathrm{for}i\ne j.$$

This S model is the most commonly used model for analyzing square contingency tables [2,3,4].

The point symmetry (PS) model proposed by Wall and Lienert [5] is defined by

$${\pi }_{ij}={\pi }_{i^{*}{j}^{*}}\mathrm{for}i,j=1,\dots ,r,$$

where $i^{*}=r+1-i$ and $j^{*}=r+1-j$. This PS model assumes the point of symmetry as a center of the square contingency table.

The double symmetry (DS) model proposed by Tomizawa [6] is defined by

$${\pi }_{ij}={\pi }_{ji}={\pi }_{i^{*}{j}^{*}}={\pi }_{j^{*}{i}^{*}}\mathrm{for}i,j=1,\dots ,r.$$

This DS model indicates that both the S and PS model hold.

When a model does not hold, we may be interested in measuring the degree of deviation from the model. For square contingency tables with nominal categories, Tomizawa et al. [7] proposed an index ${\Phi }_S^{\left(\lambda \right)}$ that represents the degree of deviation from the S model, Tomizawa et al. [8] proposed an index ${\Phi }_{PS}^{\left(\lambda \right)}$ that represents the degree of deviation from the PS model, and Yamamoto et al. [9] proposed an index ${\Phi }_{DS}^{\left(\lambda \right)}$ that represents the degree of deviation from the DS model.

This study focuses on the index that represents the degree of deviation from the DS model. Although the DS model satisfies both the S and PS models simultaneously, the above index ${\Phi }_{DS}^{\left(\lambda \right)}$ cannot concurrently measure the degree of deviation from S and PS. To address this gap, Ando et al. [10] proposed a two-dimensional index that can concurrently measure those. This two-dimensional index was constructed by combining existing indexes ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$. Ando et al. [10] points out that it is necessary to construct as a two-dimensional index rather than a univariate index because existing indexes ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$ are not independent. Ando et al. [10] considered three datasets: (1) the degree of deviation from the S model is large but the degree of deviation from the PS model is small, (2) the degree of deviation from the S model is small but the degree of deviation from the PS model is large, and (3) both the degree of deviation from the S model and the PS model are large. By using these datasets which have a different structure with respect to the deviation from the DS model, Ando et al. [10] showed that the all values of the index ${\Phi }_{DS}^{\left(\lambda \right)}$ applied to these datasets are the same, whereas all the values of the two-dimensional index are different. Thus, this two-dimensional index gives more detailed results than the index ${\Phi }_{DS}^{\left(\lambda \right)}$.

On the other hand, existing indexes ${\Phi }_S^{\left(\lambda \right)}$, ${\Phi }_{PS}^{\left(\lambda \right)}$ and ${\Phi }_{DS}^{\left(\lambda \right)}$ are constructed using power divergence, while the two-dimensional index is constructed using only Kullback-Leibler information, which is a special case of power divergence. Moreover, the power divergence includes several divergences, for example, the power divergence with $\lambda =-0.5$ is equivalent to the Freeman-Tukey type divergence, the power divergence with $\lambda =1$ is equivalent to the Pearson chi-squared type divergence. For details on power divergence, see Cressie and Read [11], Read and Cressie [12]. Previous studies (e.g., [7,8]) pointed out that it is important to use several indexes of divergence to accurately measure the degree of deviation from a model. This study proposes a two-dimensional index that is constructed by combining existing indexes ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$ based on power divergence.

The rest of this paper is organized as follows. In Section 2, we propose a generalized two-dimensional index for measuring the degree of deviation from DS. In Section 3, we develop an approximate confidence region for the proposed two-dimensional index. We then use numerical examples to show the utility of the proposed two-dimensional index in Section 4. We also present results obtained by applying the proposed two-dimensional index to real data. We close with concluding remarks in Section 5.

2. Two-Dimensional Index to Measure Deviation from DS

We propose a generalized two-dimensional index for measuring deviation from DS in square contingency tables. The proposed two-dimensional index can concurrently measure the degree of deviation from S and PS. The proposed two-dimensional index is based on power divergence.

Assume that ${\pi }_{ij}+{\pi }_{ji}>0$ for all $i\ne j$, and ${\pi }_{ij}+{\pi }_{i^{*}{j}^{*}}>0$ for all $(i,j)\in E$, where

$$\begin{array}{c}\hfill E=\left\{\begin{array}{cc}\left\{\right(i,j\left)\right|i,j=1,\dots ,r;(i,j)\ne \left(\right(r+1)/2,(r+1)/2)\}\hfill & \left(r\mathrm{is}\mathrm{odd}\right),\hfill \\ \left\{\right(i,j\left)\right|i,j=1,\dots ,r\}\hfill & \left(r\mathrm{is}\mathrm{even}\right).\hfill \end{array}\right.\end{array}$$

In order to measure the degree of deviation from DS, we consider the following two-dimensional index:

$$\begin{array}{c}\hfill {\mathsf{\Psi }}^{\left(\lambda \right)}=\left(\begin{array}{c}\hfill {\Phi }_S^{\left(\lambda \right)}\\ \hfill {\Phi }_{PS}^{\left(\lambda \right)}\end{array}\right)(\lambda >-1),\end{array}$$

where indexes ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$ are those considered by Tomizawa et al. [7] and Tomizawa et al. [8], respectively (see the Appendix A and Appendix B for the details of these indexes). Note that the $\lambda$ is a real value and is chosen by the user. We recommend choosing the $\lambda$ (e.g., $-0.5,0,1$) corresponding to the famous divergence. When $\lambda =0$, the proposed two-dimensional index is equivalent to the index by Ando et al. [10]. Thus, ${\mathsf{\Psi }}^{\left(\lambda \right)}$ is a generalization of the index by Ando et al. [10]. The two-dimensional index ${\mathsf{\Psi }}^{\left(\lambda \right)}$ has the following characteristics: (i) ${\mathsf{\Psi }}^{\left(\lambda \right)}={(0,0)}^{\prime }$ if and only if the DS model holds; (ii) ${\mathsf{\Psi }}^{\left(\lambda \right)}={(1,1)}^{\prime }$ if and only if the degree of deviation from DS is maximum, in the sense that ${\pi }_{ij}={\pi }_{j^{*}{i}^{*}}=0$ (then ${\pi }_{ji}>0$ and ${\pi }_{i^{*}{j}^{*}}>0$) or ${\pi }_{ji}={\pi }_{i^{*}{j}^{*}}=0$ (then ${\pi }_{ij}>0$ and ${\pi }_{j^{*}{i}^{*}}>0$) for all $i\ne j$, and either ${\pi }_{ii}=0$ or ${\pi }_{i^{*}{i}^{*}}=0$ for $i=1,\dots ,r/2$ (when r is even) or $i=1,\dots ,(r-1)/2$ (when r is odd); (iii) ${\mathsf{\Psi }}^{\left(\lambda \right)}={(1,\ast \ast )}^{\prime }$ if and only if the degree of deviation from S is maximum and the degree of deviation from PS is not maximum, in the sense that ${\pi }_{ij}=0$ (then ${\pi }_{ji}>0$) for all $i\ne j$; and (iv) ${\mathsf{\Psi }}^{\left(\lambda \right)}={(\ast \ast ,1)}^{\prime }$ if and only if the degree of deviation from PS is maximum and the degree of deviation from S is not maximum, in the sense that ${\pi }_{ij}=0$ (then ${\pi }_{i^{*}{j}^{*}}>0$) for all $(i,j)\in E$.

3. Approximate Confidence Region for the Proposed Two-Dimensional Index

Let

$$\mathit{n}={(n_{11},n_{12},\dots ,n_{1r},n_{21},n_{22},\dots ,n_{2r},\dots ,n_{r1},n_{r2},\dots ,n_{rr})}^{\prime },$$

$$\mathbf{\pi }={({\pi }_{11},{\pi }_{12},\dots ,{\pi }_{1r},{\pi }_{21},{\pi }_{22},\dots ,{\pi }_{2r},\dots ,{\pi }_{r1},{\pi }_{r2},\dots ,{\pi }_{rr})}^{\prime }.$$

Assume that $\mathit{n}$ has a multinomial distribution with sample size N and probability vector $\mathbf{\pi }$. The $\sqrt{N}(\mathit{p}-\mathbf{\pi })$ has an asymptotically Gaussian distribution with mean zero and covariance matrix $\mathit{D}\left(\mathbf{\pi }\right)-\mathbf{\pi }{\mathbf{\pi }}^{\prime }$, where $\mathit{p}$ = $\mathit{n}/N$ and $\mathit{D}\left(\mathbf{\pi }\right)$ is a diagonal matrix with the elements of $\mathbf{\pi }$ on the main diagonal (see, e.g., Agresti [13]). We estimate ${\mathsf{\Psi }}^{\left(\lambda \right)}$ by ${\widehat{\mathsf{\Psi }}}^{\left(\lambda \right)}={({\widehat{\Phi }}_S^{\left(\lambda \right)},{\widehat{\Phi }}_{PS}^{\left(\lambda \right)})}^{\prime }$, where ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ and ${\widehat{\Phi }}_{PS}^{\left(\lambda \right)}$ are given by ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$ with ${\pi }_{ij}$ replaced by $p_{ij}$, respectively. Using the delta method (see Agresti [13]), $\sqrt{N}({\widehat{\mathsf{\Psi }}}^{\left(\lambda \right)}-{\mathsf{\Psi }}^{\left(\lambda \right)})$ has an asymptotically bivariate Gaussian distribution with mean zero and covariance matrix

$$\begin{array}{ccc}\hfill {\mathsf{\Sigma }}^{\left(\lambda \right)}& =& \left(\frac{\partial {\mathsf{\Psi }}^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)\left(\mathit{D}\left(\mathbf{\pi }\right)-\mathbf{\pi }{\mathbf{\pi }}^{\prime }\right){\left(\frac{\partial {\mathsf{\Psi }}^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)}^{\prime }\hfill \\ & =& \left(\begin{array}{cc}{\sigma }_{11}^{\left(\lambda \right)}\hfill & {\sigma }_{12}^{\left(\lambda \right)}\hfill \\ {\sigma }_{21}^{\left(\lambda \right)}\hfill & {\sigma }_{22}^{\left(\lambda \right)}\hfill \end{array}\right),\hfill \end{array}$$

with ${\sigma }_{12}^{\left(\lambda \right)}$ = ${\sigma }_{21}^{\left(\lambda \right)}$. Let

$$\begin{array}{c}\hfill \delta =\underset{i\ne j}{\sum \sum }{\pi }_{ij},\Delta =\underset{(i,j)\in E}{\sum \sum }{\pi }_{ij}.\end{array}$$

The elements ${\sigma }_{11}^{\left(\lambda \right)}$, ${\sigma }_{12}^{\left(\lambda \right)}$, and ${\sigma }_{22}^{\left(\lambda \right)}$ are expressed as follows:

$$\begin{array}{ccc}\hfill {\sigma }_{11}^{\left(\lambda \right)}& =& \left(\frac{\partial {\Phi }_S^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)\left(\mathit{D}\left(\mathbf{\pi }\right)-\mathbf{\pi }{\mathbf{\pi }}^{\prime }\right){\left(\frac{\partial {\Phi }_S^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)}^{\prime }\hfill \\ & =& \frac{1}{{\delta }^2}\left[\underset{i\ne j}{\sum \sum }{\pi }_{ij}{\left({\Omega }_{ij}^{\left(\lambda \right)}\right)}^2-\delta {\left({\Phi }_S^{\left(\lambda \right)}\right)}^2\right],\hfill \\ \hfill {\sigma }_{12}^{\left(\lambda \right)}& =& \left(\frac{\partial {\Phi }_S^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)\left(\mathit{D}\left(\mathbf{\pi }\right)-\mathbf{\pi }{\mathbf{\pi }}^{\prime }\right){\left(\frac{\partial {\Phi }_{PS}^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)}^{\prime }\hfill \\ & =& \frac{1}{\delta \Delta }\left[\underset{i\ne j}{\sum \sum }{\pi }_{ij}\left({\Omega }_{ij}^{\left(\lambda \right)}-{\Phi }_S^{\left(\lambda \right)}\right)\left(W_{ij}^{\left(\lambda \right)}-{\Phi }_{PS}^{\left(\lambda \right)}\right)\right],\hfill \end{array}$$

$$\begin{array}{ccc}\hfill {\sigma }_{22}^{\left(\lambda \right)}& =& \left(\frac{\partial {\Phi }_{PS}^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)\left(\mathit{D}\left(\mathbf{\pi }\right)-\mathbf{\pi }{\mathbf{\pi }}^{\prime }\right){\left(\frac{\partial {\Phi }_{PS}^{\left(\lambda \right)}}{\partial {\mathbf{\pi }}^{\prime }}\right)}^{\prime }\hfill \\ & =& \frac{1}{{\Delta }^2}\left[\underset{(i,j)\in E}{\sum \sum }{\pi }_{ij}{\left(W_{ij}^{\left(\lambda \right)}\right)}^2-\Delta {\left({\Phi }_{PS}^{\left(\lambda \right)}\right)}^2\right],\hfill \end{array}$$

where for $\lambda >-1$

$$\begin{array}{ccc}\hfill {\Omega }_{ij}^{\left(\lambda \right)}& =& \left\{\begin{array}{cc}{\displaystyle \frac{1}{log2}log\left(2a_{ij}^c\right)}\hfill & (\lambda =0),\hfill \\ \\ {\displaystyle \frac{1}{2^{\lambda }-1}\left[{(2a_{ij}^c)}^{\lambda }-1+\lambda a_{ji}^c\{{(2a_{ij}^c)}^{\lambda }-{(2a_{ji}^c)}^{\lambda }\}\right]}\hfill & (\lambda \ne 0),\hfill \end{array}\right.\hfill \\ \\ \hfill W_{ij}^{\left(\lambda \right)}& =& \left\{\begin{array}{cc}{\displaystyle \frac{1}{log2}log\left(2c_{ij}^c\right)}\hfill & (\lambda =0),\hfill \\ \\ {\displaystyle \frac{1}{2^{\lambda }-1}\left[{(2c_{ij}^c)}^{\lambda }-1+\lambda c_{i^{*}{j}^{*}}^c\{{(2c_{ij}^c)}^{\lambda }-{(2c_{i^{*}{j}^{*}}^c)}^{\lambda }\}\right]}\hfill & (\lambda \ne 0),\hfill \end{array}\right.\hfill \end{array}$$

with

$$\begin{array}{c}\hfill a_{ij}^c=\frac{{\pi }_{ij}}{{\pi }_{ij}+{\pi }_{ji}},c_{ij}^c=\frac{{\pi }_{ij}}{{\pi }_{ij}+{\pi }_{i^{*}{j}^{*}}}.\end{array}$$

Note that the asymptotic variances ${\sigma }_{11}^{\left(\lambda \right)}$ and ${\sigma }_{22}^{\left(\lambda \right)}$ of ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$, respectively, have been given by Tomizawa et al. [7] and Tomizawa et al. [8], however, the asymptotic covariance ${\sigma }_{12}^{\left(\lambda \right)}$ of ${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$ is first derived in this study. An approximate bivariate $100\left(1-\alpha \right)\%$ confidence region for the index ${\mathsf{\Psi }}^{\left(\lambda \right)}$ is given by

$$N{({\widehat{\mathsf{\Psi }}}^{\left(\lambda \right)}-{\mathsf{\Psi }}^{\left(\lambda \right)})}^{\prime }{\widehat{{\mathsf{\Sigma }}^{\left(\lambda \right)}}}^{-1}({\widehat{\mathsf{\Psi }}}^{\left(\lambda \right)}-{\mathsf{\Psi }}^{\left(\lambda \right)})\le {\chi }_{\left(1-\alpha ;2\right)}^2,$$

where ${\chi }_{\left(1-\alpha ;2\right)}^2$ is the upper $1-\alpha$ percentile of the central chi-square distribution with two degrees of freedom and $\widehat{{\mathsf{\Sigma }}^{\left(\lambda \right)}}$ is given by ${\mathsf{\Sigma }}^{\left(\lambda \right)}$ with ${\pi }_{ij}$ replaced by $p_{ij}$.

4. Examples

4.1. Utility of the Proposed Two-Dimensional Index

In this section, we demonstrate the usefulness employing several divergences to compare the degrees of deviation from DS in several datasets. We consider the two artificial datasets in

Table 1. Two artificial datasets.

(a)					(b)
137	71	948	986		801	247	132	104
291	605	400	997		964	973	56	406
1	450	268	361		85	952	333	393
22	645	639	124		809	697	625	727

. We compare the degrees of deviation from DS for Table 1a,b using the confidence region for ${\mathsf{\Psi }}^{\left(\lambda \right)}$.

Table 2. Estimated indexes ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ and ${\widehat{\Phi }}_{PS}^{\left(\lambda \right)}$ and estimated covariance matrix of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ applied to the data in Table 1a,b.

(a) For Table 1a
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{\mathit{P}\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
0	0.346	0.341	0.471	0.278	0.417
1	0.371	0.370	0.472	0.267	0.416
(b) For Table 1b
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{\mathit{P}\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
0	0.287	0.259	0.853	0.488	0.538
1	0.348	0.320	1.006	0.557	0.682

gives the estimated values of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ and ${\mathsf{\Sigma }}^{\left(\lambda \right)}$ for Table 1a,b.

From Figure 1, we see that the confidence regions for ${\mathsf{\Psi }}^{\left(\lambda \right)}$ do not overlap for the data in Table 1a,b. We can conclude that Table 1a,b has a different structure in the degree of deviation from DS. That is, Table 1a,b has a different structure with regard to the degree of deviation from S or PS. From Figure 1, when $\lambda =0$, we can conclude that the degree of deviation from DS for Table 1a is greater than that for Table 1b, but when $\lambda =1$, we cannot conclude this. We should, therefore, examine the value of the two-dimensional index using several $\lambda$ to compare the degrees of deviation from DS for several datasets.

4.2. Example with Real Data

Consider the data in

Table 3. The two tables below show the three-year production and price forecasts, given by experts in July 1956, and the actual production and price figures from May 1959 for a sample of about 4000 Danish factories; from Andersen [14].

(a) For prices
	Actual
Forecast	Higher	No Change	Lower	Total
Higher	209	169	6	384
No change	190	3073	184	3447
Lower	3	62	81	146
Total	402	3304	271	3977
(b) For production
	Actual
Forecast	Higher	No Change	Lower	Total
Higher	532	394	69	995
No change	447	1727	334	2508
Lower	39	230	231	500
Total	1018	2351	634	4003

, which are taken from Anderson [14].

We are interested in the DS model for these data. We define, for example, the probability that the forecast and actual figures are “No change”and “Higher”, respectively, as ${\pi }_{N,H}$, and the probability that they are “Lower”and “No change”, respectively, by ${\pi }_{L,N}$. For Table 3, we are interested in whether the forecast accuracy changes depending on the category. When the forecast accuracy does not depend on these categories, the following holds: (1) the probabilities that the categories of the forecast and the actual are the same and are equal to one another (${\pi }_{H,H}={\pi }_{L,L}$); (2) the probabilities that the difference between the categories of the forecast and the actual is one are also equal (${\pi }_{N,H}={\pi }_{H,N}={\pi }_{N,L}={\pi }_{L,N}$); and (3) the probabilities that the difference between the categories of the forecast and the actual is two are also equal (${\pi }_{H,L}={\pi }_{L,H}$). The above probability structure indicates the DS model. Moreover, we are interested in whether the degree to which the forecast accuracy depends on the categories is greater for prices than for production, or vice versa.

Table 4. Estimated indexes ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ and ${\widehat{\Phi }}_{PS}^{\left(\lambda \right)}$ and estimated covariance matrix of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ applied to Table 3a,b.

(a) For Table 3a
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{\mathit{P}\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
−0.5	0.047	0.054	0.535	0.139	0.411
0	0.077	0.089	1.305	0.315	1.029
1	0.102	0.119	2.105	0.478	1.707
(b) For Table 3b
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{\mathit{P}\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
−0.5	0.009	0.036	0.040	0.010	0.110
0	0.015	0.060	0.111	0.027	0.290
1	0.020	0.082	0.208	0.048	0.513

shows the value of ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ and ${\widehat{\Phi }}_{PS}^{\left(\lambda \right)}$. We shall compare the degrees of deviation from DS for Table 3a,b using the confidence region for ${\mathsf{\Psi }}^{\left(\lambda \right)}$. The estimates of ${\mathsf{\Sigma }}^{\left(\lambda \right)}$, applied to the data in Table 3a,b, are shown in Table 4.

Figure 2 shows the confidence regions of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ applied to the data in Table 3a,b. We see that the confidence regions of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ do not overlap with regard to several values of $\lambda$. Therefore, it may be concluded that Table 3a,b has a different structure with regard to the degree of deviation from DS, in the sense that Table 3a,b has a different structure with regard to the degree of deviation from S. However, we cannot conclude whether the degree of deviation from DS is greater for Table 3a than for Table 3b. This is because, when both the degrees of deviation from S and PS are greater for Table 3a than for Table 3b, we can conclude that the degree of deviation from DS is greater for Table 3a than for Table 3b.

Next, consider the data in

Table 5. The two tables below show the decayed teeth of 349 female 18–39-year-old patients of a dental clinic in Sapporo City, for the period 2001–2005; from Tomizawa et al. [15].

(a) For female with left and right decayed teeth
	Right
Left	0 to 4	5 to 8	9 and above	Total
0 to 4	103	45	1	149
5 to 8	35	84	33	152
9 and above	3	17	42	62
Total	141	146	76	363
(b) For female with lower and upper decayed data
	Upper
Lower	0 to 4	5 to 8	9 and above	Total
0 to 4	97	62	15	174
5 to 8	20	63	75	158
9 and above	2	6	23	31
Total	119	131	113	363

, which are taken from Tomizawa et al. [15].

We shall compare the degrees of deviation from DS for Table 5a,b using the confidence region for ${\mathsf{\Psi }}^{\left(\lambda \right)}$. The estimates of ${\mathsf{\Sigma }}^{\left(\lambda \right)}$, applied to the data in Table 5a,b, are shown in

Table 6. Estimated indexes ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ and ${\widehat{\Phi }}_{PS}^{\left(\lambda \right)}$ and estimated covariance matrix of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ applied to Table 5a,b.

(a) For Table 5a
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{{P}{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
−0.5	0.024	0.064	0.113	0.043	0.142
0	0.040	0.105	0.302	0.127	0.361
1	0.055	0.141	0.540	0.231	0.605
(b) For Table 5b
	Index		Covariate Matrix
$\mathbf{\lambda }$	${\widehat{\mathsf{\Phi }}}_{\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathsf{\Phi }}}_{\mathit{P}\mathit{S}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{11}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{12}}^{\left(\mathbf{\lambda }\right)}$	${\widehat{\mathbf{\sigma }}}_{\mathbf{22}}^{\left(\mathbf{\lambda }\right)}$
−0.5	0.281	0.233	0.962	0.541	0.472
0	0.414	0.356	1.526	0.890	0.884
1	0.501	0.445	1.718	1.067	1.127

.

Figure 3 shows the confidence regions of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ applied to the data in Table 5a,b. We see that the confidence regions of ${\mathsf{\Psi }}^{\left(\lambda \right)}$ do not overlap in both horizontal and vertical axes with regard to several values of $\lambda$. Therefore, we can conclude that the degree of deviation from DS is greater for Table 5b than for Table 5a.

5. Concluding Remarks

This study proposed a generalized two-dimensional index that concurrently measures the degree of deviation from S and PS. Since the two indexes (${\Phi }_S^{\left(\lambda \right)}$ and ${\Phi }_{PS}^{\left(\lambda \right)}$) were used to measure the degree of deviation from S and PS are not independent (${\sigma }_{12}^{\left(\lambda \right)}\ne 0$), it is necessary to concurrently measure the degree of deviation from S and PS when we measure the degree of deviation from DS. To compare degrees of deviation from DS in several datasets using the proposed two-dimensional index, we should use several $\lambda$ rather than one specified $\lambda$. Therefore, we recommend to choose the several $\lambda$ (e.g., $-0.5,0,1$) corresponding to the famous divergence.

The estimator of the proposed two-dimensional index is the unbiased estimator when the sample size is large. When the sample size is small, however, the estimator of the proposed two-dimensional index may be the biased estimator. Through simulation study, Tomizawa et al. [16] investigated the performance of the estimator ${\widehat{\Phi }}_S^{\left(\lambda \right)}$. Tomizawa et al. [16] showed that (1) when the sample size was less than 300, the estimator ${\widehat{\Phi }}_S^{\left(\lambda \right)}$ had a bias, (2) when the sample size was above 300, it had a slight bias, and (3) when the sample size was above 1000, it had almost no bias. We believe that the proposed two-dimensional estimator ${\widehat{\mathsf{\Psi }}}^{\left(\lambda \right)}$ may be similar results to the estimator ${\widehat{\Phi }}_S^{\left(\lambda \right)}$, although it is necessary to verify by simulation study. In future research, the above concern will be investigated.