SEASONAL ADJUSTMENT VERSUS SEASONALITY MODELING: Effect on Tourism Demand Forecasting

Similar documents
Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.

Time series Decomposition method

LONG MEMORY AT THE LONG-RUN AND THE SEASONAL MONTHLY FREQUENCIES IN THE US MONEY STOCK. Guglielmo Maria Caporale. Brunel University, London

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

How to Deal with Structural Breaks in Practical Cointegration Analysis

Summer Term Albert-Ludwigs-Universität Freiburg Empirische Forschung und Okonometrie. Time Series Analysis

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

Ready for euro? Empirical study of the actual monetary policy independence in Poland VECM modelling

Types of Exponential Smoothing Methods. Simple Exponential Smoothing. Simple Exponential Smoothing

OBJECTIVES OF TIME SERIES ANALYSIS

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Quarterly ice cream sales are high each summer, and the series tends to repeat itself each year, so that the seasonal period is 4.

Chapter 16. Regression with Time Series Data

Section 4 NABE ASTEF 232

Cointegration and Implications for Forecasting

Forecasting optimally

Box-Jenkins Modelling of Nigerian Stock Prices Data

Lecture 3: Exponential Smoothing

Exponential Smoothing

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

Tourism forecasting using conditional volatility models

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN

DYNAMIC ECONOMETRIC MODELS vol NICHOLAS COPERNICUS UNIVERSITY - TORUŃ Józef Stawicki and Joanna Górka Nicholas Copernicus University

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

The electromagnetic interference in case of onboard navy ships computers - a new approach

NCSS Statistical Software. , contains a periodic (cyclic) component. A natural model of the periodic component would be

Nonstationarity-Integrated Models. Time Series Analysis Dr. Sevtap Kestel 1

Department of Economics East Carolina University Greenville, NC Phone: Fax:

Exercise: Building an Error Correction Model of Private Consumption. Part II Testing for Cointegration 1

3.1 More on model selection

A unit root test based on smooth transitions and nonlinear adjustment

Time Series Test of Nonlinear Convergence and Transitional Dynamics. Terence Tai-Leung Chong

Dynamic Econometric Models: Y t = + 0 X t + 1 X t X t k X t-k + e t. A. Autoregressive Model:

Properties of Autocorrelated Processes Economics 30331

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

THE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Nonlinearity Test on Time Series Data

Chapter 15. Time Series: Descriptive Analyses, Models, and Forecasting

DEPARTMENT OF STATISTICS

Lecture Notes 2. The Hilbert Space Approach to Time Series

Testing for a Single Factor Model in the Multivariate State Space Framework

ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

Stationary Time Series

14 Autoregressive Moving Average Models

EUROINDICATORS WORKING GROUP. A new method for assessing direct versus indirect adjustment

Vehicle Arrival Models : Headway

Dynamic Models, Autocorrelation and Forecasting

STATE-SPACE MODELLING. A mass balance across the tank gives:

Chapter 3, Part IV: The Box-Jenkins Approach to Model Building

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

Unit Root Time Series. Univariate random walk

Professorial Chair Lecture. Don Santiago Syjuco Distinguished Professorial Chair in Economics

Frequency independent automatic input variable selection for neural networks for forecasting

Morning Time: 1 hour 30 minutes Additional materials (enclosed):

Why is Chinese Provincial Output Diverging? Joakim Westerlund, University of Gothenburg David Edgerton, Lund University Sonja Opper, Lund University

International Parity Relations between Poland and Germany: A Cointegrated VAR Approach

Robust estimation based on the first- and third-moment restrictions of the power transformation model

A Point Optimal Test for the Null of Near Integration. A. Aznar and M. I. Ayuda 1. University of Zaragoza

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

New innovative 3-way ANOVA a-priori test for direct vs. indirect approach in Seasonal Adjustment

Forecast of Adult Literacy in Sudan

Wavelet Variance, Covariance and Correlation Analysis of BSE and NSE Indexes Financial Time Series

A multivariate labour market model in the Czech Republic 1. Jana Hanclová Faculty of Economics, VŠB-Technical University Ostrava

Modeling Rainfall in Dhaka Division of Bangladesh Using Time Series Analysis.

Distribution of Least Squares

Exponentially Weighted Moving Average (EWMA) Chart Based on Six Delta Initiatives

Use of Unobserved Components Model for Forecasting Non-stationary Time Series: A Case of Annual National Coconut Production in Sri Lanka

Fourier Transformation on Model Fitting for Pakistan Inflation Rate

Generalized Least Squares

Arima Fit to Nigerian Unemployment Data

Long-Term Demand Prediction using Long-Run Equilibrium Relationship of Intrinsic Time-Scale Decomposition Components

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

Distribution of Estimates

Mean Reversion of Balance of Payments GEvidence from Sequential Trend Break Unit Root Tests. Abstract

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Problem Set 5. Graduate Macro II, Spring 2017 The University of Notre Dame Professor Sims

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA

