'Return the difference between observed and fitted values.
Returns
resid : array (nobs x neqs)
The residuals.'

@cache_readonly
def resid(self):

    return (self.y_all.T[self.k_ar:] - self.fittedvalues)
