library(readr) library(AER) Titanic_2 <- read_csv("https://stanford.io/2O9RUCF")[,-3] colnames(Titanic_2) = c("Survived", "Class", "Sex", "Age", "Siblings", "Parents", "Fare") Titanic_2$Class <- as.factor(Titanic_2$Class) surv_mod <- lm(Survived ~ Class, data = Titanic_2) # calcular las estimaciones específicas de la clase # ajuste el LMP aumentado y asígnelo a `LPM_mod` # obtener un resumen robusto de los coeficientes del modelo # calcular las estimaciones específicas de la clase surv_prob_c1 <- surv_mod$coefficients[1] surv_prob_c2 <- surv_prob_c1 + surv_mod$coefficients[2] surv_prob_c3 <- surv_prob_c1 + surv_mod$coefficients[3] # ajuste el LMP aumentado y asígnelo a `LPM_mod` LPM_mod <- lm(Survived ~ ., data = Titanic_2) # obtener un resumen robusto de los coeficientes del modelo coeftest(LPM_mod, vcovHC) test_predefined_objects("Titanic_2") test_predefined_objects("surv_mod") test_object("surv_prob_c1") test_object("surv_prob_c2") test_object("surv_prob_c3") test_function("coeftest", args = c("x", "vcov.")) success_msg("Lindo. El LMP aumentado predice que la probabilidad de supervivencia aumenta con la clase de pasajeros. Los coeficientes estimados en las variables ficticias de la clase de pasajeros cambiaron solo marginalmente y siguen siendo muy significativos. Se concluye que no hay un OVB sustancial en el LPM simple del Ejercicio 5.")