This function fits an autoregression (AR) model to quantile series (QSER) separately for each quantile level using stats::ar().
qser2ar(y.qser, p = NULL, order.max = NULL, method = c("none", "gamm", "sp"))a list with the following elements:
matrix or array of AR coefficients
vector or matrix of residual covariance
order of AR model
length of time series
matrix or array of residuals
matrix or array of pre-calculated QSER, e.g., using qser()
order of AR model (default = NULL: selected by AIC)
maximum order for AIC if p = NULL (default = NULL: determined by stats::ar())
quantile smoothing method: "gamm", "sp", or "NA" (default)