Question

我在R中有一个数据框，其中列为a-g，其中cols a和b为非数字，其余为数字。

当我在控制台中运行以下行时，它按预期工作 - 给我标准偏差，n和每个变量的平均值：

df %>% 
select(a, b, c, d, e) %>%
aggregate(.~a+b, data = ., FUN = function(x) c(avg = mean(x), std = sd(x, na.rm = TRUE), n = length(x)))

但是，当我尝试将输出分配给数据帧时，它只运行均值函数，并且不会为标准偏差或n创建列。为什么会这样？

Answer 1

当我们使用dplyr时，group_by和summarise/mutate可以获得预期的输出

library(dplyr)
df %>% 
   select(a, b, c, d, e) %>%
   group_by(a, b) %>%
   mutate(n = n()) %>%
   group_by(n, add = TRUE) %>%
   summarise_all(funs(mean, sd))

关于为什么aggregate表现不同，我们连接两个或多个函数的输出，它返回一个列为'c'，'d'和'e'的matrix输出。

str(res)
#'data.frame':   5 obs. of  5 variables:
# $ a: Factor w/ 3 levels "A","B","C": 1 3 1 2 3
# $ b: Factor w/ 2 levels "a","b": 1 1 2 2 2
# $ c: num [1:5, 1:3] -0.495 0.131 0.448 -0.495 -0.3 ...
#  ..- attr(*, "dimnames")=List of 2
#  .. ..$ : NULL
#  .. ..$ : chr  "avg" "std" "n"
# $ d: num [1:5, 1:3] -0.713 1.868 -0.71 -0.508 -0.545 ...
#  ..- attr(*, "dimnames")=List of 2
#  .. ..$ : NULL
#  .. ..$ : chr  "avg" "std" "n"
# $ e: num [1:5, 1:3] -0.893 -0.546 -0.421 1.572 -0.867 ...
#  ..- attr(*, "dimnames")=List of 2
#  .. ..$ : NULL
#  .. ..$ : chr  "avg" "std" "n"

其中res是OP代码的输出

要将其转换为普通data.frame列，请使用

res1 <- do.call(data.frame, res)
str(res1)
#'data.frame':   5 obs. of  11 variables:
# $ a    : Factor w/ 3 levels "A","B","C": 1 3 1 2 3
# $ b    : Factor w/ 2 levels "a","b": 1 1 2 2 2
# $ c.avg: num  -0.495 0.131 0.448 -0.495 -0.3
# $ c.std: num  0.233 NA NA 1.589 1.116
# $ c.n  : num  2 1 1 3 2
# $ d.avg: num  -0.713 1.868 -0.71 -0.508 -0.545
# $ d.std: num  1.365 NA NA 0.727 0.322
# $ d.n  : num  2 1 1 3 2
# $ e.avg: num  -0.893 -0.546 -0.421 1.572 -0.867
# $ e.std: num  0.771 NA NA 1.371 0.255
# $ e.n  : num  2 1 1 3 2

数据

set.seed(24)
df <- data.frame(a = rep(LETTERS[1:3], each = 3), 
   b = sample(letters[1:2], 9, replace = TRUE), 
   c = rnorm(9), d = rnorm(9), e = rnorm(9))

R中的聚合函数未正确分配给df

1 个答案:

数据