Package TimeSeries.OBeu

Similar documents
Unit root problem, solution of difference equations Simple deterministic model, question of unit root

Package TSPred. April 5, 2017

Forecasting using R. Rob J Hyndman. 2.4 Non-seasonal ARIMA models. Forecasting using R 1

Univariate, Nonstationary Processes

Forecasting using R. Rob J Hyndman. 2.3 Stationarity and differencing. Forecasting using R 1

Package polypoly. R topics documented: May 27, 2017

Package atsa. July 8, 2015

Modeling and forecasting global mean temperature time series

Package thief. R topics documented: January 24, Version 0.3 Title Temporal Hierarchical Forecasting

Empirical Approach to Modelling and Forecasting Inflation in Ghana

Package Tcomp. June 6, 2018

STAT 436 / Lecture 16: Key

Package SK. May 16, 2018

Univariate ARIMA Models

0.1 ARIMA: ARIMA Models for Time Series Data

Forecasting Bangladesh's Inflation through Econometric Models

Package gtheory. October 30, 2016

Forecasting: principles and practice 1

Lab: Box-Jenkins Methodology - US Wholesale Price Indicator

Empirical Market Microstructure Analysis (EMMA)

FORECASTING THE INVENTORY LEVEL OF MAGNETIC CARDS IN TOLLING SYSTEM

TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA

Firstly, the dataset is cleaned and the years and months are separated to provide better distinction (sample below).

Package invgamma. May 7, 2017

Autoregressive Integrated Moving Average Model to Predict Graduate Unemployment in Indonesia

Automatic seasonal auto regressive moving average models and unit root test detection

Package NlinTS. September 10, 2018

Package MultisiteMediation

Forecasting using R. Rob J Hyndman. 3.2 Dynamic regression. Forecasting using R 1

Package EMVS. April 24, 2018

Seasonality. Matthieu Stigler January 8, Version 1.1

Stat 153 Time Series. Problem Set 4

Estimating trends using filters

Package RTransferEntropy

Report for Fitting a Time Series Model

Package nardl. R topics documented: May 7, Type Package

FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL

Package bhrcr. November 12, 2018

Package reinforcedpred

Package WaveletArima

BCT Lecture 3. Lukas Vacha.

Package jmcm. November 25, 2017

Testing for non-stationarity

Case Study: Modelling Industrial Dryer Temperature Arun K. Tangirala 11/19/2016

Package bigtime. November 9, 2017

STAT Financial Time Series

Economics 618B: Time Series Analysis Department of Economics State University of New York at Binghamton

Package varband. November 7, 2016

Stochastic Analysis and Forecasts of the Patterns of Speed, Acceleration, and Levels of Material Stock Accumulation in Society

Lecture 5: Estimation of time series

ARDL Cointegration Tests for Beginner

Author: Yesuf M. Awel 1c. Affiliation: 1 PhD, Economist-Consultant; P.O Box , Addis Ababa, Ethiopia. c.

Package TSF. July 15, 2017

Lecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem

Time Series Forecasting: A Tool for Out - Sample Model Selection and Evaluation

Advanced Econometrics

Package ARCensReg. September 11, 2016

On Consistency of Tests for Stationarity in Autoregressive and Moving Average Models of Different Orders

Response surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test

ECONOMETRIA II. CURSO 2009/2010 LAB # 3

Suan Sunandha Rajabhat University

11/18/2008. So run regression in first differences to examine association. 18 November November November 2008

Application of ARIMA Models in Forecasting Monthly Total Rainfall of Rangamati, Bangladesh

Time Series Analysis of United States of America Crude Oil and Petroleum Products Importations from Saudi Arabia

Exercises - Time series analysis

477/577 In-class Exercise 5 : Fitting Wine Sales

Estimation and application of best ARIMA model for forecasting the uranium price.

Scenario 5: Internet Usage Solution. θ j

Package OUwie. August 29, 2013

Package AID. R topics documented: November 10, Type Package Title Box-Cox Power Transformation Version 2.3 Date Depends R (>= 2.15.

MODELING MAXIMUM MONTHLY TEMPERATURE IN KATUNAYAKE REGION, SRI LANKA: A SARIMA APPROACH

Package RSwissMaps. October 3, 2017

STAT 520 FORECASTING AND TIME SERIES 2013 FALL Homework 05

Dynamic Time Series Regression: A Panacea for Spurious Correlations

Univariate linear models

Package flexcwm. R topics documented: July 2, Type Package. Title Flexible Cluster-Weighted Modeling. Version 1.1.

Analysis. Components of a Time Series

Package FinCovRegularization

Forecasting with R A practical workshop

Package cvequality. March 31, 2018

7. Integrated Processes

Box-Jenkins ARIMA Advanced Time Series

Response surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test

FE570 Financial Markets and Trading. Stevens Institute of Technology

Package tspi. August 7, 2017

Package exuber. June 17, 2018

Package ForecastCombinations

Classical Decomposition Model Revisited: I

7. Integrated Processes

Chapter 3: Regression Methods for Trends

Transformations for variance stabilization

Package nowcasting. April 25, Type Package

Leftovers. Morris. University Farm. University Farm. Morris. yield

Package REGARMA. R topics documented: October 1, Type Package. Title Regularized estimation of Lasso, adaptive lasso, elastic net and REGARMA

Package skellam. R topics documented: December 15, Version Date

Package SimCorMultRes

Ch 6. Model Specification. Time Series Analysis

Package sym.arma. September 30, 2018

Homework 5, Problem 1 Andrii Baryshpolets 6 April 2017

Package testassay. November 29, 2016

Transcription:

Type Package Package TimeSeries.OBeu Title Time Series Analysis 'OpenBudgets' Version 1.2.2 Date 2018-01-20 January 22, 2018 Estimate and return the needed parameters for visualisations designed for 'OpenBudgets' <http://openbudgets.eu/> time series data. Calculate time series model and forecast parameters in Budget time series data of municipalities across Europe, according to the 'OpenBudgets' data model. There are functions for measuring deterministic and stochastic trend of the input time series data, decomposing with local regression models or seasonal trend decomposition, modelling the appropriate auto regressive integrated moving average model and provide forecasts for the input 'OpenBudgets' time series fiscal data. Also, can be used generally to extract visualisation parameters convert them to 'JSON' format and use them as input in a different graphical interface. Maintainer Kleanthis Koupidis <koupidis@okfn.gr> URL https://github.com/okgreece/timeseries.obeu BugReports https://github.com/okgreece/timeseries.obeu/issues License GPL-2 file LICENSE LazyData TRUE Imports forecast, locfit, jsonlite, stats, trend, tseries RoxygenNote 6.0.1 Suggests knitr, rmarkdown VignetteBuilder knitr NeedsCompilation no Author Kleanthis Koupidis [aut, cre], Charalampos Bratsas [aut] Repository CRAN Date/Publication 2018-01-22 09:41:59 UTC 1

2 Athens_approved_ts R topics documented: Athens_approved_ts..................................... 2 Athens_draft_ts....................................... 3 Athens_executed_ts..................................... 3 Athens_reserved_ts..................................... 4 Athens_revised_ts...................................... 4 open_spending.ts...................................... 5 ts.acf............................................. 6 ts.analysis.......................................... 7 ts.forecast.......................................... 10 ts.non.seas.decomp..................................... 11 ts.non.seas.model...................................... 12 ts.seasonal.decomp..................................... 14 ts.seasonal.model...................................... 15 ts.stationary.test....................................... 17 Index 19 Athens_approved_ts Time series of Approved Expenditure Budget Phase of Municipality of Athens Time series data with the Approved Budget phase expenditure amounts of Municipality of Athens from 2004-2015 The years of the recorded approved budget phase amounts The approved budget phase amounts of this time range Format A ts object with 12 approved amounts from 2004-2015 Source The ttl and rdf expenditure data are stored in: https://github.com/openbudgets/datasets/tree/master/greekmunicipalities/municipality-of-athens/dataset

Athens_draft_ts 3 Athens_draft_ts Time series of Draft Expenditure Budget Phase of Municipality of Athens Time series data with the Draft Budget phase expenditure amounts of Municipality of Athens from 2004-2015 The years of the recorded draft budget phase amounts. The draft budget phase amounts of this time range. Format A ts object with 12 draft amounts from 2004-2015 Source The ttl and rdf expenditure data are stored in: https://github.com/openbudgets/datasets/tree/master/greekmunicipalities/municipality-of-athens/dataset Athens_executed_ts Time series of Executed Expenditure Budget Phase of Municipality of Athens Time series data with the Executed Budget phase expenditure amounts of Municipality of Athens from 2004-2015 The years of the recorded executed budget phase amounts. The executed budget phase amounts of this time range. Format A ts object with 12 draft amounts from 2004-2015 Source The ttl and rdf expenditure data are stored in: https://github.com/openbudgets/datasets/tree/master/greekmunicipalities/municipality-of-athens/dataset

4 Athens_revised_ts Athens_reserved_ts Time series of Reserved Expenditure Budget Phase of Municipality of Athens Time series data with the Reserved Budget phase expenditure amounts of Municipality of Athens from 2004-2015 The years of the recorded reserved budget phase amounts. The reserved budget phase amounts of this time range. Format A ts object with 12 reserved amounts from 2004-2015 Source The ttl and rdf expenditure data are stored in: https://github.com/openbudgets/datasets/tree/master/greekmunicipalities/municipality-of-athens/dataset Athens_revised_ts Time series of Revised Expenditure Budget Phase of Municipality of Athens Time series data with the Revised Budget phase expenditure amounts of Municipality of Athens from 2004-2015 The years of the recorded revised budget phase amounts. The revised budget phase amounts of this time range. Format A ts object with 12 revised amounts from 2004-2015 Source The ttl and rdf expenditure data are stored in: https://github.com/openbudgets/datasets/tree/master/greekmunicipalities/municipality-of-athens/dataset

open_spending.ts 5 open_spending.ts Read and analyze univariate time series data from Open Spending API Extract and analyze univariate time series data from Open Spending API, using the ts.analysis function. Usage open_spending.ts(json_data,time,amount,order=null,prediction_steps=1) Arguments json_data time amount order The json string, URL or file from Open Spending API Specify the time label of the json time series data Specify the amount label of the json time series data An integer vector of length 3 specifying the order of the Arima model prediction_steps The number of prediction steps Details This function extracts the time series data provided by the Open Spending API, in order to return the results from the ts.analysis function. Value A json string with the resulted parameters of the ts.analysis function. Author(s) Kleanthis Koupidis See Also ts.analysis

6 ts.acf ts.acf Extract the ACF and PACF parameters of time series and their model residuals Usage This function is included in ts.analysis function and aims to extract the ACF and PACF details of the input time series data and the ACF, PACF of the residuals after fitting an Arima model. ts.acf(tsdata, model_residuals, a=0.95, tojson=f) Arguments tsdata The input univariate time series data model_residuals The model s residuals after fitting a model to the time series a Details Value tojson The significant level (default a=0.95) If TRUE the results are returned in json format, default returns a list This function is used internally in ts.analysis function and the output is a list with grouped ACF and PACF parameters of the input time series data, as well as the ACF and PACF parameters of the residuals needed for the graphical purposes in OBEU. A list with the parameters: acf.parameters: acf The estimated acf values of the input time series acf.lag The lags at which the acf is estimated confidence.interval.up The upper limit of the confidence interval confidence.interval.low The lower limit of the confidence interval pacf.parameters: pacf The estimated pacf values of the input time series pacf.lag The lags at which the pacf is estimated confidence.interval.up The upper limit of the confidence interval confidence.interval.low The lower limit of the confidence interval acf.residuals.parameters: acf.res The estimated acf values of the model residuals acf.res.lag The lags at which the acf is estimated of the model residuals

ts.analysis 7 confidence.interval.up The upper limit of the confidence interval confidence.interval.low The lower limit of the confidence interval pacf.residuals.parameters: pacf.res The estimated pacf values of the model residuals pacf.res.lag The lags at which the pacf is estimated of the model residuals confidence.interval.up The upper limit of the confidence interval confidence.interval.low The lower limit of the confidence interval Author(s) Kleanthis Koupidis See Also ts.analysis, Acf, Pacf Examples ts.acf(athens_draft_ts) ts.analysis Time series analysis results for OBEU Time series Univariate time series analysis for short and long time series data using the appropriate model. Usage ts.analysis(tsdata, x.order=null, prediction.steps=1, tojson=t) Arguments tsdata x.order The input univariate time series data An integer vector of length 3 specifying the order of the Arima model prediction.steps The number of prediction steps tojson If TRUE the results are returned in json format, default returns a list

8 ts.analysis Details Value This function automatically tests for stationarity of the input time series data using ts.stationary.test function. Depending the nature of the time series data and the stationary tests there are four branches: a.)short and non seasonal, b.)short and seasonal, c.)long and non seasonal and d.)long and seasonal. For branches a and c ts.non.seas.model is used and for b and d ts.seasonal.model is used. This function also decomposes both seasonal and non seasonal time series through ts.non.seas.decomp and ts.seasonal.decomp and forecasts h steps ahead the user selected(default h=1) using ts.forecast. A json string with the parameters: acf.param acf.parameters: * acf The estimated acf values of the input time series * acf.lag The lags at which the acf is estimated * confidence.interval.up The upper limit of the confidence interval * confidence.interval.low The lower limit of the confidence interval pacf.parameters: * pacf The estimated pacf values of the input time series * pacf.lag The lags at which the pacf is estimated * confidence.interval.up The upper limit of the confidence interval * confidence.interval.low The lower limit of the confidence interval acf.residuals.parameters: * acf.res The estimated acf values of the model residuals * acf.res.lag The lags at which the acf is estimated of the model residuals * confidence.interval.up The upper limit of the confidence interval * confidence.interval.low The lower limit of the confidence interval pacf.residuals.parameters: * pacf.res The estimated pacf values of the model residuals * pacf.res.lag The lags at which the pacf is estimated of the model residuals * confidence.interval.up The upper limit of the confidence interval * confidence.interval.low The lower limit of the confidence interval param stl.plot: * trend The estimated trend component * trend.ci.up The estimated up limit for trend component (for non seasonal time series) * trend.ci.low The estimated low limit for trend component (for non seasonal time series) * seasonal The estimated seasonal component * remainder The estimated remainder component * time The time of the series was sampled

