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 <[EMAIL PROTECTED]>
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
