コクラン・アーミテージ検定
例題:
「表 1 のようなデータにおいて,5 群の比率が独立変数と線形傾向があるか検定しなさい。」
独立変数 | ケース数 | 陽性数 | 比率 | (比率の推定値) |
---|---|---|---|---|
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 群の比率が独立変数と線形傾向があるか検定しなさい。」
独立変数 | ケース数 | 陽性数 | 比率 | (比率の推定値) |
---|---|---|---|---|
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 |
検定手順:
- 表 2 のように,k 個の群があるとき,第 i 群に対応する外的基準を Xi,ケース数を ni,そのうちで対象とする特性を持つケース数(ここでは陽性数と呼ぶことにする)を ri とする。
表 2.コクラン・アーミテージ検定を適用するデータの形式 群 第 1 群 第 2 群 … 第 k 群 ケース数 n1 n2 … nk 陽性数 r1 r2 … rk 比率 p1 p2 … pk 外的基準値 X1 X2 … Xk
- 各群の比率 pi = ri / ni が外的基準 Xi と直線的な関係があると仮定すると,比率の予測値は と表せる。
a,b は次式で推定できる。
例題では,a = -0.2624224,b = 0.02158385 となる。比率の予測値は,表 1 に示すようになる。
- 直線の傾き(b)が 0 であるかどうかの検定(傾きの有意性の検定)は,次式の χ2T が自由度 1 の χ2 分布に従うことを用いて検定できる。
- 各群の比率に直線では表しきれない傾向があるかどうか(直線性の検定)は,次式の χ2Q が自由度 k - 2 の χ2 分布に従うことを用いて検定できる。
- 各群の比率に差があるかどうかを検定するのは k×2 分割表の χ2 検定で行える。式の χ2H は,自由度 k - 1 の χ2 分布に従う。
- 比率の一様性の検定では,各群に順序関係があるときにもその情報を利用していない。
コクラン・アーミテージ検定は次式のように,比率の一様性の検定を目的とする χ2H 統計量を直線式で説明できる部分 χ2T と直線式からの乖離の部分 χ2Q に分割する。
実際の計算では,χ2H と χ2T を先に計算し,χ2Q は上の式から求めれば計算が簡単である。
- 結果の表示は表 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
- 検定結果の解釈は以下のようになる。
- χ2T が大きな値であれば(有意確率が小さければ)傾きが 0 ではないことを意味する。すなわち,群の比率は,外的基準値により増減することを意味する。
- χ2Q が大きな値であれば(有意確率が小さければ)直線性の仮定は受け入れられない。すなわち,直線では表しきれないことを意味する。
- 両者をあわせて考えると,χ2T が大きく,かつ,χ2Q が小さければ,各群の比率は直線式で表現できることになる。
例題では,各群の比率は,独立変数値と直線的な傾向を持つといえる。 - χ2T が大きな値であれば(有意確率が小さければ)傾きが 0 ではないことを意味する。すなわち,群の比率は,外的基準値により増減することを意味する。
コクラン・アーミテージ検定と同じ種類の言葉
検定に関連する言葉 | クラスカルウォリス検定 フリードマン検定 コクランアーミテージ検定 スミルノフグラブス検定 コクランのQ検定 |
- コクラン・アーミテージ検定のページへのリンク