CumHDD, which we then use to estimate the density, as follows. First, we simulate 250 151-day
realizations of the temperature shock εr by drawing with replacement from the empirical distribution of
estimated temperature shocks (εf). Second, we run the 250 151-day realizations of temperature shocks
through the estimated model (1) to obtain 250 simulated 151-day realizations of daily average temperature.
Third, we convert the 250 simulated 151-day realizations of daily average temperature into 250 simulated
151-day realizations of HDD, which we cumulate over the November-March heating season,
CumHDDs = 5~^11 HDDty, s = 1, 2, ..., 250. Finally, we form the empirical distribution function of
(JumelDIp, based upon Cu=HDDs, s = 1, ..., 250.
After passing through the entire sample, we have 41 assessed distribution functions, Fy(∙∖
y o I960,..., 2000, one governing each of CumHDDy, y = I960,..., 2000. We assess the conditional
calibration of those distributional forecasts via the probability integral transform, as suggested and
extended by Rosenblatt (1952) and extended by Dawid (1984), Diebold, Gunther and Tay (1998), and
Diebold, Hahn and Tay (1999). In particular, if the estimated distribution and true distribution coincide
year-by-year, then the series of assessed distribution functions Fy(∙) evaluated at the corresponding series
of realized values of CumHDDy should be approximately iid and uniformly distributed on the unit interval.
. Hd
Formally, z ≡ Fj,(C∞wffi)Zζ,) ~ t7(0, 1). For each city, we check uniformity by examining histograms
of z, and we check independence by examining correlograms of the first four powers of z. The sample of
size 41 is of course small, but the framework has previously been applied successfully in small samples, as
for example by Diebold, Tay and Wallis (1999).
First consider assessing uniformity. We estimate the density of z using simple four-bin histograms,
which we present in the leftmost column of Figure 9, accompanied by 95% pointwise error bands under the
iid U(0, 1) null hypothesis. Interestingly, the z series differ rather noticeably from uniformity, and
moreover, they display a common pattern: too many large CumHDD realizations occur relative to the
assessed distributions, as evidenced by the increase in the histograms when moving from left to right. The
common nature of uniformity violations may indicate a neglected common temperature component, due for
-14-