'(cached attribute) estimate of variance of the residuals
sigma2 = sigma2_OLS * (1 - hii)
where hii is the diagonal of the hat matrix'

@cache_readonly
def resid_var(self):

    return (self.results.mse_resid * (1 - self.hat_matrix_diag))
