Cross Validated Asked by Bromideh on December 3, 2020

I have shared the main data d2015.txt, includes 5 columns. The first column is the $y_t$ the time series observation. The 2nd to the 4th column is the covariates/regressors and the fifth column is the mean of the three covariates. The main file shared in my google drive.

I used this code to fit an ARIMA model in which "ur" is univariate time series of observation and the "gt" is a multivariate covariates dataset with 3 columns:

```
> fit2 = auto.arima(ur,xreg=gt)
```

the output is:

```
Series: ur
Regression with ARIMA(1,1,0) errors
Coefficients:
ar1 drift Apply office resume
0.5109 -0.1562 0.0098 0.0001 -5e-04
s.e. 0.1141 0.0277 0.0013 0.0017 9e-04
sigma^2 estimated as 0.01239: log likelihood=50.62
AIC=-89.24 AICc=-87.71 BIC=-76.48
Training set error measures:
ME RMSE MAE MPE MAPE MASE ACF1
Training set 0.002114705 0.10588 0.08902272 0.04035552 0.5226742 0.497244 -0.06112877
```

My first question is: how to formulate this ARIMA into mathematical writing?

$Y_t$= ???

I used `pred2 = forecast(fit2,h=1, xreg=gt)`

to forecast the next month, but it predicts for more months!!! Why? Can you correct me, pls?

Note that

`auto.arima()`

fits a regression with ARIMA errors, not an ARIMAX model. There is a difference. Thus, your model is$$ y_t=X_tbeta+epsilon_t,$$

where $X_t$ contains your regressors at time $t$, and

$$ epsilon_tsimtext{ARIMA(1,1,0) with drift,} $$

or

$$ (epsilon_t-epsilon_{t-1}) = gamma t+phi(epsilon_{t-1}-epsilon_{t-2})+eta_t.$$

$gamma$ is estimated by your

`drift`

parameter and $phi$ by the`ar1`

one.You can find what is happening in

`?forecast.Arima`

(note the capital`A`

), under the entry for`h`

in the`Arguments`

section:If ‘xreg’ is used, ‘h’ is ignored and the number of forecast periods is set to the number of rows of ‘xreg’.

Of course, you need to supply the new regressors when forecasting, and then the number of rows in this new regressor matrix is used to determine the forecasting horizon. However, you are using the regressors

`gt`

you used in fitting for forecasting. You probably don't want to do this.

Answered by Stephan Kolassa on December 3, 2020

Get help from others!

Recent Questions

- How can I transform graph image into a tikzpicture LaTeX code?
- How Do I Get The Ifruit App Off Of Gta 5 / Grand Theft Auto 5
- Iv’e designed a space elevator using a series of lasers. do you know anybody i could submit the designs too that could manufacture the concept and put it to use
- Need help finding a book. Female OP protagonist, magic
- Why is the WWF pending games (“Your turn”) area replaced w/ a column of “Bonus & Reward”gift boxes?

Recent Answers

- Peter Machado on Why fry rice before boiling?
- Jon Church on Why fry rice before boiling?
- Lex on Does Google Analytics track 404 page responses as valid page views?
- haakon.io on Why fry rice before boiling?
- Joshua Engel on Why fry rice before boiling?

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP