対応のあるt検定 計算ツール|前後比較・治療効果・効果量 Cohen's dz 対応
対応のあるt検定(paired t-test / dependent t-test)を、差の平均・差の標準偏差・サンプル数から瞬時計算。t値・自由度・臨界値・両側p値の目安・95%信頼区間・効果量(Cohen's dz)まで一括算出。治療前後の血圧・体重・血糖値、ダイエット・トレーニング効果、広告接触前後の認知度、視力矯正の効果、A/Bテストの被験者内前後比較で必須の統計手法を、独立2群t検定・Wilcoxon符号順位検定との使い分けとともに、ICH E9統計原則・CONSORT声明・PMDA臨床試験ガイダンスの要点に沿って解説します。
対応のあるt検定 計算機
計算式と考え方
基本式
t = d̄ ÷ (s_d / √n)
自由度 df = n − 1
d̄ = 各ペアの差の平均、s_d = 差の標準偏差
対応のあるt検定は、同一被験者から得られた前後2測定値、あるいはマッチングされたペアの差を検定します。「差の平均が0である」という帰無仮説を検定し、有意に0から離れていれば「効果あり」と判定します。
手順
- 各ペアの差 d_i = 前 − 後 を計算
- 差の平均 d̄ と差の標準偏差 s_d を求める
- t = d̄ / (s_d / √n) を計算
- 自由度 df = n − 1 のt分布で臨界値と比較 / p値を算出
- |t| > 臨界値なら有意(帰無仮説を棄却)
信頼区間
95%CI = d̄ ± t*(df, 0.025) × (s_d / √n)
信頼区間には差の実質的な大きさとその精度が示されるため、p値だけでなく必ず報告することが推奨されます(CONSORT声明・ICH E9)。
対応あり/なしの違い
| 項目 | 対応のあるt検定 | 独立2群t検定 |
|---|---|---|
| データ構造 | 同一被験者の前後 or マッチペア | 異なる2群の独立サンプル |
| 自由度 | n − 1 | n1 + n2 − 2 |
| 個人差の扱い | 差を取ることで相殺 | 群間ばらつきに含まれる |
| 検出力 | ◎ 高(同じnで有意差を検出しやすい) | △ 中 |
| 典型例 | 治療前後の血圧、双子比較 | 薬剤群vsプラセボ群 |
| 前提 | 差が正規分布 | 両群が正規分布・等分散 |
被験者内デザイン(within-subject design)は個人差を排除できるため、被験者間デザインより少ないサンプルで同じ有意差を検出できます。ただし順序効果(先の測定が後に影響)や学習効果に注意が必要です。
適用条件と前提
- ペア構造:同一個体の前後、または一対一対応するマッチデータ。
- 差の正規性:差 d_i が正規分布に従うことを仮定(元データの正規性ではなく差の正規性)。Shapiro-Wilk検定・Kolmogorov-Smirnov検定・QQプロットで確認。
- 独立性:ペア同士は独立。同じ被験者が複数ペアに含まれる場合は不適用(混合効果モデル・反復測定ANOVAが必要)。
- 連続量:測定値は間隔尺度以上(血圧・体重・スコアなど)。名義尺度・順序尺度はMcNemar検定・Wilcoxon符号順位検定を使用。
- 外れ値の影響:n≧30なら中心極限定理でロバストになるが、n<30では極端な外れ値の影響大。
差の正規性が疑わしい場合、または順序尺度データの場合はWilcoxon符号順位検定(ノンパラ版)を使用します。
効果量 Cohen's dz の解釈
p値だけでは「差の実質的な大きさ」が伝わらないため、効果量(effect size)の併記が必須です。対応のあるt検定ではCohen's dzを使用。
Cohen's dz = d̄ ÷ s_d
| |dz| | 効果の大きさ | 意味 |
|---|---|---|
| 0.20 未満 | ごく小 | 実質的にほぼ差なし |
| 0.20 〜 0.50 | 小(small) | 統計的有意でも実質差は小 |
| 0.50 〜 0.80 | 中(medium) | 目視で分かる程度の差 |
| 0.80 以上 | 大(large) | 明確な効果 |
大サンプル(n=1,000等)ではごく小さな効果量でもp<0.05になるため、臨床的意義を判断するには効果量と信頼区間で評価することが重要です。
臨界値 t*(両側5%)早見表
| 自由度 df | t*(両側5%) | t*(両側1%) |
|---|---|---|
| 5 | 2.571 | 4.032 |
| 10 | 2.228 | 3.169 |
| 15 | 2.131 | 2.947 |
| 20 | 2.086 | 2.845 |
| 29 | 2.045 | 2.756 |
| 50 | 2.009 | 2.678 |
| 100 | 1.984 | 2.626 |
| 200 | 1.972 | 2.601 |
| ∞(正規分布) | 1.960 | 2.576 |
|t| > t*(df, α/2) なら「有意差あり」。df が大きくなるほど t* は正規分布のz値(1.96 / 2.576)に近づきます。
計算例
降圧薬の臨床試験
高血圧患者30名に降圧薬を8週間投与。投与前後の収縮期血圧の差の平均 d̄=12 mmHg、SD s_d=15 mmHg。
t = 12 ÷ (15 / √30) = 12 ÷ 2.74 ≒ 4.38、df=29、t*(29, 0.025)=2.045 → 有意差あり(p<0.001)。
Cohen's dz = 12/15 = 0.80(大きな効果)。95%CI: 6.4 〜 17.6 mmHg。
ダイエットプログラムの効果
被験者20名の3ヶ月ダイエット。体重差 d̄=3.5 kg、SD s_d=4.2 kg。
t = 3.5 ÷ (4.2/√20) = 3.5 ÷ 0.94 ≒ 3.72、df=19、t*=2.093 → 有意差あり。Cohen's dz = 0.83(大)。95%CI: 1.5 〜 5.5 kg。
広告接触前後の認知度
被験者100名の広告視聴前後の認知スコア(0-10点)。差 d̄=0.6 点、SD s_d=2.8 点。
t = 0.6 ÷ (2.8/√100) = 0.6 ÷ 0.28 ≒ 2.14、df=99、t*=1.984 → 有意差あり(p<0.05)。
Cohen's dz = 0.21(小)。統計的有意でも実質効果は小さいケース。
業界での使い方
医薬品・臨床試験
降圧薬・血糖降下薬・鎮痛剤等の効果検証で頻用。被験者内比較は個人差を排除できるため少数例(n=20〜50)でも有意差を検出可能。ICH E9・CONSORT声明で対応t検定の使用条件と報告方法が規定されています。
ヘルスケア・フィットネス
ダイエット・トレーニングプログラムの効果測定。体重・体脂肪率・筋力・持久力の前後比較で標準的な統計手法。効果量Cohen's dz≧0.5を「臨床的に意義ある変化(MCID)」の目安に。
教育・研修効果測定
研修前後のテストスコア、e-Learning受講前後の理解度、資格試験対策講座の効果検証で使用。ペアデータのため受講者の個人差を排除して純粋な学習効果を評価できます。
マーケティング・広告効果
広告接触前後のブランド認知・購入意向・NPS等の変化検証。同一パネル被験者を追跡調査し、対応t検定で純粋な広告効果を測定します。効果量が小さい場合は施策の投資対効果の再検討が必要。
UX・ユーザビリティテスト
UI改修前後のタスク完了時間、エラー数、SUSスコア(システムユーザビリティ尺度)の比較で使用。少数被験者(n=10〜20)でも有意差が出ることが多い被験者内デザインの利点を活かします。
製造業・工程改善
改善施策前後の生産性・不良率・作業時間の比較。同じライン・同じ製品での前後比較なので個体差の影響がなく、少ないサンプルでも改善効果を統計的に立証できます。
ノンパラ手法との使い分け
差の正規性が満たされない場合や順序尺度データでは、パラメトリックな対応t検定は使えません。代替手法として以下を使用します。
| 手法 | データ型 | 特徴 | 使用場面 |
|---|---|---|---|
| 対応のあるt検定 | 連続・差が正規分布 | 検出力最高 | 血圧・体重・スコア |
| Wilcoxon符号順位検定 | 連続・順序 | ノンパラ・順位ベース | アンケート順位・非正規 |
| 符号検定(sign test) | 連続・順序 | ノンパラ・+/−のみ | 方向性のみ検証 |
| McNemar検定 | 二値(0/1) | χ²ベース | 治療前後の陽性/陰性 |
| ブートストラップt | 連続 | 再標本化による | 小標本・非正規 |
よくある間違い・注意点
- 元データの正規性を確認して差の正規性を確認せず — 対応t検定の前提は「差 d_i の正規分布」であって、元データの正規分布ではありません。差のヒストグラム・QQプロット・Shapiro-Wilk検定で確認します。
- 独立2群t検定と混用 — 同一被験者の前後を独立2群として扱うと、個人差が誤差項に含まれて検出力が低下。「同じ人(あるいはマッチペア)か?」を必ず確認します。
- 効果量を報告しない — CONSORT声明・APAガイドラインではp値だけでなく効果量(Cohen's dz等)と信頼区間の報告が推奨。大サンプルでは効果量が小さくてもp<0.05になり実質的意味が伝わりません。
- 多重比較補正を忘れる — 複数変数を同時に検定する場合、有意水準を Bonferroni補正(α/k)またはHolm法で調整。5変数を5%水準で検定すると、少なくとも1つが偶然有意になる確率は約23%です。
- 順序効果を無視 — 「前」測定が「後」測定に影響する可能性(学習・練習・慣れ・キャリーオーバー)。可能ならクロスオーバーデザイン・ランダム化順序で対応します。
- 脱落者(drop-out)を安易に除外 — 前後両方の測定が揃った人のみ分析すると生存者バイアス。ITT解析(intention-to-treat)や欠測値補完(LOCF・MI等)を検討します。
- 両側検定と片側検定の使い分けミス — 明確な仮説方向(「下がる」と事前予測)がある場合のみ片側。事前登録なしに事後で片側にすると p値操作(p-hacking)と見なされます。
関連する規格・統計手法
- ICH E9 ― Statistical Principles for Clinical Trials. 医薬品臨床試験の統計原則。対応t検定を含む検定手順の要件を規定。
- CONSORT声明(CONSORT 2010) ― Consolidated Standards of Reporting Trials. ランダム化比較試験の報告ガイドライン。効果量・信頼区間の報告を必須化。
- JIS Z 8402-1〜6 ― 測定方法及び測定結果の精度(真度及び精度)。統計解析全般の日本国内共通ルール。
- JIS Z 9041シリーズ ― データの統計的な解釈方法。平均値と分散に関する検定手順。
- ISO 5479 ― Statistical interpretation of data - Tests for departure from the normal distribution. 正規性検定の国際規格。
- APA Publication Manual ― American Psychological Association 論文執筆規約。効果量報告のガイドライン。
- PMDA臨床試験統計ガイダンス ― 医薬品医療機器総合機構による臨床試験統計解析の指針。
参考文献・公的資料
より詳細な統計理論・臨床試験プロトコルは、以下の公的機関・学会の資料を参照してください。
- 医薬品医療機器総合機構(PMDA) ― 医薬品承認審査・臨床試験ガイダンス。ICH E9の日本語版・関連通知を公開。
- 厚生労働省 治験 ― GCP省令・臨床試験の運用ガイドライン。
- ICH Efficacy Guidelines(E9他) ― 国際医薬品規制調和会議の臨床試験統計原則。
- CONSORT Statement ― ランダム化比較試験の報告ガイドライン公式。
- 総務省統計局 ― 国家統計における統計手法の技術資料。
- 日本統計学会 ― 統計手法の学術情報・学会誌。
- 統計数理研究所 ― 大学共同利用機関。統計手法研究の一次資料。
- 日本計量生物学会(J-STAGE) ― 生物統計・臨床統計の学術論文。
- 日本産業標準調査会(JISC) ― JIS Z 8402・JIS Z 9041等の統計規格。
- NIST Statistical Engineering Division ― 米国国立標準技術研究所の統計手法ガイド。
よくある質問(FAQ)
独立t検定との使い分けは?
同一対象の前後(または一対一対応するマッチペア)なら対応あり、異なる被験者の2群(薬剤群 vs プラセボ群など)なら独立。対応ありの方が個人差を排除できるため検出力が高く、少ないサンプル数で有意差を検出できます。「同じ人か?」がキーになります。
Wilcoxon符号順位検定との違いは?
対応t検定は差が正規分布を前提とするパラメトリック検定。Wilcoxon符号順位検定は順位ベースのノンパラメトリック版で、正規性を仮定せず順序尺度データでも使用可能。差が明らかに非正規(極端な歪み・多峰)や順序尺度データではWilcoxonを使います。
小サンプル(n=5〜10)でも使える?
使えますが差の正規性の確認が困難で、外れ値の影響も大きくなります。n=5でも効果が非常に大きければ有意差は出ますが、少数例での結論は再現性が低いため、独立した追試験を実施すべきです。可能ならWilcoxon検定やブートストラップ法を併用します。
p値だけ報告すればOK?
ダメです。CONSORT声明・ICH E9・APA Publication Manualでは、p値・効果量・95%信頼区間の3点セットでの報告が推奨。大サンプルではごく小さな効果でもp<0.05になるため、実質的意義を判断するには効果量・信頼区間が必須です。
効果量 Cohen's dz はどう解釈する?
Jacob Cohen(1988)による経験的分類で、dz=0.2(小)、0.5(中)、0.8(大)が目安。医療分野ではMCID(臨床的に意義ある最小変化)と併用し、統計的有意 + 効果量が実用閾値以上の場合に「意味のある変化」と判定します。
両側検定と片側検定の使い分けは?
基本は両側検定(差の方向を事前に決めない)。事前に「絶対に下がる」等の方向性仮説がある場合のみ片側検定。事前登録(プロトコル)に明記されていない事後の片側検定変更はp値操作(p-hacking)と見なされ、学術的に否定されます。
多重比較の補正はどうする?
複数の変数(例:5つの症状)を同時に検定する場合、Bonferroni補正(α/k、5個なら0.05/5=0.01)やHolm-Bonferroni法、Benjamini-Hochberg(FDR)法で調整します。補正なしでは第一種の過誤率が膨張します。
脱落者(drop-out)が出た場合は?
前後両方揃った人だけで分析すると生存者バイアスで効果が過大評価される恐れ。ITT解析(intention-to-treat:ランダム化割付通りに全員を解析)、LOCF(Last Observation Carried Forward)、MI(多重代入法)などで欠測値補完を検討します。
3時点以上の測定はどうする?
対応t検定は2時点のみ。3時点以上は反復測定ANOVA(repeated measures ANOVA)または混合効果モデル(LMM)を使用。時点間の比較には多重比較補正(Tukey HSD等)が必要です。
Excelでの計算方法は?
=T.TEST(範囲1, 範囲2, 尾, 検定タイプ)で、尾=2(両側)、検定タイプ=1(対応あり)を指定するとp値が直接出ます。Data Analysis Toolpakの「t検定:一対の標本による平均の検定」でも詳細な出力が得られます。
R言語での実装は?
t.test(前, 後, paired=TRUE)で、t値・df・p値・95%CI・平均差が一括出力されます。効果量はeffsizeパッケージのcohen.d(前, 後, paired=TRUE)で計算可能。CONSORT準拠の報告フォーマットを自動生成できます。
順序効果・学習効果の対策は?
「前」測定が「後」測定に影響する可能性がある場合(学習・慣れ・キャリーオーバー)、クロスオーバーデザイン(半数はABの順、半数はBAの順に測定)やウォッシュアウト期間(測定間隔を十分に空ける)で対応します。順序ランダム化を事前登録することが重要。