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

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

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

或者你可以用貝塔分布或伽瑪分布來(lái)生成狄里克萊分布

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

一般來(lái)說(shuō)，用

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

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

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

更特別的是，使用

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

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

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

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

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

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

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

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