code: def data_gamma():
      rng = np.random.default_rng(121195)
      N = 200
      a, b = 0.5, 1.1
      x = rng.uniform(-1.5, 1.5, N)
      shape = np.exp(0.3 + x * 0.5 + rng.normal(scale=0.1, size=N))
      y = rng.gamma(shape, np.exp(a + b * x) / shape, N)
      return pd.DataFrame({"x": x, "y": y})
  
file_path: bambi/tests/test_glmlss.py
imports: import numpy
  import pandas
  import pytest
  import bambi
label_desc: Use image node_id 1 for calling a global function with example usage: data_gamma() and returns: pd