使用Boost生成基於Beta分佈的隨機數

Question

我正在嘗試使用Boost來根據使用C++的Beta分佈生成隨機數。我在網上看到很多例子，根據random.hpp中的分佈生成隨機數（例如this book）。但是，我無法將它們翻譯成使用beta.hpp中的beta版本。

謝謝。

Answer 1

您首先需要從範圍（0,1）中統一繪製一個隨機數。給定任何分佈，然後可以將該數字插入分佈的「分位數函數」中，結果就好像從分佈中抽取一個隨機值。從here：{ - 1}（Y）G（Y）= F ^：

的一般方法，從具有不跳躍一個CDF是使用反函數的CDF任意分佈產生的隨機數。如果u（1），...，u（n）是（0,1）分佈上均勻隨機數，則G（u（1）），...，G（u（n））是一個隨機數用cdf F（x）分佈的樣本。

那麼我們如何獲得beta分佈的分位數函數呢？ beta.hpp的文檔是here。你應該可以使用這樣的東西：

#include <boost/math/distributions.hpp> 
using namespace boost::math; 

double alpha, beta, randFromUnif; 
//parameters and the random value on (0,1) you drew 

beta_distribution<> dist(alpha, beta); 
double randFromDist = quantile(dist, randFromUnif);