This function calculates a cumulative probability using simulated data and
Martin Maechler's ecdf
function. \(Fn\) is a step function with jumps \(i/n\) at observation
values, where \(i\) is the number of tied observations at that value. Missing values are ignored. For
observations \(y = (y1, y2, ..., yn)\), \(Fn\) is the fraction of observations less or equal to \(t\), i.e.,
\(Fn(t) = sum[yi <= t]/n\). This works for continuous, ordinal, or count variables.