I am using the code from the solution in this post:
How to pass multiple function in raster::extract()
functions_summ <- function(x, na.rm) c(mean = mea
