Question

如果我有此列表

set.seed(123)
thelist <- list(a=data.frame(x1=rnorm(10), x2=rnorm(10)),
                b=data.frame(x1=rnorm(10), x2=rnorm(10)),
                c=data.frame(x1=rnorm(10), x2=rnorm(10)))

并希望使用以下代码计算每个列表中每列的平均值。

sapply(do.call("rbind",thelist),mean)

我怎样才能为每个列表中的每一列（a：c）计算标准偏差，因为sd没有等效函数（至少据我所知）？

任何建议都将不胜感激。

Answer 1

基本的 R 解决方案将使用sapply两次。

仅限于平均值：

t(sapply(thelist, sapply, mean))

导致

            x1        x2
a  0.074625644 0.2086220
b -0.424558873 0.3220446
c -0.008715537 0.2216860

如果你们俩想要：

my_summary <- function(x){
  c(mean = mean(x), sd = sd(x))
}

as.data.frame(lapply(thelist, sapply, my_summary))

导致：

           a.x1     a.x2       b.x1      b.x2         c.x1      c.x2
mean 0.07462564 0.208622 -0.4245589 0.3220446 -0.008715537 0.2216860
sd   0.95378405 1.038073  0.9308092 0.5273024  1.082518163 0.8564451

Answer 2

首先，我通过将名称放入一列来使其可堆叠：

for (i in seq_along(thelist)) thelist[[i]]$dfname <- names(thelist)[i]

然后，堆叠并采用data.table：

的方法

require(data.table)
DT <- rbindlist(thelist)

DT[,lapply(.SD,mean),by=dfname]

给出了

   dfname           x1        x2
1:      a  0.074625644 0.2086220
2:      b -0.424558873 0.3220446
3:      c -0.008715537 0.2216860

您可能还会考虑summary功能，尽管它很笨拙：

DT[,as.list(unlist(lapply(.SD,summary))),by=dfname]
#    dfname x1.Min. x1.1st Qu. x1.Median   x1.Mean x1.3rd Qu. x1.Max. x2.Min. x2.1st Qu. x2.Median x2.Mean x2.3rd Qu. x2.Max.
# 1:      a  -1.265    -0.5318  -0.07983  0.074630    0.37800   1.715 -1.9670   -0.32690    0.3803  0.2086     0.6505  1.7870
# 2:      b  -1.687    -1.0570  -0.67700 -0.424600    0.06054   1.254 -0.3805   -0.23680    0.4902  0.3220     0.7883  0.8951
# 3:      c  -1.265    -0.6377  -0.30540 -0.008716    0.56410   2.169 -1.5490   -0.03929    0.1699  0.2217     0.5018  1.5160

最后，复制my old answer，您可以制作自己的摘要统计功能：

summaryfun <- function(x) list(mean=mean(x),sd=sd(x))
DT[,as.list(unlist(lapply(.SD,summaryfun))),by=dfname]
#    dfname      x1.mean     x1.sd   x2.mean     x2.sd
# 1:      a  0.074625644 0.9537841 0.2086220 1.0380734
# 2:      b -0.424558873 0.9308092 0.3220446 0.5273024
# 3:      c -0.008715537 1.0825182 0.2216860 0.8564451

Answer 3

您可以按照自己的建议合并数据，然后汇总如下：

thelist_named <- Map(cbind, thelist, nam = names(thelist))
thelist_binded <- do.call(rbind, thelist_named)

聚合部分：

my_summary <- function(x){
  c(mean = mean(x), sd = sd(x))
}
aggregate(.~nam, thelist_binded, my_summary)

结果：

  nam      x1.mean        x1.sd   x2.mean     x2.sd
1   a  0.074625644  0.953784051 0.2086220 1.0380734
2   b -0.424558873  0.930809213 0.3220446 0.5273024
3   c -0.008715537  1.082518163 0.2216860 0.8564451

列表中多个数据帧的摘要统计信息

3 个答案: