We have started with the practical session of time series analysis in SAS software we were dealing with the SAS data set that is ln data set which consists of two variables data and passengers. It consists of the number of passengers that it contains the number of passengers that has traveled from the month 1949 that is from the date Jan 1949 till December 1960, that is it contains the number of passengers that has traveled through Jet Airways Ellen company from the month Jan 2014, till the month December 1960. And our objective is to predict for the coming 12 months So, for that we are going to try modeling to at first we have checked stationarity of the data through scatterplot. Then we had done Arima modeling we saw the cnpc graphs, the We did a DNA test that is augmented Dickey fuller test to check non stationarity all the three tests shows that the data is very much non stationary in nature, then we then we did animal modeling with different periods of differencing.
First we did animal modeling with period of differencing one then period two and then Piero differencing as one comma one in order to gradually move towards stationarity. We have seen that our ACF graphs a PC graphs were gradually changing is verging gradually moving towards stationarity we also found a huge plunge in autocorrelation between one and 12 lakh in our first resolver which we had got through animal modeling. That is this one. See there is a huge plunge in autocorrelation in the first 1112 So, now we are going to use this one comment with as the accuracy of differencing to store the stationarity. So, let's do that. We are going to use the procedure proc Arima data equals my new corn time series underscore line identify VOD passengers.
I'm using the pH DS one comma 12. That is a appropriate period of difference in which we had caught enough first result viewer see weird Then run and then quit let's run this code see her obviously visa process completely changed the lens of autocorrelation length of come between the within the standard error region that means we have distraught stationarity that is our data has distribution at this is our autocorrelation table. Now the behavioral differences comment with no working CDC 0.13206 standard deviations with countries or an infinite number of observations from 30 that is at one comma 12 behavioral difference in 13 observations are elevated. So, the meaning observations are from 31. Now, we are going to do the mini stable Just because we have restored stationarity now we'll be using the mini table to find the optimal answer for a model and a model. So we'll be using the same procedure proc Arima later.
Equals like rename is mainly been done dark data set in a nice time series underscore and name var equal to passengers passengers is man and become available. I'm reading that period of difference he has one comment well then I'm using the Keyboard mini mini stands for minimum information criteria. So we'll be using the mini table to find optimal Elsevier and a model from and then we could execute the code Let me explain on the code. I'm using the procedure Procurement Data on the library on the data set entities underscore ln which is located inside the library mainly my analytical variable passengers. So use the statement I will say var equal to passengers one comment when we set up a period or difference in that we have got in our last reserved viewer, that is towards the beginning. This is huge planter autocorrelation at first when it was 0.9 coral over here and zero to 60 days a huge plant so we'll be using one common terms that refer to differencing.
Then we are using the key word meaning that stands for minimum information criteria. Many people is used to document lakhs of arnn a model will be using this table and then run a little bit. So now let's execute this code. See this is many, many people available information criteria table here are the minimum table values given these values are basically pac, pac stands for Bayesian information criteria, which explains the amount of unexplained variation in our model. So, lesser is the PAC value better is a model because lesser is unexplained variation better is the model over here the minimum Pac values given us hope on a 369. So, this is the minimum basic minimum value of all the values in the table.
So optimal level VR will be one and optimal level a will be zero. So optimal level VR model that is P, P will be equal to one and Q will be zero. So P and Q are lsv got this optimal level v nav model This is the optimal relationship for like no stimulus given period of differences. One comment with me no CDs, same 0.13206, which was before it, I'm a modeling with one common difference externality deviation. Again, it's the same number of observations to the same one that is 13. observations were eliminated by doing the differencing. And one common trend.
So our optimal lots is areas one p equal to one, if you will, to see those that we have to move down. So let's note it down because we have to use that for forecasting the variables to one and q1 is equal to zero. That's notice Don't be using it. We'll be forecasting by applying this optimal less desirable disarmament Using the same procedure property manager will be focusing for the coming 12 months my lab reading is against in my live one.com citizen discord bar as he just as he uses meat and unethical variables that are being used very very imagery blasting and general referencing which is document one that is one comment well then I will give document labs estimate the values as I told as we calculated is equal to one from the meaning table and Q value is equal to Zero via document every animal case or lack of a model.
We haven't forecast for the coming 12 months from reading forecast lead equals to 12. Thank you. So, let me explain the code. Here we are doing Arima modeling on the data set MCs underscore nn located in segment one, we are using passages variable which is an analytical variable our auto references can comment where there is a lot of fencing, which we got in our beginning in the first result viewer where the plunging on conditions between one comma 12 So, this was the planche we obtained in the difference in the span command Well, from the mini table we got the optimum lots of air and in a model that is from the minimum be AC very open a three string air model is run an animal to zero that is equal to one to zero. And we are forecasting for the coming 12 months to forecast legal structure.
So let's run this code should be a semicolon over here. No Let's run the score. To see this is a forecast from video. passages in our original data set we had till 144 observation that is from 2014 jan jan 1842, December 2016, Mr. Cooperman 44 observations per 141 words it was dis forecasted. So, it is forecasted from 45 onwards these are these are the forecast values for the coming 12 months this is standard error This is an indicator percent confidence limits, we got the autoregressive factors because the estimated mean the period or differencing is one comment well these are the graphs these are the seven basic graphs. This is the forecast for passengers observations with the forecast oils now, we'll be learning to hear in this video.
And with this video we have completed the scripted modeling program. I hope you all have enjoyed this session. Thank you Goodbye and All the best to you