3 个矩阵中三次测量的平均值

Question

我有 3 个矩阵，它们存储来自三次测量的值，并且想要取 3 个矩阵的平均值。

假设这三个矩阵是：

m1<-t(matrix(c("text", 1:3), ncol=2, nrow=4))
m2<-t(matrix(c("text", 1:3), ncol=2, nrow=4))
m3<-t(matrix(c("text", 1:3), ncol=2, nrow=4))

> m1
     [,1]   [,2] [,3] [,4]
[1,] "text" "1"  "2"  "3" 
[2,] "text" "1"  "2"  "3" 
> m2
     [,1]   [,2] [,3] [,4]
[1,] "text" "1"  "2"  "3" 
[2,] "text" "1"  "2"  "3" 
> m3
     [,1]   [,2] [,3] [,4]
[1,] "text" "1"  "2"  "3" 
[2,] "text" "1"  "2"  "3" 

我希望矩阵的每个位置都有这个：

mean(m1[i,j], m2[i,j], m2[i,j])

所以我尝试了 2 个 for 循环：

for(i in ncol(m1)){
  for(j in nrow(m1)){
    means[i,j]<-mean(m1[i,j], m2[i,j], m2[i,j])
  }

这显然行不通

如果返回 NA，则第一列中的文本不是问题。

有人可以帮帮我吗？ 谢谢！

Answer 1

我们可以把它放在list中，转换成数字并使用Reduce

lst <- lapply(list(m1[, -1], m2[,-1], m3[, -1]), as.numeric)
Reduce(`+`,lst)/length(lst)

如果有很多以'm'开头的矩阵，我们可以使用mget

lst <- lapply(mget(paste0("m", 1:3)), function(x) as.numeric(x[,-1]))

然后执行Reduce 步骤。