Em um quadro de dados, existe uma maneira fácil de agregar sum, mean, max et c) múltiplas variáveis simultaneament

Abaixo estão alguns dados de exemplo:

library(lubridate)
days = 365*2
date = seq(as.Date("2000-01-01"), length = days, by = "day")
year = year(date)
month = month(date)
x1 = cumsum(rnorm(days, 0.05)) 
x2 = cumsum(rnorm(days, 0.05))
df1 = data.frame(date, year, month, x1, x2)

Gostaria de agregar simultaneamente ox1 ex2ariáveis dodf2 quadro de dados por ano e mês. O código a seguir agrega ox1, mas também é possível agregar simultaneamente ox2 variável

### aggregate variables by year month
df2=aggregate(x1 ~ year+month, data=df1, sum, na.rm=TRUE)
head(df2)

Qualquer sugestão seria muito apreciad