[R] [Fwd: Trend test and test for homogeneity of odd-ratios]

From: Jacques VESLOT <jacques.veslot_at_good.ibl.fr>
Date: Sat 19 Aug 2006 - 02:13:20 EST


I partly answered my question since independence_test() function in coin package apparently do Cochran-Armitage trend test just like Eric Lecoutre's function tabletrend() - slightly modified here:

> independence_test(pheno ~ geno, data = dat2, teststat = "quad", scores = list(geno = c(0, 1, 2)))

         Asymptotic General Independence Test

data: pheno by groups 1 < 2 < 3
chi-squared = 0.2268, df = 1, p-value = 0.6339

> tabletrend(with(dat2, table(pheno, geno)))
[1] 0.6338308

Dear r-users,

I am looking for some R functions to do Cochran-Armitage trend test for 2*3 tables (binary phenotype vs. genotypes) and for testing the homogeneity of odds ratios within 2*3*k tables (binary phenotype vs. genotypes vs. strata).

In R-Help archives, I've found a 2003 script by Eric Lecoutre for Cochran-Armitage trend test and a script for Breslow-Day test for 2*2*k tables.

Could someone please tell we if there were some functions available on CRAN to do such tests ?

Thanks,

jacques



R-help@stat.math.ethz.ch 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. Received on Sat Aug 19 02:36:52 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 19 Aug 2006 - 04:21:15 EST.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.