从数据框中，是否有一种简单的方法可以同时聚合（sum、、mean等max）多个变量？

以下是一些示例数据：

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)

我想同时按年和月聚合数据框中的x1和x2变量。df2下面的代码聚合了x1变量，但是是否也可以同时聚合x2变量呢？

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