
      library(AER)
      library(MASS)
      mod <- lm(medv ~ lstat + crim + age, data = Boston)


      # construir un conjunto de confianza del 99% para los coeficientes de `crim` y` lstat`
         
      # realizar la prueba F correspondiente


      # construir un conjunto de confianza del 99% para los coeficientes de `crim` y` lstat`
      confidenceEllipse(mod, which.coef = c("crim", "lstat"), levels = 0.99)
        
      # realizar la prueba F correspondiente
      linearHypothesis(mod, c("crim = 0", "lstat = 0"))


      test_or({
        fun <- ex() %>% check_function("confidenceEllipse")
        fun %>% check_arg("model") %>% check_equal()
        fun %>% check_arg("which.coef") %>% check_equal()
        fun %>% check_arg("levels") %>% check_equal()
      }, {
        fun <- ex() %>% override_solution("confidenceEllipse(mod, which.coef = c(\"lstat\", \"crim\"), levels = 0.99)") %>% check_function("confidenceEllipse")
        fun %>% check_arg("model") %>% check_equal()
        fun %>% check_arg("which.coef") %>% check_equal()
        fun %>% check_arg("levels") %>% check_equal()
      })
      test_function_result("linearHypothesis")
      success_msg("¡Correcto! Dado que (0,0) no es un elemento del conjunto de confianza del 99%, la hipótesis nula se rechaza al nivel del 1%.")