有效地评估列中不同模拟的平均值

Question

我想评估从不同模拟得到的数据的mean（）和median（）。我们假设我有一个类似下面的矩阵，其中Run表示不同的运行。

我需要评估在此期间针对不同运行评估的Value 1和Value 2的平均值。

Time Run Value1  value 2
 1    0   2.5       3
 2    0    4       2.3
 3    0    2        1 
 4    0    2        4
 1    1    4        2
 2    1    2        5
 3    1    2        3.7
 1    2    4        2
 2    2   3.5       3

Answer 1

假设您想对运行0的值1取（2.5,4,2,2）的平均值，对于运行1等值1的平均值为（4,2,2），则可以做

meansVal1   <- vector()
meansVal2   <- vector()
mediansVal1 <- vector()
mediansVal2 <- vector()

for (r in 0:max(m[,2])) {

  meansVal1[r+1]   <-   mean(m[,3][m[,2] == r])
  meansVal2[r+1]   <-   mean(m[,4][m[,2] == r])
  mediansVal1[r+1] <- median(m[,3][m[,2] == r])
  mediansVal2[r+1] <- median(m[,4][m[,2] == r])

}