답변:
이 배포판에는 Wikipedia 페이지가 있습니다.
R로 평가하려면 조건부 확률을 정의하는 공식을 사용하면됩니다.
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
또는 패키지 dnoncenhypergeom
기능을 사용하십시오 MCMCpack
.
psi <- p1/(1-p1)/(p2/(1-p2)) # this is the odds ratio
MCMCpack::dnoncenhypergeom(x=x1, n1, n2, x1+x2, psi)
[1] 0.1818838