'predict response, sum of smoothed components
TODO: What\'s this in the case of GLM, corresponds to X*beta ?'

def predict(self, exog):

    exog_smoothed = self.smoothed(exog)
    if (exog_smoothed.shape[0] == self.k_vars):
        import warnings
        warnings.warn('old    orientation,    colvars,    will    go    away', FutureWarning)
        return (np.sum(self.smoothed(exog), axis=0) + self.alpha)
    if (exog_smoothed.shape[1] == self.k_vars):
        return (np.sum(exog_smoothed, axis=1) + self.alpha)
    else:
        raise ValueError('shape    mismatch    in    predict')
