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"))