Skip to content
Prev 169776 / 398506 Next

R equivalent of SAS Cochran-Mantel-Haenszel tests?

In SAS, for a two-way (or 3-way, stratified) table, the CMH option in 
SAS PROC FREQ gives
3 tests that take ordinality of the factors into account, for both 
variables, just the column variable
or neither.   Is there an equivalent in R?
The mantelhaen.test in stats gives something quite different (a test of 
conditional independence for
*nominal* factors in a 3-way table).

e.g. I'd like to reproduce:
*-- CMH tests;
proc freq data=sexfun order=data;
  weight count;
  tables husband * wife / cmh chisq nocol norow;
  run;

                              The FREQ Procedure
                     Summary Statistics for Husband by Wife
          Cochran-Mantel-Haenszel Statistics (Based on Table Scores)

        Statistic    Alternative Hypothesis    DF       Value      Prob

            1        Nonzero Correlation        1     10.0142    0.0016
            2        Row Mean Scores Differ     3     12.5681    0.0057
            3        General Association        9     16.7689    0.0525