# On Estimating the Linear-by-Linear Parameter for Ordinal Log-Linear Models: A Computational Study

Authors: Eric J. Beh; Thomas B. Farver
Journal: ISRN Computational Mathematics
 (2012)

---

## Abstract

Estimating linear-by-linear association has long been an important topic in the analysis of contingency tables. For ordinal variables, log-linear models may be used to detect the strength and magnitude of the association between such variables, and iterative procedures are traditionally used. Recently, studies have shown, by way of example, three non-iterative techniques can be used to quickly and accurately estimate the parameter. This paper provides a computational study of these procedures, and the results show that they are extremely accurate when compared with estimates obtained using Newton’s unidimensional method.


---

## Body

## 1. Introduction

During the 1970s and 1980s parameter estimation procedures for ordinal log-linear and ordinal association models gained considerable attention in the literature. For example, Haberman [1], Goodman [2, 3], Agresti [4], and Gilula and Haberman [5] all considered various approaches to obtain the maximum likelihood estimate (MLE) of the parameters of such models. Recently there has been a renewed interest paid to estimation methods for these ordinal models. Galindo-Garre and Vermunt [6] considered two adaptations of Goodman’s [2] approach while Ait-Sidi-Allal et al. [7] considered an iterative procedure based on Fisher’s scoring algorithm. These approaches can involve computationally intensive procedures and are therefore subject to a number of problems including the poor choice of starting values, or the slow (or even lack of) convergence of the algorithm.A different approach to estimating parameters from ordinal log-linear (or association) models is to consider the direct, iterative-free, estimation procedure originally proposed by Beh and Davy [8] and recently extended by Beh and Farver [9]. The latter authors presented three noniterative estimation procedures for estimating the association parameter for such models and empirically demonstrated (by considering 14 commonly analysed contingency tables of varying dimension and sample size) that these noniterative procedures can provide estimates that are statistically indistinguishable when compared with the estimate obtained using Newton’s unidimensional method.By performing a computational study, this paper will explore the accuracy and reliability of Newton’s method and four noniterative procedures on randomly generated contingency tables with a fixed association parameter. Such a study will be made by first briefly describing the ordinal log-linear model (Section2) and Newton’s method (Section 3.1). A brief mathematical description of the non-iterative techniques of Beh and Farver [9] will be made in Section 3.2. The calculation of the standard error of these estimates is considered in Section 3.3. Section 4.1 describes the mechanism by which contingency tables of varying dimension and association parameter are generated, while the results of the study are discussed in Section 4.2. Section 5 provides an overview of these findings.


## 2. The Ordinal Log-Linear Model

Suppose we consider anI×J two-way contingency table, N, where the (i,j)th cell entry is given by nij for i=1,…,I and j=1,…,J. Let the grand total of N be n, and let the matrix of joint relative frequencies be P, so that the (i,j)th cell entry is pij=nij/n and ∑i=1I∑j=1Jpij=1. Define the ith row marginal proportion by pi•=∑j=1Jpij, and define the jth column marginal proportion by p•j=∑i=1Ipij. When the row and column variables both consist of ordered categories, the logarithm of the (i,j)th cell may be reconstituted by considering the following ordinal log-linear model:(1)ln⁡mij=μ+αi+βj+ϕ(ui-u̅)(vj-v̅),
where mij is the expected cell frequency of the (i,j)th cell. For this model μ is the grand mean and αi and βj are the main effects associated with the ith row category and jth column category, respectively, and are constrained so that ∑i=1Iαi=∑j=1Jβj=0. These parameters may be estimated such that(2)μ̂=ln⁡(n)+1I∑i=1Iln⁡(pi•)+1J∑j=1Jln⁡(p•j),α̂i=ln⁡(pi•)-1I∑i=1Iln⁡(pi•),β̂j=ln⁡(p•j)-1J∑j=1Jln⁡(p•j),see Beh and Davy [8]. Model (1) requires the specification of ui which is the score associated with the ith row category. When the row categories are ordered, these scores can be arranged in ascending, or descending, order depending on the structure of the row variable. Similarly, vj is the score associated with the jth column category for an ordinal column variable; the scores are arranged in ascending (or descending) order. For (1), u̅ and v̅ are the weighted means of ui and vj, respectively, where the weights are the marginal relative proportions.The relationship between the two-ordered categorical variables may be determined by estimating the parameterϕ. This parameter is the measure of linear-by-linear association, and its interpretation is aided by considering that the log-odds ratio measure is(3)ln⁡(mijmi+1,j+1mi,j+1mi+1,j)=ϕ(ui+1-ui)(vj+1-vj).Therefore, if equidistant scores are chosen to reflect the ordered structure of the variables so thatu2-u1=⋯=uI-uI-1=1 and v2-v1=⋯=vJ-vJ-1=1, then ϕ is the common log-odds ratio of the contingency table.If one considers the log-linear model (1), then the (i,j)th cell can be reconstituted by considering the association model(4)mij=α̃iβ̃jexp⁡[ϕ(ui-u̅)(vj-v̅)],
where α̃i and β̃j are subject to the constraint ∏i=1Iα̃i=∏j=1Jβ̃j=1 and are treated here as nuisance parameters. The estimation of ϕ can be carried out by using Newton’s unidimensional method considered by Goodman [2] if only one set of scores is fixed. One may also refer to Haberman [1], Agresti ([4], Appendix B), Nelder and Wedderburn [10], and Galindo-Garre and Vermunt [6] for iterative algorithms. Recently, non-iterative approaches have been proposed for estimating ϕ. Beh and Davy [8] originally considered such an approach for estimating ϕ. This work has since been expanded upon by Beh and Farver [9] who proposed another three non-iterative procedures for estimating this parameter and discussed the calculation of the corresponding standard errors. By performing an empirical study, they demonstrated that such non-iterative procedures estimate ϕ with as much precision as Newton’s method. Such an approach allows one to easily, and directly, estimate the linear-by-linear parameter without resorting to iterative techniques. This paper continues the study undertaken by Beh and Farver [9] by providing an evaluation of the performance and accuracy of Newton’s method, and the four non-iterative procedures, for estimating the parameter ϕ from model (1).


## 3. On Estimatingϕ

### 3.1. Newton’s Unidimensional Method

There are a number of iterative procedures that may be used for estimatingϕ from the log-linear model (1). These include the iterative proportional fitting algorithm as described by Agresti ([4], pp. 64-66) and Nelder and Wedderburn’s [10] quasi-likelihood approach. In this paper we will consider Newton’s unidimensional method due to its relative simplicity and popularity. A brief description of Newton’s method for estimating ϕ will be provided here although Beh and Farver [9] described it in more detail.Suppose that each of the cell frequencies of a contingency tableN, nij, are independent Poisson random variables. Also denote the expected value of nij (under some model) to be mij. In this paper the expected value of nij is determined by considering (4). The (t+1)th iteration of Newton’s method involves updating the MLE of ϕ, denoted as ϕ̂N(t+1), using the estimate obtained from the tth iteration such that(5)ϕ̂N(t+1)=ϕ̂N(t)+∑i=1I∑j=1J(ui-u̅)(vj-v̅)(nij-mij(t))∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2mij(t)
until convergence is reached. This is achieved when |ϕ̂N(t+1)-ϕ̂N(t)|<ε, for some reasonably small value of ε. Here mij(t) is just (4) with ϕ̂N(t) in place of ϕ. When this stopping criterion is satisfied, the MLE of ϕ using Newton’s method, ϕ̂N, is the final iterative estimate ϕ̂N(t+1). This method is adopted in this paper for obtaining the MLE of ϕ where ε=0.0001. To initiate Newton’s method, a starting value of ϕ̂N(0)=0 is used which coincides with independence between the two categorical variables. The standard error of the MLE of ϕ using the iterative procedure described by (5) is(6)SE(ϕ̂N)=1∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2mij,
where mij is the estimated value of nij given ϕ̂N. There are some adaptations of Newton’s unidimensional method that can be included to entice it to converge faster, be more accurate, or both. Beh and Farver [9] consider some possible ways of doing this.


### 3.2. Noniterative Procedures

The computational study considered in this paper evaluates four non-iterative procedures for estimatingϕ. The first approach is(7)ϕ̂LogNI=1σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)(vj-v̅)ln⁡(pijpi•p•j),
where σI2=∑i=1Iui2pi•-(∑i=1Iuipi•)2 and σJ2=∑j=1Jvj2p•j-(∑j=1Jvjp•j)2 and is termed the log-noniterative estimate or LogNI estimate of ϕ. The LogNI estimate (7) is guaranteed to exist when pij>0, for all i and j, and so is restricted to the case where none of the cell frequencies observed in a contingency table are zero. When there exist zero cell frequencies in the contingency table, this problem may be simply overcome by substituting small values where the zeros are observed. Agresti ([11], p. 397) points out that in some situations, especially where a saturated model is considered, replacing zero cells with 0.5 is appropriate, although in other circumstances adding 0.5 may create problems.An alternative approach to approximating the natural logarithm present in theLogNI estimate (7) is to consider instead(8)ln⁡(pijpi•p•j)≈2(pij-pi•p•jpij+pi•p•j)which takes into account thatpij/(pi•p•j)≥0. Alternatively, the second-order approximation(9)ln⁡(pijpi•p•j)≈2[(pij-pi•p•jpij+pi•p•j)+13(pij-pi•p•jpij+pi•p•j)3]may also be considered. By adopting these two approximations of the logarithm function, two estimates ofϕ are obtained by considering(10)ϕ̂LogNI1=2σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)(vj-v̅)(pij-pi•p•jpij+pi•p•j),ϕ̂LogNI2=2σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)×(vj-v̅)[13(pij-pi•p•jpij+pi•p•j)3(pij-pi•p•jpij+pi•p•j)+13(pij-pi•p•jpij+pi•p•j)3],
which are referred to as the LogNI1 and LogNI2 estimates of ϕ, respectively. An alternative way to deal with the problem of zero cell frequencies is to use another approximation of the logarithm function that is present in the Log NI estimate (7). For example, by considering only the first term in the power series expansion of the function(11)ln⁡(pijpi•p•j)≈pijpi•p•j-1,
the LogNI estimate (7) can be approximated to yield an alternative estimate of ϕ(12)ϕ̂BDNI=1σIσJ∑i=1I∑j=1Jpij(ui-u̅)(vj-v̅).
This estimate was originally proposed in Beh and Davy [8] and so is termed the BD-noniterative estimate or simply the BDNI estimate of ϕ and is the fourth non-iterative estimation approach evaluated in the present simulation study.


### 3.3. Standard Error of the Noniterative Estimates

Beh and Farver [9] show, by considering the BDNI estimate (12), the standard error of ϕ̂BDNI is(13)SE(ϕ̂BDNI)=[σIσJn]-1.
However, (13) is equivalent to the standard error (6) when the row and column variables are completely independent. To show this, consider the standard error (6). When the two categorical variables are independent so that mij=npi•p•j, then ϕ=1. Therefore (6) can be alternatively expressed as(14)SE(ϕ̂=1)=[n∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2pi•p•j]-1=[n(∑i=1I(ui-u̅)2pi•)(∑j=1J(vj-v̅)2p•j)]-1=[nσI2σJ2]-1=[σIσJn]-1.This standard error provides an indication of the maximum standard error that one would obtain for a contingency table with sample sizen. For the computational study considered in this paper, the standard error of the non-iterative estimates is calculated using (6).


## 3.1. Newton’s Unidimensional Method

There are a number of iterative procedures that may be used for estimatingϕ from the log-linear model (1). These include the iterative proportional fitting algorithm as described by Agresti ([4], pp. 64-66) and Nelder and Wedderburn’s [10] quasi-likelihood approach. In this paper we will consider Newton’s unidimensional method due to its relative simplicity and popularity. A brief description of Newton’s method for estimating ϕ will be provided here although Beh and Farver [9] described it in more detail.Suppose that each of the cell frequencies of a contingency tableN, nij, are independent Poisson random variables. Also denote the expected value of nij (under some model) to be mij. In this paper the expected value of nij is determined by considering (4). The (t+1)th iteration of Newton’s method involves updating the MLE of ϕ, denoted as ϕ̂N(t+1), using the estimate obtained from the tth iteration such that(5)ϕ̂N(t+1)=ϕ̂N(t)+∑i=1I∑j=1J(ui-u̅)(vj-v̅)(nij-mij(t))∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2mij(t)
until convergence is reached. This is achieved when |ϕ̂N(t+1)-ϕ̂N(t)|<ε, for some reasonably small value of ε. Here mij(t) is just (4) with ϕ̂N(t) in place of ϕ. When this stopping criterion is satisfied, the MLE of ϕ using Newton’s method, ϕ̂N, is the final iterative estimate ϕ̂N(t+1). This method is adopted in this paper for obtaining the MLE of ϕ where ε=0.0001. To initiate Newton’s method, a starting value of ϕ̂N(0)=0 is used which coincides with independence between the two categorical variables. The standard error of the MLE of ϕ using the iterative procedure described by (5) is(6)SE(ϕ̂N)=1∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2mij,
where mij is the estimated value of nij given ϕ̂N. There are some adaptations of Newton’s unidimensional method that can be included to entice it to converge faster, be more accurate, or both. Beh and Farver [9] consider some possible ways of doing this.


## 3.2. Noniterative Procedures

The computational study considered in this paper evaluates four non-iterative procedures for estimatingϕ. The first approach is(7)ϕ̂LogNI=1σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)(vj-v̅)ln⁡(pijpi•p•j),
where σI2=∑i=1Iui2pi•-(∑i=1Iuipi•)2 and σJ2=∑j=1Jvj2p•j-(∑j=1Jvjp•j)2 and is termed the log-noniterative estimate or LogNI estimate of ϕ. The LogNI estimate (7) is guaranteed to exist when pij>0, for all i and j, and so is restricted to the case where none of the cell frequencies observed in a contingency table are zero. When there exist zero cell frequencies in the contingency table, this problem may be simply overcome by substituting small values where the zeros are observed. Agresti ([11], p. 397) points out that in some situations, especially where a saturated model is considered, replacing zero cells with 0.5 is appropriate, although in other circumstances adding 0.5 may create problems.An alternative approach to approximating the natural logarithm present in theLogNI estimate (7) is to consider instead(8)ln⁡(pijpi•p•j)≈2(pij-pi•p•jpij+pi•p•j)which takes into account thatpij/(pi•p•j)≥0. Alternatively, the second-order approximation(9)ln⁡(pijpi•p•j)≈2[(pij-pi•p•jpij+pi•p•j)+13(pij-pi•p•jpij+pi•p•j)3]may also be considered. By adopting these two approximations of the logarithm function, two estimates ofϕ are obtained by considering(10)ϕ̂LogNI1=2σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)(vj-v̅)(pij-pi•p•jpij+pi•p•j),ϕ̂LogNI2=2σIσJ∑i=1I∑j=1Jpi•p•j(ui-u̅)×(vj-v̅)[13(pij-pi•p•jpij+pi•p•j)3(pij-pi•p•jpij+pi•p•j)+13(pij-pi•p•jpij+pi•p•j)3],
which are referred to as the LogNI1 and LogNI2 estimates of ϕ, respectively. An alternative way to deal with the problem of zero cell frequencies is to use another approximation of the logarithm function that is present in the Log NI estimate (7). For example, by considering only the first term in the power series expansion of the function(11)ln⁡(pijpi•p•j)≈pijpi•p•j-1,
the LogNI estimate (7) can be approximated to yield an alternative estimate of ϕ(12)ϕ̂BDNI=1σIσJ∑i=1I∑j=1Jpij(ui-u̅)(vj-v̅).
This estimate was originally proposed in Beh and Davy [8] and so is termed the BD-noniterative estimate or simply the BDNI estimate of ϕ and is the fourth non-iterative estimation approach evaluated in the present simulation study.


## 3.3. Standard Error of the Noniterative Estimates

Beh and Farver [9] show, by considering the BDNI estimate (12), the standard error of ϕ̂BDNI is(13)SE(ϕ̂BDNI)=[σIσJn]-1.
However, (13) is equivalent to the standard error (6) when the row and column variables are completely independent. To show this, consider the standard error (6). When the two categorical variables are independent so that mij=npi•p•j, then ϕ=1. Therefore (6) can be alternatively expressed as(14)SE(ϕ̂=1)=[n∑i=1I∑j=1J(ui-u̅)2(vj-v̅)2pi•p•j]-1=[n(∑i=1I(ui-u̅)2pi•)(∑j=1J(vj-v̅)2p•j)]-1=[nσI2σJ2]-1=[σIσJn]-1.This standard error provides an indication of the maximum standard error that one would obtain for a contingency table with sample sizen. For the computational study considered in this paper, the standard error of the non-iterative estimates is calculated using (6).


## 4. The Computational Study

### 4.1. Randomly Generating a Contingency Table

In order to study the accuracy and reliability of Newton’s method and the four non-iterative procedures described above, we consider a simple procedure for generating a two-way contingency table of dimensionI×J given a value of ϕ and sample size n. The procedure involves selecting an independent and random sample of I  pseudo-marginal row proportions, denoted by p̃i•, from a uniform distribution from the interval [0,1]. A sample of J pseudo-marginal column proportions, p̃•j, is generated in a similar manner. The random marginal proportions of the ith row and jth column can therefore be obtained such that(15)p̆i•=p̃i•∑i=1Ip̃i•,p̆•j=p̃•j∑j=1Jp̃•j,respectively, ensuring that∑i=1Ip̆i•=1 and ∑j=1Jp̆•j=1. Given these randomly generated marginal proportions and the pre-defined scores used to reflect the ordered structure of the categories, the mean row and column scores are(16)u̅=∑i=1IsI(i)p̆i•,v̅=∑j=1JsJ(j)p̆•j,respectively. For this computational study, natural row scores,{sI(i)=i:i=1,2,…,I} and natural column scores {sJ(j)=j:j=1,2,…,J} were chosen to reflect the ordinal structure of the categories. However, any scoring procedure may be considered to reflect the ordinal structure of categorical variables such as midrank scores [12] and Nishisato scores [13]. The (i,j)th joint marginal cell frequency is then calculated based on the ordinal log-linear model (1) such that(17)ñij=exp⁡[μ+αi+βj+ϕ(ui-u̅)(vj-v̅)],whereñij is the frequency observed in the (i,j)th cell, given the randomly generated marginal proportions, the pre-determined ϕ value, and(18)μ=ln⁡(n)+1I∑i=1Iln⁡(p̆i•)+1J∑j=1Jln⁡(p̆•j),αi=ln⁡(p̆i•)-1I∑i=1Iln⁡(p̆i•),βj=ln⁡(p̆•j)-1J∑j=1Jln⁡(p̆•j).It is apparent that theñij values obtained will not be guaranteed to be integers. Therefore, each cell value is rounded to the nearest integer to obtain the randomly generated contingency table. Thus, the sample size of n specified in the algorithm is only an approximate sample size value.


### 4.2. Results

