Question

我正在尝试模拟一场比赛，在每轮比赛中，在13场比赛中随机选出两名球员，他们的Elo rating（相对技能水平的衡量标准）会得到更新。每个玩家都有一个真正的潜在技能等级true.rating，在比赛期间，我在每轮回合后存储他们的评分我正在尝试更新current.ratings并将每次迭代存储在ratings.matrix中。这是我的代码，用R markdown编写

true.ratings <- seq(from = -2, to = 2, length.out = 13)

current.ratings <- seq(from = -2, to = 2, length.out = 13)
ratings.matrix <- matrix(nrow = 13, ncol = 10000)


i.wins <- function(i, j, p) {
    i.wins <- rbinom(1, 1, p)
    return(i.wins)
}

update <- function (K, i, j) {
  R.i <- true.ratings[i]
  R.j <- true.ratings[j]
  probability_i.wins <- 1/(1 + exp(R.j - R.i))
  probability_j.wins <- 1 - probability_i.wins

  winner.is.i <- as.logical(i.wins(i, j, probability_i.wins))

  if (winner.is.i) {
    replace(current.ratings, i , current.ratings[i] + K*(1-probability_i.wins))
    replace(current.ratings, j,  current.ratings[j] -K*probability_j.wins)
  } 
 if (!winner.is.i) {
    current.ratings[i] <- current.ratings[i] -K*probability_i.wins
    current.ratings[j] <- current.ratings[j] + K*(1-probability_j.wins)
  }

  return(current.ratings)
}


match <- function(K) {
  ij <- sample(1:13, 2, replace = F)
  i <- ij[1]
  j <- ij[2]

  return(update(K, i, j))
}


tournament <- function(K) {

  current.ratings <- match(K)

  for (j in 1:10000) {
    ratings.matrix[, j] <- match(K)

  }
  return(ratings.matrix)
}


ratings.matrix <- tournament(K = 0.01)

评分方法在每轮结束后更新，玩家i和玩家j在播放该匹配后的当前评分。但是，当两个不同的玩家在下一轮中玩时，当前玩家i和j的等级将返回到原始true.rating值。我怀疑这与我没有正确更新我的current.rating矩阵向量有关。

这是我的第一个堆栈溢出帖子，我也是R的新手。欢迎任何其他关于我如何改进的反馈：）

Answer 1

R努力成为一种函数式语言，这意味着函数影响父环境中事物的唯一方法是通过返回值。如果要更改属于函数外部某些环境的值，则需要以R方式执行。

常用方法：

返回$ sudo passwd username并在顶层明确覆盖它。换句话说，您的函数current.ratings需要返回包含tournament()以及current.ratings的列表，然后您需要在顶级环境中明确覆盖这两个变量（在您的函数调用）。
您需要将环境本身作为参数传递给函数。这是R通过引用传递的方式。例如

ratings.matrix

请注意，在你的情况下你也可以使用厌恶的tournament <- function(K, myenv){ myenv$current.ratings <- match(K) for (j in 1:10000) { ratings.matrix[, j] <- match(K) } return(ratings.matrix) } ratings.matrix <- tournament(K = 0.01, myenv=environment())超级对齐操作符，但我没有提到它是一个实际的解决方案，因为我不想被任何腐烂的西红柿击中。

如何永久更改R中的矢量？

1 个答案: