Time Series Models for Growth of Urban Population in SAARC Countries

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

Econ Autocorrelation. Sanjaya DeSilva

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

Time series Decomposition method

How to Deal with Structural Breaks in Practical Cointegration Analysis

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

Section 4 NABE ASTEF 232

Box-Jenkins Modelling of Nigerian Stock Prices Data

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

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

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

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

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

OBJECTIVES OF TIME SERIES ANALYSIS

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

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

Predicting Money Multiplier in Pakistan

Cointegration and Implications for Forecasting

Fourier Transformation on Model Fitting for Pakistan Inflation Rate

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

GDP PER CAPITA IN EUROPE: TIME TRENDS AND PERSISTENCE

INVESTIGATING THE WEAK FORM EFFICIENCY OF AN EMERGING MARKET USING PARAMETRIC TESTS: EVIDENCE FROM KARACHI STOCK MARKET OF PAKISTAN

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

Stationary Time Series

Chapter 16. Regression with Time Series Data

Tourism forecasting using conditional volatility models

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

Computer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures

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

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

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

DEPARTMENT OF STATISTICS

Estimation Uncertainty

Distribution of Estimates

FITTING OF A PARTIALLY REPARAMETERIZED GOMPERTZ MODEL TO BROILER DATA

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

BOOTSTRAP PREDICTION INTERVALS FOR TIME SERIES MODELS WITH HETROSCEDASTIC ERRORS. Department of Statistics, Islamia College, Peshawar, KP, Pakistan 2

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

The Validity of the Tourism-Led Growth Hypothesis for Thailand

GINI MEAN DIFFERENCE AND EWMA CHARTS. Muhammad Riaz, Department of Statistics, Quaid-e-Azam University Islamabad,

Comparing Theoretical and Practical Solution of the First Order First Degree Ordinary Differential Equation of Population Model

A new flexible Weibull distribution

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

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

Nonlinearity Test on Time Series Data

14 Autoregressive Moving Average Models

Solutions to Odd Number Exercises in Chapter 6

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

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

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

Dynamic Models, Autocorrelation and Forecasting

anempirical Study Strictly as per the compliance and regulations of :

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

20. Applications of the Genetic-Drift Model

Forecasting optimally

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions

Regression with Time Series Data

DEPARTMENT OF ECONOMICS AND FINANCE COLLEGE OF BUSINESS AND ECONOMICS UNIVERSITY OF CANTERBURY CHRISTCHURCH, NEW ZEALAND

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

The Properties of Procedures Dealing with Uncertainty about Intercept and Deterministic Trend in Unit Root Testing

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Modeling Economic Time Series with Stochastic Linear Difference Equations

Wisconsin Unemployment Rate Forecast Revisited

Arima Fit to Nigerian Unemployment Data

A Hybrid Model for Improving. Malaysian Gold Forecast Accuracy

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

Lecture 15. Dummy variables, continued

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Hypothesis Testing in the Classical Normal Linear Regression Model. 1. Components of Hypothesis Tests

ECON 482 / WH Hong Time Series Data Analysis 1. The Nature of Time Series Data. Example of time series data (inflation and unemployment rates)

In this paper the innovations state space models (ETS) are used in series with:

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

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

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

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

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

A unit root test based on smooth transitions and nonlinear adjustment

Forecast of Adult Literacy in Sudan

Forecasting of boro rice production in Bangladesh: An ARIMA approach

Has the Business Cycle Changed? Evidence and Explanations. Appendix

Comparison Between Regression and Arima Models in Forecasting Traffic Volume

Air Traffic Forecast Empirical Research Based on the MCMC Method

Wednesday, November 7 Handout: Heteroskedasticity

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

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

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Lecture 3: Exponential Smoothing

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.

Properties of Autocorrelated Processes Economics 30331

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

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

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

DEPARTMENT OF ECONOMICS

Chickens vs. Eggs: Replicating Thurman and Fisher (1988) by Arianto A. Patunru Department of Economics, University of Indonesia 2004

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

Improved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method

Shiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran

A Quasi-Bayesian Analysis of Structural Breaks: China s Output and Productivity Series

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Transcription:

Advances in Managemen & Applied Economics, vol., no.1, 01, 109-119 ISSN: 179-7544 (prin version), 179-755 (online) Inernaional Scienific Press, 01 Time Series Models for Growh of Urban Populaion in SAARC Counries Md. Aikur Rahman Khan 1 and Md. Mosafizur Rahman Absrac The purpose of our presen sudy is o srike ou suiable models o explain he growh paern and o forecas for urban populaion in SAARC counries. Using he daa from UNPD for he years 1950 o 000 in five years inerval, we fied boh exponenial and models. We found he superioriy of models over exponenial models o explain he ime rend behavior of he urban populaion as a percenage of oal populaion and o forecas up o he year 05. We also found ha urbanizaion is faser in Bangladesh han any oher SAARC counries. JEL classificaion numbers: C10, R11 1 Deparmen of Populaion Science and Human Resource Developmen, Universiy of Rajshahi, Rajshahi-605, Bangladesh, e-mail: khanmar77@yahoo.com Deparmen of Saisics, Saisics and Mahemaics School, Yunnan Universiy of Finance and Economics, Kunming-6501, P.R. China, e-mail:mosafiz_bd1@yahoo.com Aricle Info: Received : December 5, 011. Revised : January 16, 01 Published online : February 8, 01

110 Urban Populaion in SAARC Counries Keywords: Urban Populaion, SAARC Counries, Roo Mean Squared Forecas Error (RMSFE), Cross Validiy Predicive Power (CVPP), Exponenial Growh and Auoregressive Moving Average () 1 Inroducion Time rend behavior of populaion parameer has been exensively sudied. To fi a ime rend model a firs graphical represenaion of he daa is examined. The line graph of populaion parameers generally posses he exponenial growh crierion (Misra, 1995; Bhende and Kanikar 1997). Thus, o explain he growh paern of populaion daa exponenial growh models were preferred by many auhors (UN, 1967 and 1997; UNFPA, 1993; BBS, 003). Bu, we can no assure ha he daa will be suiable for exponenial modeling raher here migh have some oher modeling echniques ha can be more suiable for. Time series daa are ime dependen and are auocorrelaed (Pankraz, 1991; Gujarai, 1995; Cleary and Hay, 1980). So, we can apply he ime series modeling echniques based on he effecs of auocorrelaion. In his paper, an aemp is made o examine he suiabiliy of boh exponenial and ime series models for he urban populaion in SAARC counries expressed as a percenage of oal populaion. Daa and Mehods The daa have been colleced from UNPD (1999) and Earhrend (004). The daa have been shown graphically using line plo (Figure1). We have analyzed he daa using some sophisicaed models including he ime rend behavior of he daa. As he ime frame was in five years inerval from 1950 o 000 so here migh have some impac over he passages of ime. In his recogniion we need o

Md. Aikur Rahman Khan and Md. Mosafizur Rahman 111 examine wheher he daa possess saionary crierion or no. We used he uni roo es of MacKinnon (1996) o es he saionariy of he daa. An auoregressive model of order p, AR(p) (Cleary and Hay, 1980; Pankraz, 1991, Gujarai, 1995; and Hamilon, 1994) is of he form Y α p i1 β Y u (1) i i where, Y is he value of he variable Y (urban populaion as a percenage of oal populaion) a ime, α is he inercep erm, β i (i 1,,..., p) are he parameers, Y is he i-h lagged variable, -i is he ime rend, γ is he coefficien of he ime variable, and u is he error erm which is a whie noise indeed. If we include he effec of q moving average erms, inercep, and rend componen along wih he AR(p) model hen an auoregressive moving average model of order p and q, ha is, (p,q/c,t) model (Gujarai, 1995 and Pankraz, 1991) can be formed as Y α p i1 β Y i i q j1 u j j u () where α is he inercep erm and is he ime rend. (p,q/c, T) refers o auoregressive moving average model wih p auoregressive erms and q moving average erms including inercep erm (C) and rend componen (T). For esing he significance of he fied model saionariy crierion of residuals and oulier deecion echniques (Pankraz, 1991) are applied. The saionariy of residuals are examined using he uni roo es (MacKinnon, 1996). We have used he mehod of examining sandardized residuals (absolue value of sandardized residuals over 3.0 implies he presence of an oulier) o examine he presence of oulier (Pankraz, 1991). We have also used he recen developed cross validiy predicive power (CVPP) and resriced cross validiy predicive power

