Desde un marco de datos, ¿hay una manera fácil de agregar sum, mean, max et c) múltiples variables simultáneamente?

continuación se muestran algunos datos de muestra:

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)

Me gustaría agregar simultáneamente elx1 yx2 variables de ladf2 marco de datos por año y mes. El siguiente código agrega elx1 variable, pero también es posible agregar simultáneamente lax2 ¿variable

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

Cualquier sugerencia sería muy apreciada.