Added optimal_sampling() to estimate variance in phi as a function of sample allocation (from Smouse et al. 2001)