r计算同名列的均方差

Question

我有一个简短的问题

我有一个数据帧，有许多测量列。我想计算具有相同（标题）名称的列的平均值。我使用下面的代码（在stackoverflow中找到）。。

如何计算数据帧中具有相同列名的列的平均值

作为示例数据...

df <- data.frame(c(1, 2, 3, 4,5),
                 c(2, 3, 4,NA,2),
                 c(3, 4, 5,3,6),
                 c(3, 7, NA,3,6))
names(df) <- c("a", "b", "a", "b")

df <- sapply(split.default(df, names(df)), rowMeans, na.rm = TRUE) 

结果是这样的。。。

a    b
2    2.5
3    5
4    4
3.5  3
5.5  4

这段代码告诉我具有相同(标题)名称的列的含义。

但是我也想要均方差。我试着用rowSds替换rowMeans，但是不起作用。

知道如何使用相同的代码来计算均方差吗？？

Answer 1

这是一个用户定义的函数，可能会有用。你可能想看看:

rowVars

Answer 2

这应该行得通：

df <- data.frame(c(1, 2, 3),
                 c(2, 3, 4),
                 c(3, 4, 5))
names(df) <- c("a", "b", "a")


sapply(split.default(df, names(df)), function(smaller_df) {
  sapply(smaller_df, function(col) c(mean(col), sd(col)))
})

第一个应用程序适用于通过拆分生成的每个 data.frame，每个应用程序将对应于具有相同名称的一组列。第二个应用适用于每列。

如果您想获得具有给定名称的列中所有测量值的均值和均方差，而不是作为单独的样本，您可以将内部sapplication更改为：

sapply(list(unlist(smaller_df)), function(col) c(mean(col), sd(col)))

Answer 3

基于你之前的方法的一个想法是做以下事情

sapply(split.default(df, names(df)), function(x) apply(x, 1, sd, na.rm=TRUE))
#              a         b
# [1,] 1.4142136 0.7071068
# [2,] 1.4142136 2.8284271
# [3,] 1.4142136        NA
# [4,] 0.7071068        NA
# [5,] 0.7071068 2.8284271

请记住，之所以返回 NA，是因为不应在大小为 1 的样本上评估 sd。