r - 在R中生成一个随机矩阵，每行之和等于1

Question

我想在R软件中生成一个矩阵(4行30列)，其中元素随机，元素范围在0到1之间，每行之和等于1。

Answer 1

这是一个基于 softmax(多项式 logit)变换的解决方案。

m <- matrix(rnorm(4 * 30), nrow=30)
prob <- exp(m)/rowSums(exp(m))

rowSums(prob)
#[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

all(prob > 0 & prob < 1)
#[1] TRUE