回复：从矩阵中随机选择值，添加 for 循环（在 R 中）

Question

可能的重复：
从现有矩阵中随机选择值添加向量后的矩阵（在 R 中）

这是我上周问题的后续，可以找到此处。我不确定将这个问题发布在同一个地方是否合适，或者将其作为一个新问题发布。

好的，上次我询问如何在将新向量绑定到矩阵后从矩阵中随机删除值。答案非常有用，但是当我使用非方阵时我发现了一个错误。我一直在循环中运行代码并每次对矩阵求和以确保其正常工作，但我发现总和有所不同，这意味着代码有时会在矩阵中选择错误的值（我希望它只选择和替换）。

这是代码：

mat1<-matrix(c(1,0,1,0, 0,1,1,1, 1,0,0,0, 1,0,0,1, 1,1,1,1, 0,0,0,1),byrow=F, nrow=4)
I.vec<-c(0,1,1,1,0,0)


foo <- function(mat, vec) {
nr <- nrow(mat)
nc <- ncol(mat)
cols <- which(vec == 1L)
rows <- sapply(seq_along(cols), 
    function(x, mat, cols) {
        ones <- which(mat[,cols[x]] == 1L)
        sample(ones, 1)
        }, mat = mat, cols = cols)
ind <- (nr*(cols-1)) + rows
mat[ind] <- 0
mat <- rbind(mat, vec)
rownames(mat) <- NULL
mat
}

set.seed(2)

for (j in 1:1000){                             #run this vector through the simulations
     I.vec2=sample(I.vec,replace=FALSE)       #randomize interactions
     temp=foo(mat1,I.vec2)                    #run foo function
     prop=sum(temp)
     print.table(prop)
     }

在这种情况下，有时矩阵的总和是 13，有时是 14，而它应该始终是 = sum(mat1) = 13。

我试图拆开代码，我认为一切除了 rows 函数之外，它工作正常，不可否认，我并不完全理解。

Answer 1

该问题是 sample() 的一个功能。我将更新原始 Q，但问题是由于在候选矩阵的列中观察到单个 1 造成的。生成 rows 的代码试图从一组 1 中进行采样。不幸的是，我忘记了 sample() 有一个功能，当第一个参数是长度向量时1、sample() 将其视为您想要从集合 1, ..., n 中采样，其中 n< /code> 是您真正想要从中采样的集合中单个元素的值。

一个简单的示例说明了当 sample() 的参数 x 是长度为 1 的向量时的行为：

> set.seed(1)
> replicate(10, sample(4, 1))
 [1] 2 2 3 4 1 4 4 3 3 1

本能地，这些都应该是 4，但是它们不是因为这个已记录且众所周知的功能。