Current version works with 2D data, the 2nd dim stores samples/scenarios. To be extended towards multiple dimension, as implemented in `prob()`