new function for frequency calculation based on data.table: sffcNA (very fast for big data sets)