library(MASS)
A <- matrix(c(1,1,0,0,0,0,1,1,1,0,1,0,0,1,0,1), 4, 4, byrow=TRUE)
b <- c(4, 4, 4, 4)

v <- Null(t(A))

dhyper(4, 4, 4, 4)

## [1] 0.01428571

"[1] 0.01428571"

## [1] "[1] 0.01428571"

dhyper(3, 4, 4, 4)

## [1] 0.2285714

"[1] 0.2285714"

## [1] "[1] 0.2285714"

Tea <- matrix(c(3, 1, 1, 3), 2, 2, byrow = TRUE)
fisher.test(Tea, alternative = "greater")

## 
##  Fisher's Exact Test for Count Data
## 
## data:  Tea
## p-value = 0.2429
## alternative hypothesis: true odds ratio is greater than 1
## 95 percent confidence interval:
##  0.3135693       Inf
## sample estimates:
## odds ratio 
##   6.408309