Hi (again), Please find here a correction of the function I send previously --------------------------------------------------------- tableChisqMH=function(x) { n=sum(x) G2=(n-1)*(tablepearson(x)$estimate^2) dll=1 p.value=1-pchisq(G2,dll) out=list(estimate=G2,dll=dll,p.value=p.value,dim=dim(x),name="Mantel-Haenszel Chi-square") return(out) } --------------------------------------------------------- Eric __________________________________________________ Eric Lecoutre Informaticien/Statisticien Institut de Statistique UCL (+32) (0)10 47 30 50 lecoutre at stat.ucl.ac.be http://www.stat.ucl.ac.be/ISpersonnel/lecoutre __________________________________________________ Le vrai danger, ce n''est pas quand les ordinateurs penseront comme des hommes, c''est quand les hommes penseront comme des ordinateurs. Sydney Harris