Diagonals–Parameter Symmetry Model and Its Property for Square Contingency Tables with Ordinal Categories

Abstract

The diagonals–parameter symmetry (DPS) model is a proposed method for analyzing square contingency tables with ordinal categories. Previously, it was stated that the generalized DPS (DPS[f]) model was equivalent to the DPS model for any function f, but the proof was not provided. This paper presents the derivation of the DPS[f] model and the proof of the relationship between the two models. The findings offer various interpretations of the DPS model. Additionally, a new model is considered, and it is shown that the proposed model and the DPS[f] model are separable.

1. Introduction

A contingency table with identical categories for rows and columns can be produced when a categorical variable is repeatedly measured. Observations in this type of table tend to concentrate on the cells along the main diagonal. Our research focuses on applying symmetry instead of assuming independence between row and column categories. Several studies have addressed symmetry issues, such as [1,2,3,4,5,6,7,8,9].

Let X and Y represent the row and column variables for an $r\times r$ contingency table with ordinal categories. Additionally, let ${\pi }_{ij}$ represent the probability of an observation falling into the $(i,j)$th cell, where $i=1,\dots ,r$ and $j=1,\dots ,r$. The diagonals–parameter symmetry (DPS) model proposed by Goodman [10] is defined as follows.

$$\begin{array}{c}{\pi }_{ij}=\left\{\begin{array}{cc}{d}_k{\psi }_{ij}\hfill & (i<j),\hfill \\ {\psi }_{ij}\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

(1)

where ${\psi }_{ij}={\psi }_{ji}$ and $k=j-i$. The parameter $d_k$ in the DPS model represents the odds of an observation falling into cells $(i,j)$ where $j-i=k$ and $i<j$, rather than cells $(j,i)$ for $k=1,\dots ,r-1$. Moreover, the ratio between ${\pi }_{ij}$ and ${\pi }_{ji}$ can be expressed as the constant $d_k$ for $j-i=k$ and $i<j$. This ratio depends solely on the distance from the main diagonal cells.

When $d_1=d_2=\cdots =d_{r-1}=1$ in Equation (1), the DPS model reduces to the symmetry (S) model proposed by Bowker [1]. When $d_k$ is independent of i and j in Equation (1), with $d_1=\cdots =d_{r-1}$, the DPS model reduces to the conditional symmetry (CS) model proposed by McCullagh [11].

Using the f-divergence, Kateri and Papaioannou [2] proposed the generalized DPS (DPS[f]) model, defined as

$$\begin{array}{c}{\pi }_{ij}={\pi }_{ij}^S{F}^{-1}({\Delta }_k+{\zeta }_{ij})(i=1,\dots ,r;j=1,\dots ,r),\hfill \end{array}$$

(2)

where $k=i-j$, ${\pi }_{ij}^S=({\pi }_{ij}+{\pi }_{ji})/2$, ${\zeta }_{ij}={\zeta }_{ji}$ and ${\Delta }_k+{\Delta }_{-k}=0$. It should be noted that the function f is twice-differentiable and strictly convex. Additionally, $F\left(t\right)=f^{\prime }\left(t\right)$, $f\left(1\right)=0$, $f\left(0\right)={lim}_{t\to 0}f\left(t\right)$, $0·f(0/0)=0$, and $0·f(a/0)=a{\lim }_{t\to \infty }[f\left(t\right)/t]$. The model derivation is not included in their paper. They did mention that the DPS model is the closest to symmetry regarding the Kullback–Leibler distance under some conditions and that the DPS[f] model is equivalent to the DPS model. In this study, we will derive the DPS[f] model and provide proof of the relation between the two models. We can obtain various interpretations of the DPS model from the result. We discuss the necessary and sufficient condition for the S model and the property between test statistics for goodness of fit.

The paper is organized as follows: Section 2 derives Equation (2) and interprets the model from an information theory viewpoint. The proof is given that the DPS[f] model is equivalent to the DPS model regardless of the function f. Section 3 considers a new model and proves that the proposed model and the DPS[f] model are separable. A numerical example is provided in Section 4. Finally, Section 5 summarizes the paper.

2. Properties of the DPS[f] Model

Kateri and Papaioannou [2] noted that the DPS[f] model is the closest model to the S model in terms of the f-divergence under the conditions where ${\sum \sum }_{j-i=k}{\pi }_{ij}$ (and ${\sum \sum }_{i-j=k}{\pi }_{ij}$) for k=$1,\dots ,r-1$ and the sums ${\pi }_{ij}+{\pi }_{ji}$ for $i,j=1,\dots ,r$ are given. Similar research has been conducted in, for example, Ireland et al. [12], Kateri and Agresti [3], and Tahata [5]. This section derives the DPS[f] model and describes its properties.

We can obtain the following theorem, although the proof of Theorem 1 is given in Appendix A.1.

Theorem 1. 

In the class of models with given ${\sum \sum }_{i-j=k}{\pi }_{ij}$, $k\ne 0$, and ${\pi }_{ij}+{\pi }_{ji}$$(i=1,\dots ,r;j=1,\dots ,r)$, the model

$$\begin{array}{c}{\pi }_{ij}={\pi }_{ij}^S{F}^{-1}({\Delta }_k+{\zeta }_{ij})(i=1,\dots ,r;j=1,\dots ,r)\hfill \end{array}$$

with $k=i-j$, ${\zeta }_{ij}={\zeta }_{ji}$ and ${\Delta }_k+{\Delta }_{-k}=0$, is the model closest to the complete symmetry model in terms of the f-divergence.

The DPS$\left[f\right]$ model can be expressed as

$$\begin{array}{c}F\left(2{\pi }_{ij}^c\right)=\left\{\begin{array}{cc}{\gamma }_{ij}+a_k\hfill & (i<j),\hfill \\ {\gamma }_{ij}\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

(3)

where $k=j-i$, ${\gamma }_{ij}={\gamma }_{ji}$ and ${\pi }_{ij}^c={\pi }_{ij}/({\pi }_{ij}+{\pi }_{ji})$. It should be noted that ${\pi }_{ij}^c$ represents the conditional probability of an observation falling in the $(i,j)$ cell, given that it falls in either the $(i,j)$ cell or the $(j,i)$ cell. Namely, the DPS[f] model indicates that

$$\begin{array}{c}F\left(2{\pi }_{ij}^c\right)-F\left(2{\pi }_{ji}^c\right)=a_k(i<j).\hfill \end{array}$$

(4)

When $a_1=\cdots =a_{r-1}=0$, the DPS[f] model is reduced to the S model.

If $f\left(x\right)=xlog\left(x\right)$, $x>0$, then the f-divergence is reduced to the KL divergence. When we set $f\left(x\right)=xlog\left(x\right)$, Equation (3) is reduced to

$$\begin{array}{c}{\pi }_{ij}=\left\{\begin{array}{cc}{\pi }_{ij}^{S}exp({\gamma }_{ij}+a_k-1)\hfill & (i<j),\hfill \\ {\pi }_{ij}^{S}exp({\gamma }_{ij}-1)\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

where $k=j-i$ and ${\gamma }_{ij}={\gamma }_{ji}$. We shall refer to this model as the DPS_KL model. Under the DPS_KL model, the ratios of ${\pi }_{ij}$ and ${\pi }_{ji}$ for $i<j$ are expressed as

$$\begin{array}{c}\frac{{\pi }_{ij}}{{\pi }_{ji}}=d_k^{\mathrm{KL}}(i<j),\hfill \end{array}$$

(5)

where $d_k^{\mathrm{KL}}=exp\left(a_k\right)$ and $k=j-i$. Since Equation (5) indicates that the ratio of ${\pi }_{ij}$ and ${\pi }_{ji}$ depends on the distance of $k=j-i$, the DPS_KL model is equivalent to the DPS model proposed by Goodman [10]. Namely, the DPS model is the closest model to the S model in terms of the KL divergence under the conditions where ${\sum \sum }_{i-j=k}{\pi }_{ij}$, $k\ne 0$, and the sums ${\pi }_{ij}+{\pi }_{ji}$ for $i=1,\dots ,r;j=1,\dots ,r$ are given. This is a special case of Theorem 1.

If $f\left(x\right)=-log\left(x\right)$, $x>0$, then the f-divergence is reduced to the reverse KL divergence. Then, the DPS[f] model is reduced to

$$\begin{array}{c}{\pi }_{ij}=\left\{\begin{array}{cc}{\pi }_{ij}^S\left(-\frac{1}{{\gamma }_{ij}+a_k}\right)\hfill & (i<j),\hfill \\ {\pi }_{ij}^S\left(-\frac{1}{{\gamma }_{ij}}\right)\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

where $k=j-i$ and ${\gamma }_{ij}={\gamma }_{ji}$. We shall refer to this model as the DPS_RKL model. This model is the closest to the S model when the divergence is measured by the reverse KL divergence and can be expressed as

$$\begin{array}{c}\frac{1}{{\pi }_{ij}^c}-\frac{1}{{\pi }_{ji}^c}=d_k^{\mathrm{RKL}}(i<j),\hfill \end{array}$$

where $d_k^{\mathrm{RKL}}=-2a_k$ and $k=j-i$. This model indicates that the difference between inverse probabilities $1/{\pi }_{ij}^c$ and $1/{\pi }_{ji}^c$ depends on the distance of $k=j-i$.

If $f\left(x\right)={(1-x)}^2$, then the f-divergence is reduced to the ${\chi }^2$-divergence (Pearsonian distance). Then, the DPS[f] model is reduced to

$$\begin{array}{c}{\pi }_{ij}=\left\{\begin{array}{cc}{\pi }_{ij}^S\left(\frac{{\gamma }_{ij}+a_k}{2}+1\right)\hfill & (i<j),\hfill \\ {\pi }_{ij}^S\left(\frac{{\gamma }_{ij}}{2}+1\right)\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

where $k=j-i$ and ${\gamma }_{ij}={\gamma }_{ji}$. We shall refer to this model as the DPS_P model. This model is the closest to the S model when the divergence is measured by the ${\chi }^2$-divergence and can be expressed as

$$\begin{array}{c}{\pi }_{ij}^c-{\pi }_{ji}^c=d_k^{\mathrm{P}}(i<j),\hfill \end{array}$$

where $d_k^{\mathrm{P}}=a_k/4$ and $k=j-i$. This model indicates that the difference between ${\pi }_{ij}^c$ and ${\pi }_{ji}^c$ depends on the distance of $k=j-i$.

Moreover, if $f\left(x\right)={\left(\lambda (\lambda +1)\right)}^{-1}(x^{\lambda +1}-x)$, $x>0$, where $\lambda$ is a real-valued parameter, then the f-divergence is reduced to the power-divergence [13]. Then, the DPS$\left[f\right]$ model is reduced to

$$\begin{array}{c}{\pi }_{ij}=\left\{\begin{array}{cc}{\pi }_{ij}^S{\left(\lambda \left({\gamma }_{ij}+a_k\right)+\frac{1}{\lambda +1}\right)}^{\frac{1}{\lambda }}\hfill & (i<j),\hfill \\ {\pi }_{ij}^S{\left(\lambda {\gamma }_{ij}+\frac{1}{\lambda +1}\right)}^{\frac{1}{\lambda }}\hfill & (i\ge j),\hfill \end{array}\right.\hfill \end{array}$$

where $k=j-i$ and ${\gamma }_{ij}={\gamma }_{ji}$. We shall refer to this model as the DPS_PD(λ) model. This model is the closest to the S model when the power-divergence measures the divergence and can be expressed as

$$\begin{array}{c}{\left({\pi }_{ij}^c\right)}^{\lambda }-{\left({\pi }_{ji}^c\right)}^{\lambda }=d_k^{\mathrm{PD}\left(\lambda \right)}(i<j),\hfill \end{array}$$

where $d_k^{\mathrm{PD}\left(\lambda \right)}=\left(\lambda a_k\right)/2^{\lambda }$ and $k=j-i$. This model indicates that the difference between the symmetric conditional probabilities to the power of $\lambda$ depends on the distance of $k=j-i$. When we apply the DPS_PD(λ) model, we should set the value of $\lambda$.

Kateri and Papaioannou [2] reported that the DPS[f] model is equivalent to the DPS model regardless of f. That is, all the models described above (i.e., DPS_KL, DPS_RKL, DPS_P, and DPS_PD(λ)) are equivalent to the DPS model surprisingly. However, the proof was not given. We prove the following theorem.

Theorem 2. 

The DPS[f] model is equivalent to the DPS model regardless of f.

The poof is given in Appendix A.2. Theorem 2 states that the DPS model holds if and only if the DPS[f] model holds. If the DPS model fits the given dataset, we obtain various interpretations for the data.

When $a_1=\cdots =a_{r-1}$, the DPS[f] model is reduced to the conditional symmetry model based on the f-divergence (CS[f]) model. The CS[f] model is described previously Kateri and Papaioannou [2]. Additionally, Fujisawa and Tahata [14] proposed the generalization of the CS[f] model. Similarly, when $d_1=\cdots =d_{r-1}$, the DPS model is reduced to the CS model proposed by McCullagh [11]. The CS[f] model is equivalent to the CS model regardless of f (Kateri and Papaioannou [2]). Hence, Theorem 2 leads to the following result.

Corollary 1. 

The CS[f] model is equivalent to the CS model regardless of f.

3. Equivalence Conditions for Symmetry

Here, the equivalence conditions of the S model are discussed. If the S model holds, then the DPS[f] model with $a_1=\cdots =a_{r-1}=0$ holds. Conversely, if the DPS[f] model holds, then the S model does not hold generally. Therefore, we are interested in considering an additional condition to obtain the S model when the DPS[f] model holds. Other studies have discussed such conditions; see Read [15] and Tahata et al. [16].

We consider the distance global symmetry (DGS) model defined as

$$\begin{array}{c}{\delta }_k^U={\delta }_k^L(k=1,\dots ,r-1),\hfill \end{array}$$

where ${\delta }_k^U={\sum \sum }_{j-i=k}{\pi }_{ij}$, ${\delta }_k^L={\sum \sum }_{i-j=k}{\pi }_{ij}$. For $k=1,\dots ,r-1$, this model indicates that the sum of probabilities which are apart distance $k=j-i$ from main diagonal cells is equal to the sum of probabilities which are apart distance $k=i-j$ from main diagonal cells. We obtain the following theorem. (The proof is given in Appendix A.3.)

Theorem 3. 

The S model holds if and only if both the DPS[f] and DGS models hold.

Next, we consider the global symmetry (GS) model, which is defined as

$$\begin{array}{c}\underset{i<j}{\sum \sum }{\pi }_{ij}=\underset{i<j}{\sum \sum }{\pi }_{ji}.\hfill \end{array}$$

It should be noted that the DGS model implies the GS model. Read [15] noted that the S model holds if and only if both the CS and GS models hold. Fujisawa and Tahata [14] proved that the S model holds if and only if the CS[f] and GS models hold. These statements are the same as those from Corollary 1. In addition, a refined estimator for measures associated with the S, CS, and GS models was introduced by [17]. The result has a significant connection to decomposing the S model and separating the goodness-of-fit test statistic of the S model. According to Corollary 1, the refined estimator for the measure of CS can be utilized to gauge the extent of deviation from the CS[f] model.

This section proves the separation of the test statistics for the S model into those for the DPS[f] model and the DGS model. Consider a square contingency table of size $r\times r$ where $n_{ij}$ denotes the observed frequency in the cell located at the $(i,j)$ position. Assume this contingency table adheres to a multinomial distribution. In this context, let $m_{ij}$ represent the expected frequency in the $(i,j)$ cell, and ${\widehat{m}}_{ij}$ be its corresponding maximum likelihood estimate under a specified model. To test each model’s goodness of fit, we can employ the likelihood ratio chi-square statistic, denoted by $G^2\left(M\right)$. This statistic is computed using the following formula:

$$\begin{array}{c}{G}^2\left(M\right)=2\sum _{i=1}^r\sum _{j=1}^r{n}_{ij}\log \left(\frac{n_{ij}}{{\widehat{m}}_{ij}}\right).\hfill \end{array}$$

This statistic follows a chi-square distribution with the corresponding degrees of freedom (df).

It is supposed that model M₃ holds if and only if both models M₁ and M₂ hold. In this case, if the analyst has found hypothesis M3 unacceptable, their attention will move to examining components M1 and M2. For these three models, Aitchison [18] discussed the properties of the Wald test statistics, and Darroch and Silvey [19] described the properties of the likelihood ratio chi-square statistics. Assume that the following equivalence holds:

$$T\left({\mathrm{M}}_3\right)=T\left({\mathrm{M}}_1\right)+T\left({\mathrm{M}}_2\right),$$

(6)

where T is the goodness of fit test statistic and the number of df for M₃ is equal to the sum of numbers of df for M₁ and M₂. If both M₁ and M₂ are accepted with a high probability (at the $\alpha$ significance level), then M₃ is accepted. However, when (6) does not hold, an incompatible situation where both M₁ and M₂ are accepted with a high probability but M₃ is rejected may arise. In fact, Darroch and Silvey [19] showed such an interesting example. The partitions of chi-squared test statistics are also discussed in, for example, [20,21].

From Theorem 3, the S model holds if and only if the DPS[f] model and the DGS model hold. In addition, df for the DPS[f] model is $(r-1)(r-2)/2$ and that for DGS model is $(r-1)$. The df for the S model can be obtained by adding the degrees of freedom for the DPS[f] model and the DGS model. Thus, we consider partitioning test statistics.

Theorem 2 confirms that the DPS[f] model is equivalent to the DPS model. Therefore, the maximum likelihood estimates (MLEs) under the DPS[f] model are given by

$$\begin{array}{c}\left\{\begin{array}{ccc}{\widehat{m}}_{ij}\hfill & =\frac{n_k^U}{n_k^U+n_k^L}\left(n_{ij}+n_{ji}\right)\hfill & (i<j),\hfill \\ {\widehat{m}}_{ij}\hfill & =n_{ij}\hfill & (i=j),\hfill \\ {\widehat{m}}_{ij}\hfill & =\frac{n_k^L}{n_k^U+n_k^L}\left(n_{ij}+n_{ji}\right)\hfill & (i>j),\hfill \end{array}\right.\hfill \end{array}$$

(7)

where $k=|j-i|$, $n_k^U=\sum {\sum }_{k=j-i}{n}_{ij}$, and $n_k^L=\sum {\sum }_{k=j-i}{n}_{ji}$ (Goodman [10]).

Next, we consider the MLEs under the DGS model using the Lagrange function. Since the kernel of the log likelihood is ${\sum }_{i=1}^r{\sum }_{j=1}^r{n}_{ij}log{\pi }_{ij}$, Lagrange function L is written as

$$\begin{array}{c}L=\sum _{i=1}^r\sum _{j=1}^r{n}_{ij}\log {\pi }_{ij}+\lambda \left(\sum _{i=1}^r\sum _{j=1}^r{\pi }_{ij}-1\right)+\sum _{k=1}^{r-1}{\lambda }_k\left(\underset{k=j-i}{\sum \sum }\left({\pi }_{ij}-{\pi }_{ji}\right)\right).\hfill \end{array}$$

Equating the derivation of L to 0 with respect to ${\pi }_{ij}$, $\lambda$, and ${\lambda }_k$ gives

$$\begin{array}{c}\left\{\begin{array}{ccc}{\widehat{m}}_{ij}\hfill & =\frac{(n_k^U+n_k^L)n_{ij}}{2n_k^U}\hfill & (i<j),\hfill \\ {\widehat{m}}_{ij}\hfill & =n_{ij}\hfill & (i=j),\hfill \\ {\widehat{m}}_{ij}\hfill & =\frac{(n_k^U+n_k^L)n_{ij}}{2n_k^L}\hfill & (i>j),\hfill \end{array}\right.\hfill \end{array}$$

(8)

where $k=|j-i|$. It is important to note that the DPS and DGS models do not remain the same when the row and column categories are permuted. Therefore, these models should be used with data from an ordinal category.

We obtain the following equivalence from Equations (7) and (8):

$$\begin{array}{c}{G}^2\left(S\right)=G^2\left(DPS\left[f\right]\right)+G^2\left(DGS\right),\hfill \end{array}$$

because the MLEs under the S model are ${\widehat{m}}_{ij}=(n_{ij}+n_{ji})/2$. Therefore, the DPS[f] model and the DGS model are separable and exhibit independence.

Let $W\left(M\right)$ denote the Wald statistic for model M. We obtain the following theorem and prove it in Appendix A.4.

Theorem 4. 

$W\left(S\right)$ is equal to the sum of $W\left(DPS\right[f\left]\right)$ and $W\left(DGS\right)$.

4. Numerical Example

Table 1. How much influence should religious leaders and medical leaders have in government funding for decisions on stem cell research? [22].

Religious	Medical Leaders
Leaders	Great (1)	Fair (2)	Little (3)	None (4)	Total
Great (1)	36	16	7	7	66
	(36.00) a	(11.96) a	(6.22) a	(7.00) a
	(36.00) b	(60.19) b	(70.95) b	(67.00) b
Fair (2)	74	96	22	4	196
	(78.04) a	(96.00) a	(26.05) a	(4.78) a
	(42.67) b	(96.00) b	(82.76) b	(40.55) b
Little (3)	119	174	48	4	345
	(119.78) a	(169.95) a	(48.00) a	(3.99) a
	(62.59) b	(100.34) b	(48.00) b	(15.05) b
None (4)	127	93	26	18	264
	(127.00) a	(92.22) a	(26.01) a	(18.00) a
	(67.00) b	(48.91) b	(14.99) b	(18.00) b
Total	356	379	103	33	871

^a MLEs under the DPS model; ^b MLEs under the DGS model.

, which is taken from Smith et al. [22], describes the amount of influence religious leaders and medical leaders should have in government funding for decisions on stem cell research when surveying 871 people. The influence levels are divided into four categories: (1) Great influence, (2) Some influence, (3) A little influence, and (4) No influence.

The values of the likelihood ratio chi-square statistics $G^2$ and the corresponding p values for the models applied to these data are shown in

Table 2. Likelihood ratio chi-square values $G^2$ for the models applied to Table 1.

Models	df	${\mathit{G}}^2$	p-Value
S	6	545.15	<$0.0001$
DPS	3	2.45	0.4847
DGS	3	542.70	<$0.0001$

. Table 2 indicates that the sum of the test statistics DPS (i.e., DPS[f]) model and DGS model is equal to that of the S model. The S model fits the data very poorly. We can infer that the marginal distribution for religious leaders is not equal to that for medical leaders. On the other hand, the DPS model fits the data very well. The likelihood-ratio test for the null hypothesis $H_0$:$d_1=d_2=d_3=1$ uses a test statistic which is the difference between $G^2$ for the S model and the DPS model. The resulting test statistic is $545.15-2.45=542.70$ with three degrees of freedom. This indicates strong evidence of at least one difference from 1. Additionally, the DGS model fits the data poorly. From Theorem 3, the reason of the poor fit of S model is caused by the poor fit of the DGS model rather than the DPS model.

The values of MLEs of $(d_1,d_2,d_3)$ in Equation (1) are $(0.15,0.05,0.06)$. It should be noted that $(d_1,d_2,d_3)$ is equal to $(d_1^{\mathrm{KL}},d_2^{\mathrm{KL}},d_3^{\mathrm{KL}})$ in the DPS_KL model. Let $(i,j)$ denote the pair that the amount of influence religious leaders is ith level and that of medical leaders is jth level. When $k=j-i$ ($k=1,2,3$), a pair $(i,j)$ is ${\widehat{d}}_k$ times as likely as a pair $(j,i)$ on condition that a pair is $(i,j)$ or $(j,i)$. From ${\widehat{d}}_k<1$ ($k=1,2,3$), the probability distribution for religious leaders is stochastically higher than the probability distribution of medical leaders. That is, the medical leaders rather than the religious leaders should have influence in government funding for decisions on stem cell research.

Moreover, from Theorem 2, we can obtain various interpretations. Since the DPS model holds, the DPS_RKL, DPS_P, and DPS_PD(λ) models also hold. For example, we obtain

$$\begin{array}{cc}({\widehat{d}}_1^{\mathrm{RKL}},{\widehat{d}}_2^{\mathrm{RKL}},{\widehat{d}}_3^{\mathrm{RKL}})\hfill & =(6.37,19.22,18.09),\hfill \\ \hfill ({\widehat{d}}_1^{\mathrm{P}},{\widehat{d}}_2^{\mathrm{P}},{\widehat{d}}_3^{\mathrm{P}})& =(-0.73,-0.90,-0.90),\hfill \end{array}$$

and for $\lambda =3$,

$$\begin{array}{c}({\widehat{d}}_1^{\mathrm{PD}\left(3\right)},{\widehat{d}}_2^{\mathrm{PD}\left(3\right)},{\widehat{d}}_3^{\mathrm{PD}\left(3\right)})=(-0.65,-0.86,-0.85).\hfill \end{array}$$

When $k=j-i$ ($k=1,2,3$), we can infer that (i) the difference between the reciprocal of conditional probability that a pair is $(i,j)$ and the reciprocal of conditional probability that a pair is $(j,i)$ is ${\widehat{d}}_k^{\mathrm{RKL}}$ on condition that the pair is $(i,j)$ or $(j,i)$ from the DPS_RKL model, (ii) the difference between the conditional probability that a pair is $(i,j)$ and the conditional probability that a pair is $(j,i)$ is ${\widehat{d}}_k^{\mathrm{P}}$ under the same condition from the DPS_P model, and (iii) the difference between the conditional probability that a pair is $(i,j)$ to the third power and the conditional probability that a pair is $(j,i)$ to the third power is ${\widehat{d}}_k^{\mathrm{PD}\left(3\right)}$ under the same condition from the DPS_PD(3) model.

5. Concluding Remarks

This paper proves that the DPS[f] model is equivalent to the DPS model proposed by Goodman [10]. This result provides various interpretations of the DPS model. The separation of the test statistic for the S model is considered. The DPS[f] and DGS models are separable and exhibit independence. Kateri and Papaioannou [2], Kateri and Agresti [3], Tahata [5] and Fujisawa and Tahata [14] considered models based on the f-divergence for the analysis of square contingency tables with ordinal categories. In the future, it should be studied whether the model based on the f-divergence is equivalent to the conventional model.