code: def test_normal_with_splines(data_normal, data_new_normal):
      knots = np.linspace(0, 23, 8)[1:-1]
      knots_s = np.linspace(0, 23, 5)[1:-1]
      formula = bmb.Formula(
          "count_normalized ~ 0 + bs(hour, knots=knots, intercept=True)",
          "sigma ~ 0 + bs(hour, knots=knots_s, intercept=True)",
      )
      model = bmb.Model(formula, data_normal)
      idata = model.fit(tune=100, draws=100)
      model.predict(idata, kind="pps")
      model.predict(idata, kind="pps", data=data_new_normal)
  
file_path: bambi/tests/test_glmlss.py
imports: import numpy
  import pandas
  import pytest
  import bambi
label_desc: Use image node_id 5 for calling a global function with example usage: test_normal_with_splines(data_normal, data_new_normal) without return types