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")