sample_mode function

Compute a sample mode