将`cbind（）`作为函数参数传递

Question

假设我有一个很好的小数据框

df <- data.frame(x=seq(1,5),y=seq(5,1),z=c(1,2,3,2,1),a=c(1,1,1,2,2))
df
##  x y z a
## 1 1 5 1 1
## 2 2 4 2 1
## 3 3 3 3 1
## 4 4 2 2 2
## 5 5 1 1 2

我希望aggregate成为其中的一部分：

aggregate(cbind(x,z)~a,FUN=sum,data=df)
##  a x z
## 1 1 6 6
## 2 2 9 3

如何使其成为程序化的？我想通过：

1. 要汇总的变量列表cbind(x,z)
2. 分组变量a（我将在程序的其他几个部分使用它，因此传递整个事件cbind(x,z)~a没有帮助）
3. 事情发生的环境

我的出发点是

blah <- function(varlist,groupvar,df) {
# I kinda like to see what I am doing here 
cat(paste0(deparse(substitute(varlist)),"~",deparse(substitute(groupvar))),"\n")
cat(is.data.frame(df),"\n")
cat(dim(df),"\n")
# but I really need to aggregate this
return( aggregate(eval(deparse(substitute(varlist))~deparse(substitute(groupvar)),df),
FUN=sum,data=df) )
}

它中途运作：

blah(cbind(x,z),a,df)
## [1] "cbind(x, z)~a"
## TRUE 
## 5 4
## Error in FUN(X[[i]], ...) : invalid 'type' (character) of argument 

所以我能够构建我需要的公式的字符表示，但是将它放入aggregate()会失败。