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) # obtain the matrix and save it to `cov_matrix` # compute the standard errors and assign them to the vector `SEs` # obtain the matrix and save it to `cov_matrix` cov_matrix <- summary(mod)$cov.unscaled # compute the standard errors and assign them to the vector `SEs` SEs <- sqrt(diag(cov_matrix)) test_object("cov_matrix") test_object("SEs")