The computational study performed in this paper is based on random contingency tables (using the procedure described in Section4.1) of different dimensions with a sample size of approximately 1000. Tables were generated with the following dimensions: 2 × 2, 2 × 3, 2 × 4, 2 × 5, 3 × 3, 3 × 4, 3 × 5, 4 × 4, 4 × 5, and 5 × 5. For each dimension, 200 random tables were generated each with a fixed value of ϕ; the ϕ values considered in this paper were 0, 0.01, 0.05, and 0.1 to 1.5 at increments of 0.1. In terms of the exponent of ϕ (i.e., the common log-odds ratio of the table), these are equivalent to 1, 1.01, 1.05, and 1.11 to 4.48, respectively. In practice, the common log-odds ratio may be negative; however, for the purposes of exploring the behaviour of Newton’s unidimensional method and the four non-iterative procedures (the LogNI, LogNI1, LogNI2, and BDNI estimates) against the true value of ϕ, we restricted our attention to positive ϕ values. A preliminary investigation revealed that randomly generated contingency tables with a negative ϕ behaved in a similar manner to tables generated with a positive ϕ value.Tables1, 2, 3, 4, and 5 provide a summary of the mean estimate and standard error (in parentheses) of ϕ from the 200 randomly generated contingency tables. Table 1 summarises the results obtained using Newton’s method, and Table 2 provides a summary of the accuracy of ϕ̂LogNI (7) for estimating ϕ. Table 3 summarises the results obtained by considering ϕ̂LogNI1 (10) to estimate the parameter while Tables 4 and 5 summarise output when using ϕ̂LogNI2 (7) and ϕ̂BDNI (12), respectively. We also took into consideration the number of iterations needed for Newton’s method to converge to four decimal places. Thus, Table 6 summarises the median number of iterations needed for convergence to occur. Also summarised in this table is the range (maximum and minimum) of the number of iterations observed for the each set of 200 simulated tables.Table 1
Mean estimate (and standard error) ofϕ using Newton’s unidimensional method, ϕ̂N, based on 200 simulated contingency tables with a specified value of ϕ and n=1000.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)0.00020 (0.1160)0.00235 (0.0829)0.00042 (0.0583)−0.00023 (0.0611)−0.00029 (0.0488)−0.00007 (0.0346)0.00008 (0.0312)0.00004 (0.0239)0.00005 (0.0178)0.010.01617 (0.2765)0.00740 (0.1237)0.01503 (0.0792)0.00868 (0.0631)0.01040 (0.0631)0.01016 (0.0446)0.01012 (0.0344)0.01020 (0.0307)0.01002 (0.0241)0.01000 (0.0179)0.050.05454 (0.2188)0.05134 (0.1194)0.05046 (0.0824)0.04911 (0.0621)0.04964 (0.0615)0.05009 (0.0424)0.05021 (0.0338)0.04997 (0.0301)0.05001 (0.0236)0.05009 (0.0183)0.100.10278 (0.2165)0.09998 (0.1089)0.10087 (0.0775)0.10061 (0.0608)0.09994 (0.0690)0.09990 (0.0436)0.09980 (0.0341)0.10006 (0.0301)0.10024 (0.0234)0.10001 (0.0181)0.200.20473 (0.2048)0.20097 (0.1149)0.19999 (0.0770)0.20010 (0.0580)0.19969 (0.0605)0.19970 (0.0431)0.20021 (0.0317)0.19995 (0.0290)0.20008 (0.0223)0.19995 (0.0165)0.300.30392 (0.2016)0.30622 (0.1130)0.30312 (0.0770)0.30040 (0.0558)0.30019 (0.0616)0.30022 (0.0404)0.30008 (0.0314)0.30007 (0.0278)0.29995 (0.0202)0.30001 (0.0149)0.400.40491 (0.2295)0.39656 (0.1170)0.40184 (0.0771)0.40082 (0.0629)0.40079 (0.0649)0.40080 (0.0412)0.40030 (0.0295)0.40001 (0.0263)0.40015 (0.0187)0.40036 (0.0129)0.500.50762 (0.2353)0.50533 (0.1164)0.50139 (0.0869)0.50326 (0.0564)0.49981 (0.0591)0.49788 (0.0422)0.50018 (0.0284)0.50030 (0.0249)0.50022 (0.0167)0.50025 (0.0113)0.600.61482 (0.2229)0.60065 (0.1102)0.60149 (0.0751)0.60076 (0.0561)0.59980 (0.0577)0.59945 (0.0373)0.59976 (0.0272)0.60012 (0.0230)0.59999 (0.0156)0.60015 (0.0096)0.700.71150 (0.2131)0.70272 (0.1175)0.69969 (0.0690)0.70033 (0.0524)0.69953 (0.0547)0.70031 (0.0359)0.70049 (0.0250)0.70026 (0.0209)0.70015 (0.0139)0.69996 (0.0078)0.800.80482 (0.2773)0.80215 (0.1081)0.79924 (0.0785)0.79978 (0.0485)0.79985 (0.0524)0.80013 (0.0345)0.80002 (0.0226)0.80009 (0.0197)0.79995 (0.0117)0.80038 (0.0066)0.900.89570 (0.2133)0.89081 (0.1033)0.90338 (0.0687)0.90103 (0.0513)0.89998 (0.0519)0.90018 (0.0314)0.90058 (0.0214)0.89980 (0.0174)0.90041 (0.0102)0.90025 (0.0055)1.001.01856 (0.2358)0.99917 (0.1074)0.99251 (0.0804)1.00302 (0.0504)0.99954 (0.0517)0.99992 (0.0301)1.00014 (0.0188)0.99983 (0.0167)0.99987 (0.0091)0.99659 (0.0045)1.101.10279 (0.2169)1.10185 (0.1187)1.10167 (0.0646)1.10041 (0.0472)1.10064 (0.0475)1.10020 (0.0282)1.09998 (0.0174)1.10030 (0.0146)1.10028 (0.0076)1.09676 (0.0035)1.201.20637 (0.2047)1.20018 (0.1057)1.20161 (0.0626)1.20078 (0.0430)1.20001 (0.0444)1.20077 (0.0260)1.19968 (0.0160)1.19960 (0.0132)1.19975 (0.0066)1.19886 (0.0030)1.301.32157 (0.2372)1.30192 (0.1100)1.30332 (0.0551)1.30048 (0.0396)1.30092 (0.0452)1.30039 (0.0241)1.29982 (0.0143)1.30007 (0.0121)1.30016 (0.0057)1.28541 (0.0023)1.401.41818 (0.2433)1.40873 (0.0898)1.40110 (0.0564)1.39902 (0.0395)1.40079 (0.0423)1.40040 (0.0224)1.40003 (0.0132)1.40001 (0.0106)1.39921 (0.0051)1.38124 (0.0021)1.501.49777 (0.1907)1.50079 (0.0939)1.50105 (0.0531)1.50184 (0.0384)1.49974 (0.0427)1.49900 (0.0223)1.50036 (0.0118)1.50045 (0.0091)1.49371 (0.0045)1.47004 (0.0017)Table 2
Mean estimate (and standard error) ofϕ using ϕ̂LogNI based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00026 (0.1161)0.00350 (0.0828)0.00019 (0.0583)−0.00023 (0.0611)−0.00030 (0.0448)−0.00003 (0.0346)0.00002 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01618 (0.2760)0.00562 (0.1233)0.01509 (0.0783)0.00810 (0.0629)0.01028 (0.0631)0.01024 (0.0446)0.01012 (0.0344)0.01021 (0.0307)0.01005 (0.0241)0.01001 (0.0179)0.050.05454 (0.2187)0.05164 (0.1193)0.05092 (0.0824)0.04926 (0.0620)0.04962 (0.0615)0.05012 (0.0423)0.05022 (0.0338)0.05003 (0.0300)0.05000 (0.0236)0.05016 (0.0182)0.100.10278 (0.2163)0.10001 (0.1088)0.10092 (0.0773)0.10061 (0.0607)0.10010 (0.0689)0.09997 (0.0434)0.09989 (0.0339)0.10008 (0.0299)0.10034 (0.0231)0.10018 (0.0177)0.200.20474 (0.2045)0.20054 (0.1145)0.20012 (0.0765)0.20043 (0.0574)0.19977 (0.0600)0.19992 (0.0424)0.20054 (0.0308)0.19996 (0.0280)0.20032 (0.0212)0.20024 (0.0151)0.300.30392 (0.2007)0.30788 (0.1103)0.30350 (0.0750)0.30083 (0.0545)0.30022 (0.0604)0.30042 (0.0389)0.30045 (0.0296)0.30035 (0.0259)0.30004 (0.0179)0.30088 (0.0124)0.400.40491 (0.2281)0.39618 (0.1155)0.40358 (0.0748)0.40355 (0.0598)0.40126 (0.0629)0.40262 (0.0386)0.40080 (0.0265)0.40078 (0.0232)0.40167 (0.0151)0.40312 (0.0092)0.500.50762 (0.2308)0.50795 (0.1130)0.49748 (0.0845)0.50464 (0.0518)0.50083 (0.0561)0.49881 (0.0385)0.50105 (0.0240)0.50318 (0.0205)0.50290 (0.0122)0.50552 (0.0069)0.600.61482 (0.2188)0.60143 (0.1069)0.60722 (0.0700)0.60415 (0.0506)0.59982 (0.0535)0.59997 (0.0324)0.60070 (0.0218)0.60368 (0.0175)0.60521 (0.0102)0.61267 (0.0048)0.700.71150 (0.2063)0.70642 (0.1128)0.70040 (0.0639)0.70262 (0.0463)0.70007 (0.0493)0.70457 (0.0298)0.70614 (0.0185)0.70520 (0.0145)0.70790 (0.0079)0.71631 (0.0029)0.800.80483 (0.2665)0.80608 (0.1018)0.80105 (0.0717)0.79954 (0.0419)0.74372 (0.0469)0.80350 (0.0275)0.80484 (0.0154)0.80512 (0.0124)0.81298 (0.0055)0.82603 (0.0018)0.900.89570 (0.2076)0.90016 (0.0967)0.90264 (0.0605)0.90583 (0.0421)0.90236 (0.0442)0.90229 (0.0234)0.91092 (0.0135)0.90512 (0.0100)0.91617 (0.0039)0.93313 (0.0011)1.001.01857 (0.2247)0.99812 (0.0995)0.99674 (0.0698)1.01016 (0.0394)1.00154 (0.0429)1.00645 (0.0209)1.01148 (0.0105)1.01483 (0.0085)1.02763 (0.0028)1.03001 (0.0006)1.101.10280 (0.2085)1.09827 (0.1085)1.10866 (0.0533)1.10855 (0.0361)1.10441 (0.0376)1.11082 (0.0186)1.11704 (0.0088)1.12472 (0.0063)1.13542 (0.0018)1.12646 (0.0003)1.201.20638 (0.1958)1.20485 (0.0944)1.21104 (0.0507)1.20885 (0.0319)1.20120 (0.0338)1.21299 (0.0157)1.21861 (0.0070)1.22115 (0.0049)1.23445 (0.0012)1.23221 (0.0002)1.301.32158 (0.2227)1.30226 (0.0970)1.31646 (0.0494)1.30883 (0.0286)1.30870 (0.0334)1.31479 (0.0137)1.32367 (0.0053)1.32639 (0.0039)1.33670 (0.0008)1.31161 (0.0001)1.401.41819 (0.2224)1.40275 (0.0774)1.40607 (0.0435)1.41200 (0.0270)1.41515 (0.0299)1.41974 (0.0115)1.43487 (0.0043)1.43576 (0.0028)1.42865 (0.0006)1.36380 (0.0008)1.501.49778 (0.1785)1.50663 (0.0795)1.50531 (0.0402)1.51515 (0.0250)1.50567 (0.0290)1.50459 (0.0108)1.53475 (0.0032)1.53676 (0.0018)1.52404 (0.0004)1.43543 (0.0004)Table 3
Mean estimate (and standard error) ofϕ using ϕ̂LogNI1 based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00008 (0.1160)0.00299 (0.0829)0.00032 (0.0583)−0.00024 (0.0611)−0.00030 (0.0448)−0.00004 (0.0346)0.00005 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01610 (0.2759)0.00633 (0.1234)0.01424 (0.0785)0.00836 (0.0630)0.01029 (0.0631)0.01020 (0.0446)0.01011 (0.0344)0.01021 (0.0307)0.01004 (0.0241)0.01001 (0.0179)0.050.05446 (0.2187)0.05152 (0.1193)0.05075 (0.0824)0.04922 (0.0620)0.04962 (0.0615)0.05008 (0.0423)0.05016 (0.0338)0.04996 (0.0300)0.04992 (0.0236)0.05000 (0.0182)0.100.10231 (0.2154)0.09998 (0.1099)0.10078 (0.0773)0.10046 (0.0607)0.09993 (0.0689)0.09975 (0.0434)0.09954 (0.0339)0.09969 (0.0299)0.09964 (0.0231)0.09911 (0.0177)0.200.20451 (0.2045)0.20025 (0.1145)0.19952 (0.0765)0.19932 (0.0574)0.19893 (0.0600)0.19818 (0.0424)0.19766 (0.0308)0.19676 (0.0281)0.19492 (0.0213)0.19101 (0.0153)0.300.30337 (0.2007)0.30434 (0.1110)0.29977 (0.0753)0.29721 (0.0546)0.29743 (0.0604)0.29481 (0.0390)0.29092 (0.0297)0.28925 (0.0261)0.28149 (0.0182)0.26891 (0.0129)0.400.40352 (0.2275)0.39385 (0.1155)0.39791 (0.0749)0.39275 (0.0602)0.39458 (0.0629)0.38791 (0.0388)0.37732 (0.0269)0.37233 (0.0236)0.35294 (0.0159)0.32652 (0.0103)0.500.50350 (0.2310)0.49954 (0.1134)0.48823 (0.0846)0.48801 (0.0526)0.48693 (0.0563)0.47029 (0.0389)0.45542 (0.0248)0.44396 (0.0214)0.40896 (0.0135)0.36225 (0.0086)0.600.60742 (0.2165)0.59252 (0.1070)0.58724 (0.0709)0.57402 (0.0516)0.57580 (0.0539)0.55214 (0.0331)0.52330 (0.0231)0.50859 (0.0189)0.45332 (0.0121)0.38300 (0.0070)0.700.70298 (0.2073)0.69001 (0.1127)0.67442 (0.0644)0.65328 (0.0473)0.66050 (0.0499)0.62525 (0.0311)0.57783 (0.0204)0.55366 (0.0166)0.48359(0.0104)0.38935 (0.0054)0.800.79400 (0.2685)0.78320 (0.1027)0.75942 (0.0726)0.73482 (0.0431)0.74372 (0.0469)0.68976 (0.0291)0.62194 (0.0178)0.59053 (0.0151)0.49281 (0.0085)0.39437 (0.0043)0.900.88478 (0.2076)0.86906 (0.0972)0.84284 (0.0619)0.80984 (0.0446)0.81907 (0.0456)0.74370 (0.0256)0.65549 (0.0165)0.61944 (0.0130)0.50525 (0.0071)0.39864(0.0034)1.000.99208 (0.2242)0.95982 (0.1003)0.91896 (0.0711)0.87180 (0.0428)0.88914 (0.0447)0.78265 (0.0240)0.68186 (0.0140)0.64302 (0.0122)0.51164 (0.0061)0.39090(0.0026)1.101.07792 (0.2070)1.04284 (0.1102)0.99425 (0.0563)0.93162 (0.0397)0.94721 (0.0399)0.82129 (0.0221)0.69747 (0.0127)0.65111 (0.0103)0.50923 (0.0048)0.38387 (0.0019)1.201.17866 (0.1918)1.12713 (0.0961)1.06491 (0.0539)0.98924 (0.0355)1.00068 (0.0366)0.85517 (0.0198)0.71206 (0.0114)0.66053 (0.0091)0.51300 (0.0040)0.38100 (0.0016)1.301.27633 (0.2260)1.20889 (0.0981)1.12928 (0.0542)1.05000 (0.0322)1.05912 (0.0370)0.88164 (0.0181)0.71797 (0.0099)0.67223 (0.0081)0.51028 (0.0033)0.37375 (0.011)1.401.37037 (0.2244)1.29294 (0.0806)1.19599 (0.0470)1.07872 (0.0317)1.09808 (0.0340)0.90062 (0.0146)0.73551 (0.0089)0.66730 (0.0068)0.50864 (0.0029)0.37198 (0.0011)1.501.44834 (0.1785)1.37143 (0.0823)1.24810 (0.0439)1.12323 (0.0305)1.14909 (0.0337)0.92568 (0.0161)0.73350 (0.0077)0.66406 (0.0056)0.50923 (0.0025)0.36991 (0.0009)Table 4
Mean estimate (and standard error) ofϕ using ϕ̂LogNI2 based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00020 (0.1161)0.00331 (0.0828)0.00026 (0.0583)−0.00023 (0.0611)−0.00030 (0.0448)−0.00003 (0.0346)0.00003 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01617 (0.2760)0.00587 (0.1233)0.01482 (0.0783)0.00818 (0.0629)0.01028 (0.0631)0.01022 (0.0446)0.01012 (0.0344)0.01021 (0.0307)0.01005 (0.0241)0.01001 (0.0179)0.050.05454 (0.2187)0.05161 (0.1193)0.05075 (0.0824)0.04925 (0.0620)0.04962 (0.0615)0.05011 (0.0423)0.05021 (0.0338)0.05002 (0.0300)0.05000 (0.0236)0.05014 (0.0182)0.100.10271 (0.2162)0.10001 (0.1088)0.10089 (0.0773)0.10061 (0.0607)0.10007 (0.0689)0.09996 (0.0434)0.09987 (0.0339)0.10007 (0.0299)0.10029 (0.0231)0.10012 (0.0177)0.200.20473 (0.2045)0.20054 (0.1145)0.20012 (0.0765)0.20036 (0.0574)0.19974 (0.0600)0.19982 (0.0424)0.20031 (0.0308)0.19976 (0.0280)0.19977 (0.0212)0.19881 (0.0152)0.300.30391 (0.2007)0.30672 (0.1106)0.30252 (0.0752)0.30051 (0.0545)0.30013 (0.0604)0.30002 (0.0389)0.29933 (0.0296)0.29895 (0.0259)0.29652 (0.0179)0.29131 (0.0125)0.400.40486 (0.2281)0.39614 (0.1155)0.40305 (0.0748)0.40154 (0.0599)0.40072 (0.0629)0.40021 (0.0386)0.39647 (0.0266)0.39431 (0.0233)0.38530 (0.0154)0.36861 (0.0097)0.500.50682 (0.2309)0.50599 (0.1131)0.49682 (0.0845)0.50235 (0.0521)0.49920 (0.0561)0.49339 (0.0386)0.48907 (0.0242)0.48343 (0.0208)0.46052 (0.0128)0.42230 (0.0079)0.600.61277 (0.2177)0.60066 (0.1069)0.60309 (0.0704)0.59827 (0.0509)0.59648 (0.0540)0.58824 (0.0326)0.57337 (0.0223)0.56490 (0.0181)0.52303 (0.0112)0.45747 (0.0062)0.700.70991 (0.2067)0.70403 (0.1126)0.69708 (0.0640)0.69114 (0.0466)0.69288 (0.0494)0.67813 (0.0302)0.64850 (0.0193)0.63058 (0.0155)0.56907 (0.0094)0.47393 (0.0046)0.800.80340 (0.2670)0.80302 (0.1021)0.79317 (0.0719)0.78409 (0.0422)0.78794 (0.0462)0.76014 (0.0281)0.71204 (0.0166)0.68449 (0.0139)0.59183 (0.0075)0.48567 (0.0036)0.900.89526 (0.2076)0.89665 (0.0968)0.88885 (0.0610)0.87568 (0.0431)0.87881 (0.0447)0.83495 (0.0243)0.76251 (0.0152)0.73009 (0.0118)0.61468 (0.0061)0.49547 (0.0027)1.001.01160 (0.2247)0.99328 (0.0996)0.97792 (0.0701)0.95673 (0.0410)0.96469 (0.0435)0.89309 (0.0225)0.80583 (0.0126)0.76616 (0.0108)0.62836 (0.0051)0.48938 (0.0020)1.101.09777 (0.2076)1.08835 (0.1090)1.07277 (0.0546)1.03428 (0.0379)1.04616 (0.0384)0.94742 (0.0205)0.83311 (0.0113)0.78539 (0.0090)0.63087 (0.0040)0.48327 (0.0015)1.201.20275 (0.1938)1.18780 (0.0948)1.16055 (0.0520)1.10968 (0.0337)1.11730 (0.0350)1.00053 (0.0181)0.85947 (0.0099)0.80448 (0.0078)0.63668 (0.0032)0.48174 (0.0012)1.301.31224 (0.2236)1.28023 (0.0971)1.24276 (0.0517)1.18733 (0.0305)1.19350 (0.0351)1.04284 (0.0164)0.87401 (0.0085)0.82342 (0.0068)0.63865 (0.0026)0.47549 (0.0008)1.401.41026 (0.2230)1.37640 (0.0786)1.32384 (0.0451)1.23586 (0.0297)1.25046 (0.0320)1.07434 (0.0146)0.90116 (0.0075)0.82417 (0.0057)0.63906 (0.0023)0.47335 (0.0008)1.501.49244 (0.1784)1.46870 (0.0806)1.39468 (0.0420)1.29823 (0.0282)1.31951 (0.0314)1.10946 (0.0142)0.90476 (0.0064)0.82555 (0.0045)0.64242 (0.0019)0.47224 (0.0006)Table 5
Mean estimate (and standard error) ofϕ using ϕ̂BDNI based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00064 (0.2014)0.00017 (0.0060)0.00229 (0.0829)0.00042 (0.0583)−0.00023 (0.0611)−0.00029 (0.0448)−0.00007 (0.0346)0.00008 (0.0312)0.00004 (0.0239)0.00005 (0.0178)0.010.01578 (0.2745)0.00740 (0.1236)0.01351 (0.0787)0.00868 (0.0630)0.01038 (0.0631)0.01016 (0.0446)0.01012 (0.0344)0.01020 (0.0307)0.01002 (0.0241)0.01000 (0.0179)0.050.05537 (0.2182)0.05137 (0.1193)0.05045 (0.0823)0.04906 (0.0620)0.04959 (0.0615)0.05004 (0.0423)0.05008 (0.0338)0.04987 (0.0300)0.04982 (0.0236)0.04983 (0.0182)0.100.10060 (0.2125)0.09994 (0.1087)0.10065 (0.0773)0.10031 (0.0606)0.09990 (0.0688)0.09950 (0.0434)0.09900 (0.0339)0.09925 (0.0299)0.09871 (0.0231)0.09771 (0.0177)0.200.20554 (0.2037)0.20065 (0.1143)0.19853 (0.0765)0.19771 (0.0574)0.19736 (0.0600)0.19566 (0.0424)0.19334 (0.0308)0.19195 (0.0280)0.18773 (0.0214)0.18065 (0.0154)0.300.30442 (0.1998)0.30099 (0.1113)0.29682 (0.0754)0.29088 (0.0546)0.29373 (0.0604)0.28762 (0.0390)0.27933 (0.0299)0.27595 (0.0262)0.26320 (0.0184)0.24319 (0.0132)0.400.39883 (0.2243)0.38977 (0.1152)0.38999 (0.0749)0.38243 (0.0602)0.38611 (0.0628)0.37154 (0.0389)0.35374 (0.0272)0.34521 (0.0240)0.31750 (0.0164)0.28401 (0.0109)0.500.49941 (0.2300)0.48981 (0.1128)0.48156 (0.0833)0.46832 (0.0528)0.46943 (0.0564)0.44038 (0.0393)0.41736 (0.0253)0.39981 (0.0220)0.35530 (0.0142)0.30366 (0.0093)0.600.60336 (0.2113)0.58180 (0.1068)0.56500 (0.0712)0.54913 (0.0516)0.54854 (0.0541)0.50988 (0.0337)0.46819 (0.0238)0.44776 (0.0198)0.38563 (0.0130)0.31530 (0.0077)0.700.69972 (0.2058)0.66768 (0.1121)0.64706 (0.0645)0.60455 (0.0479)0.62207 (0.0503)0.56595 (0.0319)0.50408 (0.0214)0.47675 (0.0176)0.40331 (0.0113)0.31913 (0.0061)0.800.79011 (0.2568)0.75919 (0.1026)0.71197 (0.0728)0.67057 (0.0437)0.69012 (0.0476)0.61207 (0.0302)0.53082 (0.0190)0.50154 (0.0162)0.40365 (0.0094)0.32837 (0.0049)0.900.86600 (0.2049)0.83271 (0.0969)0.78390 (0.0627)0.73735 (0.0455)0.75035 (0.0464)0.64749 (0.0269)0.55024 (0.0178)0.51415 (0.0142)0.41282 (0.0079)0.34591 (0.0038)1.000.95778 (0.2192)0.92254 (0.0996)0.84033 (0.0739)0.77967 (0.0437)0.80391 (0.0458)0.67170 (0.0254)0.56617 (0.0153)0.53203 (0.0134)0.42429 (0.0069)0.36422 (0.0029)1.101.05198 (0.2011)0.99021 (0.1101)0.90467 (0.0571)0.81670 (0.0410)0.83903 (0.0413)0.69141 (0.0236)0.57224 (0.0140)0.53513 (0.0115)0.43699 (0.0055)0.38981 (0.0021)1.201.14028 (0.1859)1.05917 (0.0963)0.95848 (0.0551)0.85077 (0.0370)0.86977 (0.0383)0.71457 (0.0213)0.58876 (0.0126)0.54391 (0.0102)0.45929 (0.0045)0.43908 (0.0017)1.301.24366 (0.2000)1.12319 (0.0973)0.99862 (0.0553)0.88461 (0.0340)0.91884 (0.0386)0.72751 (0.0197)0.59232 (0.0111)0.55751 (0.0091)0.48817 (0.0037)0.60876 (0.0011)1.401.34106 (0.2180)1.19429 (0.0788)1.04569 (0.0489)0.90048 (0.0337)0.93586 (0.0359)0.74161 (0.0180)0.61512 (0.0100)0.58435 (0.0076)0.51219 (0.0032)0.65896 (0.0010)1.501.37638 (0.1752)1.25152 (0.0836)1.06724(0.0460)0.92840 (0.0324)0.98432 (0.0356)0.76259 (0.0177)0.63542 (0.0086)0.59487 (0.0063)0.59138 (0.0027)0.65896 (0.0007)Table 6
Median number of iterations needed for Newton’s unidimensional method to converge to four decimal places. In parenthesis is the range of iterations for the 200 simulated contingency tables.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.002222222222(1–7)(1–6)(1–9)(1–4)(1–3)(1–4)(1–3)(2-3)(1–3)(1–3)0.013333333333(2–14)(2–30)(2–95)(2–14)(2–5)(2-3)(2-3)(2-3)(2-3)(2-3)0.054444544455(3–13)(1–7)(2–6)(4–6)(4–6)(4-5)(4-5)(4-5)(4-5)(4–6)0.105555566667(3–29)(4–7)(5–8)(5–8)(4–6)(5–7)(5–7)(5–7)(6–8)(6–10)0.207778889101113(5–13)(7–12)(7–11)(7–11)(7–10)(7–11)(7–15)(8–13)(9–17)(8–21)0.3099910101213141723(3–23)(9–219)(9–128)(9–15)(9–15)(9–17)(10–22)(10–19)(12–27)(15–35)0.4011111213131619212839(7–25)(6–16)(11–25)(11–23)(11–22)(12–26)(13–33)(14–47)(16–49)(20–65)0.5013141516.5172127304468(4–261)(12–56)(13–28)(14–103)(13–30)(6–47)(18–53)(18–57)(23–80)(24–110)0.60151618212128374266113(11–260)(15–32)(16–33)(17–66)(16–47)(18–61)(21–81)(23–98)(30–139)(41–261)0.701719232827375364101211(14–737)(17–39)(19–52)(21–67)(19–43)(20–71)(26–150)(30–126)(25–195)(59–560)0.802024283434497589170359(8–407)(20–40)(23–58)(25–64)(22–49)(24–140)(32–160)(38–188)(61–352)(102–1026)0.90242934424367107.5128.5255605(9–39)(24–46)(27–716)(29–114)(26–101)(32–116)(39–207)(58–332)(62–653)(147–2129)1.00283242555390145.51723951028(16–987)(27–83)(17–101)(37–286)(31–122)(42–281)(67–349)(66–480)(99–1173)(3–3243)1.1032395268.569116211275656.51875.5(29–863)(4–1983)(38–187)(48–285)(39–181)(45–383)(57–580)(68–590)(151–2370)(8–6656)1.203647648888157.52853929992725(28–208)(36–94)(44–180)(51–201)(46–210)(69–579)(94–889)(129–907)(199–3120)(589–12917)1.30425577109.5109207420.55561571.54730(29–1239)(33–267)(53–275)(58–265)(52–222)(69–441)(135–1104)(154–1606)(260–4833)(11–23803)1.40486494143.5135.5284562822.522496224(23–2288)(45–2781)(58–262)(78–372)(63–431)(88–557)(129–1675)(187–3370)(24–7041)(3–25352)1.505675114177170369.582712333102.59769(38–92)(57–216)(74–302)(104–914)(72–624)(103–1071)(189–2641)(305–4608)(22–11773)(9–34759)Table1 and Figure 1 show that Newton’s method is very accurate for estimating ϕ. The Wald statistic of the mean parameter estimate indicates that such a procedure gives estimates of ϕ that are statistically indistinguishable from the true estimate for all sized contingency tables and ϕ values. However, the number of iterations needed for Newton’s method to estimate the parameter varies greatly depending on the size of the contingency table and the magnitude of ϕ (Table 6). Contingency tables of small dimension and small ϕ required no more than about 10 iterations for convergence to four decimal places to occur. Due to the accuracy of Newton’s method, the legends attached to Figure 1 highlighting the varying table sizes have been omitted.Figure 1
Graphical view of the performance of Newton’s unidimensional method for estimatingϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Table6 shows that for contingency tables of all dimensions considered, if ϕ is larger than 0.8, then the number of iterations for convergence was nearly always larger than 10 with a median number of iterations ranging from 20 (for a 2 × 2 contingency table) to 9769 (for a 5 × 5 table). For the 200 randomly generated 5 × 5 contingency tables, the number of iterations until convergence was achieved ranging from an economical 9 iterations to nearly 35000.Table2 and Figure 2 show that of the four non-iterative procedures, using ϕ̂LogNI to estimate ϕ came closest to matching the performance of Newton’s method. There were 15 cases out of 180 total cases where the mean estimate of ϕ obtained using ϕ̂LogNI was different from the true ϕ at the 1% level of significance. Ten of these 15 cases occurred when estimating a ϕ>1.00 in a contingency table having a dimension of 4 × 5 or 5 × 5. Due to the accuracy of the ϕ̂LogNI estimate, the legends attached to Figure 2 highlighting the varying table sizes have been omitted.Figure 2
Graphical view of the performance ofϕ̂LogNI (7) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Tables3, 4, and 5 and the corresponding figures show that the estimators ϕ̂LogNI1, ϕ̂LogNI2, and ϕ̂BDNI are stable and accurate for relatively small contingency tables and where ϕ is also small (not exceeding approximately 0.5). However they performed poorly for relatively large tables with a large parameter value, continually underestimating the parameter. Figures 3, 4 and 5 graphically verify these results. Using ϕ̂LogNI2 resulted in 30 cases (out of 180 total cases) where the mean estimate of ϕ was different from the true ϕ value (at the 1% level of significance) with 23 of these 30 cases occurring when estimating a ϕ≥1.00 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5, or 5 × 5. Using ϕ̂LogNI1 resulted in 35 cases where the mean estimate of ϕ was statistically different from the true ϕ value with 27 of these 35 cases occurring when estimating a ϕ≥0.90 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5 or 5 × 5. Using ϕ̂BDNI resulted in 43 cases where the mean estimate of ϕ was statistically different from the true ϕ value with 31 of these 43 cases occurring when estimating a ϕ≥0.80 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5 or 5 × 5.Figure 3
Graphical view of the performance ofϕ̂LogNI1 (10) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Figure 4
Graphical view of the performance ofϕ̂LogNI2 (7) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Figure 5
Graphical view of the performance ofϕ̂BDNI (12) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Two hundred (200) contingency tables of dimension 4 × 5 were randomly generated (using the algorithm described in Section4.1) with ϕ=0.4 and sample size approximately 1000. For these tables, Figure 6 illustrates the stability of each of the five estimation procedures and summarises the P value obtained from the Wald test of the estimates ϕ̂N, ϕ̂LogNI, ϕ̂LogNI1, ϕ̂LogN2, and ϕ̂BDNI when compared with ϕ=0.4. All five methods provide reliable and accurate estimates of the parameter and, for all the tables generated, are statistically indistinguishable when compared with the true value of the parameter. These figures show that Newton’s method (Figure 6(a)) provides the most reliable estimate since the P value calculated from each of the 200 tables is consistently near 1.0. Of the non-iterative procedures Figure 6(b) shows that ϕ̂LogNI is nearly as stable (although slightly more variable) when compared with Newton’s method.TheP-values obtained from the Wald Test comparing the true ϕ value against (a) ϕ̂N (b) ϕ̂LogNI, (c) ϕ̂LogNI2, (d) ϕ̂LogNI1 and (e) ϕ̂BDNI. These P-values are based on the estimates obtained from randomly generating 200 contingency tables of dimension 4 × 5 where ϕ = 0.4 and n≈1000.
(a)(b)(c)(d)(e)The three remaining non-iterative procedures also behave well, but provide poorer estimates ofϕ than Newton’s method or ϕ̂LogNI. Figure 6(c) shows that ϕ̂LogNI2 provides estimates that are closer to ϕ than ϕ̂LogNI1 Figure 6(d) shows that the P value of ϕ̂LogNI1 is centred around 0.7 while Figure 6(c) shows that the P value of ϕ̂LogNI2 is centred around 0.9. This is to be expected since ϕ̂LogNI1 is a first-order approximation of the LogNI estimate (7) while ϕ̂LogNI2 is its second-order approximation. For the 200 4 × 5 tables, these two non-iterative estimation procedures provide more accurate estimates of ϕ than ϕ̂BDNI. Figure 6(e) shows that the P value from the Wald test of the estimate is centred around 0.5.


