for循环列表中的每个矩阵

Question

我正在尝试将以下 for 循环应用于列表per_d中的每个矩阵，并创建一个名为per_hole的新列表。 我不确定该怎么做，我应该使用 lapply 吗？

非常感谢您的帮助！

per_hole <- per_d
  for (i in 1:S) {
    for (j in 1:t){
  if (per_hole [i,j] > CS) {
    per_hole [i,j] <- per_hole [i,j] - rnorm (1, mean = 1, sd = 0.5)
  } else {
    per_hole [i,j] <- per_hole [i,j] + rnorm (1, mean = 1, sd = 0.5)
  }}}

复制代码

N <- 1
CS <- 10.141 
S <- seq (7.72,13,0.807) 
t <- 15
l <- length (S)
m0 <- 100 
exps <- c(0.2, 0.5, 0.9, 1.5, 2) 
sd_per <- c(0.2, 0.5, 0.8, 1.3, 1.8) 
sd_noise <- 3 


per <- lapply(sd_per, function(x){ 
  per <- matrix(nrow = length(S)*N, ncol = t+1)
  for (i in 1:dim(per)[1]) {
    for (j in 1:t+1){
      per [,1] <- replicate (n = N, S)
      per [i,j] <- round (abs (rnorm (1, mean = per[i,1], sd =x)),digits=3)
      colnames(per) <- c('physical',paste('t', 1:15, sep = ""))
      per <- as.data.frame (per)
    }
  }
  per <- per [,-1]
  return(per)
}
)  
names(per) <- paste("per", seq_along(sd_per), sep = "")

per_d <- lapply(per, function(x){ 
  per_d <- abs (x - 10.141) 
}
)  
names(per_d) <- paste("per_d", seq_along(sd_per), sep = "")

Answer 1

你可以试试

per_hole <- lapply(per_d,function(x) x + ifelse(x>CS,-1,1)*rnorm(prod(dim(x)),1,0.5))

要么

per_hole <- lapply(per_d, function(x) x + rnorm(prod(dim(x)), 1-2*(x > CS), 0.5))