[R] R equivalent of SAS Cochran-Mantel-Haenszel tests?

[David Freedman]

How about 'cmh_test' in the coin package?

>From the PDF: The null hypothesis of the independence of y and x is tested,
block defines an optional factor for stratification. chisq_test implements
Pearson’s chi-squared test, cmh_test the Cochran-Mantel-Haenzsel test and
lbl_test the linear-by-linear association test for ordered data.

David Freedman



Michael Friendly wrote:
> 
> 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
> 
> -- 
> Michael Friendly     Email: friendly AT yorku DOT ca 
> Professor, Psychology Dept.
> York University      Voice: 416 736-5115 x66249 Fax: 416 736-5814
> 4700 Keele Street    http://www.math.yorku.ca/SCS/friendly.html
> Toronto, ONT  M3J 1P3 CANADA
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
> 
> 

-- 
View this message in context: http://www.nabble.com/R-equivalent-of-SAS-Cochran-Mantel-Haenszel-tests--tp21914012p21918306.html
Sent from the R help mailing list archive at Nabble.com.