FUNCTIONAL COEFFICIENT AUTOREGRESSIVE NONLINEAR TIME-SERIES MODEL FOR DESCRIBING INDIA S LAC EXPORT DATA USING SAS VERSION 9.3
|
|
- Jeffery Gallagher
- 5 years ago
- Views:
Transcription
1 FUNCTIONAL COEFFICIENT AUTOREGRESSIVE NONLINEAR TIME-SERIES MODEL FOR DESCRIBING INDIA S LAC EXPORT DATA USING SAS VERSION 9.3 Ranit Kumar Paul and Himadri Ghosh I.A.S.R.I., Library Avenue, Pusa, New Delhi ranitstat@gmail.com, him_adri@iasri.res.in INTRODUCTION Box Jenkins linear autoregressive integrated moving average (ARIMA) methodology is widely used for analyzing time-series data. Beyond linear domain, there are many nonlinear forms to be explored. In fact, nonlinear time-series analysis has been one of the maor areas of research in Time-series analysis for more than two decades now. These models are generally more appropriate than linear models for accurately describing dynamics of the series and for making multistep-ahead forecasts. Early development of nonlinear time-series analysis focused on various parametric forms. Engle (1982), in a path breaking work, proposed Autoregressive conditional heteroscedastic (ARCH) model for modelling volatility present in a data set. However, the conditional variance of ARCH(q) model, where q indicates the order of lag, has the property that the unconditional autocorrelation function of squared residuals, if it exists, decays very rapidly, unless q is large. To overcome this limitation of ARCH model, Bollerslev (1986) proposed the Generalized ARCH (GARCH) model, in which the unconditional autocorrelation function of squared residuals has slow decay rate. Unlike ARIMA model, these models are able to capture the presence of heteroscedasticity of conditional error variances. Another important family is bilinear timeseries model, proposed by C. W. G. Granger in 1978, which is capable of modelling data sets in which outliers appear in random epochs (Ghosh et al., 2006b). A heartening feature of the third important family, viz. Self exciting threshold autoregressive (SETAR) family, proposed by H. Tong, is that it is able to describe cyclical data quite efficiently (Tong, 1995). Quite often it is not possible to postulate appropriate parametric form for the underlying phenomenon and, in such cases; Nonparametric approach is called for. In the versatile model, viz. Functional-coefficient autoregressive (FCAR) model, introduced by Chen and Tsay (1993), the coefficient function changes gradually rather than abruptly. Fan and Zhang (2000) applied this model for the analysis of longitudinal data, in which different coefficient functions share a single smoothing variable time t, and the two-step estimates of coefficient functions in the model were given. In the present study, an attempt is made to apply FCAR model to country s yearly export data of lac from 1900 to Comparison of the performance of FCAR model with the SETAR and ARIMA models for modelling as well as forecasting is also carried out. DESCRIPTION OF FCAR MODEL Functional-coefficient autoregressive (FCAR) nonparametric nonlinear time-series model, introduced by Chen and Tsay (1993), admits the form X a ( X ) X... a ( X ) X ε (1) t 1 t-d t-1 where { t } is a sequence of independent random variables with zero mean and unit variance, and t is independent of X t-1, X t-2,. The coefficient- functions a 1 (.), a 2 (.),, a p (.) p t-d t-p t
2 are unknown and change gradually rather than abruptly. It is a direct extension of linear AR model, but allows the coefficients to vary according to a threshold variable X t-d. Nonparametric procedures are used to estimate the functions in the model, hence allowing data to speak for themselves regarding the model to be used. FCAR model can be regarded as a Stochastic regression model by introducing dependent variable Y as current observation X t, the i-th independent variable X i as lag i variable X t-i, and U as lag d variable X t-d. With induced variables, FCAR model (Cai et al., 2000) can be written as Y t = a 1 (U) X t-1 + a 2 (U) X t-2 + +a p (U) X t-p + t. (2) It is noted that many of the successful parametric nonlinear models belong to FCAR family. For example, if the functions (X ) a ( X ) I X c, it a i in (1) are step functions 2 x reduces to Threshold autoregressive (TAR) model. When ai ( X ) i i e, the model becomes an Exponential AR (EXPAR) model. SETAR and many other models also belong to this class. Hence, nonparametric determination of functional forms in the model may provide obective guidelines on choosing an appropriate parametric model. It also allows researchers to develop new models that are useful in their applications by specifying a parametric form for the coefficient functions based on nonparametric estimates. Estimation of parameters In order to apply FCAR model to data, the coefficient-functions may be estimated by using a local linear regression technique. These coefficient-functions are expanded by Taylor s series expansion in which unknown coefficients are estimated by the Method of weighted least squares, weights being the kernel density function (Fan and Yao, 2003). For any given u 0 and u in a neighborhood of u 0 and using Taylor s series expansion: a u) a ( u ) a ( u )( u u ) a b ( u ) (3) ( u0 where a and b are local intercept and slope. Using the data with U i around u 0 and local model (3), the following expression is minimized: n p Yi a b U i u0 X i K h U i u0 (4) i1 1 where K h (.) = h -1 K(./h), K(.) is a kernel function and h is the bandwidth. Then, the local linear regression estimator is simply aˆ ( u ) aˆ 0. The local linear regression estimator a and b can be easily obtained. Let e, 2 p be the 2p 1 unit vector with 1 at the th position, X ~ denote an n 2p matrix with X X( u ) as its i th row, and W diag K minimizing h i, i U i 0 T Y ( Y 1,..., Y n ), where superscript T indicates transpose. Set ( U1 u0 ),..., K h ( U n u0 ). Then the local regression problem reduces to ~ T ~ Y Xβ WY Xβ, where β a..., a, b,..., T 1. The local least square estimator is β XW X XW Y 1, 1 p b p 2 i ˆ ~ i i 1 ~ 1 ~ 2
3 Selection of bandwidth and model dependent variable As discussed in Cai et al. (2000), optimal bandwidth h in the local linear regression methodology is selected by Modified multifold cross-validation criterion. Taking m and Q as two positive integers such that n > mq, the parameters are estimated using various bandwidth values h and the Q-fitted models are used for carrying out one-step forecasting error of the next section of the time-series of length m based on the estimated models. Let a ˆ,q. be the estimated coefficients using q th, q = 1, 2,, Q subseries {(U i, X i, Y i ), 1 i n-qm}. The average prediction error using q th subseries is given by 1 APE q m The overall average prediction error is given by APE 2 nqmm p i, q i X i, 1 1 i n qm h Y aˆ U h Q 1 Q q1 APE q h. (5) The proposed data driven bandwidth is the one that minimizes APE(h). In practice, generally m = [0.1n] and Q = 4 are considered. The selected bandwidth does not depend critically on the choice of m and Q so long as mq is reasonably large, thereby ensuring that evaluation of prediction errors is stable. Choosing an appropriate model dependent variable U is also very important. Knowledge of physical background of data may be very helpful. Without any prior information, it is pertinent to choose U in terms of some data driven methods, such as Akaike information criterion (AIC), cross - validation and other criteria. Let APE(h,d) be the average prediction error defined by (2.5) using lagged variable U = X t-d. A simple and practical approach is to minimize APE(h,d) simultaneously for h in a certain range and d over the set {1, 2,, p}. The order p can also be chosen to minimize the APE. An Illustration (Ghosh et al, 2010) Lac is a resinous protective secretion of the tiny lac insect. Conventionally, it is used for ewellery, scaling wax, wood polish, mirror coating, leather finishing, etc. It is also used for coating of citrus fruits, chocolate chewing gum, printing ink, nail polish, solar cell, etc. Climatic factor have a profound influence on the development, life-cycle, larval emergence, etc. of this insect. Naturally, therefore, crop production is determined to a large extent by temperature, humidity and rainfall. These factors along with presence of interactions among three-species system involving host trees biomass and populations of lac insect and its predators are mainly responsible for cyclical variations in production as well as in export quantity of lac. The graph of country s lac export data during the period 1900 to 2000 (Fig. 1), obtained from Shellac and Forest Products Export Promotion Council, Kolkata, India, depicts this type of pattern. In general setup of time series observations the above data can be taken to be observed against chronological evolution of time like 1,2,3, Further, the ascent periods tend to exceed descent periods by approximately 50%, thereby indicating asymmetricity of the data set. 3
4 Fitting of ARIMA model by SAS Partial sas output Fig. 1. Indian Lac export data data lac; input y; /* y denotes the lac export*/ cards; ; proc arima data= lac; identify var=y(2) ; /* differencing the series twice*/ estimate p=2 q=4; forecast lead=4 out=flac ; run; proc print data=flac; run; Name of Variable = export Period(s) of Differencing 2 Mean of Working Series Standard Deviation Number of Observations 99 Observation(s) eliminated by differencing 2 Autocorrelation Check for White Noise To Chi- Pr > Lag Square DF ChiSq Autocorrelations
5 Conditional Least Squares Estimation Standard Approx Parameter Estimate Error t Value Pr > t Lag MU MA1, < MA1, MA1, MA1, AR1, AR1, < A perusal of Fig. 1 indicates that the country s lac export data is not stationary. Therefore, the original series is differenced repeatedly unless it becomes stationary. It is found that this is achieved after second differencing. On investigating the autocorrelation functions (acf) and partial autocorrelation functions (pacf) of the resultant series, as well as minimum AIC and BIC values computed respectively by (6) and (7), ARIMA(2,2,4) model is selected. AIC = - 2 log(l) + 2p (6) BIC = - 2 log(l) + p log(n) (7) where L is the likelihood function, p denotes the number of parameters, and N is the sample size. The parameters of this model, which is a particular case of Cyclical trend model given in Harvey (1996), are estimated by using SPSS, Ver software package and the same are reported along with other related statistics in Table 1. For selected model, the AIC and BIC values are respectively computed as and Table 1: Parameter estimates along with other related statistics Parameters Estimate Standard t-statistic Probability Error Constant AR <.0001 AR MA MA MA MA <.0001 Fitting of FCAR model Optimum FCAR model is selected by the methodology discussed in Section 2.2. Estimates of parameters of this model are obtained by using computer program written in SAS- IML and appended as an Annexure. To ease the calculation, logarithm of the data is considered for analysis. The optimum values for p and d are found as p = 4 and d = 3. The estimated coefficient functions are plotted in Fig. 2, which clearly shows that the coefficient functions a i (u), i = 1, 2, 3, 4 are no longer constants but are changing as u varies. The function APE against the bandwidth over a grid of point h = 0.1 ( = 1, 2,, 10) is computed. The selected bandwidth along with autoregressive order p and delay parameter d which minimizes the APE defined in (5) is found as h =
6 Finding optimum values of (p, d, h) on the basis of APE criterion in SAS IML data lac; input export; cards; ; Proc iml; use lac; read all into Y;*/specify the time-series vector; pdh=(1,4,999); */Initializing the vector; APEh=(4,1,0); u0=sum(y)/nrow(y); do p=1 to 8; do d= 1 to p; x0=(nrow(y),p,0); do i = p+1 to nrow(y); x0[i,]=y[i-1:i-p]`; do h=0.1 to 1 by 0.1; do =1 to 4; t=nrow(y)-7*; n=t-p; x=(t,p,0); /*forming the matrix of the lagged values*/ u=(t,1,0); do i = p+1 to t; x[i,]=y[i-1:i-p]`; u[i]=y[i-d]; x1=x[p+1:t,]; U1=u[p+1:t,]; w=(n,n,0); w1=(n,1,0); /* kernel function is k(x)=1/h*k(x/h)*/ /*kx=(1/(h*sqrt(2*3.14))*exp(-x*x/(2*h*h)));*/ /*calculating the weight matrix by using kernel function*/ x2=(n,p,1); do i = 1 to n; ind=(u1[i]-u0)/h; /* ind is the indicator variable*/ if ind<=1 then in=1;else in=0; w[i,i]=(0.75/h)*(1-((u1[i]-u0)*(u1[i]-u0))/(h*h))*in; w1[i,]=w[i,i]; x2[i,]=x1[i,]*(u1[i]-u0); Xcurl=(n,2*p,0); xcurl=x1 x2; y1=y[p+1:t,]; beta=inv(xcurl`*w*xcurl)*(xcurl`*w*y1); beta1=beta[1:p,]; 6
7 yhat=xcurl*beta; error=y1-yhat; apeerror=(7,1,0); apeerror=error[n-7*+1:n-7*+7,]; APEh[]=(1/7)*apeerror`*apeerror; APE=sum(APEh)/4; pdh=pdh//(p d h APE); pdh=pdh[2:nrow(pdh),]; minape=min(pdh[,4]); do i=1 to nrow(pdh); if (pdh[i,4]=minape) then rw=i; result=pdh[rw,]; print beta1; run; Partial SAS Output p d h APE
8
9 FCAR nonlinear time-series model for describing Indian lac export data minape p d h Beta p and d are the order of FCAR model where p denotes AR lag and d is the delay parameter. h is the optimal bandwidth. minape denotes the minimum Average Prediction Error a 1 (u) a 2 (u) a 3 (u) a 4 (u) Fig. 2. The functional forms of a 1 (u), a 2 (u), a 3 (u) and a 4 (u) 9
10 The fitted FCAR (4,3) model for log data is obtained by using equations 2 and 3, and using the parameter estimates Beta1 (SAS output) as follows: log X t = { (log X t )} log X t-1 + { (log X t )} log X t-2 +{ (log X t )} log X t-3 + { (log X t )} log X t-4 Here, X t-3 is the threshold value and is the mean of the data By taking antilog of fitted values by above model, performance of the fitted model is evaluated by computing the AIC and BIC values: AIC = N log(rss/n) + 2p (8) BIC = N log(rss/n) + p log(n) (9) where RSS is the residual sum of squares by minimizing (2.4), p is the total number of parameters to be estimated and N is the sample size. For the present example, the AIC and BIC values computed form (8) and (9) are and respectively. Ghosh et al. (2006a) had earlier fitted the following SETAR (2;1,2) model to the same data: X t X t-1 = X t-1, if X t-1 < = X t X t-2, if X t-1 > with AIC value computed as , which is slightly more than the one for FCAR model. Thus, for describing the data set under consideration, FCAR model is slightly better than SETAR model and both of these are superior to ARIMA model. To get a visual idea, graph of fitted FCAR model along with data points is exhibited in Fig. 3. Fig. 3. The fitted FCAR(4, 3) model along with the data points Forecasting One-step ahead forecasts and iterative two-step ahead forecasts of lac export during the years 2001 to 2004 for ARIMA (2,2,4), SETAR(2;1,2) and FCAR(4,3) models are carried out and reported in Table 2 and Table 3 respectively. For one-step ahead forecast, let T=2000 be the origin of forecast. Let ŷt i1 T i be the one step ahead forecast of yt i1 given information upto time T+i, i=0,1,2,3. Similarly, for two-step ahead forecast, let T=1999 be the origin of forecast. Let ŷt i2 Ti be the one step ahead forecast of yt i2 given information upto time T+i, i=0,1,2,3. The models are compared on the basis of Mean absolute prediction error (MAPE), Mean square prediction error (MSPE) and Relative mean absolute prediction error (RMAPE) given by 10
11 4 MAPE = 1 / 4y t ŷ i1 4 i ti MSPE = 1/ 4 y yˆ i1 4 ti 2 ti RMAPE = 1/ 4y ˆ ti yti / yti100 i1 Computed values of MAPE, MSPE and RMAPE for the above mentioned three models are reported in Table 4. lead time of forecast Table 2. One-step ahead Forecast of lac export Actual Forecast by ARIMA model SETAR model FCAR model lead time of forecast Table 3. Two-step ahead Forecast of lac export Actual Forecast by ARIMA model SETAR model FCAR model Table 4. Comparison of forecast performance One-step ahead Two-step ahead Models MAPE MSPE RMAPE MAPE MSPE RMAPE ARIMA(2,2,4) SETAR(2,1,2) FCAR(4,3) A perusal of Table 4 indicates that FCAR model has performed better than both the SETAR as well as ARIMA models for forecasting. To sum up, for the data under consideration, FCAR model is found to be the best for modelling as well as forecasting. Concluding Remarks FCAR nonparametric nonlinear time-series model is applied for modelling and forecasting of Indian lac export data by using SAS. Superiority of FCAR model over ARIMA and SETAR models for the data under consideration is clearly demonstrated. One-step ahead forecasts and iterative two-step ahead forecasts of lac export during the years 2001 to 2004 for ARIMA (2,2,4), SETAR(2;1,2) and FCAR(4,3) models are carried out. On the basis of three statistics 11
12 namely MSPE, MAPE and RMAPE, it is concluded that FCAR model performs better than the SETAR model and ARIMA model for the data under consideration. References Bollerslev, T. (1986). Generalized autoregressive conditional heteroscedasticity. J. Econ. 31, Cai, Z., Fan, J. and Yao, Q. (2000). Functional-coefficient regression models for nonlinear time-series. J. Amer. Statist. Assoc. 95, Chen, R. and Tsay, R. S. (1993). Functional coefficient autoregressive models. J. Amer. Statist. Assoc. 88, Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of U.K. inflation. Econometrica 50, Fan, J., Zhang, J. (2000). Two-step estimation of functional linear models with applications to longitudinal data. J. Roy. Statist. Soc. Ser. B 62, Fan, J. and Yao, Q. (2003). Nonlinear Time-Series: Nonparametric and Parametric Methods. Springer, U.S.A. Ghosh, H., Sunilkumar, G. and Praneshu. (2006a). Self exciting threshold autoregressive models for describing cyclical data. Cal. Statist. Assoc. Bull. 58, Ghosh, H., Sunilkumar, G. and Praneshu. (2006b). Modelling and forecasting of bilinear time-series using frequency domain approach. J. Combi. Info. Sys. Sci. 31, Ghosh, H., Paul, R. K. and Praneshu. (2010). Functional coefficient autoregressive model for forecasting Indian lac export data. Mod. Assist. Statist. Appl. 5, Granger, C. and Terasvirta, T. (1993) Modelling Nonlinear Economic Relationships. Oxford University Press, Oxford Harvey, A. C. (1996). Forecasting Structural Time-Series Models and the Kalman Filter Cambridge University Press, U. K. Tong, H. (1995). Nonlinear Time-series Analysis: A Dynamic Approach. 2 nd Ed. Oxford University Press, Oxford. 12
FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL
FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL B. N. MANDAL Abstract: Yearly sugarcane production data for the period of - to - of India were analyzed by time-series methods. Autocorrelation
More informationTIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA
CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis
More informationDesign of Time Series Model for Road Accident Fatal Death in Tamilnadu
Volume 109 No. 8 2016, 225-232 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Design of Time Series Model for Road Accident Fatal Death in Tamilnadu
More informationDiagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations
Diagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations Farhat Iqbal Department of Statistics, University of Balochistan Quetta-Pakistan farhatiqb@gmail.com Abstract In this paper
More informationRevisiting linear and non-linear methodologies for time series prediction - application to ESTSP 08 competition data
Revisiting linear and non-linear methodologies for time series - application to ESTSP 08 competition data Madalina Olteanu Universite Paris 1 - SAMOS CES 90 Rue de Tolbiac, 75013 Paris - France Abstract.
More informationUnivariate ARIMA Models
Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.
More informationSTRUCTURAL TIME-SERIES MODELLING
1: Structural Time-Series Modelling STRUCTURAL TIME-SERIES MODELLING Prajneshu Indian Agricultural Statistics Research Institute, New Delhi-11001 1. Introduction. ARIMA time-series methodology is widely
More informationNonparametric Estimation of Functional-Coefficient Autoregressive Models
Nonparametric Estimation of Functional-Coefficient Autoregressive Models PEDRO A. MORETTIN and CHANG CHIANN Department of Statistics, University of São Paulo Introduction Nonlinear Models: - Exponential
More informationINTRODUCTION TO TIME SERIES ANALYSIS. The Simple Moving Average Model
INTRODUCTION TO TIME SERIES ANALYSIS The Simple Moving Average Model The Simple Moving Average Model The simple moving average (MA) model: More formally: where t is mean zero white noise (WN). Three parameters:
More informationSugarcane Productivity in Bihar- A Forecast through ARIMA Model
Available online at www.ijpab.com Kumar et al Int. J. Pure App. Biosci. 5 (6): 1042-1051 (2017) ISSN: 2320 7051 DOI: http://dx.doi.org/10.18782/2320-7051.5838 ISSN: 2320 7051 Int. J. Pure App. Biosci.
More informationEconometric Forecasting
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend
More informationForecasting using R. Rob J Hyndman. 2.4 Non-seasonal ARIMA models. Forecasting using R 1
Forecasting using R Rob J Hyndman 2.4 Non-seasonal ARIMA models Forecasting using R 1 Outline 1 Autoregressive models 2 Moving average models 3 Non-seasonal ARIMA models 4 Partial autocorrelations 5 Estimation
More informationIssues on quantile autoregression
Issues on quantile autoregression Jianqing Fan and Yingying Fan We congratulate Koenker and Xiao on their interesting and important contribution to the quantile autoregression (QAR). The paper provides
More informationEconometric modeling of the relationship among macroeconomic variables of Thailand: Smooth transition autoregressive regression model
The Empirical Econometrics and Quantitative Economics Letters ISSN 2286 7147 EEQEL all rights reserved Volume 1, Number 4 (December 2012), pp. 21 38. Econometric modeling of the relationship among macroeconomic
More informationDynamic Time Series Regression: A Panacea for Spurious Correlations
International Journal of Scientific and Research Publications, Volume 6, Issue 10, October 2016 337 Dynamic Time Series Regression: A Panacea for Spurious Correlations Emmanuel Alphonsus Akpan *, Imoh
More informationArma-Arch Modeling Of The Returns Of First Bank Of Nigeria
Arma-Arch Modeling Of The Returns Of First Bank Of Nigeria Emmanuel Alphonsus Akpan Imoh Udo Moffat Department of Mathematics and Statistics University of Uyo, Nigeria Ntiedo Bassey Ekpo Department of
More informationNANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION MAS451/MTH451 Time Series Analysis TIME ALLOWED: 2 HOURS
NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER II EXAMINATION 2012-2013 MAS451/MTH451 Time Series Analysis May 2013 TIME ALLOWED: 2 HOURS INSTRUCTIONS TO CANDIDATES 1. This examination paper contains FOUR (4)
More informationSeasonal Autoregressive Integrated Moving Average Model for Precipitation Time Series
Journal of Mathematics and Statistics 8 (4): 500-505, 2012 ISSN 1549-3644 2012 doi:10.3844/jmssp.2012.500.505 Published Online 8 (4) 2012 (http://www.thescipub.com/jmss.toc) Seasonal Autoregressive Integrated
More informationLecture 16: ARIMA / GARCH Models Steven Skiena. skiena
Lecture 16: ARIMA / GARCH Models Steven Skiena Department of Computer Science State University of New York Stony Brook, NY 11794 4400 http://www.cs.sunysb.edu/ skiena Moving Average Models A time series
More informationFunctional-Coefficient Dynamic Nelson-Siegel Model
Functional-Coefficient Dynamic Nelson-Siegel Model Yuan Yang March 12, 2014 Abstract I propose an Functional-coefficient Dynamic Nelson-Siegel (FDNS) model to forecast the yield curve. The model fits each
More informationEcon 423 Lecture Notes: Additional Topics in Time Series 1
Econ 423 Lecture Notes: Additional Topics in Time Series 1 John C. Chao April 25, 2017 1 These notes are based in large part on Chapter 16 of Stock and Watson (2011). They are for instructional purposes
More informationTime Series Analysis -- An Introduction -- AMS 586
Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data
More informationA stochastic modeling for paddy production in Tamilnadu
2017; 2(5): 14-21 ISSN: 2456-1452 Maths 2017; 2(5): 14-21 2017 Stats & Maths www.mathsjournal.com Received: 04-07-2017 Accepted: 05-08-2017 M Saranyadevi Assistant Professor (GUEST), Department of Statistics,
More informationECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications
ECONOMICS 7200 MODERN TIME SERIES ANALYSIS Econometric Theory and Applications Yongmiao Hong Department of Economics & Department of Statistical Sciences Cornell University Spring 2019 Time and uncertainty
More informationARIMA modeling to forecast area and production of rice in West Bengal
Journal of Crop and Weed, 9(2):26-31(2013) ARIMA modeling to forecast area and production of rice in West Bengal R. BISWAS AND B. BHATTACHARYYA Department of Agricultural Statistics Bidhan Chandra Krishi
More informationForecasting Area, Production and Yield of Cotton in India using ARIMA Model
Forecasting Area, Production and Yield of Cotton in India using ARIMA Model M. K. Debnath 1, Kartic Bera 2 *, P. Mishra 1 1 Department of Agricultural Statistics, Bidhan Chanda Krishi Vishwavidyalaya,
More informationEmpirical Market Microstructure Analysis (EMMA)
Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg
More informationMODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH. I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo
Vol.4, No.2, pp.2-27, April 216 MODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo ABSTRACT: This study
More informationAnalysis of Violent Crime in Los Angeles County
Analysis of Violent Crime in Los Angeles County Xiaohong Huang UID: 004693375 March 20, 2017 Abstract Violent crime can have a negative impact to the victims and the neighborhoods. It can affect people
More informationGoodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach
Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach By Shiqing Ling Department of Mathematics Hong Kong University of Science and Technology Let {y t : t = 0, ±1, ±2,
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationMODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY
MODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY The simple ARCH Model Eva Rubliková Ekonomická univerzita Bratislava Manuela Magalhães Hill Department of Quantitative Methods, INSTITUTO SUPERIOR
More informationFORECASTING YIELD PER HECTARE OF RICE IN ANDHRA PRADESH
International Journal of Mathematics and Computer Applications Research (IJMCAR) ISSN 49-6955 Vol. 3, Issue 1, Mar 013, 9-14 TJPRC Pvt. Ltd. FORECASTING YIELD PER HECTARE OF RICE IN ANDHRA PRADESH R. RAMAKRISHNA
More informationPrashant Pant 1, Achal Garg 2 1,2 Engineer, Keppel Offshore and Marine Engineering India Pvt. Ltd, Mumbai. IJRASET 2013: All Rights are Reserved 356
Forecasting Of Short Term Wind Power Using ARIMA Method Prashant Pant 1, Achal Garg 2 1,2 Engineer, Keppel Offshore and Marine Engineering India Pvt. Ltd, Mumbai Abstract- Wind power, i.e., electrical
More informationStudy on Modeling and Forecasting of the GDP of Manufacturing Industries in Bangladesh
CHIANG MAI UNIVERSITY JOURNAL OF SOCIAL SCIENCE AND HUMANITIES M. N. A. Bhuiyan 1*, Kazi Saleh Ahmed 2 and Roushan Jahan 1 Study on Modeling and Forecasting of the GDP of Manufacturing Industries in Bangladesh
More informationRomanian Economic and Business Review Vol. 3, No. 3 THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS
THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS Marian Zaharia, Ioana Zaheu, and Elena Roxana Stan Abstract Stock exchange market is one of the most dynamic and unpredictable
More informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationFIN822 project 2 Project 2 contains part I and part II. (Due on November 10, 2008)
FIN822 project 2 Project 2 contains part I and part II. (Due on November 10, 2008) Part I Logit Model in Bankruptcy Prediction You do not believe in Altman and you decide to estimate the bankruptcy prediction
More informationOn Autoregressive Order Selection Criteria
On Autoregressive Order Selection Criteria Venus Khim-Sen Liew Faculty of Economics and Management, Universiti Putra Malaysia, 43400 UPM, Serdang, Malaysia This version: 1 March 2004. Abstract This study
More informationLab: Box-Jenkins Methodology - US Wholesale Price Indicator
Lab: Box-Jenkins Methodology - US Wholesale Price Indicator In this lab we explore the Box-Jenkins methodology by applying it to a time-series data set comprising quarterly observations of the US Wholesale
More informationSuan Sunandha Rajabhat University
Forecasting Exchange Rate between Thai Baht and the US Dollar Using Time Series Analysis Kunya Bowornchockchai Suan Sunandha Rajabhat University INTRODUCTION The objective of this research is to forecast
More informationImplementation of ARIMA Model for Ghee Production in Tamilnadu
Inter national Journal of Pure and Applied Mathematics Volume 113 No. 6 2017, 56 64 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Implementation
More informationBox-Jenkins ARIMA Advanced Time Series
Box-Jenkins ARIMA Advanced Time Series www.realoptionsvaluation.com ROV Technical Papers Series: Volume 25 Theory In This Issue 1. Learn about Risk Simulator s ARIMA and Auto ARIMA modules. 2. Find out
More informationat least 50 and preferably 100 observations should be available to build a proper model
III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or
More informationCh 6. Model Specification. Time Series Analysis
We start to build ARIMA(p,d,q) models. The subjects include: 1 how to determine p, d, q for a given series (Chapter 6); 2 how to estimate the parameters (φ s and θ s) of a specific ARIMA(p,d,q) model (Chapter
More informationForecasting of Soybean Yield in India through ARIMA Model
Available online at www.ijpab.com Kumar et al Int. J. Pure App. Biosci. 5 (5): 1538-1546 (2017) ISSN: 2320 7051 DOI: http://dx.doi.org/10.18782/2320-7051.5834 ISSN: 2320 7051 Int. J. Pure App. Biosci.
More informationEcon 623 Econometrics II Topic 2: Stationary Time Series
1 Introduction Econ 623 Econometrics II Topic 2: Stationary Time Series In the regression model we can model the error term as an autoregression AR(1) process. That is, we can use the past value of the
More informationProblem Set 2 Solution Sketches Time Series Analysis Spring 2010
Problem Set 2 Solution Sketches Time Series Analysis Spring 2010 Forecasting 1. Let X and Y be two random variables such that E(X 2 ) < and E(Y 2 )
More informationDecision 411: Class 9. HW#3 issues
Decision 411: Class 9 Presentation/discussion of HW#3 Introduction to ARIMA models Rules for fitting nonseasonal models Differencing and stationarity Reading the tea leaves : : ACF and PACF plots Unit
More informationEmpirical Approach to Modelling and Forecasting Inflation in Ghana
Current Research Journal of Economic Theory 4(3): 83-87, 2012 ISSN: 2042-485X Maxwell Scientific Organization, 2012 Submitted: April 13, 2012 Accepted: May 06, 2012 Published: June 30, 2012 Empirical Approach
More informationAutomatic seasonal auto regressive moving average models and unit root test detection
ISSN 1750-9653, England, UK International Journal of Management Science and Engineering Management Vol. 3 (2008) No. 4, pp. 266-274 Automatic seasonal auto regressive moving average models and unit root
More informationForecasting the unemployment rate when the forecast loss function is asymmetric. Jing Tian
Forecasting the unemployment rate when the forecast loss function is asymmetric Jing Tian This version: 27 May 2009 Abstract This paper studies forecasts when the forecast loss function is asymmetric,
More informationUnivariate linear models
Univariate linear models The specification process of an univariate ARIMA model is based on the theoretical properties of the different processes and it is also important the observation and interpretation
More informationMinitab Project Report - Assignment 6
.. Sunspot data Minitab Project Report - Assignment Time Series Plot of y Time Series Plot of X y X 7 9 7 9 The data have a wavy pattern. However, they do not show any seasonality. There seem to be an
More informationOil price volatility in the Philippines using generalized autoregressive conditional heteroscedasticity
Oil price volatility in the Philippines using generalized autoregressive conditional heteroscedasticity Carl Ceasar F. Talungon University of Southern Mindanao, Cotabato Province, Philippines Email: carlceasar04@gmail.com
More informationProblem Set 2: Box-Jenkins methodology
Problem Set : Box-Jenkins methodology 1) For an AR1) process we have: γ0) = σ ε 1 φ σ ε γ0) = 1 φ Hence, For a MA1) process, p lim R = φ γ0) = 1 + θ )σ ε σ ε 1 = γ0) 1 + θ Therefore, p lim R = 1 1 1 +
More informationAPPLIED ECONOMETRIC TIME SERIES 4TH EDITION
APPLIED ECONOMETRIC TIME SERIES 4TH EDITION Chapter 2: STATIONARY TIME-SERIES MODELS WALTER ENDERS, UNIVERSITY OF ALABAMA Copyright 2015 John Wiley & Sons, Inc. Section 1 STOCHASTIC DIFFERENCE EQUATION
More informationTime Series I Time Domain Methods
Astrostatistics Summer School Penn State University University Park, PA 16802 May 21, 2007 Overview Filtering and the Likelihood Function Time series is the study of data consisting of a sequence of DEPENDENT
More informationAR, MA and ARMA models
AR, MA and AR by Hedibert Lopes P Based on Tsay s Analysis of Financial Time Series (3rd edition) P 1 Stationarity 2 3 4 5 6 7 P 8 9 10 11 Outline P Linear Time Series Analysis and Its Applications For
More information5 Transfer function modelling
MSc Further Time Series Analysis 5 Transfer function modelling 5.1 The model Consider the construction of a model for a time series (Y t ) whose values are influenced by the earlier values of a series
More informationNonlinear Time Series
Nonlinear Time Series Recall that a linear time series {X t } is one that follows the relation, X t = µ + i=0 ψ i A t i, where {A t } is iid with mean 0 and finite variance. A linear time series is stationary
More informationLesson 13: Box-Jenkins Modeling Strategy for building ARMA models
Lesson 13: Box-Jenkins Modeling Strategy for building ARMA models Facoltà di Economia Università dell Aquila umberto.triacca@gmail.com Introduction In this lesson we present a method to construct an ARMA(p,
More informationOn Forecast Strength of Some Linear and Non Linear Time Series Models for Stationary Data Structure
American Journal of Mathematics and Statistics 2015, 5(4): 163-177 DOI: 10.5923/j.ajms.20150504.01 On Forecast Strength of Some Linear and Non Linear Imam Akeyede 1,*, Babatunde Lateef Adeleke 2, Waheed
More informationStudy of Time Series and Development of System Identification Model for Agarwada Raingauge Station
Study of Time Series and Development of System Identification Model for Agarwada Raingauge Station N.A. Bhatia 1 and T.M.V.Suryanarayana 2 1 Teaching Assistant, 2 Assistant Professor, Water Resources Engineering
More informationTMA4285 December 2015 Time series models, solution.
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 5 TMA4285 December 205 Time series models, solution. Problem a) (i) The slow decay of the ACF of z t suggest that
More informationIDENTIFICATION OF ARMA MODELS
IDENTIFICATION OF ARMA MODELS A stationary stochastic process can be characterised, equivalently, by its autocovariance function or its partial autocovariance function. It can also be characterised by
More informationFunctional Coefficient Regression Models for Nonlinear Time Series: A Polynomial Spline Approach
Functional Coefficient Regression Models for Nonlinear Time Series: A Polynomial Spline Approach JIANHUA Z. HUANG University of Pennsylvania HAIPENG SHEN University of North Carolina at Chapel Hill ABSTRACT.
More informationFinQuiz Notes
Reading 9 A time series is any series of data that varies over time e.g. the quarterly sales for a company during the past five years or daily returns of a security. When assumptions of the regression
More informationModeling climate variables using time series analysis in arid and semi arid regions
Vol. 9(26), pp. 2018-2027, 26 June, 2014 DOI: 10.5897/AJAR11.1128 Article Number: 285C77845733 ISSN 1991-637X Copyright 2014 Author(s) retain the copyright of this article http://www.academicjournals.org/ajar
More informationEvaluation of Some Techniques for Forecasting of Electricity Demand in Sri Lanka
Appeared in Sri Lankan Journal of Applied Statistics (Volume 3) 00 Evaluation of Some echniques for Forecasting of Electricity Demand in Sri Lanka.M.J. A. Cooray and M.Indralingam Department of Mathematics
More informationTIME SERIES DATA PREDICTION OF NATURAL GAS CONSUMPTION USING ARIMA MODEL
International Journal of Information Technology & Management Information System (IJITMIS) Volume 7, Issue 3, Sep-Dec-2016, pp. 01 07, Article ID: IJITMIS_07_03_001 Available online at http://www.iaeme.com/ijitmis/issues.asp?jtype=ijitmis&vtype=7&itype=3
More informationarxiv: v1 [stat.me] 5 Nov 2008
arxiv:0811.0659v1 [stat.me] 5 Nov 2008 Estimation of missing data by using the filtering process in a time series modeling Ahmad Mahir R. and Al-khazaleh A. M. H. School of Mathematical Sciences Faculty
More informationThe log transformation produces a time series whose variance can be treated as constant over time.
TAT 520 Homework 6 Fall 2017 Note: Problem 5 is mandatory for graduate students and extra credit for undergraduates. 1) The quarterly earnings per share for 1960-1980 are in the object in the TA package.
More informationRoss Bettinger, Analytical Consultant, Seattle, WA
ABSTRACT USING PROC ARIMA TO MODEL TRENDS IN US HOME PRICES Ross Bettinger, Analytical Consultant, Seattle, WA We demonstrate the use of the Box-Jenkins time series modeling methodology to analyze US home
More informationForecasting of meteorological drought using ARIMA model
Indian J. Agric. Res., 51 (2) 2017 : 103-111 Print ISSN:0367-8245 / Online ISSN:0976-058X AGRICULTURAL RESEARCH COMMUNICATION CENTRE www.arccjournals.com/www.ijarjournal.com Forecasting of meteorological
More informationAnalysis. Components of a Time Series
Module 8: Time Series Analysis 8.2 Components of a Time Series, Detection of Change Points and Trends, Time Series Models Components of a Time Series There can be several things happening simultaneously
More information5 Autoregressive-Moving-Average Modeling
5 Autoregressive-Moving-Average Modeling 5. Purpose. Autoregressive-moving-average (ARMA models are mathematical models of the persistence, or autocorrelation, in a time series. ARMA models are widely
More informationFORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL
FORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL S.Poyyamozhi 1, Dr. A. Kachi Mohideen 2. 1 Assistant Professor and Head, Department of Statistics, Government Arts College (Autonomous), Kumbakonam
More informationIntroduction to ARMA and GARCH processes
Introduction to ARMA and GARCH processes Fulvio Corsi SNS Pisa 3 March 2010 Fulvio Corsi Introduction to ARMA () and GARCH processes SNS Pisa 3 March 2010 1 / 24 Stationarity Strict stationarity: (X 1,
More informationPrediction of Grain Products in Turkey
Journal of Mathematics and Statistics Original Research Paper Prediction of Grain Products in Turkey Özlem Akay, Gökmen Bozkurt and Güzin Yüksel Department of Statistics, The Faculty of Science and Letters,
More informationCircle the single best answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 6, 2017 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 32 multiple choice
More informationFive Decades of Research in Statistical Modelling: An Overview
Five Decades of Research in Statistical Modelling: An Overview Prajneshu and Himadri Ghosh 1. Introduction Statistical modelling plays a very important role in comprehending underlying relationships among
More informationShort-Term Load Forecasting Using ARIMA Model For Karnataka State Electrical Load
International Journal of Engineering Research and Development e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Volume 13, Issue 7 (July 217), PP.75-79 Short-Term Load Forecasting Using ARIMA Model For
More informationDEPARTMENT OF ECONOMICS
ISSN 0819-64 ISBN 0 7340 616 1 THE UNIVERSITY OF MELBOURNE DEPARTMENT OF ECONOMICS RESEARCH PAPER NUMBER 959 FEBRUARY 006 TESTING FOR RATE-DEPENDENCE AND ASYMMETRY IN INFLATION UNCERTAINTY: EVIDENCE FROM
More informationTheoretical and Simulation-guided Exploration of the AR(1) Model
Theoretical and Simulation-guided Exploration of the AR() Model Overview: Section : Motivation Section : Expectation A: Theory B: Simulation Section : Variance A: Theory B: Simulation Section : ACF A:
More informationThe autocorrelation and autocovariance functions - helpful tools in the modelling problem
The autocorrelation and autocovariance functions - helpful tools in the modelling problem J. Nowicka-Zagrajek A. Wy lomańska Institute of Mathematics and Computer Science Wroc law University of Technology,
More informationLATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION
LATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION BEZRUCKO Aleksandrs, (LV) Abstract: The target goal of this work is to develop a methodology of forecasting Latvian GDP using ARMA (AutoRegressive-Moving-Average)
More informationTime Series Analysis of United States of America Crude Oil and Petroleum Products Importations from Saudi Arabia
International Journal of Applied Science and Technology Vol. 5, No. 5; October 2015 Time Series Analysis of United States of America Crude Oil and Petroleum Products Importations from Saudi Arabia Olayan
More informationA Data-Driven Model for Software Reliability Prediction
A Data-Driven Model for Software Reliability Prediction Author: Jung-Hua Lo IEEE International Conference on Granular Computing (2012) Young Taek Kim KAIST SE Lab. 9/4/2013 Contents Introduction Background
More informationTesting for Regime Switching in Singaporean Business Cycles
Testing for Regime Switching in Singaporean Business Cycles Robert Breunig School of Economics Faculty of Economics and Commerce Australian National University and Alison Stegman Research School of Pacific
More informationRoss Bettinger, Analytical Consultant, Seattle, WA
ABSTRACT DYNAMIC REGRESSION IN ARIMA MODELING Ross Bettinger, Analytical Consultant, Seattle, WA Box-Jenkins time series models that contain exogenous predictor variables are called dynamic regression
More informationExample: Baseball winning % for NL West pennant winner
(38) Likelihood based statistics ñ ln( _) = log of likelihood: big is good. ñ AIC = -2 ln( _) + 2(p+q) small is good Akaike's Information Criterion ñ SBC = -2 ln( _) + (p+q) ln(n) Schwartz's Bayesian Criterion
More informationTransformations for variance stabilization
FORECASTING USING R Transformations for variance stabilization Rob Hyndman Author, forecast Variance stabilization If the data show increasing variation as the level of the series increases, then a transformation
More informationChapter 3: Regression Methods for Trends
Chapter 3: Regression Methods for Trends Time series exhibiting trends over time have a mean function that is some simple function (not necessarily constant) of time. The example random walk graph from
More informationA Course in Time Series Analysis
A Course in Time Series Analysis Edited by DANIEL PENA Universidad Carlos III de Madrid GEORGE C. TIAO University of Chicago RUEY S. TSAY University of Chicago A Wiley-Interscience Publication JOHN WILEY
More informationDo Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods
Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods Robert V. Breunig Centre for Economic Policy Research, Research School of Social Sciences and School of
More informationForecasting. Simon Shaw 2005/06 Semester II
Forecasting Simon Shaw s.c.shaw@maths.bath.ac.uk 2005/06 Semester II 1 Introduction A critical aspect of managing any business is planning for the future. events is called forecasting. Predicting future
More informationVolatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
More informationTechnical Appendix-3-Regime asymmetric STAR modeling and exchange rate reversion
Technical Appendix-3-Regime asymmetric STAR modeling and exchange rate reversion Mario Cerrato*, Hyunsok Kim* and Ronald MacDonald** 1 University of Glasgow, Department of Economics, Adam Smith building.
More informationThe Identification of ARIMA Models
APPENDIX 4 The Identification of ARIMA Models As we have established in a previous lecture, there is a one-to-one correspondence between the parameters of an ARMA(p, q) model, including the variance of
More information