title: Simulating a basic sarima model in R
body: <p>I'm looking for something like </p>
  
  <pre><code>   arima.sim()
  </code></pre>
  
  <p>but for sarima models.  I've looked at simulate.Arima() in the forecast package, but it seems to require an input dataset parsed by Arima(), which I don't want to do.  I also looked at the gsarima library, but it seems to be only able to simulate seasonal AR models.  Is there any way to simulate a sarima model if you only want to provide the following information:</p>
  
  <ol>
  <li><p>The values of all pure ARMA and seasonal ARMA terms for different lags.</p></li>
  <li><p>The number of differences for both the seasonal and non-seasonal integrated terms.</p></li>
  <li><p>The length of a season.</p></li>
  <li><p>The number of terms I want to simulate.</p></li>
  </ol>