ts.analysis 9 Author(s) See Also stl.general: * stl.degree The degree of fit * degfr The effective degrees of freedom for non seasonal time series * degfr.fitted The fitted degrees of freedom for non seasonal time series * fitted The model s fitted values residuals The residuals of the model (fitted innovations) compare: * arima.order The Arima order for seasonal time series * arima.coef A vector of AR, MA and regression coefficients * arima.coef.se The standard error of the coefficients * covariance.coef The matrix of the estimated variance of the coefficients * resid.variance The residuals variance * not.used.obs The number of not used observations for the fitting * used.obs The used observations for the fitting * loglik The maximized log-likelihood (of the differenced data), or the approximation to it used * aic The AIC value corresponding to the log-likelihood * bic The BIC value corresponding to the log-likelihood * gcv The generalized cross-validation statistic time series * aicc The second-order Akaike Information Criterion corresponding to the log-likelihood for seasonal time series forecasts ts.model a string indicating the arima orders data_year The time that time series data were sampled data The time series values predict_time The time that defined by the prediction_steps parameter predict_values The predicted values that defined by the prediction_steps parameter up80 The upper limit of the 80% predicted confidence interval low80 The lower limit of the 80% predicted confidence interval up95 The upper limit of the 95% predicted confidence interval low95 The lower limit of the 95% predicted confidence interval Kleanthis Koupidis ts.stationary.test, ts.acf, ts.seasonal.model, ts.seasonal.decomp, ts.non.seas.model, ts.non.seas.decomp, ts.forecast Examples ts.analysis(athens_draft_ts, prediction.steps=3)

