Generar un número aleatorio basado en la distribución Beta utilizando Boost

Question

Estoy tratando de usar Boost para generar números aleatorios según la distribución beta usando C++. He visto muchos ejemplos en línea para generar números aleatorios según las distribuciones en random.hpp (por ejemplo, this book). Sin embargo, no veo traducirlos para usar la distribución beta encontrada en beta.hpp.

Gracias.

Answer 1

En primer lugar, querrá dibujar un número aleatorio del rango (0,1). Dada cualquier distribución, puede conectar ese número en la "función de cuantiles" de la distribución, y el resultado es como si se extrajera un valor aleatorio de la distribución. De here:

Un método general para generar números aleatorios de una distribución arbitraria que tiene una cdf sin saltos es utilizar la función inversa a la cdf: G (y) = F^{- 1} (y) . Si u (1), ..., u (n) son números aleatorios de la distribución uniforme en (0,1), entonces G (u (1)), ..., G (u (n)) es un azar muestra de la distribución con cdf F (x).

Entonces, ¿cómo obtenemos una función cuantil para una distribución beta? La documentación para beta.hpp es here. Debería poder usar algo como esto:

#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);