code: def test_simulate_vs_R(setup_model):
      model, params, results_R = setup_model
  
      fit = fit_austourists_with_R_params(
          model, results_R, set_state=True
      )
  
      innov = np.asarray(
          [[1.76405235, 0.40015721, 0.97873798, 2.2408932]]
      ).T
      sim = fit.simulate(
          4, anchor="end", repetitions=1, random_errors=innov
      )
      expected = np.asarray(results_R["simulation"])
  
      assert_allclose(expected, sim.values, rtol=1e-5, atol=1e-5)
  
imports: from statsmodels.compat.pandas import QUARTER_END
  from statsmodels.compat.platform import PLATFORM_LINUX32, PLATFORM_WIN
  from itertools import product
  import json
  import pathlib
  import numpy
  from numpy.testing import assert_allclose, assert_almost_equal
  import pandas
  import pytest
  import scipy.stats
  from statsmodels.tsa.exponential_smoothing.ets import ETSModel
  import statsmodels.tsa.holtwinters
  import statsmodels.tsa.statespace.exponential_smoothing
label_desc: Use image node_id 20 for calling a global function with example usage: test_simulate_vs_R(setup_model) without return types