10 ts.forecast ts.forecast Time series forecast results of OBEU Time Series Univariate time series forecasts for short and long time series data using the appropriate model. Usage ts.forecast(ts_modelx, h=1, tojson=f) Arguments ts_modelx h tojson The input univariate time series data The number of prediction steps If TRUE the results are returned in json format, default returns a list Details This function is used internally in ts.analysis and forecasts the model that fits the input data using the auto.arima function(see forecast package). The model selection depends on the results of some diagnostic tests (acf,pacf,pp adf and kpss). For short time series the selected arima model is among various orders of the AR part using the first differences and the first order moving average component, with the lower AIC value. Value A list with the parameters: ts.model a string indicating the arima orders data_year The time that time series data were sampled data The time series values predict_time The time that defined by the prediction_steps parameter predict_values The predicted values that defined by the prediction_steps parameter up80: The upper limit of the 80% predicted confidence interval low80: The lower limit of the 80% predicted confidence interval up95: The upper limit of the 95% predicted confidence interval low95: The lower limit of the 95% predicted confidence interval Author(s) Kleanthis Koupidis

ts.non.seas.decomp 11 See Also ts.analysis, forecast Examples Athens_draft <- ts.non.seas.model(athens_draft_ts) #Hold the model object of non seasonal modeling draft<-athens_draft$model.summary ts.forecast(draft) ts.non.seas.decomp Non seasonal decomposition Decomposition of time series with no seasonal component using local regression models Usage ts.non.seas.decomp(tsdata, tojson=f) Arguments tsdata tojson The input univariate non seasonal time series data If TRUE the results are returned in json format, default returns a list Details Value For non-seasonal time series there is no seasonal component. Local regression and likelihood models (locfit package) are used in order to extract the trend and remaider components. A list with the following components: stl.plot: trend The estimated trend component trend.ci.up The estimated up limit for trend component trend.ci.low The estimated low limit for trend component seasonal The estimated seasonal component remainder The estimated remainder component time The time of the series was sampled stl.general: stl.degree The degree of fit degfr The effective degrees of freedom

12 ts.non.seas.model degfr.fitted The fitted degrees of freedom residuals_fitted: residuals The residuals of the model (fitted innovations) fitted The model s fitted values time the time of tsdata line The y=0 line compare: resid.variance The residuals variance used.obs The used observations for the fitting loglik The maximized log-likelihood (of the differenced data), or the approximation to it used aic The AIC value corresponding to the log-likelihood bic The BIC value corresponding to the log-likelihood gcv The generalized cross-validation statistic Author(s) Kleanthis Koupidis References add See Also ts.analysis, locfit, predict.locfit Examples ts.non.seas.decomp(athens_draft_ts) ts.non.seas.model Model fit of non seasonal time series Model fit of non seasonal time series Usage ts.non.seas.model(tsdata, x.ord=null, tojson=f)

ts.non.seas.model 13 Arguments tsdata x.ord tojson The input univariate non seasonal time series data An integer vector of length 3 specifying the order of the Arima model If TRUE the results are returned in json format, default returns a list Details Value Model fit of non seasonal time series using arima models of non seasonal time series data. The model with the lowest AIC value is selected for forecasts. A list with the following components: Author(s) See Also model.summary: ts_model The summary model details returned as Arima object for internal use in ts.analysis function model: ts_model: arima.order The Arima order arima.coef A vector of AR, MA and regression coefficients arima.coef.se The standard error of the coefficients residuals: The residuals of the model (fitted innovations) compare: variance.coef The matrix of the estimated variance of the coefficients resid.variance The MLE of the innovations variance not.used.obs The number of not used observations for the fitting used.obs the number of used observations for the fitting loglik The maximized log-likelihood (of the differenced data), or the approximation to it used aic The AIC value corresponding to the log-likelihood bic The BIC value corresponding to the log-likelihood aicc The second-order Akaike Information Criterion corresponding to the log-likelihood Kleanthis Koupidis ts.analysis, Arima Examples ts.non.seas.model(athens_draft_ts)

14 ts.seasonal.decomp ts.seasonal.decomp Decomposition of seasonal time series Usage Decomposition of seasonal time series data using stlm from forecast package. This function is used internally in ts.analysis. ts.seasonal.decomp(tsdata, tojson=f) Arguments tsdata tojson The input univariate seasonal time series data If TRUE the results are returned in json format, default returns a list Details Value Decomposition of seasonal time series data through arima models is based on stlm from forecast package and returns a list with useful parameters for OBEU. A list with the following components: stl.plot: trend: The estimated trend component seasonal: The estimated seasonal component remainder: The estimated remainder component time: The time of the series was sampled stl.general: model.summary The summary object of the arima model to use in forecast if needed stl.win: An integer vector of length 3 indicating the spans used for the "s", "t", and "l" smoothers stl.degree: An integer vector of length 3 indicating the polynomial degrees for these smoothers residuals_fitted: residuals: The residuals of the model (fitted innovations) fitted: The model s fitted values time the time of tsdata line The y=0 line compare: arima.order: The Arima order

