Использование r: Как сгенерировать экспоненциальную дистанцию ​​со скоростью 1, используя только функцию runif

Question

Используя R, как я смогу сгенерировать экспоненциально распределенную случайную переменную со скоростью 1, только сгенерировав равномерно распределенную случайную переменную (функция runif в r).

Answer 1

Найдено в https://stephens999.github.io/fiveMinuteStats/inverse_transform_sampling.html:

# inverse transfrom sampling
num.samples <-  1000
U           <-  runif(num.samples)
X           <- -log(1-U)/2

# plot
hist(X, freq=F, xlab='X', main='Generating Exponential R.V.')
curve(dexp(x, rate=2) , 0, 3, lwd=2, xlab = "", ylab = "", add = T)