Add mfv() function from modeest package for calculating mode