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@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. ______________________________________________ R-help@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.
