[R] Fast version of Fisher's Exact Test
Steven McKinney
smckinney at bccrc.ca
Sat Apr 9 00:44:53 CEST 2011
Do you mean a test something such as this?
> fisher.test(matrix(c(502,498,490, 510), nrow = 2))
Fisher's Exact Test for Count Data
data: matrix(c(502, 498, 490, 510), nrow = 2)
p-value = 0.6228
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.8770113 1.2550998
sample estimates:
odds ratio
1.049119
This runs quickly on my machine.
> system.time(fisher.test(matrix(c(502,498,490, 510), nrow = 2)))
user system elapsed
0.008 0.001 0.010
> sessionInfo()
R version 2.12.2 (2011-02-25)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_CA.UTF-8/en_CA.UTF-8/C/C/en_CA.UTF-8/en_CA.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
loaded via a namespace (and not attached):
[1] tools_2.12.2
>
Can you provide an example that is running slowly for you?
Steven McKinney
________________________________________
From: r-help-bounces at r-project.org [r-help-bounces at r-project.org] On Behalf Of Jim Silverton [jim.silverton at gmail.com]
Sent: April 8, 2011 9:43 AM
To: r-help at r-project.org
Subject: Re: [R] Fast version of Fisher's Exact Test
Is anyone aware of a fast way of doing fisher's exact test for a series of 2
x 2 tables in R? The fisher.test is really slow if n1=1000 and n2 = 1000.
--
Thanks,
Jim.
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org 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.
More information about the R-help
mailing list