Джулия: генерирует нормально распределенное случайное число с ограниченным диапазоном

Question

Вопрос: Как я могу сгенерировать случайное число в интервале [0,1] из гауссовского распределения в Юлии?

Я собираю randn - это способ генерирования нормально распределенных случайных чисел, но описание того, как задать диапазон, в документации довольно непрозрачно.

Answer 1

Используйте пакет Distributions .Если у вас его еще нет:

using Pkg ; Pkg.add("Distributions")

, тогда:

using Distributions
mu = 0    #The mean of the truncated Normal
sigma = 1 #The standard deviation of the truncated Normal
lb = 0    #The truncation lower bound
ub = 1    #The truncation upper bound
d = Truncated(Normal(mu, sigma), lb, ub)  #Construct the distribution type
x = rand(d, 100) #Simulate 100 obs from the truncated Normal

или все в одной строке:

x = rand(Truncated(Normal(0, 1), 0, 1), 100)