<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJS</journal-id><journal-title-group><journal-title>Open Journal of Statistics</journal-title></journal-title-group><issn pub-type="epub">2161-718X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojs.2015.56054</article-id><article-id pub-id-type="publisher-id">OJS-60338</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Decomposition of Independence Using the Logit Uniform Association Model and Equality of Concordance and Discordance for Two-Way Classifications
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ouji</surname><given-names>Tahata</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nobuko</surname><given-names>Miyamoto</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Sadao</surname><given-names>Tomizawa</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Information Sciences, Tokyo University of Science, Chiba, Japan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>kouji_tahata@is.noda.tus.ac.jp(OT)</email>;<email>miyamoto@is.noda.tus.ac.jp(NM)</email>;<email>tomizawa@is.noda.tus.ac.jp(ST)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>13</day><month>10</month><year>2015</year></pub-date><volume>05</volume><issue>06</issue><fpage>514</fpage><lpage>518</lpage><history><date date-type="received"><day>4</day>	<month>September</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>11</month>	<year>October</year>	</date><date date-type="accepted"><day>16</day>	<month>October</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   For two-way contingency tables with ordered categories, the present paper gives a theorem that the independence model holds if and only if the logit uniform association model holds and equality of concordance and discordance for all pairs of adjacent rows and all dichotomous collapsing of the columns holds. Using the theorem, we analyze the cross-classification of duodenal ulcer patients according to operation and dumping severity. 
 
