コクラン・アーミテージ検定     Last modified: Nov 07, 2002

外的基準が連続変数であるときに,各群の比率が外的基準と線形傾向があるかどうかを検定する。
注:元々は線形傾向の検定であったが,傾向検定(単調増加の検定)として使われることも多い。そのような場合には,外的基準は必ずしも必要ない(外的基準がない場合には,等間隔な整数値を仮定する)。


例題

 「表 1 のようなデータにおいて,5 群の比率が独立変数と線形傾向があるか検定しなさい。」

表 1.数値例
独立変数 ケース数 陽性数 比率 比率の推定値
10 30 2 0.067 -0.047  
20 35 4 0.114 0.169  
30 47 14 0.298 0.385  
40 21 13 0.619 0.601  
50 45 39 0.867 0.817  


検定手順:

  1. 表 2 のように,k 個の群があるとき,第 i 群に対応する外的基準を Xi,ケース数を ni,そのうちで対象とする特性を持つケース数(ここでは陽性数と呼ぶことにする)を ri とする。

    表 2.コクラン・アーミテージ検定を適用するデータの形式
    第 1 群 第 2 群    k
    ケース数 n1 n2    nk
    陽性数 r1 r2    rk
    比率 p1 p2    pk
    外的基準値 X1 X2    Xk
    n=i=1kni,   r=i=1kri,   p¯=rn,   X¯=i=1kni Xin

  2. 各群の比率 pi=rini が外的基準 Xi と直線的な関係があると仮定すると,比率の予測値は p^i=a+b Xi と表せる。

    ab は次式で推定できる。 b=i=1kni (pip¯) (XiX¯)i=1kni (XiX¯)2=i=1kni pi Xin p¯ X¯i=1kni Xi2n X¯2a=p¯b X¯ 例題では,a=0.2624224b=0.02158385 となる。比率の予測値は,表 1 に示すようになる。

  3. 直線の傾き(b)が 0 であるかどうかの検定(傾きの有意性の検定)は,次式の χT2 が自由度 1χ2 分布に従うことを用いて検定できる。 χT2=b2p¯ (1p¯) i=1kni (XiX¯)2
  4. 各群の比率に直線では表しきれない傾向があるかどうか(直線性の検定)は,次式の χQ2 が自由度 k2χ2 分布に従うことを用いて検定できる。 χQ2=i=1kni (pip^i)2p¯ (1p¯)
  5. 各群の比率に差があるかどうかを検定するのは k×2 分割表の χ2 検定で行える。式の χH2 は,自由度 k1χ2 分布に従う。 χH2=n2r (nr) (i=1kri2nir2n)
  6. 比率の一様性の検定では,各群に順序関係があるときにもその情報を利用していない。
    コクラン・アーミテージ検定は次式のように,比率の一様性の検定を目的とする χH2 統計量を直線式で説明できる部分 χT2 と直線式からの乖離の部分 χQ2 に分割する。 (1)χH2=χT2+χQ2 実際の計算では,χH2χT2 を先に計算し,χQ2(1) 式から求めれば計算が簡単である。

  7. 結果の表示は表 3 のように行う。これらはちょうど一元配置分散分析において,全分散を級内分散と級間分散に分解することや結果を分散分析表の形で表現することに似ている。

    表 3.コクラン・アーミテージ検定の結果
    要因 χ2 自由度
    傾き χT2 1
    直線からの乖離 χQ2 k2
    合計(非一様性) χH2 k1

    例題では,表 4 のようにまとめられる。

    表 4.コクラン・アーミテージ検定の結果
    要因 χ2 自由度 有意確率
    傾き 68.57273 1 < 0.001
    直線からの乖離 4.014484 3 0.25990
    合計(非一様性) 72.58722 4 < 0.001

  8. 検定結果の解釈は以下のようになる。

    例題では,各群の比率は,独立変数値と直線的な傾向を持つといえる。

・ R で計算してみる


演習問題


応用問題


・ 計算プログラム [R]
・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI