simulates bootstrap samples for the stochastic process y, using a stationary
auto-regressive model of order "pmax", AR(pmax). If
pmax = NULL (default), the function estimates the process maximum
lags using an AIC as a model selection criteria.
References
Bulmann, P. (1997). Sieve Bootstrap for time series. Bernoulli.
3(2), 123 -148.