I wrote code in positional_utils.R for positional methods with k candidates. The code choked trying to make the simplices for integration, though -- errors from qhull.

Could do for simulation pretty easily I think