def reg_predict(model, X):
    """
        Shim for making predictions using a fitted model and a (full) set of
        predictor data. Re-orders columns and adds an intercept as needed.
    """
    model_vars = list(model.params.index)
    if 'Intercept' in model_vars:
        X = X.assign(Intercept=1)

    X = X[model_vars]
    return model.predict(X)