Solutions to Odd Number Exercises in Chapter 6

Økonomisk Kandidateksamen 2005(II) Econometrics 2. Solution

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

"FORECASTING OF THE RAINFALL AND THE DISCHARGE OF THE NAMORONA RIVER IN VOHIPARARA AND FFT ANALYSES OF THESE DATA

A New Unit Root Test against Asymmetric ESTAR Nonlinearity with Smooth Breaks

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T

ACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin

PROC NLP Approach for Optimal Exponential Smoothing Srihari Jaganathan, Cognizant Technology Solutions, Newbury Park, CA.

Comparing Means: t-tests for One Sample & Two Related Samples

STAD57 Time Series Analysis. Lecture 17

STAD57 Time Series Analysis. Lecture 17

Linear Dynamic Models

12: AUTOREGRESSIVE AND MOVING AVERAGE PROCESSES IN DISCRETE TIME. Σ j =

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

Choice of Spectral Density Estimator in Ng-Perron Test: A Comparative Analysis

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

Transcription:

Seasonal Asian-African Adjusmen Journal Versus of Economics Seasonaliy and Modeling: Economerics, Effec Vol. on Tourism 13, No. 1, Demand 2013: 71-84 Forecasing 71 SEASONAL ADJUSTMENT VERSUS SEASONALITY MODELING: Effec on Tourism Demand Forecasing Amira GASMI * ABSTRACT In his sudy, we rea he seasonal variaion in monhly ime series in he conex of he Wesern- European ourism demand for Tunisia, by presening differen echniques of deecion of seasonaliy and he parameric and non-parameric approaches of seasonal adjusmen. Then, we compare he forecasing performance of hese mehods. The empirical resuls miliae in favour of he TRAMO-SEATS mehod. In fac, his approach provides he bes forecas. In erms of forecasing efficiency, we noe in addiion, ha he modelling of he seasonal variaion using seasonal ARIMA model (SARIMA) may lead o beer predicive resuls compared wih oher echniques of seasonal adjusmen used in his research, namely: he X-12-ARIMA, regression on seasonal dummies and he raio-o-moving average mehods. Keywords: seasonaliy, ourism demand, forecasing performance, seasonal adjusmen, seasonal modelling. JEL Classificaion: C22, C52, C53, L83. 1. INTRODUCTION Seasonaliy is a major characerisic of he ourism aciviy. I reveals he influence of he seasons on he ourism demand. This phenomenon is relaed o weaher changes as well as insiuional facors (school holidays, professional vacaion, public (Chrismas or Easer), religious and commemoraive fesivals). The calendar can also generae a seasonal movemen in monhly ime series since he number of working days varies from one monh o anoher. I is also relaed o cerain socio-culural characerisics (spor pracices; social or religious habis). Taking ha ino accoun, he phenomenon of high and low seasons consiues a problem of size which worries he acors of he ourism field. Time series analysis aims o separae he shor-erm behaviour from ha of long-erm of an economic daa series and o give reliable forecass for hese separae componens and for he oaliy of he series. The seasonal variaions explain mos of he variaion in he growh raes of he majoriy of he economic ime series. In order o draw conclusions on he naure of he business cycles and * Laboraoire d Economie e Finance Appliquees IHEC Carhage, Tunisia, E-mail: amira.gasmi@yahoo.fr

