有一种方法可以将函数f
应用于矩阵的每一列:
M <- matrix(seq(1,16), 4, 4)
apply(M, 2, mean)
#[1] 2.5 6.5 10.5 14.5
但是如果我想建立一个关于矩阵的描述性统计数据,我应该使用更多的indeces。例如,max, min, mean
等。
但是R不允许做这样的事情:
apply(M, 2, c(mean, max))
获得此输出:
# [,1] [,2] [,3] [,4]
#mean 2.5 6.5 10.5 14.5
#max 4 8 12 16
你能告诉我如何处理这个问题吗?
答案 0 :(得分:5)
apply(M, 2, function(x) c(mean(x), max(x)))
# [,1] [,2] [,3] [,4]
# [1,] 2.5 6.5 10.5 14.5
# [2,] 4.0 8.0 12.0 16.0
答案 1 :(得分:2)
尝试以下方法:
f <- c("max", "min", "mean")
sapply(f, function(x) apply(M, 2, x))
max min mean
[1,] 4 1 2.5
[2,] 8 5 6.5
[3,] 12 9 10.5
[4,] 16 13 14.5