fisher.test() gives a very small p-value, which is underflow on my machine. However, the log of it should not be underflow. I'm wondering if there is a way get log() of a small p-value. An approximation is acceptable in this case. Thank you! > fisher.test(rbind(c(10000,100000),c(100000,10000000)))$p.value [1] 0