Normal Distribution
nothing is called 95% until it has been counted
Concept

Covariance estimation — where it appears

Estimating how the errors of a model vary together, from data in which the errors themselves are not observed. In a regression on dependent errors it is read off the residuals, which carry less dependence than the errors did, and the estimate then sets the weights of the fit.

Named by 4 essays across 3 fields — each of them below, with the objects they name alongside it.

What a lag of band costs, as the sample grows. The optimism one more lag of a Bartlett band adds, measured by fitting on one sample and scoring on an independent one, 2,400 pairs at each size, with two standard errors: 0.3745 ± 0.0176 at 120 rows, 0.3842 ± 0.0266 at 240 rows, 0.4523 ± 0.0392 at 480 rows, 0.4022 ± 0.0563 at 960 rows. The window's summed weights give one half a lag and the convention charges one.

The charge at four times the sample

A band of lags is charged a full unit apiece and is worth about three eighths of one at a hundred and twenty rows. The field's last essay said the error in that convention was harmless only because the error surface is flat, and that a steeper surface would give it a cost. A longer sample was the obvious place for one. It is not there: from 120 rows to 480 the error's spread across band widths falls by three quarters, and the derived charge's advantage over the convention falls from 5.5 standard errors to 1.7. What the sample does move is the charge itself, from 0.37 of a unit a lag towards the half a unit the window's weights predict.

dimension · Criterion Past the cap, long memory's surface is flat. The prediction error at every band width from four lags to eighty, above the lowest on the extended list, under long memory at d = 4/9 at three sample sizes, 600 draws each. At 120 rows the floor is at 12 lags and eighty is 0.0039 above it. At 480 the floor is at 30 lags and every width from twenty-four to eighty is within 0.0003 of it. At 960 it is at 48 and thirty to eighty lags are within 0.0000. The dashed line is the list's cap: the surface does not rise behind it, it levels off.

A cap that decides nothing

Every band-width list this field priced stopped at thirty lags, and long memory at a large sample is the one dependence whose best band wants to be wider. At 960 rows it does: the draw's own best width lies past the cap on 48.7% of draws. Lift the cap to eighty and the measured charge follows, from 28 lags to 46 — and the error moves by two ten-thousandths, in the wrong direction and inside its standard error, because the surface behind the cap is flat. The field's last conditional fails where it had its best chance.

dimension · Criterion What iterating the band is worth. How much lower the prediction error is when the Bartlett band is iterated to its fixed point than when it is read once off the least-squares residuals, paired on the draw, with two standard errors, 300 draws at 120 rows for each law and window. First-order autoregression: 4 lags 2.84 ± 0.36, 8 lags 3.43 ± 0.58, 20 lags 4.16 ± 0.80; five-period moving average: 4 lags 2.04 ± 0.28, 8 lags 2.69 ± 0.42, 20 lags 4.32 ± 0.68; long memory: 4 lags 0.83 ± 0.14, 8 lags 1.13 ± 0.25, 20 lags 1.33 ± 0.44 thousandths. Every one is below zero by more than two and a half standard errors; the gain is largest at the widest band under the two short-memory laws.

A fixed point with no likelihood under it

A tapered band of autocovariances has no likelihood to maximise, so the one rule that can fit its dependence together with the regression is to iterate: read the band off the residuals, refit, read it again. It converges on every draw, in seven to eleven fits at the median, gives back about half the memory the first fit took out of the residuals under a short-memory law, and lowers prediction error by two to four thousandths — five to eight paired standard errors. Under the law a first-order likelihood describes, its maximum still wins; under the two it does not, the band's fixed point beats it.

together · Dependence The floor a repair has to choose. Regret of the truncated covariance repaired by flooring its eigenvalues, at a persistence of 0.85 on 120 rows, 120 draws, by floor, on a log scale. Window 4: 1.5649 at 0.001, 0.0945 at 0.01, 0.0385 at 0.05, 0.0285 at 0.1, 0.0226 at 0.2, 0.0235 at 0.4, 0.0292 at 0.7; the taper at the same window 0.0275. Window 8: 1.1919 at 0.001, 0.0980 at 0.01, 0.0334 at 0.05, 0.0258 at 0.1, 0.0197 at 0.2, 0.0209 at 0.4, 0.0249 at 0.7; the taper at the same window 0.0178. Window 12: 0.8692 at 0.001, 0.0621 at 0.01, 0.0239 at 0.05, 0.0209 at 0.1, 0.0212 at 0.2, 0.0239 at 0.4, 0.0282 at 0.7; the taper at the same window 0.0150. Window 20: 1.7110 at 0.001, 0.2495 at 0.01, 0.0399 at 0.05, 0.0238 at 0.1, 0.0202 at 0.2, 0.0219 at 0.4, 0.0259 at 0.7; the taper at the same window 0.0135. Least squares gives up 0.0890 and the whitening told the truth 0.0075.

A covariance repaired

A truncated autocovariance sequence fails to be a covariance matrix on most draws of a persistent series, and this line of essays has always scored it on the draws where it survives. The textbook repair is to floor its eigenvalues. At a floor of 0.001 that moves the matrix by four per cent and gives up 1.56 of regret, seventeen times what least squares gives up, because a whitening reads an inverse and an inverse reads the smallest eigenvalues. At a floor of 0.2 it beats the taper at the shortest windows. So the repair is not a repair: it is a second dial.

banded · Dependence

Named alongside it

The objects these essays reach for when they reach for this one.

Generalised least squaresBandwidthModel selectionTaperingEffective degrees of freedomFeasible generalised least squaresInformation criterionLong memoryOptimismRegretSample sizeAutocorrelation

All concepts