コクランアーミテージ検定

統計学用語辞典

コクラン・アーミテージ検定

例題：
　「表 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

R による解析:

前もって定義されている関数を使う場合

> prop.trend.test (c(2, 4, 14, 13, 39), c(30, 35, 47, 21, 45))

	Chi-squared Test for Trend in Proportions

data:  c(2, 4, 14, 13, 39) out of c(30, 35, 47, 21, 45) ,
 using scores: 1 2 3 4 5 
X-squared = 68.5727, df = 1, p-value = < 2.2e-16

新たに定義した関数を使う場合

> Cochran.Armitage(c(2, 4, 14, 13, 39), c(30, 35, 47, 21, 45))

        Chi-sq. d.f.      P value
Trend 68.572732    1 1.110223e-16
Quad.  4.014484    3 2.599043e-01
Homo. 72.587216    4 6.439294e-15

コクラン・アーミテージ検定

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


例題：
　「表 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

   
   
2. 各群の比率 pi = ri / ni が外的基準 Xi と直線的な関係があると仮定すると，比率の予測値は と表せる。
   a，b は次式で推定できる。
   
   
   例題では，a = -0.2624224，b = 0.02158385 となる。比率の予測値は，表 1 に示すようになる。
   
3. 直線の傾き（b）が 0 であるかどうかの検定（傾きの有意性の検定）は，次式の χ²T が自由度 1 の χ² 分布に従うことを用いて検定できる。
   
   
4. 各群の比率に直線では表しきれない傾向があるかどうか（直線性の検定）は，次式の χ²Q が自由度 k - 2 の χ² 分布に従うことを用いて検定できる。
   
   
5. 各群の比率に差があるかどうかを検定するのは k×2 分割表の χ² 検定で行える。式の χ²H は，自由度 k - 1 の χ² 分布に従う。
   
   
6. 比率の一様性の検定では，各群に順序関係があるときにもその情報を利用していない。
   コクラン・アーミテージ検定は次式のように，比率の一様性の検定を目的とする χ²H 統計量を直線式で説明できる部分 χ²T と直線式からの乖離の部分 χ²Q に分割する。
   
   実際の計算では，χ²H と χ²T を先に計算し，χ²Q は上の式から求めれば計算が簡単である。
   
7. 結果の表示は表 3 のように行う。これらはちょうど一元配置分散分析において，全分散を級内分散と級間分散に分解することや結果を分散分析表の形で表現することに似ている。
   

   表 3．コクラン・アーミテージ検定の結果

   要因	χ2値	自由度
   傾き	χ2T	1
   直線からの乖離	χ2Q	k - 2
   合計（非一様性）	χ2H	k - 1

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

   表 4．コクラン・アーミテージ検定の結果

   要因	χ2値	自由度	有意確率
   傾き	68.57273	1	＜ 0.001
   直線からの乖離	4.014484	3	0.25990
   合計（非一様性）	72.58722	4	＜ 0.001

   
   
8. 検定結果の解釈は以下のようになる。
   
   • χ²T が大きな値であれば（有意確率が小さければ）傾きが 0 ではないことを意味する。すなわち，群の比率は，外的基準値により増減することを意味する。
     
   • χ²Q が大きな値であれば（有意確率が小さければ）直線性の仮定は受け入れられない。すなわち，直線では表しきれないことを意味する。
     
   • 両者をあわせて考えると，χ²T が大きく，かつ，χ²Q が小さければ，各群の比率は直線式で表現できることになる。
   
   例題では，各群の比率は，独立変数値と直線的な傾向を持つといえる。