Estou tentando transferir meu entendimento do plyr para o dplyr, mas não consigo descobrir como agrupar por várias colunas.

# make data with weird column names that can't be hard coded
data = data.frame(
  asihckhdoydkhxiydfgfTgdsx = sample(LETTERS[1:3], 100, replace=TRUE),
  a30mvxigxkghc5cdsvxvyv0ja = sample(LETTERS[1:3], 100, replace=TRUE),
  value = rnorm(100)
)

# get the columns we want to average within
columns = names(data)[-3]

# plyr - works
ddply(data, columns, summarize, value=mean(value))

# dplyr - raises error
data %.%
  group_by(columns) %.%
  summarise(Value = mean(value))
#> Error in eval(expr, envir, enclos) : index out of bounds

O que estou perdendo para traduzir o exemplo plyr em uma sintaxe dplyr-esque?

Edit 2017: O Dplyr foi atualizado, portanto, uma solução mais simples está disponível. Veja a resposta atualmente selecionada.