I'd need to know what you expect to get from such function. I'm going to assume you need to compute a level-wise summary of the data using a list of functions (here I use the mean and SD as an example).
You can try with the combination of group_by() and summarise() from the package dplyr. I don't have the code or data you are using but it would look like this:
data %>%
group_by(categorical_variable) %>% # group observations from the same level of
categorical_variable
summarise(mn1 = mean(numeric_variable1), # compute the mean of numeric_variable1 for each level of categorical_variable
std1 = sd(numeric_variable1),
mn2 = mean(numeric_variable2),
std2 = sd(numeric_variabl2))
You can also group observations from combinations of levens from two categorical variables using group_by(categorical_variable1,categorical_variable2).
Hope this addresses you problem! 
NOTE: Edited to indet the code properly.