我有一个矩阵,其中的列根据分组变量进行分组:
x <- matrix(sample(1:10,20,T),ncol=4)
[,1] [,2] [,3] [,4]
[1,] 7 8 5 3
[2,] 7 5 4 7
[3,] 7 1 9 3
[4,] 4 8 8 8
[5,] 9 9 1 5
group <- sample(1:2, 4, TRUE)
[1] 1 2 1 2
计算按列分组变量分组的每一行的平均值的最优雅方法是什么?所示示例的结果应为5 x 2矩阵:
1 2
[1,] 6.0 5.5
[2,] 5.5 6.0
[3,] 8.0 2.0
[4,] 6.0 8.0
[5,] 5.0 7.0
我查看了rowsum
命令,但是该命令通过分组变量对行求和,并且没有现成的选项来计算均值。
非常感谢您的帮助。
答案 0 :(得分:3)
我们可以按“组”在split.default
上使用data.frame
按列进行拆分,然后对数据帧rowMeans
进行list
>
sapply(split.default(as.data.frame(x), group), rowMeans)
# 1 2
#[1,] 6.0 5.5
#[2,] 5.5 6.0
#[3,] 8.0 2.0
#[4,] 6.0 8.0
#[5,] 5.0 7.0
x <- structure(c(7, 7, 7, 4, 9, 8, 5, 1, 8, 9, 5, 4, 9, 8, 1, 3, 7,
3, 8, 5), .Dim = 5:4)
group <- c(1, 2, 1, 2)