No caso , a suposição distributiva é dada por duas variáveis aleatórias binomiais independentes e . A hipótese nula é a igualdade . Mas o teste exato de Fisher é um teste condicional: ele depende da distribuição condicional de dado . Essa distribuição é uma distribuição hipergeométrica com um parâmetro desconhecido: a razão de chances e depois a hipótese nula é .2×2X1∼Bin(n1,θ1)X2∼Bin(n2,θ2)θ1=θ2X1X1+X2ψ=θ11−θ1θ21−θ2ψ=1
Esta distribuição tem sua página da Wikipedia .
Para avaliá-lo com R, você pode simplesmente usar a fórmula que define a probabilidade condicional:
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 use a dnoncenhypergeom
função do MCMCpack
pacote:
psi <- p1/(1-p1)/(p2/(1-p2)) # this is the odds ratio
MCMCpack::dnoncenhypergeom(x=x1, n1, n2, x1+x2, psi)
[1] 0.1818838