dummyinstance(distr = "rnorm", mean = 10, sd = 1)
# Make a centered (zero-mean) exponential distribution:
lambda = 4
# 10000 observations
set.seed(1234)
y <- numeric(10000)
for (i in 1:10000) y[i] <- dummyinstance(distr = "rexp", rate = lambda,
bias = -1/lambda)
mean(y)
hist(y, breaks = 50, xlim = c(-0.5, 2.5))
Run the code above in your browser using DataLab