##########################################
###  超幾何検定 (hypergeometric test)
##########################################
N <- 10000
n <- 7
X <- 1500
x <- 6
sum(dhyper(x=x:X, m=n, n=N-n, k=X))


?dhyper
x <- 6
m <- 7
n <- 9993
k <- 1500
sum(dhyper(x=x:X, m=m, n=n, k=k))

##########################################
###  カイ二乗検定 (Chi-squared Test)
##########################################
N <- 10000
n <- 7
X <- 1500
x <- 6
data <- matrix(c((n-x), (N-n)-(X-x), x, (X-x)), ncol=2, byrow=T)
data
chisq.test(data)

##########################################
###  Distribution of p-values
##########################################
N <- 10000
n <- 7
X <- 1500

###  x=6  ###
x <- 6
data <- matrix(c((n-x), (N-n)-(X-x), x, (X-x)), ncol=2, byrow=T)
sum(dhyper(x=x:X, m=n, n=N-n, k=X))
chisq.test(data)$p.value
fisher.test(data)$p.value

###  x=5  ###
x <- 5
data <- matrix(c((n-x), (N-n)-(X-x), x, (X-x)), ncol=2, byrow=T)
sum(dhyper(x=x:X, m=n, n=N-n, k=X))
chisq.test(data)$p.value
fisher.test(data)$p.value