Forecasting Framework for Open Access Time Series in Energy

06/02/2016 ∙ by Gergo Barta, et al. ∙ Budapest University of Technology and Economics Dmlab 0

In this paper we propose a framework for automated forecasting of energy-related time series using open access data from European Network of Transmission System Operators for Electricity (ENTSO-E). The framework provides forecasts for various European countries using publicly available historical data only. Our solution was benchmarked using the actual load data and the country provided estimates (where available). We conclude that the proposed system can produce timely forecasts with comparable prediction accuracy in a number of cases. We also investigate the probabilistic case of forecasting - that is, providing a probability distribution rather than a simple point forecast - and incorporate it into a web based API that provides quick and easy access to reliable forecasts.



There are no comments yet.


page 6

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

I Introduction

There is a great deal of interest in forecasting energy-related time series, specifically regarding the generation of renewable energy sources, load forecasting, electricity price forecasting or gas consumption forecasting [3, 1, 2]. For transparency purposes two European platforms offer open access to such data: ENTSO-E and ENTSOG. The data from these sources are rich, some data quality issues arise, however, that need some investigation for the platform to be useful. In this paper we utilize the free data sources to show that releasing data to the public is beneficial and should be promoted. By building a framework for the automatic prediction of energy-related time series published on an open portal we wish to add our part of encouragement.

Furthermore we incorporate the idea promoted in the Global Energy Forecasting Competition [4]; a new and unique approach to forecasting hourly renewable energy production, electric load and price was introduced: the forecasting of a probability distribution of the target variable in contrast to forecasting a single value. Probabilistic forecasts provide more detailed insight for stakeholders offering immediate business value. Our solution achieved good results in terms of pinball loss in the competition, and we assess if the framework is suitable for reliable probabilistic forecasts in a rolling fashion for a number of energy-related time series.

The document is organized as follows. After some basic information on the used data source in Section II, we present in Section III our comprehensive report on data quality of the platform. A brief methodological overview can be found in Section IV and the extensive tour of the internal framework structure in Section V. We perform an in-depth analysis of real-world data, then present and discuss detailed experimental results (Section VI). Conclusions and further ideas are provided in the final Section VII.

Ii Data source

ENTSO-E, the European Network of Transmission System Operators, represents electricity transmission system operators (TSOs) from countries across Europe. ENTSO-E was established and given legal mandates by the EU’s Third Legislative Package for the Internal Energy Market in 2009, which aims at further liberalizing the gas and electricity markets in the EU. The main objectives of ENTSO-E centre on the integration of renewable energy sources such as wind and solar power into the power system, and the completion of the internal energy market, which is central to meeting the European Union’s energy policy objectives of affordability, sustainability and security of supply[5].

The ENTSO-E Transparency Platform was launched in early 2015. It is aimed at providing free, continuous access to pan-European electricity market data for all users, across categories such as load, generation, transmission, balancing, outages and congestion. Our work focuses on the country load data published on this new portal.

Iii Data quality report

We performed the data validation on the available load data on the ENTSO-E Transparency Platform. The time span of the validation started on January 2015 and went on until October 2015, spanning 3 seasons of the year. The portal is used to access historical load data in a timely manner. Every day we perform a screen scrape of the data source for the last day of available data. Observation of the data showed that there are major types of data reporting frequencies in the publicly available load data: hourly, half-hourly and quarter hourly data reporting.

Data consistency was measured by

