def METHOD_NAME():
    # Generalized Linear Model (GLM)
    data = get_dataset("scotland")
    data.exog = sm.add_constant(data.exog)
    glm = sm.GLM(data.endog, data.exog, family=sm.families.Gamma())
    model = glm.fit()
    return ModelWithResults(model=model, alg=glm, inference_dataframe=data.exog)