创建概率取决于索引列的矩阵

这里有一个选项,使用 sapply

n_col <- 10
n_row <- 10
mat <- matrix(nrow = n_row,
              ncol = n_col)

set.seed(1)
sapply(1:n_col, function(x) {
  mat[, x] <- sample(x = c(1, 0),
                     size = n_row,
                     replace = TRUE,
                     prob = c(1/x, 1 - 1/x))
})
#      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
# [1,]    1    0    1    0    1    0    1    0    0     0
# [2,]    1    0    0    0    0    1    0    0    0     0
# [3,]    1    1    0    0    0    0    0    0    0     0
# [4,]    1    0    0    0    0    0    0    0    0     0
# [5,]    1    1    0    1    0    0    0    0    0     0
# [6,]    1    0    0    0    0    0    0    1    0     0
# [7,]    1    1    0    1    0    0    0    0    0     0
# [8,]    1    1    0    0    0    0    0    0    0     0
# [9,]    1    0    1    0    0    0    0    0    0     0
#[10,]    1    1    0    0    0    0    1    1    0     0

希望有帮助。