A better way to do this is to use scipy.stats.mode which accepts an axis argument.

def find_mode(arr, axis):
    m = stats.mode(arr, axis=axis)
    return m[0]