Estimated error of prediction, assuming multivariate normality
gmsep(model)
lm
GMSEP
model <- lm(mpg ~ disp + hp + wt + qsec, data = mtcars) gmsep(model)#> [1] 8.18104