moda function

Computes the mode