72 Amira Gasmi he long-erm growh, he radiional approach is o remove he seasonal componen of a series hrough he use of deerminisic seasonal dummies, seasonal differeniaion, or using he seasonal adjusmen echniques such as he X-12-ARIMA mehod. However, someimes we show ha i can be more appropriae o sudy he seasonal models hemselves (Lee and Siklos (1993), Reimers (1997)) since hey could give informaion on he behaviour of he economic agens which are exposed o changes of endencies a he momen of planning and he formaion of waiings. Thus, alhough seasonal variaions-correced daa can be useful, i is ypically recommended o use no adjused daa. Moreover, several recen empirical sudies showed ha many mehods of seasonal adjusmen lead o seriously denaured daa, in he sense ha he key properies such as he endencies, he business cycles, he nonlineariies are affeced by he seasonal adjusmen (Ghysels and Perron (1993), Ghysels, Granger and Siklos (1995), Hylleberg (1994), Miron (1996), Maravall (1995)). In conras, King and Kulendran (1997) evaluae several models, including he seasonal uni roos model, in he forecas of quarerly ouris arrivals in Ausralia coming from many counries. Their principal conclusion is ha compared o ime series models, he forecasing performance of he seasonal uni roos models is weak. This may be due o he lack of power of some uni roos ess. On he oher hand, Paap, Franses and Hoek (1997) use empirical and simulaion examples o demonsrae ha he negleced seasonal average changes can desroy considerably he forecasing performance of he univariae auoregressive processes. Thus, appropriae reamen of seasonaliy is imporan o make reliable forecass. In his aricle, we propose firsly, o sudy he seasonal aspec srongly characerizing he ourism ime series, by presening differen ess of seasonaliy deecion and various mehods of reamen of seasonaliy, in paricular seasonal adjusmen mehods versus seasonaliy modeling. Then, secondly, we will compare he forecasing performance of hese mehods. 2. ANALYSIS OF SEASONALITY 2.1. Deecion of Seasonaliy During he analysis of a ime series, i is necessary o idenify he seasonal variaion which can be probably observed. Various ypes of ess are se up o deec he presence of his componen. 2.1.1. Auocorrelaions Seasonaliy can be deeced graphically by examining auocorrelaion (ACF) and parial auocorrelaion funcions (PACF) necessary for he idenificaion of suiable ARIMA models. Indeed, he correlogram of a seasonal series ofen akes a sinusoidal form (see able 1). 2.1.2. Tradiional Tess of Presence of Seasonaliy To es he presence of seasonal variaion, a mulipliciy of ess were suggesed, namely: he sable seasonaliy and moving seasonaliy ess which are Fisher ypes ess based on models of analysis of he variance o one (he monh or he quarer) and wo facors (he monh or he quarer and he year), respecively. Indeed, sable seasonaliy is a ype of seasonaliy which is repeaed a he same ime each year, and his sable aspec faciliaes he forecass. While he

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 73 moving seasonaliy is represened by a movemen effec from one monh o anoher. This lack of sabiliy makes is forecas difficul. Lasly, we disinguish he idenifiable seasonaliy es compleing he ess evoked above. I is buil saring from he values of Fisher saisics of sable and moving seasonaliy ess (Lohian and Morry, 1978). The es saisic, noed T, is expressed as follows: 1/ 2 1 2 7 T wih T1 and T2 2 FS T T 3F F If saisic T is lower han 1 so we concludes he presence of idenifiable seasonal componen and he seasonal adjusmen of he series is hen necessary. 2.1.3. Seasonal Uni Roos Tess The deecion of seasonal variaion can be done using seasonal uni roos ess. For his purpose, a cerain number of ess were implemened in he eighies and nineies, in paricular o es he seasonal variaion a order 4 and order 12. Table 1: Correlogram of Wesern-European Touris Arrivals Series M S.

