code: def test_2_sample_model():
      """
      Tests sampling function
      """
      model = ARIMAX(
          formula="y ~ x1 + x2", data=data, ar=2, ma=2, family=Skewt()
      )
      x = model.fit("BBVI", iterations=100, quiet_progress=True)
      sample = model.sample(nsims=100)
      assert sample.shape[0] == 100
      assert sample.shape[1] == len(data) - 2
  
imports: import numpy
  import pandas
  from pyflux.arma import ARIMAX
  from pyflux.families import Skewt
label_desc: Use image node_id 20 for calling a global function with example usage: test_2_sample_model() without return types