A probabilistic derivation of van der Waals' equation is given where the non-ideal contributions to the entropy coincide with the log-likelihood ratios, at their maxima, of the binomial distribution and the upper or lower bounds on the hypergeometric and Polya distributions. The hypergeometric distribution accounts for the repulsive interactions of a hard sphere gas while the Polya distribution accounts for the attractive interactions.