## 4.1. Randomly Generating a Contingency Table

In order to study the accuracy and reliability of Newton’s method and the four non-iterative procedures described above, we consider a simple procedure for generating a two-way contingency table of dimensionI×J given a value of ϕ and sample size n. The procedure involves selecting an independent and random sample of I  pseudo-marginal row proportions, denoted by p̃i•, from a uniform distribution from the interval [0,1]. A sample of J pseudo-marginal column proportions, p̃•j, is generated in a similar manner. The random marginal proportions of the ith row and jth column can therefore be obtained such that(15)p̆i•=p̃i•∑i=1Ip̃i•,p̆•j=p̃•j∑j=1Jp̃•j,respectively, ensuring that∑i=1Ip̆i•=1 and ∑j=1Jp̆•j=1. Given these randomly generated marginal proportions and the pre-defined scores used to reflect the ordered structure of the categories, the mean row and column scores are(16)u̅=∑i=1IsI(i)p̆i•,v̅=∑j=1JsJ(j)p̆•j,respectively. For this computational study, natural row scores,{sI(i)=i:i=1,2,…,I} and natural column scores {sJ(j)=j:j=1,2,…,J} were chosen to reflect the ordinal structure of the categories. However, any scoring procedure may be considered to reflect the ordinal structure of categorical variables such as midrank scores [12] and Nishisato scores [13]. The (i,j)th joint marginal cell frequency is then calculated based on the ordinal log-linear model (1) such that(17)ñij=exp⁡[μ+αi+βj+ϕ(ui-u̅)(vj-v̅)],whereñij is the frequency observed in the (i,j)th cell, given the randomly generated marginal proportions, the pre-determined ϕ value, and(18)μ=ln⁡(n)+1I∑i=1Iln⁡(p̆i•)+1J∑j=1Jln⁡(p̆•j),αi=ln⁡(p̆i•)-1I∑i=1Iln⁡(p̆i•),βj=ln⁡(p̆•j)-1J∑j=1Jln⁡(p̆•j).It is apparent that theñij values obtained will not be guaranteed to be integers. Therefore, each cell value is rounded to the nearest integer to obtain the randomly generated contingency table. Thus, the sample size of n specified in the algorithm is only an approximate sample size value.


