例3.7 Cochran-Mantel-Haenszel統計量

データセットMigraineには、偏頭痛治療の臨床試験に関する 仮説データが含まれています。この試験では、男女両方からなる複数の被験者が、新薬かプラシーボのいずれかによる処置を受けます。処置に対する被験者の反応は、'Better'または'Same'のいずれかでコード化されます。これらのデータはセルカウントとして記録され、処置と反応の組み合わせごとの被験者数が変数Countに記録されます。

data Migraine;
   input Gender $ Treatment $ Response $ Count @@;
   datalines;
female Active  Better 16   female Active  Same 11
female Placebo Better  5   female Placebo Same 20
male   Active  Better 12   male   Active  Same 16
male   Placebo Better  7   male   Placebo Same 19
;

次のPROC FREQステートメントは、Genderにより層化された多次元クロス表を作成します。この表は、Treatmentが行でResponseが列になります。TABLESステートメントのRELRISKオプションは、行がTreatmentで列がResponseの2元表に関するオッズ比と相対リスクを要求します。PLOTS=オプションは、Genderの各レベルの相対リスクとその信頼限界を示す相対リスクプロットを要求します。 CMHオプションは、この多次元クロス表に関するCochran-Mantel-Haenszel統計量を要求します。この層化された表の場合、CMHオプションは、共通相対リスクの推定値と、オッズ比の等質性に関するBreslow-Day検定も提供します。NOPRINTオプションが指定されているため、クロス集計表は表示されません。

ods graphics on;
proc freq data=Migraine;
   tables Gender*Treatment*Response / 
          relrisk plots(only)=relriskplot(stats) cmh noprint; 
   weight Count;
   title 'Clinical Trial for Treatment of Migraine Headaches';
run;
ods graphics off;

出力3.7.1出力3.7.4に、分析結果を示します。相対リスクプロット(図3.7.1)には、Genderの2つのレベルに関する相対リスクと信頼限界が示されています。出力3.7.2にCMH統計量を示します。層化された表の場合、3つのCMH統計量により同じ仮説の検定が行われます。有意な値(0.004)は、性別に関する補正を行った後に、処置と応答間に連関性が残ることを示しています。

CMHオプションは、出力3.7.3に示すような全体的な相対リスクに関する表も作成します。これは予備試験であるため、相対リスクの推定値により新薬の効果が評価されます。"Cohort (Col1 Risk)"の値は、最初の列に対応する推定値(改善のリスク)となります。新薬により偏頭痛が改善される確率は、プラシーボにより偏頭痛が改善される確率よりも2倍以上高くなります。

Breslow-Day検定での大きな値(0.2218)は、オッズ比に関して有意な性差が存在しないことを示しています(出力3.7.4を参照)。

出力3.7.1 相対リスク図
 相対リスク図

出力3.7.2 Cochran-Mantel-Haenszel統計量
Cochran-Mantel-Haenszel Statistics (Based on Table Scores)
Statistic Alternative Hypothesis DF Value Prob
1 Nonzero Correlation 1 8.3052 0.0040
2 Row Mean Scores Differ 1 8.3052 0.0040
3 General Association 1 8.3052 0.0040

出力3.7.3 CMHオプション:共通相対リスク
Estimates of the Common Relative Risk (Row1/Row2)
Type of Study Method Value 95% Confidence Limits
Case-Control Mantel-Haenszel 3.3132 1.4456 7.5934
(Odds Ratio) Logit 3.2941 1.4182 7.6515
Cohort Mantel-Haenszel 2.1636 1.2336 3.7948
(Col1 Risk) Logit 2.1059 1.1951 3.7108
Cohort Mantel-Haenszel 0.6420 0.4705 0.8761
(Col2 Risk) Logit 0.6613 0.4852 0.9013

出力3.7.4 CMHオプション: Breslow-Day検定
Breslow-Day Test for
Homogeneity of the Odds Ratios
Chi-Square 1.4929
DF 1
Pr > ChiSq 0.2218