code: def test_sample_model():
      """
      Tests sampling function
      """
      model = ARIMA(data=data, ar=2, ma=2, family=Skewt())
      x = model.fit("BBVI", iterations=100, quiet_progress=True)
      sample = model.sample(nsims=40)
      assert sample.shape[0] == 40
      assert sample.shape[1] == len(data) - 2
  
file_path: pyflux/pyflux/arma/tests/test_arima_skewt.py
imports: import numpy
  from pyflux.arma import ARIMA
  from pyflux.families import Skewt
label_desc: Use image node_id 10 for calling a global function with example usage: test_sample_model() without return types