## 4.2. Results

The computational study performed in this paper is based on random contingency tables (using the procedure described in Section4.1) of different dimensions with a sample size of approximately 1000. Tables were generated with the following dimensions: 2 × 2, 2 × 3, 2 × 4, 2 × 5, 3 × 3, 3 × 4, 3 × 5, 4 × 4, 4 × 5, and 5 × 5. For each dimension, 200 random tables were generated each with a fixed value of ϕ; the ϕ values considered in this paper were 0, 0.01, 0.05, and 0.1 to 1.5 at increments of 0.1. In terms of the exponent of ϕ (i.e., the common log-odds ratio of the table), these are equivalent to 1, 1.01, 1.05, and 1.11 to 4.48, respectively. In practice, the common log-odds ratio may be negative; however, for the purposes of exploring the behaviour of Newton’s unidimensional method and the four non-iterative procedures (the LogNI, LogNI1, LogNI2, and BDNI estimates) against the true value of ϕ, we restricted our attention to positive ϕ values. A preliminary investigation revealed that randomly generated contingency tables with a negative ϕ behaved in a similar manner to tables generated with a positive ϕ value.Tables1, 2, 3, 4, and 5 provide a summary of the mean estimate and standard error (in parentheses) of ϕ from the 200 randomly generated contingency tables. Table 1 summarises the results obtained using Newton’s method, and Table 2 provides a summary of the accuracy of ϕ̂LogNI (7) for estimating ϕ. Table 3 summarises the results obtained by considering ϕ̂LogNI1 (10) to estimate the parameter while Tables 4 and 5 summarise output when using ϕ̂LogNI2 (7) and ϕ̂BDNI (12), respectively. We also took into consideration the number of iterations needed for Newton’s method to converge to four decimal places. Thus, Table 6 summarises the median number of iterations needed for convergence to occur. Also summarised in this table is the range (maximum and minimum) of the number of iterations observed for the each set of 200 simulated tables.Table 1
Mean estimate (and standard error) ofϕ using Newton’s unidimensional method, ϕ̂N, based on 200 simulated contingency tables with a specified value of ϕ and n=1000.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)0.00020 (0.1160)0.00235 (0.0829)0.00042 (0.0583)−0.00023 (0.0611)−0.00029 (0.0488)−0.00007 (0.0346)0.00008 (0.0312)0.00004 (0.0239)0.00005 (0.0178)0.010.01617 (0.2765)0.00740 (0.1237)0.01503 (0.0792)0.00868 (0.0631)0.01040 (0.0631)0.01016 (0.0446)0.01012 (0.0344)0.01020 (0.0307)0.01002 (0.0241)0.01000 (0.0179)0.050.05454 (0.2188)0.05134 (0.1194)0.05046 (0.0824)0.04911 (0.0621)0.04964 (0.0615)0.05009 (0.0424)0.05021 (0.0338)0.04997 (0.0301)0.05001 (0.0236)0.05009 (0.0183)0.100.10278 (0.2165)0.09998 (0.1089)0.10087 (0.0775)0.10061 (0.0608)0.09994 (0.0690)0.09990 (0.0436)0.09980 (0.0341)0.10006 (0.0301)0.10024 (0.0234)0.10001 (0.0181)0.200.20473 (0.2048)0.20097 (0.1149)0.19999 (0.0770)0.20010 (0.0580)0.19969 (0.0605)0.19970 (0.0431)0.20021 (0.0317)0.19995 (0.0290)0.20008 (0.0223)0.19995 (0.0165)0.300.30392 (0.2016)0.30622 (0.1130)0.30312 (0.0770)0.30040 (0.0558)0.30019 (0.0616)0.30022 (0.0404)0.30008 (0.0314)0.30007 (0.0278)0.29995 (0.0202)0.30001 (0.0149)0.400.40491 (0.2295)0.39656 (0.1170)0.40184 (0.0771)0.40082 (0.0629)0.40079 (0.0649)0.40080 (0.0412)0.40030 (0.0295)0.40001 (0.0263)0.40015 (0.0187)0.40036 (0.0129)0.500.50762 (0.2353)0.50533 (0.1164)0.50139 (0.0869)0.50326 (0.0564)0.49981 (0.0591)0.49788 (0.0422)0.50018 (0.0284)0.50030 (0.0249)0.50022 (0.0167)0.50025 (0.0113)0.600.61482 (0.2229)0.60065 (0.1102)0.60149 (0.0751)0.60076 (0.0561)0.59980 (0.0577)0.59945 (0.0373)0.59976 (0.0272)0.60012 (0.0230)0.59999 (0.0156)0.60015 (0.0096)0.700.71150 (0.2131)0.70272 (0.1175)0.69969 (0.0690)0.70033 (0.0524)0.69953 (0.0547)0.70031 (0.0359)0.70049 (0.0250)0.70026 (0.0209)0.70015 (0.0139)0.69996 (0.0078)0.800.80482 (0.2773)0.80215 (0.1081)0.79924 (0.0785)0.79978 (0.0485)0.79985 (0.0524)0.80013 (0.0345)0.80002 (0.0226)0.80009 (0.0197)0.79995 (0.0117)0.80038 (0.0066)0.900.89570 (0.2133)0.89081 (0.1033)0.90338 (0.0687)0.90103 (0.0513)0.89998 (0.0519)0.90018 (0.0314)0.90058 (0.0214)0.89980 (0.0174)0.90041 (0.0102)0.90025 (0.0055)1.001.01856 (0.2358)0.99917 (0.1074)0.99251 (0.0804)1.00302 (0.0504)0.99954 (0.0517)0.99992 (0.0301)1.00014 (0.0188)0.99983 (0.0167)0.99987 (0.0091)0.99659 (0.0045)1.101.10279 (0.2169)1.10185 (0.1187)1.10167 (0.0646)1.10041 (0.0472)1.10064 (0.0475)1.10020 (0.0282)1.09998 (0.0174)1.10030 (0.0146)1.10028 (0.0076)1.09676 (0.0035)1.201.20637 (0.2047)1.20018 (0.1057)1.20161 (0.0626)1.20078 (0.0430)1.20001 (0.0444)1.20077 (0.0260)1.19968 (0.0160)1.19960 (0.0132)1.19975 (0.0066)1.19886 (0.0030)1.301.32157 (0.2372)1.30192 (0.1100)1.30332 (0.0551)1.30048 (0.0396)1.30092 (0.0452)1.30039 (0.0241)1.29982 (0.0143)1.30007 (0.0121)1.30016 (0.0057)1.28541 (0.0023)1.401.41818 (0.2433)1.40873 (0.0898)1.40110 (0.0564)1.39902 (0.0395)1.40079 (0.0423)1.40040 (0.0224)1.40003 (0.0132)1.40001 (0.0106)1.39921 (0.0051)1.38124 (0.0021)1.501.49777 (0.1907)1.50079 (0.0939)1.50105 (0.0531)1.50184 (0.0384)1.49974 (0.0427)1.49900 (0.0223)1.50036 (0.0118)1.50045 (0.0091)1.49371 (0.0045)1.47004 (0.0017)Table 2
Mean estimate (and standard error) ofϕ using ϕ̂LogNI based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00026 (0.1161)0.00350 (0.0828)0.00019 (0.0583)−0.00023 (0.0611)−0.00030 (0.0448)−0.00003 (0.0346)0.00002 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01618 (0.2760)0.00562 (0.1233)0.01509 (0.0783)0.00810 (0.0629)0.01028 (0.0631)0.01024 (0.0446)0.01012 (0.0344)0.01021 (0.0307)0.01005 (0.0241)0.01001 (0.0179)0.050.05454 (0.2187)0.05164 (0.1193)0.05092 (0.0824)0.04926 (0.0620)0.04962 (0.0615)0.05012 (0.0423)0.05022 (0.0338)0.05003 (0.0300)0.05000 (0.0236)0.05016 (0.0182)0.100.10278 (0.2163)0.10001 (0.1088)0.10092 (0.0773)0.10061 (0.0607)0.10010 (0.0689)0.09997 (0.0434)0.09989 (0.0339)0.10008 (0.0299)0.10034 (0.0231)0.10018 (0.0177)0.200.20474 (0.2045)0.20054 (0.1145)0.20012 (0.0765)0.20043 (0.0574)0.19977 (0.0600)0.19992 (0.0424)0.20054 (0.0308)0.19996 (0.0280)0.20032 (0.0212)0.20024 (0.0151)0.300.30392 (0.2007)0.30788 (0.1103)0.30350 (0.0750)0.30083 (0.0545)0.30022 (0.0604)0.30042 (0.0389)0.30045 (0.0296)0.30035 (0.0259)0.30004 (0.0179)0.30088 (0.0124)0.400.40491 (0.2281)0.39618 (0.1155)0.40358 (0.0748)0.40355 (0.0598)0.40126 (0.0629)0.40262 (0.0386)0.40080 (0.0265)0.40078 (0.0232)0.40167 (0.0151)0.40312 (0.0092)0.500.50762 (0.2308)0.50795 (0.1130)0.49748 (0.0845)0.50464 (0.0518)0.50083 (0.0561)0.49881 (0.0385)0.50105 (0.0240)0.50318 (0.0205)0.50290 (0.0122)0.50552 (0.0069)0.600.61482 (0.2188)0.60143 (0.1069)0.60722 (0.0700)0.60415 (0.0506)0.59982 (0.0535)0.59997 (0.0324)0.60070 (0.0218)0.60368 (0.0175)0.60521 (0.0102)0.61267 (0.0048)0.700.71150 (0.2063)0.70642 (0.1128)0.70040 (0.0639)0.70262 (0.0463)0.70007 (0.0493)0.70457 (0.0298)0.70614 (0.0185)0.70520 (0.0145)0.70790 (0.0079)0.71631 (0.0029)0.800.80483 (0.2665)0.80608 (0.1018)0.80105 (0.0717)0.79954 (0.0419)0.74372 (0.0469)0.80350 (0.0275)0.80484 (0.0154)0.80512 (0.0124)0.81298 (0.0055)0.82603 (0.0018)0.900.89570 (0.2076)0.90016 (0.0967)0.90264 (0.0605)0.90583 (0.0421)0.90236 (0.0442)0.90229 (0.0234)0.91092 (0.0135)0.90512 (0.0100)0.91617 (0.0039)0.93313 (0.0011)1.001.01857 (0.2247)0.99812 (0.0995)0.99674 (0.0698)1.01016 (0.0394)1.00154 (0.0429)1.00645 (0.0209)1.01148 (0.0105)1.01483 (0.0085)1.02763 (0.0028)1.03001 (0.0006)1.101.10280 (0.2085)1.09827 (0.1085)1.10866 (0.0533)1.10855 (0.0361)1.10441 (0.0376)1.11082 (0.0186)1.11704 (0.0088)1.12472 (0.0063)1.13542 (0.0018)1.12646 (0.0003)1.201.20638 (0.1958)1.20485 (0.0944)1.21104 (0.0507)1.20885 (0.0319)1.20120 (0.0338)1.21299 (0.0157)1.21861 (0.0070)1.22115 (0.0049)1.23445 (0.0012)1.23221 (0.0002)1.301.32158 (0.2227)1.30226 (0.0970)1.31646 (0.0494)1.30883 (0.0286)1.30870 (0.0334)1.31479 (0.0137)1.32367 (0.0053)1.32639 (0.0039)1.33670 (0.0008)1.31161 (0.0001)1.401.41819 (0.2224)1.40275 (0.0774)1.40607 (0.0435)1.41200 (0.0270)1.41515 (0.0299)1.41974 (0.0115)1.43487 (0.0043)1.43576 (0.0028)1.42865 (0.0006)1.36380 (0.0008)1.501.49778 (0.1785)1.50663 (0.0795)1.50531 (0.0402)1.51515 (0.0250)1.50567 (0.0290)1.50459 (0.0108)1.53475 (0.0032)1.53676 (0.0018)1.52404 (0.0004)1.43543 (0.0004)Table 3
Mean estimate (and standard error) ofϕ using ϕ̂LogNI1 based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00008 (0.1160)0.00299 (0.0829)0.00032 (0.0583)−0.00024 (0.0611)−0.00030 (0.0448)−0.00004 (0.0346)0.00005 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01610 (0.2759)0.00633 (0.1234)0.01424 (0.0785)0.00836 (0.0630)0.01029 (0.0631)0.01020 (0.0446)0.01011 (0.0344)0.01021 (0.0307)0.01004 (0.0241)0.01001 (0.0179)0.050.05446 (0.2187)0.05152 (0.1193)0.05075 (0.0824)0.04922 (0.0620)0.04962 (0.0615)0.05008 (0.0423)0.05016 (0.0338)0.04996 (0.0300)0.04992 (0.0236)0.05000 (0.0182)0.100.10231 (0.2154)0.09998 (0.1099)0.10078 (0.0773)0.10046 (0.0607)0.09993 (0.0689)0.09975 (0.0434)0.09954 (0.0339)0.09969 (0.0299)0.09964 (0.0231)0.09911 (0.0177)0.200.20451 (0.2045)0.20025 (0.1145)0.19952 (0.0765)0.19932 (0.0574)0.19893 (0.0600)0.19818 (0.0424)0.19766 (0.0308)0.19676 (0.0281)0.19492 (0.0213)0.19101 (0.0153)0.300.30337 (0.2007)0.30434 (0.1110)0.29977 (0.0753)0.29721 (0.0546)0.29743 (0.0604)0.29481 (0.0390)0.29092 (0.0297)0.28925 (0.0261)0.28149 (0.0182)0.26891 (0.0129)0.400.40352 (0.2275)0.39385 (0.1155)0.39791 (0.0749)0.39275 (0.0602)0.39458 (0.0629)0.38791 (0.0388)0.37732 (0.0269)0.37233 (0.0236)0.35294 (0.0159)0.32652 (0.0103)0.500.50350 (0.2310)0.49954 (0.1134)0.48823 (0.0846)0.48801 (0.0526)0.48693 (0.0563)0.47029 (0.0389)0.45542 (0.0248)0.44396 (0.0214)0.40896 (0.0135)0.36225 (0.0086)0.600.60742 (0.2165)0.59252 (0.1070)0.58724 (0.0709)0.57402 (0.0516)0.57580 (0.0539)0.55214 (0.0331)0.52330 (0.0231)0.50859 (0.0189)0.45332 (0.0121)0.38300 (0.0070)0.700.70298 (0.2073)0.69001 (0.1127)0.67442 (0.0644)0.65328 (0.0473)0.66050 (0.0499)0.62525 (0.0311)0.57783 (0.0204)0.55366 (0.0166)0.48359(0.0104)0.38935 (0.0054)0.800.79400 (0.2685)0.78320 (0.1027)0.75942 (0.0726)0.73482 (0.0431)0.74372 (0.0469)0.68976 (0.0291)0.62194 (0.0178)0.59053 (0.0151)0.49281 (0.0085)0.39437 (0.0043)0.900.88478 (0.2076)0.86906 (0.0972)0.84284 (0.0619)0.80984 (0.0446)0.81907 (0.0456)0.74370 (0.0256)0.65549 (0.0165)0.61944 (0.0130)0.50525 (0.0071)0.39864(0.0034)1.000.99208 (0.2242)0.95982 (0.1003)0.91896 (0.0711)0.87180 (0.0428)0.88914 (0.0447)0.78265 (0.0240)0.68186 (0.0140)0.64302 (0.0122)0.51164 (0.0061)0.39090(0.0026)1.101.07792 (0.2070)1.04284 (0.1102)0.99425 (0.0563)0.93162 (0.0397)0.94721 (0.0399)0.82129 (0.0221)0.69747 (0.0127)0.65111 (0.0103)0.50923 (0.0048)0.38387 (0.0019)1.201.17866 (0.1918)1.12713 (0.0961)1.06491 (0.0539)0.98924 (0.0355)1.00068 (0.0366)0.85517 (0.0198)0.71206 (0.0114)0.66053 (0.0091)0.51300 (0.0040)0.38100 (0.0016)1.301.27633 (0.2260)1.20889 (0.0981)1.12928 (0.0542)1.05000 (0.0322)1.05912 (0.0370)0.88164 (0.0181)0.71797 (0.0099)0.67223 (0.0081)0.51028 (0.0033)0.37375 (0.011)1.401.37037 (0.2244)1.29294 (0.0806)1.19599 (0.0470)1.07872 (0.0317)1.09808 (0.0340)0.90062 (0.0146)0.73551 (0.0089)0.66730 (0.0068)0.50864 (0.0029)0.37198 (0.0011)1.501.44834 (0.1785)1.37143 (0.0823)1.24810 (0.0439)1.12323 (0.0305)1.14909 (0.0337)0.92568 (0.0161)0.73350 (0.0077)0.66406 (0.0056)0.50923 (0.0025)0.36991 (0.0009)Table 4
Mean estimate (and standard error) ofϕ using ϕ̂LogNI2 based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00068 (0.2014)−0.00020 (0.1161)0.00331 (0.0828)0.00026 (0.0583)−0.00023 (0.0611)−0.00030 (0.0448)−0.00003 (0.0346)0.00003 (0.0312)0.00006 (0.0239)0.00005 (0.0178)0.010.01617 (0.2760)0.00587 (0.1233)0.01482 (0.0783)0.00818 (0.0629)0.01028 (0.0631)0.01022 (0.0446)0.01012 (0.0344)0.01021 (0.0307)0.01005 (0.0241)0.01001 (0.0179)0.050.05454 (0.2187)0.05161 (0.1193)0.05075 (0.0824)0.04925 (0.0620)0.04962 (0.0615)0.05011 (0.0423)0.05021 (0.0338)0.05002 (0.0300)0.05000 (0.0236)0.05014 (0.0182)0.100.10271 (0.2162)0.10001 (0.1088)0.10089 (0.0773)0.10061 (0.0607)0.10007 (0.0689)0.09996 (0.0434)0.09987 (0.0339)0.10007 (0.0299)0.10029 (0.0231)0.10012 (0.0177)0.200.20473 (0.2045)0.20054 (0.1145)0.20012 (0.0765)0.20036 (0.0574)0.19974 (0.0600)0.19982 (0.0424)0.20031 (0.0308)0.19976 (0.0280)0.19977 (0.0212)0.19881 (0.0152)0.300.30391 (0.2007)0.30672 (0.1106)0.30252 (0.0752)0.30051 (0.0545)0.30013 (0.0604)0.30002 (0.0389)0.29933 (0.0296)0.29895 (0.0259)0.29652 (0.0179)0.29131 (0.0125)0.400.40486 (0.2281)0.39614 (0.1155)0.40305 (0.0748)0.40154 (0.0599)0.40072 (0.0629)0.40021 (0.0386)0.39647 (0.0266)0.39431 (0.0233)0.38530 (0.0154)0.36861 (0.0097)0.500.50682 (0.2309)0.50599 (0.1131)0.49682 (0.0845)0.50235 (0.0521)0.49920 (0.0561)0.49339 (0.0386)0.48907 (0.0242)0.48343 (0.0208)0.46052 (0.0128)0.42230 (0.0079)0.600.61277 (0.2177)0.60066 (0.1069)0.60309 (0.0704)0.59827 (0.0509)0.59648 (0.0540)0.58824 (0.0326)0.57337 (0.0223)0.56490 (0.0181)0.52303 (0.0112)0.45747 (0.0062)0.700.70991 (0.2067)0.70403 (0.1126)0.69708 (0.0640)0.69114 (0.0466)0.69288 (0.0494)0.67813 (0.0302)0.64850 (0.0193)0.63058 (0.0155)0.56907 (0.0094)0.47393 (0.0046)0.800.80340 (0.2670)0.80302 (0.1021)0.79317 (0.0719)0.78409 (0.0422)0.78794 (0.0462)0.76014 (0.0281)0.71204 (0.0166)0.68449 (0.0139)0.59183 (0.0075)0.48567 (0.0036)0.900.89526 (0.2076)0.89665 (0.0968)0.88885 (0.0610)0.87568 (0.0431)0.87881 (0.0447)0.83495 (0.0243)0.76251 (0.0152)0.73009 (0.0118)0.61468 (0.0061)0.49547 (0.0027)1.001.01160 (0.2247)0.99328 (0.0996)0.97792 (0.0701)0.95673 (0.0410)0.96469 (0.0435)0.89309 (0.0225)0.80583 (0.0126)0.76616 (0.0108)0.62836 (0.0051)0.48938 (0.0020)1.101.09777 (0.2076)1.08835 (0.1090)1.07277 (0.0546)1.03428 (0.0379)1.04616 (0.0384)0.94742 (0.0205)0.83311 (0.0113)0.78539 (0.0090)0.63087 (0.0040)0.48327 (0.0015)1.201.20275 (0.1938)1.18780 (0.0948)1.16055 (0.0520)1.10968 (0.0337)1.11730 (0.0350)1.00053 (0.0181)0.85947 (0.0099)0.80448 (0.0078)0.63668 (0.0032)0.48174 (0.0012)1.301.31224 (0.2236)1.28023 (0.0971)1.24276 (0.0517)1.18733 (0.0305)1.19350 (0.0351)1.04284 (0.0164)0.87401 (0.0085)0.82342 (0.0068)0.63865 (0.0026)0.47549 (0.0008)1.401.41026 (0.2230)1.37640 (0.0786)1.32384 (0.0451)1.23586 (0.0297)1.25046 (0.0320)1.07434 (0.0146)0.90116 (0.0075)0.82417 (0.0057)0.63906 (0.0023)0.47335 (0.0008)1.501.49244 (0.1784)1.46870 (0.0806)1.39468 (0.0420)1.29823 (0.0282)1.31951 (0.0314)1.10946 (0.0142)0.90476 (0.0064)0.82555 (0.0045)0.64242 (0.0019)0.47224 (0.0006)Table 5
Mean estimate (and standard error) ofϕ using ϕ̂BDNI based on 200 simulated contingency tables with a specified value of ϕ and n≈1000. The values in bold are where the mean of the 200 P values obtained from the Wald test (of the difference between the estimated and true parameter) is significant at the 1% level. The italicised values are significant at the 5% level.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.00−0.00064 (0.2014)0.00017 (0.0060)0.00229 (0.0829)0.00042 (0.0583)−0.00023 (0.0611)−0.00029 (0.0448)−0.00007 (0.0346)0.00008 (0.0312)0.00004 (0.0239)0.00005 (0.0178)0.010.01578 (0.2745)0.00740 (0.1236)0.01351 (0.0787)0.00868 (0.0630)0.01038 (0.0631)0.01016 (0.0446)0.01012 (0.0344)0.01020 (0.0307)0.01002 (0.0241)0.01000 (0.0179)0.050.05537 (0.2182)0.05137 (0.1193)0.05045 (0.0823)0.04906 (0.0620)0.04959 (0.0615)0.05004 (0.0423)0.05008 (0.0338)0.04987 (0.0300)0.04982 (0.0236)0.04983 (0.0182)0.100.10060 (0.2125)0.09994 (0.1087)0.10065 (0.0773)0.10031 (0.0606)0.09990 (0.0688)0.09950 (0.0434)0.09900 (0.0339)0.09925 (0.0299)0.09871 (0.0231)0.09771 (0.0177)0.200.20554 (0.2037)0.20065 (0.1143)0.19853 (0.0765)0.19771 (0.0574)0.19736 (0.0600)0.19566 (0.0424)0.19334 (0.0308)0.19195 (0.0280)0.18773 (0.0214)0.18065 (0.0154)0.300.30442 (0.1998)0.30099 (0.1113)0.29682 (0.0754)0.29088 (0.0546)0.29373 (0.0604)0.28762 (0.0390)0.27933 (0.0299)0.27595 (0.0262)0.26320 (0.0184)0.24319 (0.0132)0.400.39883 (0.2243)0.38977 (0.1152)0.38999 (0.0749)0.38243 (0.0602)0.38611 (0.0628)0.37154 (0.0389)0.35374 (0.0272)0.34521 (0.0240)0.31750 (0.0164)0.28401 (0.0109)0.500.49941 (0.2300)0.48981 (0.1128)0.48156 (0.0833)0.46832 (0.0528)0.46943 (0.0564)0.44038 (0.0393)0.41736 (0.0253)0.39981 (0.0220)0.35530 (0.0142)0.30366 (0.0093)0.600.60336 (0.2113)0.58180 (0.1068)0.56500 (0.0712)0.54913 (0.0516)0.54854 (0.0541)0.50988 (0.0337)0.46819 (0.0238)0.44776 (0.0198)0.38563 (0.0130)0.31530 (0.0077)0.700.69972 (0.2058)0.66768 (0.1121)0.64706 (0.0645)0.60455 (0.0479)0.62207 (0.0503)0.56595 (0.0319)0.50408 (0.0214)0.47675 (0.0176)0.40331 (0.0113)0.31913 (0.0061)0.800.79011 (0.2568)0.75919 (0.1026)0.71197 (0.0728)0.67057 (0.0437)0.69012 (0.0476)0.61207 (0.0302)0.53082 (0.0190)0.50154 (0.0162)0.40365 (0.0094)0.32837 (0.0049)0.900.86600 (0.2049)0.83271 (0.0969)0.78390 (0.0627)0.73735 (0.0455)0.75035 (0.0464)0.64749 (0.0269)0.55024 (0.0178)0.51415 (0.0142)0.41282 (0.0079)0.34591 (0.0038)1.000.95778 (0.2192)0.92254 (0.0996)0.84033 (0.0739)0.77967 (0.0437)0.80391 (0.0458)0.67170 (0.0254)0.56617 (0.0153)0.53203 (0.0134)0.42429 (0.0069)0.36422 (0.0029)1.101.05198 (0.2011)0.99021 (0.1101)0.90467 (0.0571)0.81670 (0.0410)0.83903 (0.0413)0.69141 (0.0236)0.57224 (0.0140)0.53513 (0.0115)0.43699 (0.0055)0.38981 (0.0021)1.201.14028 (0.1859)1.05917 (0.0963)0.95848 (0.0551)0.85077 (0.0370)0.86977 (0.0383)0.71457 (0.0213)0.58876 (0.0126)0.54391 (0.0102)0.45929 (0.0045)0.43908 (0.0017)1.301.24366 (0.2000)1.12319 (0.0973)0.99862 (0.0553)0.88461 (0.0340)0.91884 (0.0386)0.72751 (0.0197)0.59232 (0.0111)0.55751 (0.0091)0.48817 (0.0037)0.60876 (0.0011)1.401.34106 (0.2180)1.19429 (0.0788)1.04569 (0.0489)0.90048 (0.0337)0.93586 (0.0359)0.74161 (0.0180)0.61512 (0.0100)0.58435 (0.0076)0.51219 (0.0032)0.65896 (0.0010)1.501.37638 (0.1752)1.25152 (0.0836)1.06724(0.0460)0.92840 (0.0324)0.98432 (0.0356)0.76259 (0.0177)0.63542 (0.0086)0.59487 (0.0063)0.59138 (0.0027)0.65896 (0.0007)Table 6
Median number of iterations needed for Newton’s unidimensional method to converge to four decimal places. In parenthesis is the range of iterations for the 200 simulated contingency tables.
Size of contingency tableTrueϕ2 × 22 × 32 × 42 × 53 × 33 × 43 × 54 × 44 × 55 × 50.002222222222(1–7)(1–6)(1–9)(1–4)(1–3)(1–4)(1–3)(2-3)(1–3)(1–3)0.013333333333(2–14)(2–30)(2–95)(2–14)(2–5)(2-3)(2-3)(2-3)(2-3)(2-3)0.054444544455(3–13)(1–7)(2–6)(4–6)(4–6)(4-5)(4-5)(4-5)(4-5)(4–6)0.105555566667(3–29)(4–7)(5–8)(5–8)(4–6)(5–7)(5–7)(5–7)(6–8)(6–10)0.207778889101113(5–13)(7–12)(7–11)(7–11)(7–10)(7–11)(7–15)(8–13)(9–17)(8–21)0.3099910101213141723(3–23)(9–219)(9–128)(9–15)(9–15)(9–17)(10–22)(10–19)(12–27)(15–35)0.4011111213131619212839(7–25)(6–16)(11–25)(11–23)(11–22)(12–26)(13–33)(14–47)(16–49)(20–65)0.5013141516.5172127304468(4–261)(12–56)(13–28)(14–103)(13–30)(6–47)(18–53)(18–57)(23–80)(24–110)0.60151618212128374266113(11–260)(15–32)(16–33)(17–66)(16–47)(18–61)(21–81)(23–98)(30–139)(41–261)0.701719232827375364101211(14–737)(17–39)(19–52)(21–67)(19–43)(20–71)(26–150)(30–126)(25–195)(59–560)0.802024283434497589170359(8–407)(20–40)(23–58)(25–64)(22–49)(24–140)(32–160)(38–188)(61–352)(102–1026)0.90242934424367107.5128.5255605(9–39)(24–46)(27–716)(29–114)(26–101)(32–116)(39–207)(58–332)(62–653)(147–2129)1.00283242555390145.51723951028(16–987)(27–83)(17–101)(37–286)(31–122)(42–281)(67–349)(66–480)(99–1173)(3–3243)1.1032395268.569116211275656.51875.5(29–863)(4–1983)(38–187)(48–285)(39–181)(45–383)(57–580)(68–590)(151–2370)(8–6656)1.203647648888157.52853929992725(28–208)(36–94)(44–180)(51–201)(46–210)(69–579)(94–889)(129–907)(199–3120)(589–12917)1.30425577109.5109207420.55561571.54730(29–1239)(33–267)(53–275)(58–265)(52–222)(69–441)(135–1104)(154–1606)(260–4833)(11–23803)1.40486494143.5135.5284562822.522496224(23–2288)(45–2781)(58–262)(78–372)(63–431)(88–557)(129–1675)(187–3370)(24–7041)(3–25352)1.505675114177170369.582712333102.59769(38–92)(57–216)(74–302)(104–914)(72–624)(103–1071)(189–2641)(305–4608)(22–11773)(9–34759)Table1 and Figure 1 show that Newton’s method is very accurate for estimating ϕ. The Wald statistic of the mean parameter estimate indicates that such a procedure gives estimates of ϕ that are statistically indistinguishable from the true estimate for all sized contingency tables and ϕ values. However, the number of iterations needed for Newton’s method to estimate the parameter varies greatly depending on the size of the contingency table and the magnitude of ϕ (Table 6). Contingency tables of small dimension and small ϕ required no more than about 10 iterations for convergence to four decimal places to occur. Due to the accuracy of Newton’s method, the legends attached to Figure 1 highlighting the varying table sizes have been omitted.Figure 1
Graphical view of the performance of Newton’s unidimensional method for estimatingϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Table6 shows that for contingency tables of all dimensions considered, if ϕ is larger than 0.8, then the number of iterations for convergence was nearly always larger than 10 with a median number of iterations ranging from 20 (for a 2 × 2 contingency table) to 9769 (for a 5 × 5 table). For the 200 randomly generated 5 × 5 contingency tables, the number of iterations until convergence was achieved ranging from an economical 9 iterations to nearly 35000.Table2 and Figure 2 show that of the four non-iterative procedures, using ϕ̂LogNI to estimate ϕ came closest to matching the performance of Newton’s method. There were 15 cases out of 180 total cases where the mean estimate of ϕ obtained using ϕ̂LogNI was different from the true ϕ at the 1% level of significance. Ten of these 15 cases occurred when estimating a ϕ>1.00 in a contingency table having a dimension of 4 × 5 or 5 × 5. Due to the accuracy of the ϕ̂LogNI estimate, the legends attached to Figure 2 highlighting the varying table sizes have been omitted.Figure 2
Graphical view of the performance ofϕ̂LogNI (7) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Tables3, 4, and 5 and the corresponding figures show that the estimators ϕ̂LogNI1, ϕ̂LogNI2, and ϕ̂BDNI are stable and accurate for relatively small contingency tables and where ϕ is also small (not exceeding approximately 0.5). However they performed poorly for relatively large tables with a large parameter value, continually underestimating the parameter. Figures 3, 4 and 5 graphically verify these results. Using ϕ̂LogNI2 resulted in 30 cases (out of 180 total cases) where the mean estimate of ϕ was different from the true ϕ value (at the 1% level of significance) with 23 of these 30 cases occurring when estimating a ϕ≥1.00 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5, or 5 × 5. Using ϕ̂LogNI1 resulted in 35 cases where the mean estimate of ϕ was statistically different from the true ϕ value with 27 of these 35 cases occurring when estimating a ϕ≥0.90 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5 or 5 × 5. Using ϕ̂BDNI resulted in 43 cases where the mean estimate of ϕ was statistically different from the true ϕ value with 31 of these 43 cases occurring when estimating a ϕ≥0.80 in a contingency table having a dimension of 3 × 5, 4 × 4, 4 × 5 or 5 × 5.Figure 3
Graphical view of the performance ofϕ̂LogNI1 (10) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Figure 4
Graphical view of the performance ofϕ̂LogNI2 (7) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Figure 5
Graphical view of the performance ofϕ̂BDNI (12) for estimating ϕ. Mean estimate of parameter versus true parameter for the different values of ϕ and contingency tables of varying dimension.Two hundred (200) contingency tables of dimension 4 × 5 were randomly generated (using the algorithm described in Section4.1) with ϕ=0.4 and sample size approximately 1000. For these tables, Figure 6 illustrates the stability of each of the five estimation procedures and summarises the P value obtained from the Wald test of the estimates ϕ̂N, ϕ̂LogNI, ϕ̂LogNI1, ϕ̂LogN2, and ϕ̂BDNI when compared with ϕ=0.4. All five methods provide reliable and accurate estimates of the parameter and, for all the tables generated, are statistically indistinguishable when compared with the true value of the parameter. These figures show that Newton’s method (Figure 6(a)) provides the most reliable estimate since the P value calculated from each of the 200 tables is consistently near 1.0. Of the non-iterative procedures Figure 6(b) shows that ϕ̂LogNI is nearly as stable (although slightly more variable) when compared with Newton’s method.TheP-values obtained from the Wald Test comparing the true ϕ value against (a) ϕ̂N (b) ϕ̂LogNI, (c) ϕ̂LogNI2, (d) ϕ̂LogNI1 and (e) ϕ̂BDNI. These P-values are based on the estimates obtained from randomly generating 200 contingency tables of dimension 4 × 5 where ϕ = 0.4 and n≈1000.
(a)(b)(c)(d)(e)The three remaining non-iterative procedures also behave well, but provide poorer estimates ofϕ than Newton’s method or ϕ̂LogNI. Figure 6(c) shows that ϕ̂LogNI2 provides estimates that are closer to ϕ than ϕ̂LogNI1 Figure 6(d) shows that the P value of ϕ̂LogNI1 is centred around 0.7 while Figure 6(c) shows that the P value of ϕ̂LogNI2 is centred around 0.9. This is to be expected since ϕ̂LogNI1 is a first-order approximation of the LogNI estimate (7) while ϕ̂LogNI2 is its second-order approximation. For the 200 4 × 5 tables, these two non-iterative estimation procedures provide more accurate estimates of ϕ than ϕ̂BDNI. Figure 6(e) shows that the P value from the Wald test of the estimate is centred around 0.5.


