code: def test_gamma(data_gamma, data_new_gamma):
      formula = bmb.Formula("y ~ x", "alpha ~ x")
  
      # Default links
      model = bmb.Model(formula, data_gamma, family="gamma")
      idata = model.fit(tune=100, draws=100, random_seed=1234)
      model.predict(idata, kind="pps")
      model.predict(idata, kind="pps", data=data_new_gamma)
  
      # Custom links
      model = bmb.Model(
          formula,
          data_gamma,
          family="gamma",
          link={"mu": "log", "alpha": "log"},
      )
      idata = model.fit(tune=100, draws=100, random_seed=1234)
      model.predict(idata, kind="pps")
      model.predict(idata, kind="pps", data=data_new_gamma)
  
file_path: bambi/tests/test_glmlss.py
imports: import numpy
  import pandas
  import pytest
  import bambi
label_desc: Use image node_id 6 for calling a global function with example usage: test_gamma(data_gamma, data_new_gamma) without return types