cs <- c(23, 19, 30, 22, 23, 29, 35, 36, 33, 25) ts <- c(430, 430, 333, 410, 390, 377, 325, 310, 328, 375) mod <- lm(ts ~ cs) # obtener la matriz y guardarla en `cov_matrix` # calcular los errores estándar y asignarlos al vector `SEs` # obtener la matriz y guardarla en `cov_matrix` cov_matrix <- summary(mod)$cov.unscaled # calcular los errores estándar y asignarlos al vector `SEs` SEs <- sqrt(diag(cov_matrix)) test_object("cov_matrix") test_object("SEs")