74 Amira Gasmi Tes DHF (Dickey, Hasza and Fuller, 1984): i allows o es he null assumpion = 1 in he model x = x S +. Under H 0 rue, he series is seasonal and he filer S = (1 L S ) suggesed by Box & Jenkins (1970) is appropriae o adjus i (S being he period of seasonaliy). Tes HEGY (Hylleberg, Engle, Granger and Yoo, 1990): he lieraure on he concep of unis roos (e.g., Dickey, Bell and Miller (1986)) shows ha he assumpion of exisence of cerain filers of differeniaion amouns o emi he assumpion of presence of a cerain number of seasonal and non-seasonal unis roos in a ime series. This can be easily seen by wriing: S = (1 L S ), and by solving he equaion: (1 z S ) = 0. The general soluion o his equaion is: {1, cos(2k / S) + i sin(2k / S); wih he erm (2k / s) for k = 1, 2,, represens he corresponding seasonal frequency, giving S differen soluions which all of hem are on he circle uni. The HEGY mehod, mainly conceived for he quarerly series, was adaped o he monhly case hanks o Franses (1990) and of Beaulieu and Miron (1993) works. In fac, if S = 12, soluions of he equaion (1 z 12 ) = 0 are: 1 for he non-seasonal uni roo corresponding o frequency 0; and 11 seasonal uni roos 1 1 1 1 1, i, (1 3 i), (1 3 i), ( 3 i), ( 3 i) corresponding respecively o he 2 2 2 2 2 5 following frequencies:,,,,, and o he operaors of differeniaion: 2 3 3 6 6 (1+ L), (1+ L 2 ), (1+ L+ L 2 ), (1 L+ L 2 ), (1+ 3 L+ L 2 ) and (1 3 L+ L 2 ). S 1 Thus, a filer of differeniaion ( S ) can be wrien as: s (1 L)(1 L L ), and can hus be decomposed ino a par wih non-seasonal uni roo and a par wih (S 1) seasonal uni roos. In his es, we resor o he decomposiion of he polynomial (1-L 12 ), wih 12 roos unis and we consider he following form: ( L) z z z z z z z z z 8 1 1, 1 2 2, 1 3 3, 1 4 3, 2 5 4, 1 6 4, 2 7 5, 1 8 5, 2 z z z z 9 6, 1 10 6, 2 11 7, 1 12 7, 2 The variables z i are in such a way ha: z i = P i (L) y, where polynomials P i are defined as follows : 2 4 8 P1 ( L) (1 L)(1 L )(1 L L ) 2 4 8 P2 ( L) (1 L)(1 L )(1 L L ) 2 4 8 P3 ( L) (1 L )(1 L L ) P L L L L L L P L L L L L L 4 2 4 2 P6 ( L) (1 L )(1 L L )(1 L L ) 4 2 4 2 P7 ( L) (1 L )(1 L L )(1 L L ) 12 P8 ( L) (1 L ) 4 2 2 4 4( ) (1 )(1 3 )(1 ) 4 2 2 4 5( ) (1 )(1 3 )(1 ) (1)

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 75 Wih also: (L) is an auoregressive polynomial in L, and µ may conain a consan, 11 seasonal dummies and/or a rend. The variables z i are hen associaed o he differen roos of he polynomial. The equaion (1) is esimaed using leas squares ordinary mehod. We may carry ou ess for he parameers 1 and 2, and F ess associaed o he couples ( 3, 4 ), ( 5, 6 ), ( 7, 8 ), ( 9, 10 ) and ( 11, 12 ) : i is a quesion of esing he joined significance of he coefficiens. This amouns o es he assumpion of exisence of uni roos a he differen frequencies. For his purpose, we mus compare he es saisics relaed o he esimaed parameers wih he criical values provided by Franses (1990) and Beaulieu and Miron (1993). To check he exisence of he roos 1 and -1 corresponding o frequencies 0 and respecively, we carry ou wo individual ess on parameers 1 and 2. As for he oher seasonal uni roos, we can perform eiher joined ess whose null assumpion akes he form k = k 1 = 0 and his for he even values of k, from 4 o 12; or quie simply, individual ess, suggesed in Franses (1990), allowing o verify he non-saionariy of he ime series a all he seasonal frequencies and his by esing he null assumpion according o which here is a seasonal uni roo ( k = 0, k [3, 12]). However, i should be noed ha he applicaion of he OLS o he regression (1) is made where he order of (L) is given in such a way ha he errors are roughly whie noises, or a leas, non-auocorrelaed residuals. For his purpose, Hylleberg and al. (1990) and Engle and al. (1993) propose o inroduce addiional lags of he variable unil we obain non-auocorrelaed residuals. 2.2. Seasonal Adjusmen Mehods Seasonal adjusmen mehods can be classified in wo caegories, namely: parameric approach and nonparameric approach. 2.2.1. Non-parameric Approaches The X-12-ARIMA mehod: when he seasonal variaion is very apparen in he ime series, a firs approach consiss in removing such seasonal flucuaions by using a seasonal adjusmen programs. They are echniques allowing he idenificaion of he differen componens of he iniial series (rend-cycle, seasonaliy, irregular) by applying linear filers, which cancels or preserves a well defined componen (endency-cycle or seasonal variaion). The irregular one is represened hereafer by he residual of he decomposiion. These linear filers are moving averages which consiue he principal ool of Census X- 11 mehod buil from successive ieraions of moving averages of differen naures for beer esimaing he series componens. However, his echnique leads o a loss of informaion in he final end of he series. This gap is filled by he forecas of fuure values of he ime series before is seasonal adjusmen, and his using ARIMA models. I is wha made i possible o exend he X-11 echnique (Census Bureau, 1967) o X-11-ARIMA (Dagum, 1988) and hen o X-12-ARIMA (Findlay and Al, 1998). The laer conains he RegARIMA module which allows o deec and o remove any undesirable effec of he series (ouliers, calendar effecs ).

76 Amira Gasmi The raio o moving average mehod: Monhly values of he sudied series (X ) are divided by he moving average figure corresponding for each monh (MA ), and expressed in % o X generae he raio-o-moving average: Mraio 100. The moving average is calculaed MA as follows: 1 M X 2 X 1... 2X 1 2X 2 X 1... 2 X 1 X, m 2 p m 2m p p p p 1 MA M ( X ) X 2X 2 X... 2X X 24 12 6 5 4 5 6 A. Carpener (2003) reveals ha his moving average eliminaes seasonal variaions from monhly series, preserves linear rends and reduces of more han 90% he variance of a whie noise. These raios are weighed by he monh and hereafer will separae he seasonal and cyclic componens. 2.2.2. Parameric Approaches The regression mehod: his approach is based on he Buys-Ballo model (1847) which consiss in carrying ou he regression below, using seasonal dummies (S,i ) in such a way ha S,i akes value 1 if T corresponds o he seasonal period, and 0 if no. The model is wrien as follows: T 1 X S. 0 1 i, i i1 Wih: T being he period of seasonaliy (T = 4 for a quarerly series, T = 12 for monhly daa). The use of only (T 1) dummies makes i possible o avoid he problem of colineariy which could exis wih he vecor uni relaing o he consan. We esimae hus (T 1) seasonal coefficiens and we check he T h using he principle of conservaion of he surfaces : i 0. Seasonal adjusmen by mehod TRAMO-SEATS: TRAMO-SEATS program (Gomez & Maravall, 1996) belongs o he parameric seasonal adjusmen mehods based on he signal exracion. I is composed of wo independen subrouines bu which are complemenary since hey are generally used ogeher: TRAMO program (Time series Regression wih ARIMA noise, Missing observaions and Ouliers) falls under he same opic as ARIMA modelling, or more exacly, i is abou an exension o hese models. Is principle is in fac o model he iniial series using he univariae approach of Box & Jenkins via ARIMA or seasonal ARIMA (SARIMA) models, while deecing, esimaing and correcing as a preliminary he ouliers, he missing values, he calendar effecs (holidays, public holidays ) as well as srucural changes, likely o disurb he esimaion of he model coefficiens. T i1

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 77 SEATS program (Signal Exracion ARIMA Time Series) comes o complee TRAMO procedure by decomposing he iniial series hus modelled in is componens (rend, cycle, irregular and seasonaliy) by signal exracion, using he specral analysis of he iniial series. 2.3. Seasonal Differeniaion and Seasonaliy Modelling The verificaion of he exisence of seasonal uni roos using specific ess such as DHF (1984) and HEGY (1990) requires special reamen of seasonaliy. The use of he filer (1 L s ), suggesed by Box and Jenkins (1970), o differeniae he seasonal series, depends on he fac ha he variable is non-saionary a frequency 0 and a all he seasonal frequencies (Pichery and Ouerfelli, 1998). The exisence of seasonal uni roos leads o model he seasonal variaion insead of correcing or removing i using seasonal adjusmen mehods. The mos largely used seasonal model is he muliplicaive seasonal ARIMA model or SARIMA(p,d,q)(P,D,Q)S proposed by Box & Jenkins (1970) as a generalizaion of ARIMA(p,d,q) models conaining a seasonal par and which is wrien in his form: ( L) ( L s ) d D y ( L) ( L s ) (2) p P s q Q Where: S is he period of seasonaliy (S = 12 for monhly daa, S = 4 for quarerly daa); = 1 L, S = 1 L S, p, P, q, Q are polynomials of degrees: p, P, q, Q and he roos are of module higher han 1; ( ) is a whie noise; d and D are respecively he orders of non seasonal and seasonal differeniaion. 3. THE DATA The Wesern-European marke being he principal marke ransmiing ouriss owards Tunisia, he empirical applicaion is carried ou based on he series of he Wesern- European ouris arrivals in Tunisia ransformed o logarihm and subsequenly noed LTOEU. The sample covers he period from January 1997 o December 2009. For he esimaion, we use he daa beween January 1997 and June 2009, he six remaining observaions are used for he ex-pos forecas and for he predicive performance evaluaion of he various mehods. Daa are provided by he Naional office of Tunisian Tourism relaing o he minisry of ourism and rade. Table 2 Resuls of Seasonaliy Tess F S F M T Value 304,959 9,797 0,2439 Decision Presence of sable Absence of moving Presence of idenifiable seasonaliy seasonaliy seasonaliy

78 Amira Gasmi Table 3 Resul of Tes DHF (1984) Wesern-European ouris arrivals in Tunisia Tes saisic Level 5% Decision Filer 1,166-5,84 Accepe H 0 (1- L 12 )] Table 4 Resul of Tes HEGY (1990) Wesern-European ouris arrivals in Tunisia Frequency Tes saisic Level 5% Decision Filer o -0,795-2,76 Accepe H 0 (1- L) / 6 3,187-1,85 Accepe H 0 (1 2 3 L L ) / 3 3,21-3,25 Accepe H 0 (1 L + L 2 ) / 2 6,080-3,25 Accepe H 0 (1 + L 2 ) 2/ 3-1,783-1,85 Accepe H 0 (1+ L + L 2 ) 5/ 6 1,61-1,85 Accepe H 0 (1 2 3 L L ) -4,027-2,76 Reje de H 0-4. EMPIRICAL RESULTS Deecion of he Seasonaliy The presence of seasonal variaion noed graphically in able 1 is confirmed hanks o he resuls of he combined es which indicaes he presence of an idenifiable seasonal variaion, since he es saisic provides a value lower han 1 (see able 2). This is marked hanks o he resuls of es HEGY presened in able 4. In fac, by using he comparison of he T-saisic calculaed in he able wih he criical values provided in Beaulieu and Miron (1993), his es reveals he presence of he non-seasonal uni roo 1 corresponding o he zero frequency. This allows us o conclude of he nonsaionariy of he variable. Hence, is differeniaion wih he filer (1 L) is required. Furhermore, he es leads o he accepance of he assumpion H 0 of presence of uni roos a all he seasonal frequencies, excep for he frequency. Consequenly, he produc of he filers indicaed in able 4 mus be applied o eliminae he seasonal and non-seasonal uni roos, ha is o say: (1 L) (1 + L 2 + L 4 + L 6 + L 8 + L 10 ). Taking ha ino accoun, we can conclude ha he suiabiliy of he applicaion of he filer (1 L 12 ) o a seasonal series, as i is recommended by Box & Jenkins (1970), depends on he fac ha he series is inegraed a he seasonal frequency zero and a all frequencies. This being, hese resuls imply ha he auomaic applicaion of he filer of seasonal differeniaion is likely o produce a specificaion error. The proof presened here indicaes ha he uni roos are someimes missing a cerain seasonal frequencies, hen heir presence have o be checked by using he es HEGY, raher han o impose hem a priori a all he frequencies.

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 79 However, and by conrariey of simplificaion, and aking ino accoun he exisence of only one seasonal frequency where he assumpion H 0 is rejeced, we have preferred he applicaion of he filer of seasonal differeniaion (1 L 12 ) suggesed by Box & Jenkins (1970) and recommended by he es of Dikey, Hasza and Fuller (1984) whose resul arises in able 3. Comparison of he Seasonal Adjusmen Mehods Figures 1, 2 and 3 presen he series of he Wesern-European ouris arrivals adjused by he differen seasonal adjusmen mehods considered in his sudy. We propose o compare he forecasing performance. For his purpose, we followed he forecas process of Box & Jenkins (1970) and he seps of idenificaions, esimaion and validaion enabled us o reain he following forecasing models: ARIMA (2,1,2), SARMA (1,1) (1,1,1) 12, ARMED (1,1), ARIMA (2,1,2), ARIMA (2,1,1) and ARMA(1,1) for each one of hese mehods of reamen of he seasonal variaion, respecively: he filer of seasonal differeniaion (1-L 12 ) suggesed by es DHF (forecass 1), he X-12-ARIMA mehod (forecass 2), he raio-o-moving average echnique (forecass 3), he regression on seasonal dummies (forecass 4) and he TRAMO-SEATS program (forecass 5). To compare he forecasing efficiency of hese models, we reained various crieria of evaluaion of he predicive precision, namely: he MAPE, he RMSE, he RMSPE and he U- Theil inequaliy coefficien. The reading of able 5 makes i possible o conclude ha overall (six-monhs-ahead horizon), he TRAMO-SEATS seasonal adjusmen mehod allows o obain he mos precise forecass since hey admi he weakes evaluaion crieria, followed by he seasonal model SARIMA (second rank) and he X-12-ARIMA mehod (hird rank). Therefore, modelling seasonaliy by he recourse o he SARIMA model (applicaion of he filer of seasonal differeniaion (1-L 12 )) is more advised in erms of forecasing efficiency han he seasonal adjusmen by he X-12-ARIMA, he raio-o-moving average and he regression on seasonal dummies mehods. Figure 1: Non-parameric Seasonal Adjusmen Approaches