11 Urban Populaion in SAARC Counries (RCVPP) (Khan and Ali, 003a). The cross validiy predicive power (CVPP) due o Sevens (1996) is (n 1)(n )(n 1) ρ cv 1- w(1- R ); w (3) n(n k 1)(n k ) and he resriced cross validiy predicive power is ρ 1- w(1 - R ); R = 0; oherwise. 1 w (n 1)(n )(n 1), w n(n k 1)(n k ), n k 1 (4) rcv where R is he coefficien of muliple deerminaion, n is he sample size, k is n 1 (n 1)(n -1)(n - ) he number of regressors used in he model, and w. n n Furher, he shrinkage of R in compuing RCVPP has been compued from he absolue difference beween RCVPP and 003a), ha is, η R (Sevens, 1996; and Khan and Ali, ρ rcv R 1 w(1 R ) R (5) Now, ρ rcv 0. 95 indicaes ha if we fi he same model o some oher daa from he same populaion hen he fied model will be able o explain 95% variaion of he dependen variable. Furher, η 0. 01 indicaes ha over he populaion he fied model is 99% sable. Thereafer, we have forecased from he fied models and have compued he roo mean squared forecas error o examine he forecasing performance of he fied model. Roo mean squared forecas error is compued (Pankraz, 1991) using he formula n r (O F ) n 1 RMSFE (6) r where O is he observed value, F is he corresponding forecased value, and r is he number of periods o forecas.

Md. Aikur Rahman Khan and Md. Mosafizur Rahman 113 Unforunaely, hese observed values, O s are no always in exisence. So, we have used he esimaed roo mean squared forecas error (ERMSFE) o examine he forecasing performance of he fied model (Khan and Ali, 003b) ha is compued as (nr)(nr k1)(nr k) SS(e) ERMSFE SS nr ()(1ρ rcv) r(n r 1)(n r )(n r 1) r (7) where, SS n r () is he sum of squared oal of n observed values and r forecased values, SS n () SS(oal) is he sum of squared oal of n observed values, SS(e) is he sum of squared residuals, and k is he number of predicors used in he model. 3 Resuls and Discussions In 1950 urban populaion of Bangladesh was 4.% of oal populaion. Bu, in 000 his percenage becomes 4.5, ha is, wihin 50 years he incremen of urban populaion is 0.3%. Wihin hese 50 years (1950-000) urban populaion in India, Pakisan, Sri Lanka, Nepal, Bhuan, and Maldives lifed up 11.1%, 19.5%, 9.%, 9.6%, 5.0%, and 15.5%, respecively. If we order he counry wih respec o he incremen of urban populaion wihin 1950 o 000 hen we can wrie BD>PK>ML>IN>NP>SL>BH. Thus, we can say ha he urbanizaion is faser in Bangladesh han any oher SAARC counries. To fi ime rend models we have examined he saionariy of he variables using he uni roo es and have found ha all he variables are saionary a level, ha is, inegraed of order zero. The esimaed resuls are given a Table 1.

114 Urban Populaion in SAARC Counries 40 35 30 5 0 15 10 5 0 1945 1955 1965 1975 1985 1995 bd in pk nepal maldiv sril bhuan Figure1: Urban populaion (in percen of oal populaion) in SAARC counries Table 1: Uni Roo Tes Variables Specificaion DF-Value MacKinnon Criical Saionary a Bangladesh None 5.764310 -.86* Level India None 9.37166 -.86* Level Pakisan None 1.3148 -.86* Level Nepal None 11.005 -.86* Level Bhuan None 14.3454 -.86* Level Maldives C, T and lag 1-4.4453-4.0815** Level Sri Lanka None.40738-1.9791** Level Here * (**) indicaes values a 1% (5%) level. We have fied models o explain he ime rend behavior of hose variables (Table ). Similarly, he fied exponenial growh models have shown in Table 3. From Table and Table 3 we can say ha he compued shrinkage of he fied exponenial growh model for urban populaion in Bangladesh, Nepal, Bhuan, and Maldives are less han heir corresponding fied models.

Md. Aikur Rahman Khan and Md. Mosafizur Rahman 115 Similarly, he compued shrinkage of he fied model for urban populaion as a percenage of oal populaion in India, Pakisan, and Srilanka are less han heir corresponding fied exponenial growh models. Table : Fied models for SAARC counries Counry & Fied Models R RCVPP Shrinkage model Bangladesh.111777T 0.70947X 0.97606 0.96614 0.009916-1 (1,0/C) Prob. (0.00000) (0.0178) Invered AR 0.71 India 1.168571.473118T 0.69199X 0.99715 0.99464 0.00516-1 (1,0/C,T) Pr ob. (0.00001) (0.00000) (0.0014) Invered AR 0.63 Pakisan.650389T 0.937693X 0.9937 0.9911 0.00599-1 (1,0/T) Prob. (0.00000) (0.0000) Invered AR 0.94 Sri Lanka 15.83038 0.771334T 0.989817û 0.9595 0.86037 0.065584-1 (0,1/C,T) Pr ob. (0.0000) (0.000) (0.0000) Invered MA -.99 Nepal 1.00508T 0.9018û 0.98516 0.9703 0.013137-1 (0,1/T) Pr ob. (0.00000) (0.0000) Invered MA -.9 Bhuan 0.65149T 0.659135X-1 0.80589 0.98768 0.96748 0.00199 (1,1/T) Pr ob. (0.00000) (0.0087) (0.0107 Invered AR 0.66, Invered MA 0. Maldives.49305T 0.77347X (1,0/T) Prob. (0.0001) ( 0.0061) Invered AR 0.73-1 0.93481 0.87707 0.057739 Here RCVPP is resriced cross validiy predicive power. I is o noe ha he residual analysis (normaliy and oulier deecion) were performed. No oulier was presen in he daa se as all he sandardized residuals

116 Urban Populaion in SAARC Counries were wihin 3. We found ha he residuals were normal wih respec o heir normal probabiliy plo. For simpliciy of conen we would like o relax hose graphical represenaions. Table 3: Fied Exponenial growh models for SAARC counries Counry & Fied Models R RCVPP Shrinkage model Bangladesh 1.61346 0.18178T e 0.97373 0.963554 0.009719 Prob. = (0.0000) (0.0000) India.74515 0.054399T e 0.987994 0.98368 0.004366 Prob. = (0.0000) (0.0000) Pakisan.85748 0.0683T e 0.99787 0.990164 0.0063 sig. = (0.0000) (0.0000) Sri Lanka.78547 0.035485T e 0.76173 0.66599 0.099573 Prob. = (0.0000) (0.0011) Nepal 0.63347 0.169864T e 0.99601 0.989910 0.00690 Prob. = (0.0000) (0.0000) Bhuan 0.498468 0.19555T e Prob. = (0.0000) (0.0000).18693 e Prob. = (0.0000) (0.0000) Maldives 0.105518T Here RCVPP is resriced cross validiy predicive power. 0.990309 0.986785 0.00354 0.879701 0.835956 0.043745 Forecased urban populaion in SAARC counries have shown in Table 4 and in Table 5. Table 5 shows he forecased urban populaion as a percenage of oal populaion from he fied exponenial models (Table 3) and Table 4 shows he forecased values from he fied models (Table ). Also, he ERMSFE have been incorporaed in hose ables. From Table 4 and Table 5 we see ha he ERMSFEs from he fied models are comparaively smaller han hose from he fied exponenial growh models. Fied exponenial model for he urban

Md. Aikur Rahman Khan and Md. Mosafizur Rahman 117 populaion of Nepal and Bhuan have greaer RCVPP as well as greaer ERMSFE han he corresponding fied models. Bu, RMSFE decreases wih he increase of resriced cross validiy predicive power (as explained by Khan and Ali, 003b). Table 4: Forecased urban populaion using he fied models of Table Year Bangladesh India Pakisan Nepal Maldives Sri Lanka Bhuan 005 6.467 9.86306 39.16154 1.51775 8.9538 5.09739 7.964004 (.5387) (0.75763) (1.4465) (0) (7.4813) (.58357) (0.4594) 010 8.0958 (.3171) 31.3985 (0.659) 41.35354 (1.493) 13.06537 (0.3368) 31.7097 (5.895) 5.85773 (.14494) 8.431600 (0.31404) 015 30.01857 (.39157) 3.79898 (0.636) 43.57411 (1.0350) 14.07040 (0.4105) 34.39317 (5.49477) 6.6906 (1.99391) 8.95901 (0.33781) 00 31.99854 (.47794) 34.6959 (0.61559) 45.8145 (1.19410) 15.07543 (0.4668) 37.050 (5.37585) 7.4004 (1.9708) 9.5096 (0.3563) 05 34.0168 (.56901) 35.74113 (0.6163) 48.09391 (1.1999) 16.08045 (0.49694) 39.6190 (5.3740) 8.17173 (1.89667) 10.08963 (0.3643381) Table 5: Forecased urban populaion using he fied models of Table 3 Year Bangladesh India Pakisan Nepal Maldives Sri Lanka Bhuan 005 31.518 9.9043 39.49 14.4505 3.61903 4.8101 7.7937 (3.433106) (1.01147) (1.391) (0.80357) (5.7090) (.96774) (0.4748) 010 37.4799 (3.74415) 31.576 (1.0514) 4.805 (1.588) 17.1561 (0.87193) 36.49 (5.987) 5.7063 (.91585) 8.870073 (0.45337) 015 44.94917 (4.1841) 33.3415 (1.07606) 45.656 (1.33041) 0.9631 (0.97633) 40.891 (6.36453) 6.6349 (3.00487) 10.097 (0.4960) 00 53.90698 (4.810053) 35.054 (1.13988) 48.4615 (1.41913) 4.05403 (1.10554) 44.76574 (6.9148) 7.597 (3.13703) 11.49363 (0.54809) 05 64.64997 (5.53798) 37.1736 (1.1163) 51.8831 (1.51933) 8.50747 (1.59875) 49.74744 (7.550968) 8.5938 (3.91503) 13.08345 (0.608558) Here, values in parenhesis are esimaed roo mean squared forecas error (ERMSFE).

118 Urban Populaion in SAARC Counries On he oher hand, he fied exponenial models for urban populaion of Bangladesh, Nepal, Bhuan, and Maldives posses less shrinkages bu more ERMSFEs compared o heir associaed fied models. Table 4 and Table 5 depic ha he ERMSFEs depend on boh he iniial populaion and he regression coefficien. If he iniial populaion is greaer bu he regression coefficien is smaller hen he ERMSFE becomes smaller and he ERMSFE becomes larger if he regression coefficien is larger. From Table 3 we find ha he regression coefficien of he fied exponenial model for Bangladesh is greaer han ha for Pakisan and India and so he ERMSFE for Bangladesh is larger. 5 Conclusions In fiing growh models for urban populaion in SAARC counries exponenial models for Bangladesh, Nepal, Bhuan, and Maldives lead less shrinkage bu more esimaed roo mean squared forecas error han he corresponding fied models. Bu, models for all he seven SAARC counries lead accepable shrinkage, smaller and more significanly acceped roo mean squared forecas error han he fied exponenial growh models. ACKNOWLEDGEMENTS. This research was suppored by Naional Naural Science Foundaion of China (YCT1017). References [1] Bangladesh Bureau of Saisics, Saisical Year Book, Various Issues. [] A.A. Bhende and T. Kanikar, Principles of Populaion Sudies, Himalay Publishing House, Mumbai, 1997.

Md. Aikur Rahman Khan and Md. Mosafizur Rahman 119 [3] Earhrend (004) hp://earhrends.wri.or [4] D.N. Gujarai, Basic Economerics, Mc Graw-Hill Inernaional, 3 rd Edi. 1995. [5] J.D. Hamilon, Time Series Analysis, Princeon Universiy Press, Princeon, New Jersey, 1994. [6] M.A.R. Khan and M.A. Ali, Resriced Cross Validiy Predicive Power, Pakisan Journal of Saisics, 19(), (003a). [7] M.A.R. Khan and M.A. Ali, VAR Modeling wih Mixed Series, Inernaional Journal of Saisical Sciences,, (003b). [8] R.G. King, C.I. Plosser, J.H. Sock and M.W. Wason, Sochasic Trends and Economic Flucuaions, Amer. Econ. Rev., 81, (Sepember, 1991), 819-840. [9] J.G. MacKinnon, Numerical Disribuion Funcion for Uni Roo and Coinegraion Tess, Journal of Applied Economerics, 11, (1996), 601-618. [10] R. Mc Cleary and R.A. Hay, Applied Time Series Analysis for he Social Science, Sega Publicaions, California, 1980. [11] B.D. Misra, An Inroducion o he Sudy of Populaion, Second Ediion, Souh Asian Publishers Pv. Ld, New Delhi, 1995. [1] A. Pankraz, Forecasing wih Dynamic Regression Models, John Wiley & Sons Inc., New York, 1991. [13] J. Sevens, Applied mulivariae saisics for he social sciences, Third Ediion, Lawrence Erlbaum Associaes, Inc., Publishers, Mahwah, New Jersey, 1996. [14] UN., The Deerminan and Consequences of Populaion Trends, Series A, Populaion Sudies, 1(4), New York, 1967. [15] UN., Feriliy Trends Among Low Feriliy Counries, Exper Group Meeing on Below-Replacemen Feriliy, Populaion Division, UN Secrearia, New York, 1997. [16] UNPD, Unied Naions Populaion Division, New York, 1999. [17] UNFP, Readings in Populaion Research Mehodology, 3 and 4, (1993).