如何运行 for loop R 程序更快？

Question

我正在使用以下 r code 计算每个 i = 1,2,...,200 左侧和右侧的对数似然。

但我想对大量生成的数据集执行此过程，例如 a = 10000 并将整个循环迭代 1000 次。我怎样才能加快以下程序的速度？我可以使用 apply 函数代替 for 函数吗？

提前致谢！

n1 = 100
n2 = 100
a = 1000 
n= n1 + n2
# number of simulated copies of y
  sim.data = matrix(NA, nrow = n, ncol = a)
  for (i in 1:a) {
    #for(j in 1:a){
    sim.data[,i] = c(rnorm(n1, 2, 1), rnorm(n-n1, 4, 1))
    #}
  }
  dim(sim.data)


  # Compute the log-likelihood
  B = ncol(sim.data)
  loglike_profb = matrix(NA, n - 1, B)
  for (j in 1:B) {
    for (i in 1:(n - 1)) {
      loglike_profb[i, j] = -0.5*(sum(((sim.data[1:i,j]) - mean(sim.data[1:i,j]))^2) + sum(((sim.data[(i + 1):n,j]) - mean(sim.data[(i +1):n,j]))^2))
    }
  }

Answer 1

你可以把loglike_profb的计算放到一个函数里，然后用mapply

loglike_profb_func <- function(i,j){
  -0.5*(sum(((sim.data[1:i,j]) - mean(sim.data[1:i,j]))^2) + sum(((sim.data[(i + 1):n,j]) - mean(sim.data[(i +1):n,j]))^2))
}
mapply(loglike_profb_func, rep(1:(n-1),B), rep(1:B,(n-1)))

如何运行 for loop R 程序更快？

How to run for loop R program faster?

for-loop

r

apply

sapply

如何 运行 for loop R 程序更快？

How to run for loop R program faster?

for-loop

r

apply

sapply

如何运行 for loop R 程序更快？