80 Amira Gasmi Figure 2: Seasonal Adjusmen wih Seasonal Dummies Figure 3: Seasonal Adjusmen wih he TRAMO-SEATS Mehod

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 81 Table 5 Forecasing Performance of Seasonal Adjusmen Mehods Horizons * MAPE RMSE RMSPE U-Theil one monh (2) 0,7226 819,4837 0,7226 0,3626% Forecass 1 2 monhs (2) 2,4325 3772,6082 2,7568 1,5085% (filer (1 L 12 )) 3 monhs (3) 2,7088 5130,5156 2,7411 1,2751% 6 monhs (2) 2,1184 9516,7616 3,0800 1,4653% one monh(3) 0,7527 853,5948 0,7527 0,3749% 2 monhs (3) 2,4037 3797,8214 2,7680 1,5184% Forecass 2 3 monhs (2) 1,6840 3501,5787 2,0260 0,8725% (X-12-ARIMA) 6 monhs (3) 4,0014 15139,9351 4,5367 2,3034% one monh(4) 2,4143 2738 2,4143 1,1928% Forecass 3 2 monhs (4) 3,6977 4593,102 3,7609 1,7802% (raio-o-moving 3 monhs (4) 8,2109 15714,649 8,2483 3,7769% average) 6 monhs (5) 12,8957 32761,35 14,8750 4,9032% Forecass 4 one monh(5) 4,6037 5220,842 4,6037 2,2500% (regression on seasonal) 2 monhs (5) 4,2033 5371,418 4,5867 2,0788% dummies mehods) 3 monhs (5) 9,2562 18323,210 9,3456 4,3768% 6 monhs (4) 11,6240 28989,436 13,418 4,3701% one monh(1) 0,05674 64,3498 0,05674 0,0284% Forecass 5 2 monhs (1) 1,6191 2219,9875 1,6724 0,8833% (Tramo-Seas) 3 monhs (1) 1,0769 1777,527 1,3353 0,4448% 6 monhs (1) 1,2054 4719,821 1,4008 0,7291% MAPE: Mean Absolue Percenage Error; RMSE: Roo Mean Square Error; RMSPE: Roo Mean Square Percenage Error. (*): Figures in blue represen he forecass order by horizon.

82 Amira Gasmi This order is mainained for one-monh and wo-monh-ahead horizons. On he conrary, for he hree-monh-ahead horizon, he forecass resuling from he X-12-ARIMA mehod become beer han hose obained using he SARIMA model. In consequence, he forecasing performance of he various mehods can vary according o he horizon of forecas, which corroboraes wih he resuls found in preceding sudies (Wong and al., 2007; Shen and al., 2009; Shen and al., 2011). By elsewhere, he empirical evidence suggess ha he echniques of reamen of he seasonal variaion affec he forecasing performance of he models, and ha differs according o sochasic or deerminisic naure of he seasonal variaion. In effec, he resuls obained in his empirical exercise reveal ha he bes forecass resul from he TRAMO-SEATS and he X-12-ARIMA mehods and also from he seasonal ARIMA model which consider he sochasic seasonal variaion (Bourbonnais and Terraza, 2008). 5. CONCLUSION In his paper, we applied four seasonal adjusmen mehods: wo parameric mehods (TRAMO- SEATS and regression on seasonal dummies) and wo non-parameric ones (he X-12-ARIMA and he raio-o-moving average), o a monhly series represening he Wesern-European ouris arrivals in Tunisia. We compared he forecasing performance of hese mehods in paricular, seasonal adjusmen versus seasonaliy modelling. The obained resuls miliae in favour of he TRAMO-SEATS mehod. In fac, his approach provides he bes forecas a all he forecas horizons. Always in erms of forecasing performance, we have been able o noe ha he seasonaliy modelling using seasonal ARIMA (SARIMA) models may lead o beer predicive resuls compared wih he oher echniques of seasonal adjusmen, namely: he X-12-ARIMA, he regression on seasonal dummies and he raio-o-moving average. Consequenly, i could be someimes more appropriae o model he seasonal variaion raher han o resor o is correcion or suppression by he means of seasonal adjusmen mehods. Anoher conclusion ha we could draw from he resuls is ha he forecasing performance is influenced by he manner wih which he seasonal variaion is reaed in he series, i.e. i differs according o sochasic or deerminisic naure of seasonaliy. Indeed, he empirical resuls reveal ha he bes predicive performance rises from he TRAMO-SEATS program, he X-12-ARIMA mehod and of he SARIMA model which consider he sochasic seasonaliy (Bourbonnais and Terraza, 2008). This is on line wih oher researches which sugges he sochasic reamen of he seasonal variaion (for example, Shen and Al, 2009). Acknowledgemens I would like o hank Mr. Regis Bourbonnais (Paris-Dauphine Universiy) for his consrucive commens and helpful suggesions hroughou his work.

