In [24]:
lambda <- .01
In [25]:
u <- runif(1)
(x <- -log(u)/lambda)
51.7508590925093
In [26]:
x <- -log(runif(1000000))/lambda
In [27]:
hist(x)
In [29]:
mean(x)
sd(x)
99.9972939338158
99.9690048361242
In [32]:
qqplot(rexp(100000,lambda), x)
In [ ]: