аюсь найти средства, не считая NA, для нескольких столбцов с кадрами данных несколькими группами.

airquality <- data.frame(City = c("CityA", "CityA","CityA",
                                  "CityB","CityB","CityB",
                                  "CityC", "CityC"),
                         year = c("1990", "2000", "2010", "1990", 
                                  "2000", "2010", "2000", "2010"),
                         month = c("June", "July", "August",
                                   "June", "July", "August",
                                   "June", "August"),
                         PM10 = c(runif(3), rnorm(5)),
                         PM25 = c(runif(3), rnorm(5)),
                         Ozone = c(runif(3), rnorm(5)),
                         CO2 = c(runif(3), rnorm(5)))
airquality

Итак, я получаю список имен с номером, чтобы я знал, какие столбцы выбрать:

nam<-names(airquality)
namelist <- data.frame(matrix(t(nam)));namelist

Я хочу рассчитать среднее значение по городу и году для PM25, озона и CO2. Это означает, что мне нужны столбцы 1,2,4,6: 7)

acast(datadf, year ~ city, mean, na.rm=TRUE)

Но это не совсем то, чего я хочу, потому что оно включает в себя значение того, что мне не нужно, и не в формате фрейма данных. Я мог бы преобразовать это и затем уронить, но это кажется очень неэффективным способом сделать это.

Есть ли способ лучше?