例題:
「表 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))
カイ二乗値 自由度 P 値
トレンド 68.572732 1 1.222834e-16
直線からの乖離 4.014484 3 2.599043e-01
非一様性 72.587216 4 6.449422e-15