## 5. Discussion

Iterative procedures can be computationally very intensive. Certainly for those contingency tables that are of large dimension, or whereϕ is large (but unknown), the number of iterations needed for convergence to occur can be excessively large. This is not necessarily a major issue for those who are experienced in statistical programming. However, for those researchers that are not, a simple direct estimation procedure can prove to be advantageous. This paper has demonstrated the benefit of considering alternative, iteration-free, procedures for estimating the linear-by-linear parameter of an ordinal log-linear, and its related association model.For small contingency tables (in terms of dimension andϕ) any of the non-iterative procedures considered in this paper will produce estimates that are statistically similar to their true value. For large contingency tables with moderately small values of ϕ (less than 0.4), any of the non-iterative procedures will work extremely well without the need for programming and executing iterative procedures.Of the four non-iterative procedures considered in this paper,ϕ̂LogNI (7) provides more accurate and reliable results than its counterparts. In fact there were only a few cases where the estimate produced was statistically different (according to the Wald test) from the true parameter value. However, an inspection of Table 2 reveals that these differences appear negligible and that such a result can be attributed to the very small standard error of the non-iterative estimate. Certainly Figure 2 shows the accurate performance of ϕ̂LogNI.The enhanced accuracy ofϕ̂LogNI, when compared with the remaining non-iterative estimates, is not unexpected since it is easily obtained by rearranging (4) and aggregating across i and j. Since ϕ̂BDNI, ϕ̂LogNI1, and ϕ̂LogNI2 are calculated using functional approximations of ϕ̂LogNI, this latter procedure will always produce more accurate estimates. Such a conclusion is observed by considering the results in Section 4.One can also observe the relative accuracy of the non-iterative methods. Since the right-hand side of (11) is less than its left-hand side, ϕ̂BDNI>ϕ̂LogNI≈ϕ for all contingency tables. This is indeed the case when comparing the results summarised in Tables 2 and 5. This inequality becomes more severe for contingency tables of relatively large dimension. For similar reasons ϕ̂LogNI1<ϕ̂LogNI2, and this is clear when comparing Table 3 with Table 4. However, such a consideration is not a trivial exercise since it allows one to better understand how many terms (from the approximation of the logarithm function) will yield an accurate and reliable estimate of ϕ.Beh and Davy [8], who originally considered the accuracy of the non-iterative estimation procedures for ordinal log-linear models, showed from a practical perspective that ϕBDNI will produce accurate estimates of the true parameter. The more extensive study by Beh and Farver [9] confirmed these findings. However, the approximation (11) is valid only if -1<pij/(pi•pj)≤1. If only the positive values of this ratio are considered, ϕBDNI will provide reasonably good estimates of ϕ only when pij/(pi•pj)≤1 for all i and j. Therefore, such a result suggests that the BDNI procedure will produce better estimates of ϕ for relatively small contingency tables than for larger tables.The computations performed here have focused on generating contingency tables with a sample size of approximately 1000. Further investigations need to be undertaken to study the behaviour of the non-iterative estimation procedures for varying sample sizes. Preliminary investigations reveal that, for contingency tables with a sample size exceeding about 200, the estimateϕ̂LogNI is similar to the estimate obtained using Newton’s unidimensional method and the true parameter. Further work on the application of such non-iterative procedures needs investigating for other types of models of association for ordinal contingency tables.



---

*Source: 340415-2012-05-28.xml*