将嵌套循环的结果写入 R 中的另一个向量

Question

我对 R 还很陌生，并且在使用它方面遇到了一些困难。我有以下代码：

repeat {
    if (t > 1000) 
        break
    else {
        y1 <- rpois(50, 15)
        y2 <- rpois(50, 15)
        y <- c(y1, y2)
        p_0y <- matrix(nrow = max(y) - min(y), ncol = 1)
        i = min(y)
        while (i <= max(y)) {
            p_0y[i - min(y), ] = (length(which(y1 == i))/50)
            i <- i + 1
        }
        p_y <- matrix(nrow = max(y) - min(y), ncol = 1)
        j = min(y)
        while (j <= max(y)) {
            p_y[j - min(y), ] = (length(which(y == j))/100)
            j <- j + 1
        }
        p_0yx <- p_0y[rowSums(p_0y == 0) == 0]
        p_yx <- p_y[rowSums(p_0y == 0) == 0]
        g = 0
        logvect <- matrix(nrow = (length(p_yx)), ncol = 1)
        while (g <= (length(p_yx))) {
            logvect[g, ] = (p_0yx[g])/(p_yx[g])
            g <- g + 1
        }
        p_0yx %*% (log2(logvect))
        print(p_0yx %*% (log2(logvect)))
        t <- t + 1
    }
}

我对最后一行之前的所有内容都很满意，但我不想将 p_0yx%*%(log2(logvect)) 的值打印到屏幕上，而是想将其存储为另一个向量。有什么想法吗？我尝试过以与嵌套循环类似的方式进行操作，但似乎不起作用。

谢谢

Answer 1

简短的答案是首先声明一个变量。将其放在您在此发布的所有内容之前。我将其称为临时。它将保存所有值。

temp <- numeric(1000)

然后，

temp[t] <- p_0yx %*% log2(logvect)

您的代码正在做一些奇怪的事情，而不是使用打印行。查看 p_0y 的第一个索引。它实际上是该矩阵中项目 0, 的索引。 R 从 1 开始索引。当您在该矩阵中创建行数时，您可以使用 max(y) - min(y)。如果最大值为 10，最小值为 1，则只有 9 行。我打赌你真的想添加一个。另外，您的代码与 R 非常不一样，带有所有不必要的 while 循环。例如，您的整个最后一个循环（以及 logvect 的初始化）可以替换为：

logvect = (p_0yx)/(p_yx)

但是回到错误..以及更多 Rness...下面的代码...可以

p_0y <- matrix(nrow = max(y) - min(y), ncol = 1)
i = min(y)
while (i <= max(y)) {
    p_0y[i - min(y), ] = (length(which(y1 == i))/50)
    i <- i + 1
    }

更正确地替换为吗？

p_0y <- numeric(max(y) - min(y) + 1)
p_0y[sort(unique(y1)) - min(y1) + 1] = table(y1)/50
p_0y <- matrix(p_0y, ncol = 1)

（对其余代码进行类似的重新思考也可以消除其余循环）