如何在R中添加总计以及group_by统计信息

Question

当计算使用任何统计summarise和group_by我们只能得到每类的汇总统计，而不是所有的人口（总）的值。 如何获得两者？

我正在寻找干净整洁的东西。 到目前为止，我只能想到：

bind_rows( 
  iris %>% group_by(Species) %>% summarise(
    "Mean" = mean(Sepal.Width), 
    "Median" = median(Sepal.Width), 
    "sd" = sd(Sepal.Width), 
    "p10" = quantile(Sepal.Width, probs = 0.1))
  , 
  iris %>% summarise(
    "Mean" = mean(Sepal.Width), 
    "Median" = median(Sepal.Width), 
    "sd" = sd(Sepal.Width), 
    "p10" = quantile(Sepal.Width, probs = 0.1)) %>% 
  mutate(Species = "Total")
  )

但是我想要更紧凑的东西。 特别是，我不想键入两次代码（用于摘要），每个组一次，总计一次。

Answer 1

如果解开要尝试的操作，则可以简化它：您拥有包含多种物种的iris数据，并且希望将其与所有物种的数据一起汇总。 绑定之前，无需计算这些摘要统计信息。 而是将iris与设置为Species = "Total"的iris版本绑定，然后进行分组和汇总。

library(tidyverse)

bind_rows(
  iris,
  iris %>% mutate(Species = "Total")
) %>%
  group_by(Species) %>%
  summarise(Mean = mean(Sepal.Width),
            Median = median(Sepal.Width),
            sd = sd(Sepal.Width),
            p10 = quantile(Sepal.Width, probs = 0.1))
#> # A tibble: 4 x 5
#>   Species     Mean Median    sd   p10
#>   <chr>      <dbl>  <dbl> <dbl> <dbl>
#> 1 setosa      3.43    3.4 0.379  3   
#> 2 Total       3.06    3   0.436  2.5 
#> 3 versicolor  2.77    2.8 0.314  2.3 
#> 4 virginica   2.97    3   0.322  2.59

我喜欢上面的评论中的注意事项，尽管我必须进行足够的这种计算才能在个人软件包中使用类似的速记功能。 对于诸如标准偏差之类的事情来说，这可能没有什么意义，但是我需要做很多事来增加人口统计总数等。（如果有用，该函数在此处 ）。

Answer 2

有点短，尽管与bind_rows非常相似

    q10 <- function(x){quantile(x , probs=0.1)}

    iris %>% 
      select(Species,Sepal.Width)%>%
      group_by(Species) %>% 
      summarise_all(c("mean", "sd", "q10")) %>% 
      t() %>% 

      cbind(c("total", iris %>% select(Sepal.Width) %>% summarise_all(c("mean", "sd", "q10")))) %>% 
      t()

可能更干净：

  bind_rows( 
    iris %>% 
      group_by(Species) %>%  
      select(Sepal.Width)%>%
      summarise_all(c("mean", "sd", "q10"))
    , 
    iris %>% 
      select(Sepal.Width)%>%
      summarise_all(c("mean", "sd", "q10")) %>% 
      mutate(Species = "Total")
  )