  In this paper we address the issue of modeling electricity loads. After
analyzing properties of the deseasonalized loads from the California power
market we fit an ARMA(1,6) model to the data. The obtained residuals seem to be
independent but with tails heavier than Gaussian. It turns out that the
hyperbolic distribution provides an excellent fit.
