r - 複数の変数を持つ summrise_each() の複数の関数

Question

以下の作品。より良い解決策があると確信しています。

library(dplyr)
library(tidyr)

iris %>%
  group_by(Species) %>%
  summarise_each(funs(mean, median)) %>%
  gather(var, val, -Species) %>%
  separate(var, c("variable", "summary"), sep = "_") %>%
  spread(summary, val)