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

【问题标题】：Efficiently evaluate the mean of different simulations in columns有效评估列中不同模拟的平均值
【发布时间】：2014-09-16 05:42:07
【问题描述】：

我想评估来自不同模拟的数据的均值（）和中位数（）。假设我有一个像下面这样的矩阵，其中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

【问题讨论】：

请在问这样的基本问题之前表现出一些努力（使用的功能和代码）。

标签： r simulation multiple-columns mean median

【解决方案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])

}

【讨论】：