ts.seasonal.model 15 arima.coef: A vector of AR, MA and regression coefficients arima.coef.se: The standard error of the coefficients covariance.coef: The matrix of the estimated variance of the coefficients resid.variance: The MLE of the innovations variance not.used.obs: The number of not used observations for the fitting used.obs: the number of used observations for the fitting loglik: The maximized log-likelihood (of the differenced data), or the approximation to it used aic: The AIC value corresponding to the log-likelihood bic: The BIC value corresponding to the log-likelihood aicc: The second-order Akaike Information Criterion corresponding to the log-likelihood Author(s) Kleanthis Koupidis References add See Also ts.analysis, stlm ts.seasonal.model Model fit of seasonal time series Model fit of seasonal time series Usage ts.seasonal.model(tsdata, x.ord=null, tojson=f) Arguments tsdata x.ord tojson The input univariate seasonal time series data An integer vector of length 3 specifying the order of the Arima model If TRUE the results are returned in json format, default returns a list Details Model fit of seasonal time series using arima models of seasonal time series data. The model with the lowest AIC value is selected for forecasts.

16 ts.seasonal.model Value A list with the following components: model.summary: ts_model The summary model details returned as Arima object for internal use in ts.analysis function model: ts_model arima.order The Arima order arima.coef A vector of AR, MA and regression coefficients arima.coef.se The standard error of the coefficients residuals_fitted: residuals The residuals of the model (fitted innovations) fitted The model s fitted values time the time of tsdata line The y=0 line compare: variance.coef The matrix of the estimated variance of the coefficients resid.variance The MLE of the innovations variance not.used.obs The number of not used observations for the fitting used.obs the number of used observations for the fitting loglik The maximized log-likelihood (of the differenced data), or the approximation to it used aic The AIC value corresponding to the log-likelihood bic The BIC value corresponding to the log-likelihood aicc The second-order Akaike Information Criterion corresponding to the log-likelihood Author(s) Kleanthis Koupidis See Also ts.analysis, Arima

ts.stationary.test 17 ts.stationary.test Stationarity testing This functions tests the stationarity of the input time series data. Usage ts.stationary.test(tsdata) Arguments tsdata The input univariate time series data Details This function tests the deterministic and stochastic trend of the input time series data. This function uses ACF and PACF functions from forecast package, Phillips Perron test, Augmented Dickey Fuller (ADF) test, Kwiatkowski Phillips Schmidt Shin (KPSS) test, from tseries package and Mann Kendall test for Monotonic Trend Cox Stuart trend test from trend package. Phillips Perron test tests the null hypothesis of whether a unit root is present in a time series sample, against a stationary alternative. The truncation lag parameter is set to trunc(4*(n/100)^0.25), where n the length of the in input time series data Augmented Dickey Fuller (ADF) test, tests the null hypothesis of whether a unit root is present in a time series sample. The truncation lag parameter is set to trunc((n-1)^(1/3))), where n the length of the input time series data Kwiatkowski Phillips Schmidt Shin (KPSS) test, tests a null hypothesis that an observable time series is stationary around a deterministic trend (i.e. trend stationary) against the alternative of a unit root. The truncation lag parameter is set to trunc(3*sqrt(n)/13), where n the length of the input time series data The non parametric Mann Kendall test is used to detect monotonic trends. The null hypothesis, H0, is that the data come from a population with independent realizations and are identically distributed. The alternative hypothesis, HA, is that the data follow a monotonic trend. The Cox Stuart test is a modified sign test. The null hypothesis, H0, is that the input time series assumed to be independent against the fact that there is a time dependent trend (monotonic trend). Value A string indicating if the time series is stationary or non stationary for internal use in ts.analysis. Author(s) Kleanthis Koupidis

18 ts.stationary.test References tseries, trend See Also ts.analysis, Acf, Pacf, pp.test, adf.test, kpss.test, mk.test, cs.test Examples ts.stationary.test(athens_approved_ts)

Index Acf, 7, 18 adf.test, 18 Arima, 13, 16 Athens_approved_ts, 2 Athens_draft_ts, 3 Athens_executed_ts, 3 Athens_reserved_ts, 4 Athens_revised_ts, 4 cs.test, 18 forecast, 11 kpss.test, 18 locfit, 12 mk.test, 18 open_spending.ts, 5 Pacf, 7, 18 pp.test, 18 predict.locfit, 12 stlm, 15 ts.acf, 6, 9 ts.analysis, 5, 7, 7, 11 13, 15, 16, 18 ts.forecast, 8, 9, 10 ts.non.seas.decomp, 8, 9, 11 ts.non.seas.model, 8, 9, 12 ts.seasonal.decomp, 8, 9, 14 ts.seasonal.model, 8, 9, 15 ts.stationary.test, 8, 9, 17 19