Working with column names as variables in formula
New argument 'pairs' and now the function returns units and ele_names. But should be tested.