# Definir una función para el estimador
# Calcular estimaciones repetidamente y almacenar los resultados en est_biased
set.seed(123)
# Graficar un histograma de est_biased
# Agregue una línea vertical roja en mu = 10
# Definir una función para el estimador
Y_tilde <- function(x){sum(x)/(length(x)-1)}
# Calcular estimaciones repetidamente y almacenar los resultados en est_biased
set.seed(123)
est_biased <- replicate(n = 10000, expr = Y_tilde(rnorm(5, 10, 5)))
# Graficar un histograma de est_biased
hist(est_biased)
# Agregue una línea vertical roja en mu = 10
abline(v = 10, col = "red")
test_function_definition("Y_tilde", function_test = {
test_expression_result("Y_tilde(1:10)")
test_expression_result("Y_tilde(1:25)")
test_expression_result("Y_tilde(seq(0, 1, 0.1))")
})
test_object("est_biased")
test_function("hist", args = "x")
test_function("abline", args = c("v", "col"))