SPB Git

spb/wp29_uqo

Public
3.2 KB · 50 lines markdown
1# F2-conf-02 — Is the memo's slope really significant?23## Question4The memo estimates y_t = α + β x_t + u_t by OLS on 400 monthly observations and, using conventional (iid) standard errors, declares β significant at 5%. The disturbances are serially correlated, so we re-estimate β with a HAC (Newey–West, Bartlett kernel) standard error and assess the claim.56## Data and exclusions7`data/series.csv`: 400 rows read, 1990-01 to 2023-04. Checks: no missing values, no duplicate months, no gaps in the monthly sequence. No exclusions; n_analysed = 400. The data are simulated (synthetic), not observed. No data issues were found.89## Method10- OLS of y on x with an intercept (statsmodels).11- Conventional iid SE reproduced as the memo's baseline.12- Residual first-order autocorrelation ρ̂ = Σ û_t û_{t-1} / Σ û_t².13- HAC SE: Newey–West with Bartlett kernel, small-sample df correction applied. **Bandwidth choice:** Andrews (1991) AR(1) plug-in for the Bartlett kernel, applied to the slope score v_t = (x_t − x̄) û_t: α̂ = 4ρ_v²/((1−ρ_v)²(1+ρ_v)²), S = 1.1447(α̂T)^{1/3}. This gives S = 8.18, rounded up to **L = 9 lags**. I preferred this data-dependent rule over the fixed Newey–West (1994) rule of thumb, floor(4(T/100)^{2/9}) = 5, because the errors are known to be persistent. The sensitivity table below shows that the conclusion does not depend on this choice.14- Inference uses the asymptotic normal distribution: two-sided p-value and 95% CI = β̂ ± 1.96·SE_HAC.1516## Results17| quantity | value |18|---|---|19| α̂ (intercept) | 2.0925 |20| β̂ (slope) | 0.1017 |21| SE iid (memo baseline) | 0.0448 (t = 2.27, p = 0.024; this reproduces the memo exactly) |22| residual ρ̂₁ | 0.609 |23| SE HAC (Bartlett, 9 lags) | 0.0654 |24| t_HAC | 1.56 |25| p_HAC (two-sided) | 0.120 |26| 95% CI (HAC) | [−0.026, 0.230] |2728I reproduced the memo's point estimates and its iid SE. After allowing for serial correlation, however, the SE is about 46% larger, the t-statistic falls to 1.56 and the 95% CI includes zero. **The memo's claim that the slope is significant at 5% is not supported.**2930## Uncertainty31Sensitivity of the HAC SE to the bandwidth. All values come from `out/diagnostics.json`.3233| lags | SE | t | p (normal) | p (t, 398 df) |34|---|---|---|---|---|35| 0 (White, heterosk. only) | 0.0490 | 2.08 | 0.038 | 0.039 |36| 5 (NW rule of thumb) | 0.0645 | 1.58 | 0.115 | 0.116 |37| 8 | 0.0653 | 1.56 | 0.119 | 0.120 |38| 9 (chosen) | 0.0654 | 1.56 | 0.120 | 0.121 |39| 12 | 0.0657 | 1.55 | 0.121 | 0.122 |40| 18 | 0.0663 | 1.54 | 0.125 | 0.125 |41| 24 | 0.0660 | 1.54 | 0.123 | 0.124 |4243The slope is significant only if serial correlation is ignored (0 lags). For every bandwidth from 5 to 24 lags, p is between 0.115 and 0.125.4445## Limitations and open issues46- Newey–West tests with Bartlett kernels tend to over-reject in finite samples when the errors are persistent. Fixed-b or other more conservative inference would probably push p even higher, so it would strengthen the conclusion rather than reverse it.47- The Andrews plug-in assumes an AR(1) approximation for the score.48- A non-rejection here does not show that β = 0. The CI is wide and includes economically non-trivial positive values.49- β is a linear-projection coefficient, not a causal effect.50
Files · main
Menu

Theme