在函数中编写一个 for 循环以获取矩阵 [关闭]

Question

我需要将delta.vec 和sigma.vec 值通过我的required.replicates 函数并将它们存储在我的practice1 矩阵中。

但我得到NULL。

sigma.vec <- c(2,4,6,8,10,12)
delta.vec <- c(1,2,5,8,10)

practice1 <- matrix(0, nrow=length(delta.vec), ncol=length(sigma.vec))

required.replicates <- function(delta, sigma, z.alpha = 1.959964, z.beta=0.8416212) {
  for(i in 1:length(delta.vec)) {
    for(j in 1:length(sigma.vec))
      practice1[i,j] <- ceiling((2*(z.alpha + z.beta)^2)* (sigma[j]/delta[i])^2)
    }
  }

practice1 <- required.replicates(delta=delta.vec, sigma=sigma.vec)
practice1

Answer 1

这样更有效率：

required.replicates <- function (delta, sigma, z.alpha, z.beta) {
 oo <- 1 / outer(delta, sigma, "/")
 ceiling(oo ^ 2 * 2 * (z.alpha + z.beta) ^ 2)
 }

practice1 <- required.replicates(delta.vec, sigma.vec, 1.959964, 0.8416212)

---

修复您的原始代码

required.replicates <- function(delta, sigma, z.alpha = 1.959964, z.beta=0.8416212) { 
 oo <- matrix(0, nrow=length(delta), ncol=length(sigma))
 for(i in 1:length(delta)) 
  for(j in 1:length(sigma))
   oo[i,j] <- ceiling((2*(z.alpha + z.beta)^2)* (sigma[j]/delta[i])^2)
 return(oo)
 }

practice1 <- required.replicates(delta.vec, sigma.vec, 1.959964, 0.8416212)

谢谢！还有一个问题，如果我希望矩阵中小于 3 的任何值的值为 3，并且任何大于 1000 的值都返回为 NA，我应该添加什么？

practice1[practice1 < 3] <- 3
practice1[practice1 > 1000] <- NA
practice1