R 我如何循环代码，使其具有4个直方图而不是1个直方图？ N

我们可以使用

lappy

在

N

值上循环并应用函数。因为只有一个参数，所以我们不需要任何lambda函数。它将返回一个

列表

N <- 10000 ; maxbound = 20 

one_side_normal <- function(N) {
  f <- numeric(N);
  for (i in 1:N) {
    X = rexp(1,rate=1);
    uniform = runif(1);
    while (uniform>exp((-X^2)+X) ){
      X = rexp(1,rate=1);
      uniform=runif(1);
    }
    f[i] <- X
  }
  return(f)
}

lst1 <- lapply(c(10, 100, 1000, 10000), one_side_normal)

nvec <- c(10, 100, 1000, 10000)
lst1 <- vector('list', length(nvec))
for(i in seq_along(nvec)) lst1[[i]] <- one_side_normal(nvec[[i]])