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


      # factor de corrección
      n <- 
      k <- 
      CF <- 

      # obtener tanto R^2 como el R^2 ajustado


      # verifique que el R^2 ajustado se pueda calcular como se indicó anteriormente


      # factor de corrección
      n <- nrow(bh_mult_mod$model)
      k <- ncol(bh_mult_mod$model)-1
      CF <- (n-1)/(n-k-1)

      # obtener tanto R^2 como el R^2 ajustado
      summary(bh_mult_mod)$r.squared
      summary(bh_mult_mod)$adj.r.squared

      # verifique que el R^2 ajustado se pueda calcular como se indicó anteriormente
      1 - (1 - summary(bh_mult_mod)$r.squared) * CF == summary(bh_mult_mod)$adj.r.squared


      test_predefined_objects("bh_mult_mod")
      test_object("n")
      test_object("k")
      test_object("CF")
      test_output_contains("summary(bh_mult_mod)$r.squared")
      test_output_contains("summary(bh_mult_mod)$adj.r.squared")
      ex() %>% check_operator("==") %>% check_result() %>% check_equal()
      success_msg("Well done!")