# How to create a loop that creates a set of ARIMA models and estimates the models to give LogL results

16 views (last 30 days)
harry gallow on 22 Jun 2018
Commented: Kevin Chng on 1 Aug 2018
Hi there,
I have manually created 16 arima models and have estimated them to produce LogL results in order for me to select the most appropriate model for my data.
See below for how I have manually created them.
mdl20=arima(1,0,1);
mdl21=arima(1,0,2);
mdl22=arima(1,0,3);
mdl23=arima(1,0,4);
% all the way to
mdl32=arima(4,0,1);
mdl33=arima(4,0,2);
mdl34=arima(4,0,3);
mdl35=arima(4,0,4);
then I have used this function to estimate them using my dataset L_Ret:
[EstMdl20,~,logL20] = estimate(Mdl20,L_Ret);
[EstMdl21,~,logL21] = estimate(Mdl21,L_Ret);
[EstMdl22,~,logL22] = estimate(Mdl22,L_Ret);
[EstMdl23,~,logL23] = estimate(Mdl23,L_Ret);
% all the way to
[EstMdl32,~,logL32] = estimate(Mdl32,L_Ret);
[EstMdl33,~,logL33] = estimate(Mdl33,L_Ret);
[EstMdl34,~,logL34] = estimate(Mdl34,L_Ret);
[EstMdl35,~,logL35] = estimate(Mdl35,L_Ret);
Any help will be immensly appreciated

#### 1 Comment

Kevin Chng on 1 Aug 2018
Hi,
I'm interested to know how you compare the model. AIC&BIC? or Residual?RMSE?
Regards, Kevin Chng

Prajit T R on 25 Jun 2018
Hi Harry
A good programming approach would be to use an array to store the corresponding values in a specific index position instead of declaring a new variable each time.
However, in this case the variable is of ARIMA type, which can not be stored in arrays. You can make a class to hold these ARIMA type variables, and create an array of objects of the class to simplify your job.
However, in your case you just want to see which model suits your workflow. You can follow this crude method to obtain the model estimates:
for i=1:4
for j=1:4
disp('Estimate for values:')
%Print i and j so that you can see which model is being processed
i
j
Mdl=arima(i,0,j)
[EstMdl,~,logL] = estimate(Mdl,L_Ret)
end
end
In each iteration, the values of EstMdl and logL will be printed and you can see which estimate suits you. Note that this after each iteration, the previous value of EstMdl and logL would be replaced. This is not an elegant solution, but in your case with just 16 models, this should be helpful.
Cheers
Prajit

#### 1 Comment

harry gallow on 25 Jun 2018
Thank you very much Prajit!!

R2017b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!