r - 使用 R 中给定概率的随机条目创建矩阵

Question

我想创建一个 100x100 矩阵 A，其中条目 a_ij 从集合 {0,1} 中随机选择，其中 P(a_ij=1)=0.2 和 P(a_ij=0)=0.8。
这是我迄今为止尝试过的:

n<-100
matrix<-matrix(0,100,100)
mynumbers<-c(1,0)
myprobs<-c(0.2,0.8)
for(i in 1:100){
  for (j in 1:100){
    matrix[i,j]<-sample(mynumbers, 1, replace=TRUE, prob=myprobs)
  }
}
matrix

我不确定样本大小是否为 1，但这种方法似乎只有在我选择 size=1 时才有效……这是正确的方法吗？先感谢您!

Answer 1

正如@akrun 指出的那样，有更简单的方法。 100 x 100 的矩阵意味着 10,000 个条目。 prob = .2是说成功 = 1 = P(a_ij=1)=0.2, size在这种情况下是指一次审判。 matrix参数应该是不言而喻的。

set.seed(2020)
trials <- rbinom(n = 10000, size = 1, prob = .2)
my.matrix <- matrix(trials, nrow = 100, ncol = 100)

或者更接近您的代码

n <- 10000
mynumbers<-c(1,0)
myprobs<-c(0.2,0.8)
trials2 <- sample(x = mynumbers, 
                  size = n, 
                  replace = TRUE, 
                  prob = myprobs)
my.matrix2 <- matrix(trials2, nrow = 100, ncol = 100)