如何生成一個隨機(jī)矩陣？你可以看到一種叫做狄里克萊分布的東西或者你可以用貝塔分布或伽瑪分布來生成狄里克萊分布]要生成m組冪，比如@charm，每組N個隨機(jī)數(shù)的和是s一般來說，用來生成m組隨機(jī)數(shù)，每組N個

如何生成一個隨機(jī)矩陣？

你可以看到一種叫做狄里克萊分布的東西

或者你可以用貝塔分布或伽瑪分布來生成狄里克萊分布

]要生成m組冪，比如@charm，每組N個隨機(jī)數(shù)的和是s

一般來說，用

來生成m組隨機(jī)數(shù)，每組N個，和是s

這里q是一個參數(shù)，數(shù)據(jù)越大，平均值越高（生成的每個數(shù)據(jù)都更接近s/N）

q越小，生成的數(shù)據(jù)越不均勻（一個數(shù)據(jù)很大，另一個數(shù)據(jù)很?。?/p>

更特別的是，使用

生成m組隨機(jī)數(shù)每組的隨機(jī)數(shù)等于alpha的長度

第i個隨機(jī)數(shù)的平均值是alpha[i]/sum（alpha）*s

這是gtools:：rdirichlet函數(shù)的算法

你可以看到，它實際上是生成gamma分布隨機(jī)數(shù)并對它們進(jìn)行歸一化

當(dāng)然，如果要生成的每個隨機(jī)數(shù)都是整數(shù)，可以用rmultinom（n，size，prob）生成

這里n是生成這類隨機(jī)數(shù)的幾組

size是每組隨機(jī)數(shù)的和

prob是一個向量，第i個元素表示每組第i個元素在總數(shù)中的比例