we still have the dangling issue of a possible "mode" function for
inclusion in the 'extra' package.
Here's the latest suggestion from the stack-replacement thread (another
function matrix mode(matrix v)
# v should be a vector
E = ecdf(vec(v))
howmuch = diff(0 | E[,2])[2:] # make sure the 1st is also diffed
where = imaxc(howmuch)
return E[where, 1] | howmuch[where]
It was noted that this is not optimal for multi-modal input, but "what
should be done is debatable". And I mentioned that "this is partly
inherited from the imax*/imin* suite of gretl functions", which also
return _something_ when in fact it is ambiguous.
My proposal is to simply document that the function only returns one
mode, not all (if several exist). But I'm happy to hear improvements.
Show replies by date