From: Thomas Lumley <tlumley_at_u.washington.edu>

Date: Sat 09 Jul 2005 - 03:34:19 EST

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 Received on Sat Jul 09 03:43:23 2005

Date: Sat 09 Jul 2005 - 03:34:19 EST

On Fri, 8 Jul 2005, Florian Hahne wrote:

> Dear listers,

*>
**> I am trying to compute the exact conditional test given strata margins
**> of a 2 by 2 by K array using the mantelhaen.test function to get a
**> common odds ratio estimate.
**> The estimate for the test on the following data is 0, which in my
**> opinion dosen't make any sense.
**>
*

Indeed it doesn't. Some probabilities are underflowing to zero. It is informative to look at the estimates from

mantelhaen.test(round(x/s), exact=TRUE) for various values of s. The estimate is stable at about 0.04 for large s until s gets to just above 1.8, then quickly decreases to zero.

I think you need exact=FALSE for data sets this large.

-thomas

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 Received on Sat Jul 09 03:43:23 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:33:26 EST
*