Included a function colMode(x) to compute mode of each column in a matrix x