我有一个矩阵,其中的列根据分组变量进行分组:
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
命令,但这是通过分组变量对行求和,并且没有现成的选项来计算均值。
我将不胜感激任何帮助。
我们可以用split.default
一个data.frame
由“集团”由列分割,然后做rowMeans
对list
的data.frames的
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)
本文收集自互联网,转载请注明来源。
如有侵权,请联系 [email protected] 删除。
我来说两句