Seasonal Adjusmen Versus Seasonaliy Modeling: Effec on Tourism Demand Forecasing 83 References Beaulieu, J. J. & Miron, J. A., (1991), The Seasonal Cycle in U.S. Manufacuring. Economics Leers, 37, 115 118. Beaulieu, J. J. & Miron, J. A., (1993), Seasonal Uni Roos in Aggregae U.S. Daa. Journal of Economerics, 55, 305 328. Bell, W. R. & Hillmer, S. C., (1984), Issues Involved wih he Seasonal Adjusmen of Economic Time Series. Journal of Business and Economic Saisics, 2, 291 320. Bourbonnais, R. & Terraza, M., (2008), Analyse des séries Temporelles. Dunod. Box, G. E. P. & Jenkins, G. M., (1970), Time Series Analysis, Forecasing and Conrol. San Francisco: Holden- Day. Charpenier, A., (2003), Inroducion à la héorie des processus en emps discre. Modèles ARIMA e Méhode Box-Jenkins. Universié Paris Dauphine. Clemens, M. P. & Hendry, D. F., (1997), An Empirical Sudy of Uni Roos in Forecasing. Inernaional Journal of Forecasing, 13, 341 355. Dagum, E. B., (1988), The X11ARIMA/88 Seasonal Adjusmen Mehod - Foundaions and User s Manual, Time Series Research and Analysis Division. Saisics Canada echnical repor. Darné, O. & Terraza, M., (2002), La désaisonnalisaion des séries emporelles viicoles. Version préliminaire. Colloque d oenomérie. Darné, O., Liago, J. & Terraza, M., (2002), Tess de racines uniaires saisonnières pour des données journalières. Revue de Saisique appliquée, 50, 71 91. Dickey, D. A., D. P. Hasza, & Fuller, W. A., (1984), Tesing for Uni Roos in Seasonal Time Series, Journal of he American Saisical Associaion, 79, 355 367. Dickey, D. A., Bell, W. R. & Miller, R. B., (1986), Uni Roos in Time Series Models : Tess and Implicaions, The American Saisician, 40, p. 12-26. Engle, R. F., Granger, C. W. J. & Hallman, J. J., (1989), Merging Shor- and Long-run Forecas. An Applicaion of Seasonal Coinegraion o Monhly Elecriciy Sales Forecasing. Journal of Economerics, 40, 45 62. Engle, R. F., Granger, C. W. J., Hylleberg, S. & Lee, H. S., (1993), Seasonal Coinegraion: he Japanese Consumpion Funcion. Journal of Economerics, 55, 275 298. Findley, D. F., B. C. Monsell, W. R. Bell, Oo, M. C. & Chen, B. C., (1998), New Capabiliies of he X-12- ARIMA Seasonal Adjusmen Program (wih Discussion). Journal of Business and Economic Saisics, 16, 127-177. Franses, P. H., (1991), Seasonaliy, Nonsaionariy and he Forecasing of Monhly Time Series. Inernaional Journal of Forecasing, 7, 199 208. Franses, P. H., (1996), Recen Advances in Modelling Seasonaliy. Journal of Economic Surveys, 10, 299 345. Franses, P. H. & Hobijn, B., (1997), Criical Values for Uni Roo Tess in Seasonal Time Series. Journal of Applied Saisics, 24, 25 47. Ghysels, E., (1994), On he Economics and Economerics of Seasonaliy. In: Sims, C.A. (Ed.), Advances in Economerics, Sixh World Congress of he Economeric Sociey. Cambridge Universiy Press, Cambridge. Ghysels, E., (1994), L analyse économérique e la saisonnalié. L acualié économique, 70, 43 62. Ghysels, E., Lee, H.S. & Noh, J., (1994), Tesing for Uni Roos in Seasonal Time Series. Journal of Economerics, 62, 415 442.

