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