Means multiple columns by multiple groups [duplicate]

Answer 1

You should try dplyr::mutate_at :

library(dplyr)
airquality %>%
  group_by(City, year) %>%
  summarise_at(.vars = c("PM10", "PM25", "Ozone", "CO2"), .funs = mean)

# A tibble: 8 x 6
# Groups:   City [?]
    City   year         PM10       PM25      Ozone         CO2
                              
1  CityA   1990  0.004087379  0.5146409 0.44393422  0.61196671
2  CityA   2000  0.039414194  0.8865582 0.06754322  0.69870187
3  CityA   2010  0.116901563  0.6608619 0.51499227  0.32952099
4  CityB   1990 -1.535888778 -0.9601897 1.17183649  0.08380664
5  CityB   2000  0.226046487  0.4037230 0.86554997 -0.05698204
6  CityB   2010 -0.824719956  0.1508471 0.32089806 -0.12871853
7  CityC   2000 -0.824509111 -0.6928741 0.85553837  0.12137923
8  CityC   2010 -1.626150294  1.5176198 0.21183149 -0.63859910