Dans le cas , l'hypothèse de distribution est donnée par deux variables aléatoires binomiales indépendantes et . L'hypothèse nulle est l'égalité . Mais le test exact de Fisher est un test conditionnel: il repose sur la distribution conditionnelle de étant donné . Cette distribution est une distribution hypergéométrique avec un paramètre inconnu: le rapport de cotes , puis l'hypothèse nulle est .2×2X1∼Bin(n1,θ1)X2∼Bin(n2,θ2)θ1=θ2X1X1+X2ψ=θ11−θ1θ21−θ2ψ=1
Cette distribution a sa page Wikipedia .
Pour l'évaluer avec R, vous pouvez simplement utiliser la formule définissant la probabilité conditionnelle:
p1 <- 7/27
p2 <- 14/70
x1 <- 7; n1 <- 27
x2 <- 14; n2 <- 56
#
m <- x1+x2
dbinom(x1, n1, p1)*dbinom(x2, n2, p2)/sum(dbinom(0:m, n1, p1)*dbinom(m-(0:m), n2, p2))
[1] 0.1818838
Ou utilisez la dnoncenhypergeom
fonction du MCMCpack
package:
psi <- p1/(1-p1)/(p2/(1-p2)) # this is the odds ratio
MCMCpack::dnoncenhypergeom(x=x1, n1, n2, x1+x2, psi)
[1] 0.1818838