如何在 R 中创建具有概率分布的矩阵

How to create a matrix with probability distribution in R

我想在 R 中创建一个矩阵,其中元素 [-1,0,1] 的概率分别为 [1/6, 2/3, 1/6]。概率可能会在运行时发生变化。对于静态概率,我得到了输出,但问题是概率的动态变化。 例如,如果我使用 [sqrt(3),0,-sqrt(3)] 创建上述概率的矩阵,则所需的输出是。

注意:概率不应如前所述是静态的。它在运行时可能会有所不同。

请帮忙解决这个问题。

假设您想要一个 2x3 矩阵:

matrix(sample(c(-1,0,1), size=6, replace=TRUE, prob=c(1/6,2/3,1/6)), nrow=2)

因此,您可以根据 prob 中定义的概率从您想要的值中抽样。这只是一个向量,但您可以在之后使用 matrix 将其变成所需形状的矩阵。用变量而不是值替换概率,以免使其成为静态的。

如果数字应该根据某种方案分配而不是根据概率随机抽取,请复制向量元素并打乱它们:

 matrix(sample(rep(c(-1,0,1), times=c(1,4,1))), nrow=2)

您可以试试这个来生成 mxn 矩阵:

sample.dynamic.matrix <- function(pop.symbols, probs, m, n) {
   samples <- sample(pop.symbols, m*n, prob = probs, replace=TRUE)
   return(matrix(samples, nrow=m))
}

set.seed(123)
sample.dynamic.matrix(-1:1, c(1/6,2/3,1/6), 2, 3)
#     [,1] [,2] [,3]
#[1,]    0    0   -1
#[2,]    1   -1    0