Statistical analysis of autoregressive fractionally integrated moving average models in R