code: def test_links():
      rng = np.random.default_rng(121195)
      data = pd.DataFrame(
          {
              "g": rng.choice([0, 1], size=100),
              "y": rng.integers(3, 10, size=100),
              "x": rng.integers(3, 10, size=100),
          }
      )
  
      FAMILIES = {
          "asymmetriclaplace": ["identity", "log", "inverse"],
          "bernoulli": ["identity", "logit", "probit", "cloglog"],
          "beta": ["logit", "probit", "cloglog"],
          "gamma": ["identity", "inverse", "log"],
          "gaussian": ["identity", "log", "inverse"],
          "negativebinomial": ["identity", "log", "cloglog"],
          "poisson": ["identity", "log"],
          "vonmises": ["identity", "tan_2"],
          "wald": ["inverse", "inverse_squared", "identity", "log"],
      }
      for family, links in FAMILIES.items():
          for link in links:
              if family == "bernoulli":
                  bmb.Model("g ~ x", data, family=family, link=link)
              else:
                  bmb.Model("y ~ x", data, family=family, link=link)
  
imports: import logging
  from functools import reduce
  from operator import add
  from os.path import dirname, join
  import arviz
  import numpy
  import pandas
  import pymc
  import pytest
  from formulae import design_matrices
  import bambi
  from bambi.terms import CommonTerm, GroupSpecificTerm
label_desc: Use image node_id 22 for calling a global function with example usage: test_links() without return types