'degrees of freedom of residuals, ddof is sum of all smoothers df'

def df_resid(self):

    return (self.results.Y.shape[0] - np.array([self.smoothers[i].df_fit() for i in range(self.exog.shape[1])]).sum())
