Y_tilde <- function(x){sum(x)/(length(x)-1)}
# compute repeatedly estimates and store the results in est_consistent
set.seed(123)
# plot a histogram of est_biased
# add a red vertical line at mu = 10
# compute repeatedly estimates and store the results in est_consistent
set.seed(123)
est_consistent <- replicate(n = 10000, expr = Y_tilde(rnorm(1000, 10, 5)))
# plot a histogram of est_consistent
hist(est_consistent)
# add a red vertical line at mu = 10
abline(v = 10, col = "red")
test_object("est_consistent")
test_function("hist", args = "x")
test_function("abline", args = c("v", "col"))