求助：R语言如何多个矩阵的对应元素的中位数？

如果有多个矩阵，矩阵的元素表示相同类型的特征，希望找到这些元素的中位数。
例如，有M1，M2，M3和M4共四个矩阵，维度为3 X 5，希望得到矩阵M，它也是3X5，它的值是M1-M4中各相同位置取值的中位数。

可以用循环或者apply来搞

循环：

1. > set.seed(1)
   
2. > m1 <- matrix(rnorm(15), 3, 5)
   
3. > m2 <- matrix(rnorm(15), 3, 5)
   
4. > m3 <- matrix(rnorm(15), 3, 5)
   
5. > m4 <- matrix(rnorm(15), 3, 5)
   
6. >
   
7. > out1 <- matrix(NA, nrow=nrow(m1), ncol=ncol(m1))
   
8. > for (i in 1:nrow(m1)) {
   
9. +   for (j in 1:ncol(m1)) {
   
10. +     out1[i, j] <- median(c(m1[i, j], m2[i, j], m3[i, j], m4[i, j]))
    
11. +   }
    
12. + }
    
13. > out1
    
14.             [,1]         [,2]        [,3]       [,4]        [,5]
    
15. [1,] -0.33569371  0.383708077  0.04656955  0.6915007 -0.83268760
    
16. [2,]  0.08372653  0.461704547  0.20784234  0.7278262  0.03925657
    
17. [3,]  0.57810227 -0.008444341 -0.27679087 -0.2045786  0.14144348

复制代码

用apply函数：

1. > all.mat <- array(c(m1, m2, m3, m4), c(3, 5, 4))
   
2. > out2 <- apply(all.mat, 1:2, median)
   
3. > out2
   
4.             [,1]         [,2]        [,3]       [,4]        [,5]
   
5. [1,] -0.33569371  0.383708077  0.04656955  0.6915007 -0.83268760
   
6. [2,]  0.08372653  0.461704547  0.20784234  0.7278262  0.03925657
   
7. [3,]  0.57810227 -0.008444341 -0.27679087 -0.2045786  0.14144348

复制代码

一样的结果