84 Amira Gasmi Gomez, V. & Maravall. A., (1996), Programs TRAMO (Time Series Regression wih Arima Noise, Missing Observaions, and Ouliers) and SEATS (Signal Exracion in Arima Time Series). Insrucions for he User. Documen de ravail 9628 du Déparemen de recherche de la Banque d Espagne. Gooijer, J. G. & Franses, P. H., (1997), Forecasing and Seasonaliy. Inernaional Journal of Forecasing, 13, p. 303-305. Hasza, D. P. & Fuller, W. A., (1982), Tesing for Non-saionary Parameer Specificaions in Seasonal Time Series Models, Annals of Saisics, 10, 1209-1216. Herwarz, H., (1997), Performance of Periodic Error Correcion Models in Forecasing Consumpion Daa. Inernaional Journal of Forecasing, 13, p. 421-431. Hylleberg, S., Engle, R. F., Granger, C. W. J. & Yoo, B. S., (1990), Seasonal Inegraion and Coinegraion. Journal of Economerics, 44, 215 238. Hylleberg, S., (1994), Modelling Seasonal Variaion. In: Hargreaves, C. P. (Ed.), Nonsaionary Time Series Analysis and Coinegraion, Oxford Universiy Press, Oxford. Hylleberg, S. & Pagan, A. R., (1997), Seasonal Inegraion and he Evolving Seasonals Model. Inernaional Journal of Forecasing, 13, 329 340. Lee, H. S., (1992), Maximum Likelihood Inference on Coinegraion and Seasonal Coinegraion. Journal of Economerics, 54, 351 365. Lee, H. S. & Siklos, P. L., (1993), The Influence of Seasonal Adjusmen on he Canadian Consumpion Funcion, 1947 1991. Journal of Economics, 26, 575 589. Lohian, J. and Morry, M., (1978), A Tes for he Presence of idenifiable Seasonaliy when using he X-11 program. Research paper 78-10-002E, Seasonal Adjusmen and Time Series Saff, Saisics Canada. Miron, J. A., (1996), The Economics of Seasonal Cycles. MIT Press, Cambridge USA. Novales, A. & Flores, R., (1997), Forecasing wih Periodic Models. A Comparison wih Time Invarian Coefficien Models. Inernaional Journal of Forecasing, 13, 393-405. Paap, R., Franses, P. H. e Hoek, H., (1997), Mean Shifs, Uni Roos, and Forecasing Seasonal Time Series. Inernaional Journal of Forecasing, 13, 357-368. Pichery, M. C. & Ouerfelli, C., Juille, (1998), La non saionnarié dans les séries saisonnières. Applicaion au ourisme unisien. Universié de Bourgogne, Dijon. Reimers, H. E., (1997), Forecasing of Seasonal Coinegraed Processes. Inernaional Journal of Forecasing, 13, 369 380. Shen, S., Li, G., & Song, H., (2009), Effec of Seasonaliy Treamen on he Forecasing Performance of Tourism Demand Models. Tourism Economics, 15 (4), 693 708. Shen, S., Li, G., & Song, H., (2011), Combinaion Forecass of Inernaional Tourism Demand. Annals of Tourism Research, 38 (1), 72 89. Shiskin, J., Young, A. H. & Musgrave, J. C., (1967), The X-11 Varian of he Census mehod II Seasonal Adjusmen. Technical Paper No. 15, Bureau of he Census, U.S. Deparmen of Commerce. Wong, K. K. F., Song, H., Wi, S. F., & Wu, D. C., (2007), Tourism Forecasing: o combine or no o combine? Tourism Managemen, 28, 1068 1078.