Hi all,
I have just uploaded a new package to the ftp site called exact2x2. The package
does exact conditional tests for 2x2 tables. Specifically, it does Fisher's
exact test and Blaker's exact test (see Blaker, 2000, Canadian Journal of
Statistics, 783-798). The motivation for the package is to have inferences from
p-values and confidence intervals match as much as is possible for the tests,
so that when p<0.05 the 95 percent confidence interval on the odds ratio does
not contain 1.
Here is an example which shows the problem with fisher.test in this regard:
> fisher.test(matrix(c(6,12,12,5),2,2))
Fisher's Exact Test for Count Data
data: matrix(c(6, 12, 12, 5), 2, 2)
p-value = 0.04371
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.03888003 1.05649145
sample estimates:
odds ratio
0.2189021
Using the exact2x2 function, there are two options for two-sided Fisher's
exact test: "minlike" (default) is the usual one that defines extreme in terms
of values smaller than the observed probability, and "central" is the one that
defines the two-sided p-value as twice the minimum of the one-sided ones. The
problem with fisher.test is that it gives confidence intervals which are an
inversion of the later but p-values from the former. Here is that example
from exact2x2:
> exact2x2(matrix(c(6,12,12,5),2,2))
Two-sided Fisher's Exact Test (usual method using minimum likelihood)
data: matrix(c(6, 12, 12, 5), 2, 2)
p-value = 0.04371
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.0435 0.9170
sample estimates:
odds ratio
0.2189021
>
> exact2x2(matrix(c(6,12,12,5),2,2),tsmethod="central")
Central Fisher's Exact Test
data: matrix(c(6, 12, 12, 5), 2, 2)
p-value = 0.06059
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
0.03888003 1.05649145
sample estimates:
odds ratio
0.2189021
>
Let me know if there are comments or problems,
Mike
******************************************************************
Michael P. Fay, PhD
Mathematical Statistician
National Institute of Allergy and Infectious Diseases
Tel: 301-451-5124 Fax:301-480-0912
(U.S. postal mail address)
6700B Rockledge Drive MSC 7609
Bethesda, MD 20892-7609
(Overnight mail address)
6700-A Rockledge Drive, Room 5133
Bethesda, MD 20817
http://www3.niaid.nih.gov/about/organization/dcr/BRB/staff/michael.htm
**********************************************************************
[[alternative HTML version deleted]]
_______________________________________________
R-packages mailing list
r-packa...@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-packages
______________________________________________
R-help@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.