</p></abstract><kwd-group><kwd>Concordance</kwd><kwd> Discordance</kwd><kwd> Independence</kwd><kwd> Logit Uniform Association Model</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Consider the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x5.png" xlink:type="simple"/></inline-formula> contingency tables with ordered categories, let X and Y denote the row and column variables, and let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x6.png" xlink:type="simple"/></inline-formula> (&gt;0) for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x7.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x8.png" xlink:type="simple"/></inline-formula>. Goodman [<xref ref-type="bibr" rid="scirp.60338-ref1">1</xref>] considered the uniform association (U) model which was defined by</p><disp-formula id="scirp.60338-formula380"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x9.png"  xlink:type="simple"/></disp-formula><p>See also Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 76). The U model may also be expressed as</p><disp-formula id="scirp.60338-formula381"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x10.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60338-formula382"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x11.png"  xlink:type="simple"/></disp-formula><p>Namely this model indicates the constant of the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x12.png" xlink:type="simple"/></inline-formula> local odds ratios <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x13.png" xlink:type="simple"/></inline-formula> defined for adjacent rows and adjacent columns. A special case of the U model obtained by putting <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x14.png" xlink:type="simple"/></inline-formula> is the independence (I) model.</p><p>If the I model holds, the correlation coefficient of X and Y equals zero; but the converse does not hold. We are interested in what structure between X and Y is necessary for obtaining the I model, in addition to the correlation coefficient being to zero.</p><p>Tomizawa, Miyamoto and Sakurai [<xref ref-type="bibr" rid="scirp.60338-ref3">3</xref>] give the theorem that the I model holds if and only if the Pearson’s correlation coefficient <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x15.png" xlink:type="simple"/></inline-formula> for X and Y equals zero and the U model holds.</p><p>Tomizawa et al. [<xref ref-type="bibr" rid="scirp.60338-ref3">3</xref>] also give the theorem that the I model holds if and only if the Kendall’s <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x16.png" xlink:type="simple"/></inline-formula> equals zero and the U model holds. For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x17.png" xlink:type="simple"/></inline-formula>, see Kendall [<xref ref-type="bibr" rid="scirp.60338-ref4">4</xref>] and Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 161).</p><p>Tahata, Miyamoto and Tomizawa [<xref ref-type="bibr" rid="scirp.60338-ref5">5</xref>] give the theorem that the I model holds if and only if the Spearman’s <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x18.png" xlink:type="simple"/></inline-formula> equals zero and the U model holds. For<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x19.png" xlink:type="simple"/></inline-formula>, see Stuart [<xref ref-type="bibr" rid="scirp.60338-ref6">6</xref>] , Kendall and Gibbons ([<xref ref-type="bibr" rid="scirp.60338-ref7">7</xref>] , p. 8), and Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 164). Also, Tahata and Tomizawa [<xref ref-type="bibr" rid="scirp.60338-ref8">8</xref>] review topics related to the quasi-uniform association model (Goodman [<xref ref-type="bibr" rid="scirp.60338-ref1">1</xref>] ), and the decomposition of symmetry into some models for the analysis of square contingency tables.</p><p>Suppose that the column variable Y is a response variable. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x20.png" xlink:type="simple"/></inline-formula> denote the jth cumulative logit within row i; i.e.,</p><disp-formula id="scirp.60338-formula383"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x21.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60338-formula384"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x22.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.60338-formula385"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x23.png"  xlink:type="simple"/></disp-formula><p>The logit uniform association (logit U) model (Agresti [<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 122) is defined by</p><disp-formula id="scirp.60338-formula386"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x24.png"  xlink:type="simple"/></disp-formula><p>namely</p><disp-formula id="scirp.60338-formula387"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x25.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.60338-formula388"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x26.png"  xlink:type="simple"/></disp-formula><p>Thus the logit U model indicates the constant of the odds ratios for the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x27.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x28.png" xlink:type="simple"/></inline-formula> tables obtained by taking all pairs of adjacent rows and all dichotomous collapsing of the response (Agresti [<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 122). A special case of the logit U model obtained by putting <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x29.png" xlink:type="simple"/></inline-formula> (i.e.,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x30.png" xlink:type="simple"/></inline-formula>) is the I model. We are now interested in what structure of probabilities <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x31.png" xlink:type="simple"/></inline-formula> is necessary for obtaining the I model, in addition to the logit U model (instead of the U model).</p><p>The purpose of the present paper is to give the decomposition of the I model by using the logit U model (in Section 2).</p></sec><sec id="s2"><title>2. Decomposition of Independence</title><p>Let</p><disp-formula id="scirp.60338-formula389"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x32.png"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.60338-formula390"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x33.png"  xlink:type="simple"/></disp-formula><p>For a randomly selected pair of observations, 1) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x34.png" xlink:type="simple"/></inline-formula>is the probability of concordance such that the</p><p>member that ranks in row <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x35.png" xlink:type="simple"/></inline-formula> rather than in row i also ranks in column <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x36.png" xlink:type="simple"/></inline-formula> or above rather than in column</p><p>j or below, and 2) <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x37.png" xlink:type="simple"/></inline-formula>is the probability of discordance such that the member that ranks in row <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x38.png" xlink:type="simple"/></inline-formula></p><p>rather than in row i ranks in column j or below rather than in column <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x39.png" xlink:type="simple"/></inline-formula> or above. Therefore <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x40.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x41.png" xlink:type="simple"/></inline-formula> indicate the sum of probabilities of such concordance and those of such discordance, respectively.</p><p>We shall consider the model of equality of concordance and discordance (say, CDE model) by</p><disp-formula id="scirp.60338-formula391"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x42.png"  xlink:type="simple"/></disp-formula><p>Then we obtain the following theorem.</p><p>Theorem 1. The I model holds if and only if both the CDE model and the logit U model hold.</p><p>Proof. If the I model holds, i.e., <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x43.png" xlink:type="simple"/></inline-formula>, then</p><disp-formula id="scirp.60338-formula392"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x44.png"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.60338-formula393"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x45.png"  xlink:type="simple"/></disp-formula><p>Thus, the CDE model holds. Also, if the I model holds, then the logit U model (with<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x46.png" xlink:type="simple"/></inline-formula>) holds.</p><p>Assuming that both the CDE model and the logit U model hold, then we shall show that the I model holds. Since the logit U model holds, we see</p><disp-formula id="scirp.60338-formula394"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x47.png"  xlink:type="simple"/></disp-formula><p>Thus</p><disp-formula id="scirp.60338-formula395"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x48.png"  xlink:type="simple"/></disp-formula><p>Since the CDE model holds, we obtain<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x49.png" xlink:type="simple"/></inline-formula>. The proof is completed. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x50.png" xlink:type="simple"/></inline-formula></p><p>Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x51.png" xlink:type="simple"/></inline-formula> denote the observed frequency in the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x52.png" xlink:type="simple"/></inline-formula> cell<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x53.png" xlink:type="simple"/></inline-formula>. Assume that a multinomial distribution applies to the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x54.png" xlink:type="simple"/></inline-formula> table. Let <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x55.png" xlink:type="simple"/></inline-formula> denote the likelihood ratio chi-squared statistic for testing goodness-of-fit of model M defined by</p><disp-formula id="scirp.60338-formula396"><graphic  xlink:href="http://html.scirp.org/file/5-1240565x56.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x57.png" xlink:type="simple"/></inline-formula> is the maximum likelihood estimate of expected frequency <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x58.png" xlink:type="simple"/></inline-formula> under the model M. The numbers of degrees of freedom (df) for testing the I, logit U, and CDE models are<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x59.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x60.png" xlink:type="simple"/></inline-formula>, and 1, respectively.</p></sec><sec id="s3"><title>3. An Example</title><p>The data in <xref ref-type="table" rid="table1">Table 1</xref> are taken directly from Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p. 12), which originally was presented by Grizzle, Starmer and Koch [<xref ref-type="bibr" rid="scirp.60338-ref9">9</xref>] . Four different operations for treating duodenal ulcer patients correspond to removal of various amounts of the stomach. Operation A is drainage and vagotomy, B is 25% resection (antrectomy) and vagotomy, C is 50% resection (hemigastrectomy) and vagotomy, and D is 75% resection. The categories of operation variable have a natural ordering. The dumping severity variable describes the extent of an undesirable potential consequence of the operation. The categories of this variable are also ordered. For these data, the I model fits well with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x61.png" xlink:type="simple"/></inline-formula> based on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x62.png" xlink:type="simple"/></inline-formula>. The logit U model also fits these data well with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x63.png" xlink:type="simple"/></inline-formula></p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Cross-classification of duodenal ulcer patients according to operation and dumping severity</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Operation</th><th align="center" valign="middle"  colspan="4"  >Dumping Severity</th></tr></thead><tr><td align="center" valign="middle" >None</td><td align="center" valign="middle" >Slight</td><td align="center" valign="middle" >Moderate</td><td align="center" valign="middle" >Total</td></tr><tr><td align="center" valign="middle" >A</td><td align="center" valign="middle" >61</td><td align="center" valign="middle" >28</td><td align="center" valign="middle" >7</td><td align="center" valign="middle" >96</td></tr><tr><td align="center" valign="middle" >B</td><td align="center" valign="middle" >68</td><td align="center" valign="middle" >23</td><td align="center" valign="middle" >13</td><td align="center" valign="middle" >104</td></tr><tr><td align="center" valign="middle" >C</td><td align="center" valign="middle" >58</td><td align="center" valign="middle" >40</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >110</td></tr><tr><td align="center" valign="middle" >D</td><td align="center" valign="middle" >53</td><td align="center" valign="middle" >38</td><td align="center" valign="middle" >16</td><td align="center" valign="middle" >107</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" >240</td><td align="center" valign="middle" >129</td><td align="center" valign="middle" >48</td><td align="center" valign="middle" >417</td></tr></tbody></table></table-wrap><p>Source: Grizzle et al. [<xref ref-type="bibr" rid="scirp.60338-ref9">9</xref>] .</p><p>based on <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x64.png" xlink:type="simple"/></inline-formula> (see Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p.123) and Tomizawa [<xref ref-type="bibr" rid="scirp.60338-ref10">10</xref>] ). Note that the U model also fits well with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x65.png" xlink:type="simple"/></inline-formula> based on <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x66.png" xlink:type="simple"/></inline-formula> (see Agresti ([<xref ref-type="bibr" rid="scirp.60338-ref2">2</xref>] , p.81) and Tomizawa [<xref ref-type="bibr" rid="scirp.60338-ref10">10</xref>] ).</p><p>For testing the hypothesis that the I model holds assuming that the logit U model holds, the difference be- tween the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x67.png" xlink:type="simple"/></inline-formula> values for the I model and the logit U model is 6.61 based on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x68.png" xlink:type="simple"/></inline-formula>. Therefore this hypothesis is rejected at the 0.05 level. Hence the logit U model is preferable to the I model for these data.</p><p>Also the CDE model fits these data poorly with <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x69.png" xlink:type="simple"/></inline-formula> based on<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x69.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x70.png" xlink:type="simple"/></inline-formula>. We see that the rejection of the hypothesis that the I model holds assuming that the logit U model holds is caused by the influence of the lack of structure of the CDE model (i.e., the lack of equality of the sum of probabilities of concordance and those of discordance), because the hypothesis that the I model holds assuming that the logit U model holds is equivalent to the CDE model from Theorem 1.</p></sec><sec id="s4"><title>4. Concluding Remarks</title><p>When the I model fits the data poorly, Theorem 1 may be useful for seeing the reason for the poor fit; namely, which of the lack of structure of the CDE model and that of the logit U model influences stronger.</p><p>From Theorem 1 we point out that the hypothesis that the I model holds under the assumption that the logit U model holds is equivalent to the hypothesis that the CDE model holds.</p><p>The U model indicates the constant of the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x71.png" xlink:type="simple"/></inline-formula> local odds ratios defined for adjacent rows and adjacent columns. On the other hand, the logit U model indicates the constant of the odds ratios for the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x72.png" xlink:type="simple"/></inline-formula> <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x73.png" xlink:type="simple"/></inline-formula> tables obtained by taking all pairs of adjacent rows and all dichotomous collapsing of the response. Thus, when the I model fits the data poorly, if the user wants to see the structure of cumulative probabilities (i.e., the structures of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x74.png" xlink:type="simple"/></inline-formula> collapsed <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-1240565x75.png" xlink:type="simple"/></inline-formula> tables), then Theorem 1 may be preferable to preceding studies which are described in Section 1.</p></sec><sec id="s5"><title>Acknowledgements</title><p>We thank the referee for comments and suggestions.</p></sec><sec id="s6"><title>Cite this paper</title><p>KoujiTahata,NobukoMiyamoto,SadaoTomizawa, (2015) Decomposition of Independence Using the Logit Uniform Association Model and Equality of Concordance and Discordance for Two-Way Classifications. Open Journal of Statistics,05,514-518. doi: 10.4236/ojs.2015.56054</p></sec></body><back><ref-list><title>References</title><ref id="scirp.60338-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Goodman, L.A. (1979) Simple Models for the Analysis of Association in Cross-Classifications Having Ordered Categories. Journal of the American Statistical Association, 74, 537-552. &lt;/br&gt;http://dx.doi.org/10.1080/01621459.1979.10481650</mixed-citation></ref><ref id="scirp.60338-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Agresti, A. (1984) Analysis of Ordinal Categorical Data. Wiley, New York.</mixed-citation></ref><ref id="scirp.60338-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Tomizawa, S., Miyamoto, N. and Sakurai, M. (2008) Decomposition of Independence Model and Separability of Its Test Statistic for Two-Way Contingency Tables with Ordered Categories. Advances and Applications in Statistics, 8, 209-218.</mixed-citation></ref><ref id="scirp.60338-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kendall, M.G. (1945) The Treatment of Ties in Ranking Problems. Biometrika, 33, 239-251.&lt;/br&gt;http://dx.doi.org/10.1093/biomet/33.3.239</mixed-citation></ref><ref id="scirp.60338-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Tahata, K., Miyamoto, N. and Tomizawa, S. (2008) Decomposition of Independence Using Pearson, Kendall and Spearman’s Correlations and Association Model for Two-Way Classifications. Far East Journal of Theoretical Statistics, 25, 273-283.</mixed-citation></ref><ref id="scirp.60338-ref6"><label>6</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Stuart</surname><given-names> A. </given-names></name>,<etal>et al</etal>. (<year>1963</year>)<article-title>Calculation of Spearman’s Rho for Ordered Two-Way Classifications</article-title><source> The American Statistician</source><volume> 17</volume>,<fpage> 23</fpage>-<lpage>24</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.60338-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Kendall, M. and Gibbons, J.D. (1990) Rank Correlation Methods. 5th Edition, Edward Arnold, London.</mixed-citation></ref><ref id="scirp.60338-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Tahata, K. and Tomizawa, S. (2014) Symmetry and Asymmetry Models and Decompositions of Models for Contingency Tables. SUT Journal of Mathematics, 50, 131-165.</mixed-citation></ref><ref id="scirp.60338-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Grizzle, J.E., Starmer, C.F. and Koch, G.G. (1969) Analysis of Categorical Data by Linear Models. Biometrics, 25, 489-504. &lt;/br&gt;http://dx.doi.org/10.2307/2528901</mixed-citation></ref><ref id="scirp.60338-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Tomizawa, S. (1992) More Parsimonious Linear-by-Linear Association Model in the Analysis of Cross-Classifications Having Ordered Categories. Biometrical Journal, 34, 129-140. &lt;/br&gt;http://dx.doi.org/10.1002/bimj.4710340202</mixed-citation></ref></ref-list></back></article>