simple heuristics; check if the country provided data was either ’N/A’ or ’0’ at the time of the download. We summarize the occurrence of the above events for both of the provided data columns: Day-ahead Total Load Forecast, Actual Total Load. To get a bigger picture and a general score of data quality for each country, we observe the frequency of errors relative to the theoretical maximum reportable values for the validation time interval (Table


We selected a number of countries for investigation (boldfaced) based on their data quality and market size; Belgium, Czech Republic, Finland, France, Germany, Hungary, Italy, Norway, Poland and Slovakia. There are several countries that do not have any data on the platform at all in the validation time span, we omit these from further analysis. These countries are the following: Bosnia and Herzegovina, Cyprus, Belarus, Iceland, Moldavia, Malta, Russia, Turkey and Ukraine.

Country code Country name Data quality Frequency
10YAL-KESH—–5 Albania 58.5% Quarter-hourly
10YAT-APG——L Austria 98.5% Quarter-hourly
10YBE———-2 Belgium 91.2% Quarter-hourly
10YCA-BULGARIA-R Bulgaria 46.7% Hourly
10YHR-HEP——M Croatia 97.4% Hourly
10YCZ-CEPS—–N Czech Republic 98.2% Hourly
10Y1001A1001A65H Denmark 91.6% Hourly
10Y1001A1001A39I Estonia 97.2% Hourly
10YFI-1——–U Finland 98.8% Hourly
10YFR-RTE——C France 97.7% Hourly
10Y1001A1001A83F Germany 83.9% Quarter-hourly
10YGR-HTSO—–Y Greece 97.6% Hourly
10YHU-MAVIR—-U Hungary 98.8% Quarter-hourly
10YIE-1001A00010 Ireland 49.8% Half-hourly
10YIT-GRTN—–B Italy 77.1% Hourly
10YLV-1001A00074 Latvia 98.5% Hourly
10YLT-1001A0008Q Lithuania 96.8% Hourly
10YLU-CEGEDEL-NQ Luxembourgh 95.7% Quarter-hourly
10YMK-MEPSO—-8 Macedonia 63.3% Hourly
10YCS-CG-TSO—S Montenegro 91.3% Hourly
10YNL———-L Netherlands 98.1% Quarter-hourly
10YNO-0——–C Norway 88.7% Hourly
10YPL-AREA—–S Poland 98.0% Hourly
10YPT-REN——W Portugal 98.2% Hourly
10YRO-TEL——P Romania 92.8% Hourly
10YCS-SERBIATSOV Serbia 85.8% Hourly
10YSK-SEPS—–K Slovakia 98.4% Hourly
10YSI-ELES—–O Slovenia 96.4% Hourly
10YES-REE——0 Spain 98.6% Hourly
10YSE-1——–K Sweden 79.4% Hourly
10YCH-SWISSGRIDZ Switzerland 82.4% Hourly
GB United Kingdom 48.1% Half-hourly
TABLE I: Overall data quality assessment score

In general we can say that in the data quality assessment period the majority 21 countries out of 32 have less than bad data frequency. After cross-referencing the missing values from our historical real-time database, with the state of the database of the Transparency Platform, we see that the majority of these values are backfilled and the errors are corrected. However if one tries to rely on the timely updates of the Transparency Platform one could get into trouble with missing load values and forecasts as well. Surprisingly the better reporting countries are from the Central-Eastern Europe, most notably Croatia, Poland, Slovakia and Hungary.

All of the above countries have missing values in both the forecasted load and the actual reported load, however only a handful have zeros (Croatia, Portugal, Lithuania, Montenegro, Sweden, Serbia, Macedonia and Ireland). It should be self-evident not to report zeros as forecasts or load values, we strongly advise a check to be made at the time of the report. Table II summarizes the two metrics as percentages for each country.

Country Target is NA Target is 0 Forecast is NA Forecast is 0
Albania 49.47% 0.00% 33.59% 0.00%
Austria 1.97% 0.00% 0.99% 0.00%
Belgium 9.24% 0.00% 8.28% 0.00%
Bulgaria 13.68% 0.00% 13.00% 0.00%
Croatia 0.42% 0.28% 0.45% 0.16%
Czech Republic 0.54% 0.00% 0.36% 0.00%
Denmark 2.36% 0.00% 1.82% 0.00%
Estonia 0.91% 0.00% 0.50% 0.00%
Finland 0.51% 0.00% 0.10% 0.00%
France 0.90% 0.00% 0.25% 0.00%
Germany 10.37% 0.00% 21.86% 0.00%
Greece 0.95% 0.00% 0.28% 0.00%
Hungary 1.39% 0.00% 1.08% 0.00%
Ireland 12.81% 0.25% 37.17% 0.00%
Italy 10.20% 0.00% 1.23% 0.00%
Latvia 0.43% 0.00% 0.34% 0.00%
Lithuania 0.79% 0.03% 0.75% 0.00%
Luxembourg 5.88% 0.00% 2.63% 0.01%
Macedonia 11.03% 0.08% 6.78% 0.47%
Montenegro 1.74% 1.02% 0.57% 1.01%
Netherlands 1.41% 0.00% 2.34% 0.00%
Norway 3.32% 0.00% 2.36% 0.00%
Poland 0.43% 0.00% 0.58% 0.00%
Portugal 0.56% 0.08% 0.25% 0.00%
Romania 2.79% 0.00% 0.82% 0.00%
Serbia 3.36% 0.01% 3.71% 0.00%
Slovakia 0.47% 0.00% 0.33% 0.00%
Slovenia 0.93% 0.03% 0.86% 0.00%
Spain 0.53% 0.00% 0.17% 0.00%
Sweden 4.86% 0.40% 5.02% 0.00%
Switzerland 8.34% 0.00% 0.44% 0.00%
United Kingdom 13.53% 0.00% 38.32% 0.00%
TABLE II: Observed data issues in countries with sufficient data

In general it is beneficial that the errors are taken care of in time, but this issue makes it harder to act upon the reported information on the portal. Also it would be beneficial to have, at least for the research community, a simple API for the data that is otherwise only accessible through screen scraping.

Iv Methodology

The aim of this paper is not only to assess the utility of the new ENTSO-E platform, but to motivate publishing open access time series data in energy by building accurate models using the state-of-the-art in machine learning.

Point forecasts for time series are often provided by reliable statistical methods like ARMA and ARIMA. Our approach to forecasting is based on a relatively new machine learning method; Gradient Boosted Regression Trees (GBRTs)

[6] were successful in a number of competitions in the past [7, 2].

Iv-a Gradient Boosting Regression Trees

Gradient boosting is an ensemble method responsible for combining weak learners for higher model accuracy, as suggested by Friedman in 2000 [8]. The predictor generated in gradient boosting is a linear combination of weak learners, here tree models are utilized for this purpose.

A sequence of models is built iteratively, starting from a constant initial prediction (in a least squares regression task, the mean of the target variable), and with each tree built on the residuals of the previous tree. To avoid overfitting, there is a shrinkage parameter (also called learning rate), which controls how much each tree contributes to the final model; this effectively slows learning down, yielding better performance at the expense of requiring more iterations.

The final prediction is calculated as follows:


Where is the number of individual trees built, is the initial prediction, is the prediction at stage , is the prediction of an individual tree learner and is the shrinkage parameter to reduce the effect of an individual tree in the sequence.

The general idea is to compute a sequence of simple trees, where each successive tree is built for the prediction residuals of the preceding tree. Each new base-learner is chosen to be maximally correlated with the negative gradient of the loss function, associated with the whole ensemble. This way the subsequent stages will focus harder on fitting these examples and the resulting predictor is a linear combination of weak learners.

Utilizing boosting has many beneficial properties; various risk functions are applicable, intrinsic variable selection is carried out, also resolves multicollinearity issues, and works well with large number of features without overfitting.

Iv-B Evaluation

Model performance is measured using the Mean Absolute Percentage Error (MAPE) for point forecasts, a well-established error measure throughout the energy industry domain. Probabilistic error rates were assessed using the Pinball loss function [9]

, a well-known tool in both statistics and machine learning and quantile time series forecasting in particular


For a quantile forecast with as the target quantile, the pinball loss score is defined as:


where is the observation used for verification, and a = 1, 2, …, 99.

The pinball loss function returns a value that can be interpreted as the accuracy of a quantile forecasting model. Being a loss function the lower rates are better.

V Framework structure

V-a Module overview

The framework consists of modules responsible for scraping the data, storing it in a NoSQL database, building GBRT models and using the models to issue timely forecasts (see Figure 1). The main output of the framework is electric load value predictions for the next hours for the selected countries. The prediction function launches subprocesses with prediction horizon set between and hours, it is also capable of parallel execution. Each subprocess consists of three main parts; data preparation, model building (if model rebuild is requested) and model application.

Fig. 1: Block diagram of framework modules

The end result of the prediction function is a dataset of predicted hourly load values ahead of the test time for the given country. This means that if the forecasting system is in continuous operation, there will be forecasts for each hour, with horizons decreasing from to .

V-B Data preparation

First, Vertical Load and total load values and the calendar of the target country are loaded from the NoSQL database. Total load means a load equal to generation and any imports deducting any exports and power used for energy storage. While Vertical Load equals the total amount of power flowing out of the transmission network to the distribution networks, to directly connected final customers or to the consuming part of generation[11]. There is an issue that total load values are not available prior to , but in previous years substituting total load with Vertical Load data improves performance. This is a curious behaviour as Vertical Load and Total Load are inherently different concepts. Needless to say, when enough Total Load data has been acquired this tweak no longer will be necessary. For some countries, Total Load values are available in quarter-hourly or half-hourly granularity, so the data preparation module aggregates them to uniform hourly values.

Next, feature vectors are created out of the load value observations. Calculating the feature vectors involves both calendar variables and lagged values. Do note, although desired at this point no weather data is available to be used. Currently two lagged values are used; the load value one week before the predicted time slot and the last known load value (adjusted by the forecast horizon). If a lagged load value happens to be missing in the training data that row is discarded entirely. It is crucial that no test row has missing values, so if a lagged load is missing the data preparation function imputes it using values from previous periods. If no lagged values are found for the past 6 weeks, the module switches to the more robust setting where no lagged load variables are used.

V-C Model building

This part is executed only if the model rebuild is requested, during normal operations model rebuild is scheduled for every midnight. There are two distinct models built: basic and advanced. Basic models are constructed on calendar-based features only, and are used as fall back method. Advanced models use lagged load features in addition meaning their performance is much better. The underlying algorithm used for both models is Gradient Boosted Regression Trees. For each country the basic version builds 50 trees with a depth of 5, while the up-scaled version has twice as many trees with larger depth. All generated models are saved in the NoSQL database along with their metadata.

V-D Model application

In this step predictions for each test row are generated by applying the most recent model in the database that matches the country, horizon and model type (basic or advanced).

Vi Results

Vi-a Model evaluation & result comparison

To justify the usage of a new approach in time series forecasting it is often benchmarked with an established methodology like ARMA on a set of time series in the respective domain [2]. In this particular application, if provided in the open access data, the country issued estimates of energy consumption are used to assess the relative accuracy of the GBRT methodology. These estimates are issued by the respective country as part of their transparency regulations. The framework provides 24 estimates for every given time point with a forecast horizon spanning between 1 and 24 hours. Day ahead forecasts are typical in energy industry so if otherwise not stated all evaluation was performed on the 24 hours horizon.

Performance is measured using the Mean Absolute Percentage Error (MAPE) and model forecasts are compared to Actual Load values directly, and error measures are compared to country benchmark model forecasts as seen on the ENTSO-E Transparency Platform. Figure 2 shows how different time series stack up in the sample period of late March.

Fig. 2: Sample forecasted load, actual load & benchmark for Germany in March 2015
Mar Apr May Jun Jul Aug Sep Oct
BE 4.481 3.474 5.232 3.347 5.309 4.817 3.846 4.335
CZ 3.409 2.958 3.395 2.275 4.587 4.242 3.582 5.564
DE 5.478 7.171 6.222 3.146 5.606 4.288 3.822 5.981
FI 4.232 2.738 3.347 3.803 3.339 2.624 5.437 5.841
FR 6.237 3.483 4.664 3.384 5.485 4.954 4.461 8.671
HU 2.603 3.278 3.270 3.102 6.107 9.859 4.577 4.053
IT 4.279 8.589 5.891 4.961 6.191
NO 5.014 3.618 9.551 4.056 5.262 4.907 7.941 5.554
PL 3.181 3.597 4.505 2.314 3.953 3.929 2.985 3.814
SK 3.567 3.610 3.030 2.281 5.069 3.711 3.407 3.570
TABLE III: Measured error by country & month

Table III summarizes measured MAPE over the test period for each country. The GBRT framework generally performs well, but accuracy varies from country to country, the best being Czech Republic, Poland, Slovakia and to some extent Hungary. The model coped very well with the different seasons; all 3 have very similar mean error rates. Do note, that there was no Actual Load data published for Italy from August and onwards.

Performance measurements are put into context when looking at the country provided estimates used as benchmark here (see Table IV). Observed MAPE rates are typically lower here, with the sole exception of Czech Republic. Cases where our GBRT framework offers comparable results include Czech Republic, Hungary and to a lesser extent Germany and Slovakia. Scandinavian and Benelux countries operate with extremely low error rates. At this point, very little is known of the models and data used by any of the countries. Though country estimates could potentially (and should) incorporate additional information besides historical and calendar-related data such as weather forecasts, longer spanning historical data, customer data, market-related data. Again, providing at least partial open-access to this data could greatly boost development in this field.

Mar Apr May Jun Jul Aug Sep Oct
BE 1.848 1.964 1.694 1.377 2.139 3.426 1.434 1.709
CZ 5.429 5.898 4.638 5.360 6.532 5.935 6.897 5.374
DE 2.886 3.744 3.299 5.958 4.711 4.042 2.962 3.617
FI 1.368 1.761 1.657 1.840 1.287 1.405 2.934 1.632
FR 1.984 1.690 1.640 1.494 1.644 1.712 1.540 1.701
HU 3.159 4.371 4.294 3.698 4.577 4.864 4.551 4.468
IT 1.617 1.980 1.633 1.683 2.043
NO 0.092 0.357 0.165 0.092 0.506 2.102 4.542 1.990
PL 1.643 1.915 1.585 1.560 2.176 2.338 2.147 1.725
SK 2.152 2.318 2.267 2.199 2.773 2.258 2.560 2.423
TABLE IV: Benchmark error by country & month

Looking at the error statistics in Table V

one can see that all error distributions are naturally right skewed suggesting a distribution close to log-normal, meaning the GBRT powered framework more frequently commits smaller errors than higher ones. The framework performance is rather robust over countries, major error percentiles are well aligned with the distribution tail having varying length. Figure

3 shows sample distributions for Germany comparing framework and benchmark errors. The GBRT errors demonstrate a much longer tail indicating that some patterns or information might be missing from this representation. On the other hand errors close to zero are also more frequent in the model output.

Fig. 3: Model output and benchmark error distributions for Germany
mean std min 25% 50% 75% max
BE 4.325 5.265 0.002 1.247 2.714 5.066 47.422
CZ 3.904 6.116 0 0.832 1.97 4.144 53.674
DE 5.17 8.015 0.001 1.267 2.747 5.307 116.77
FI 3.98 5.612 0 0.994 2.103 4.211 42.351
FR 5.321 8.482 0 1.138 2.575 5.404 76.028
HU 4.924 6.318 0.001 0.994 2.476 5.944 39.584
IT 5.924 10.815 0.001 1.154 2.588 5.724 88.475
NO 5.567 9.122 0 1.185 2.604 5.315 79.224
PL 3.571 6.15 0.002 0.754 1.724 3.793 59.342
SK 3.607 4.619 0.001 0.911 2.061 4.326 40.367
TABLE V: MAPE error statistics by country

The framework provides 24 estimates for every given time point, shorter horizons often produce better results as underlined by Table VI. Forecast errors between 1 and 24 hours length typically draw an arc, with higher error rates in the middle and the lowest on the edges. But the distribution is far from symmetrical; a 1-length horizon typically yields approximately 12% lower error rates than its 24 hour counterpart.

1 3.70 3.71 4.47 3.44 4.79 4.82 5.27 4.34 3.43 3.43
2 3.96 3.79 4.73 3.60 4.86 4.88 5.92 4.54 3.49 3.47
3 4.23 3.83 4.86 3.72 4.91 4.85 6.15 4.73 3.53 3.53
4 4.37 3.91 4.99 3.78 4.97 4.81 6.36 4.88 3.57 3.55
5 4.47 3.93 5.11 3.86 5.03 4.79 6.53 4.95 3.62 3.59
6 4.52 3.96 5.21 3.91 5.09 4.81 6.69 5.08 3.65 3.60
7 4.68 3.98 5.31 3.97 5.11 4.81 6.77 5.10 3.66 3.65
8 4.75 3.96 5.37 4.02 5.17 4.82 6.82 5.16 3.66 3.61
9 4.78 4.02 5.41 4.04 5.19 4.82 6.76 5.24 3.63 3.62
10 4.86 4.01 5.44 4.05 5.21 4.83 6.71 5.28 3.62 3.62
11 4.87 4.02 5.45 4.09 5.22 4.81 6.70 5.37 3.64 3.60
12 4.95 4.04 5.46 4.10 5.24 4.77 6.89 5.42 3.63 3.62
13 4.92 4.06 5.51 4.09 5.23 4.78 6.88 5.44 3.61 3.62
14 5.07 4.02 5.46 4.15 5.25 4.75 6.96 5.43 3.62 3.63
15 5.01 4.03 5.48 4.14 5.24 4.74 6.93 5.47 3.65 3.64
16 4.98 4.04 5.45 4.14 5.27 4.76 7.02 5.55 3.61 3.60
17 4.90 4.03 5.43 4.17 5.29 4.77 6.98 5.51 3.62 3.61
18 4.88 3.99 5.39 4.16 5.28 4.76 6.97 5.51 3.64 3.61
19 4.78 4.00 5.38 4.10 5.33 4.78 6.89 5.51 3.67 3.63
20 4.67 3.96 5.32 4.11 5.38 4.80 6.85 5.51 3.65 3.61
21 4.65 3.99 5.29 4.06 5.33 4.83 6.73 5.46 3.67 3.58
22 4.54 3.99 5.28 4.04 5.33 4.87 6.60 5.49 3.64 3.58
23 4.46 3.93 5.23 4.03 5.29 4.89 6.31 5.49 3.58 3.60
24 4.33 3.90 5.17 3.98 5.32 4.92 5.92 5.57 3.57 3.61
TABLE VI: Measured MAPE by horizon & country
Fig. 4: Probabilistic forecasts for Hungary in the first week of September

Furthermore we investigate the forecasting of a probability distribution of the target variable in contrast to forecasting a single value. Probabilistic forecasts provide more detailed insight for stakeholders offering immediate business value. For this purpose we utilize the quantile regression capabilities of the GBRT model and estimate target variable distribution over the major deciles. Probabilistic forecasts were experimented with Hungary only for the month of September, training the model on the period of January 2012 - August 2015. Figure

4 shows a sample of the estimated time series distribution along with the Actual Load. The quantile intervals here are relatively narrow on most days, with slight bias issues in some cases.

Figure 5 shows the probabilistic Pinball Loss distribution over the quantiles. The distribution is almost symmetrical with higher rates observed close to the median. The average pinball loss measured for this setting is . Unfortunately with only point forecasts published on the portal there is no direct way to compare the loss to country estimates. We look forward to probabilistic forecasts being released in the near future as they provide more detailed insight thus generating business value.

Fig. 5: Distribution of pinball loss over the quantiles
Fig. 6: Displaying forecasts on the web directly from the API

Vii Future work & Conclusions

The framework presented in this paper is functional and fully operational, yet there are still plenty of ideas to further extend it. As mentioned in the paper, at the moment no external data sources are used. Previous experience shows that incorporating country specific weather forecasts can greatly boost prediction accuracy and provide richer data sources encapsulated in the framework. ENTSOG integration, the gas counterpart of ENTSO-E, is also a reasonable next step as gas usage is even more weather dependent than electricity. To further justify the usage of GBRT in time series forecasting there is a need to benchmark it to an established methodology like ARMA on a relatively large set of time series in the energy domain in a rolling fashion. Country specific load data is particularly suitable for this purpose.

Benchmarking is not the sole purpose of this paper, as it deals with the data analysis process as a whole from ETL to data preparation (ENTSO-E integration, data cleansing and data auditing), through modeling and result presentation (see Figure 6) with an API providing parametric forecasts for a number of countries in Europe. It will hopefully contribute to the agenda of releasing more open access data to the public that can be put to good use in the future.


  • [1] E. Mangalova, E. Agafonov, Wind power forecasting using the k-nearest neighbors algorithm, International Journal of Forecasting 30 (2) (2014) 402–406.
  • [2] G. Barta, G. Borbely, G. Nagy, S. Kazi et al., GEFCOM 2014 - Probabilistic Electricity Price Forecasting Intelligent Decision Technologies, 39 (2015) 67–76.
  • [3] S. Aggarwal, L. Saini, Solar energy prediction using linear and non-linear regularization models: A study on ams (american meteorological society) 2013–14 solar energy prediction contest, Energy 78 (C) (2014) 247–256.
  • [4] T. Hong, P. Pinson et al. Probabilistic energy forecasting: State-of-the-art 2015, International Journal of Forecasting.
  • [5] Central collection and publication of electricity generation, transportation and consumption data and information for the pan-European market. [online] 2015, (Accessed: 20 November 2015).
  • [6] J. H. Friedman, Stochastic gradient boosting, Comput. Stat. Data Anal. 38 (4) (2002) 367–378.
  • [7] G. Nagy, Gefcom2012: Wind power forecasting using gradient boosting machines, in: Power and Energy Society General Meeting, Power and Energy Society, 2012, p. GM2555.
  • [8] J. H. Friedman, ”Greedy function approximation: a gradient boosting machine,” Annals of Statistics, pp. 1189-1232, 2001.
  • [9] R. Koenker and G. Bassett Jr., “Regression quantiles,” Econometrica, vol. 46, pp. 33–50, 1978
  • [10] Biau, Gérard, and Benoît Patra. ”Sequential quantile prediction of time series.” Information Theory, IEEE Transactions on 57.3 (2011): 1664-1674.
  • [11] ENTSO-E Metadata Repository Glossary. [online] 2015, (Accessed: 21 November 2015).