[R] exact conditional mantelhaen.test estimate is 0 ?!

From: Florian Hahne <f.hahne_at_dkfz-heidelberg.de>
Date: Sat 09 Jul 2005 - 01:00:21 EST


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.

x <- array(c(53, 6098, 1006, 4521, 63, 8070, 1163, 6137), dim=c(2,2,2))  x
, , 1

     [,1] [,2]
[1,] 53 1006
[2,] 6098 4521

, , 2

     [,1] [,2]
[1,] 63 1163
[2,] 8070 6137

mantelhaen.test(x, exact=TRUE)

      Exact conditional test of independence in 2 x 2 x k tables

data: x
S = 116, p-value < 2.2e-16
alternative hypothesis: true common odds ratio is not equal to 1 95 percent confidence interval:
 0.0000000000 0.0001190186
sample estimates:
common odds ratio

                0

Could this be some kind of memory overflow issue since some of the values are quite high? Or am I completely off here?

I would appreciate any help on this,

Florian

-- 
Florian Hahne <f.hahne@dkfz-heidelberg.de>

______________________________________________
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 01:05:38 2005

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