--- id: JoelLewis/finance_skills/historical-risk version: "bece6a5b" license: MIT install: manual updated: 2026-07-18 --- # historical-risk — Quantify how risky an investment has been using historical price data. This skill computes annualized volatility, maximum drawdown, recovery time, historical VaR, downside deviation, and tracking error—supporting multiple volatility estimators including Parkinson and Yang-Zhang methods. Use it to analyze peak-to-trough losses and asymmetric downside risk. Publisher: JoelLewis · Stars: 159 · Updated: 2026-07-18 Install (manual): `git clone https://github.com/JoelLewis/finance_skills` ## SKILL.md # Historical Risk Analysis ## Core Concepts ### Close-to-Close Volatility The simplest and most common volatility estimator. Compute the standard deviation of log returns and annualize. ``` sigma_annual = sigma_daily * sqrt(N) ``` where N = number of trading periods per year (typically 252 for daily, 52 for weekly, 12 for monthly). Log returns are preferred: r_t = ln(P_t / P_{t-1}). ### Parkinson (High-Low) Estimator Uses intraday high and low prices to capture intraday volatility that close-to-close misses. More efficient than close-to-close when the true process is continuous. ``` sigma^2_Park = (1 / (4 * n * ln(2))) * sum( ln(H_i / L_i)^2 ) ``` This estimator is roughly 5x more efficient than close-to-close for a diffusion process, but is biased downward when there are jumps or when the range is discretized. ### Yang-Zhang Estimator Combines overnight (close-to-open), open-to-close, and Rogers-Satchell components. It is unbiased for processes with both drift and opening jumps. ``` sigma^2_YZ = sigma^2_overnight + k * sigma^2_open-to-close + (1 - k) * sigma^2_RS ``` where k is chosen to minimize estimator variance: ``` k = 0.34 / (1.34 + (n + 1) / (n - 1)) ``` with n the number of observations, and sigma^2_RS is the Rogers-Satchell estimator that uses all four OHLC prices within each period. ### Drawdown Analysis Drawdown at time t measures the decline from the running peak: ``` DD_t = (Peak_t - Value_t) / Peak_t ``` where Peak_t = max(Value_s) for all s <= t. - **Maximum Drawdown (MDD):** MDD = max(DD_t) over the evaluation period. - **Drawdown Duration:** The number of periods from a peak until a new peak is reached. - **Recovery Time:** The number of periods from the trough back to the prior peak level. ### Historical VaR The non-parametric (empirical) Value-at-Risk is simply the alpha-percentile of the historical return distribution. No distributional assumptions are made. ``` VaR_alpha = -Percentile(R, alpha) ``` For example, 95% VaR uses the 5th percentile of returns. The negative sign is a convention so that VaR is expressed as a positive loss number. ### Downside Deviation Measures dispersion of returns below a Minimum Acceptable Return (MAR): ``` sigma_d = sqrt( (1/n) * sum( min(R_i - MAR, 0)^2 ) ) ``` Common choices for MAR: 0%, the risk-free rate, or the mean return. ### Tracking Error Standard deviation of the difference between portfolio and benchmark returns, annualized: ``` TE = std(R_p - R_b) * sqrt(N) ``` This measures how consistently the portfolio tracks (or deviates from) its benchmark. ### Semi-Variance Variance computed using only returns below the mean (or below a threshold): ``` SV = (1/n) * sum( min(R_i - mean(R), 0)^2 ) ``` Semi-variance isolates downside risk and is the foundation for the Sortino ratio (see performance-metrics). ## Key Formulas | Formula | Expression | Use Case | |---------|-----------|----------| | Annualized Volatility | sigma_ann = sigma_period * sqrt(N) | Convert period vol to annual vol | | Log Return | r_t = ln(P_t / P_{t-1}) | Compute continuously compounded returns | | Parkinson Variance | sigma^2 = (1 / (4n ln2)) * sum(ln(H/L)^2) | Volatility from high-low data | | Drawdown | DD_t = (Peak_t - Value_t) / Peak_t | Measure peak-to-trough decline | | Max Drawdown | MDD = max(DD_t) | Worst historical decline | | Historical VaR (95%) | 5th percentile of return series | Non-parametric loss estimate | | Downside Deviation | sigma_d = sqrt((1/n) * sum(min(R_i - MAR, 0)^2)) | Asymmetric risk below MAR | | Tracking Error | TE = std(R_p - R_b) * sqrt(N) | Portfolio vs benchmark deviation | | Semi-Variance | (1/n) * sum(min(R_i - mean(R), 0)^2) | Below-mean variance | ## Worked Examples ### Example 1: Annualized Volatility from Daily Returns **Given:** A stock has daily log returns with a sample standard deviation of 1.2%. Assume 252 trading days per year. **Calculate:** Annualized volatility. **Solution:** ``` sigma_annual = 0.012 * sqrt(252) = 0.012 * 15.875 = 0.1905 ~ 19.05% ``` The stock's annualized volatility is approximately 19%. ### Example 2: Maximum Drawdown from a Price Series **Given:** A fund's NAV follows this path over six months: $120, $135, $150, $130, $105, $125. **Calculate:** Maximum drawdown and identify the peak and trough. **Solution:** Running peaks: $120, $135, $150, $150, $150, $150. Drawdowns at each point: - $120: (120-120)/120 = 0% - $135: (135-135)/135 = 0% - $150: (150-150)/150 = 0% - $130: (150-130)/150 = 13.3% - $105: (150-105)/150 = 30.0% - $125: (150-125)/150 = 16.7% **Maximum Drawdown = 30.0%**, occurring from the peak of $150 to the trough of $105. As of the last observation ($125), the drawdown has not yet fully recovered. ### Example 3: Historical VaR **Given:** 500 daily returns sorted from worst to best. The 25th-worst return is -2.8% and the 26th-worst is -2.6%. **Calculate:** 95% 1-day historical VaR. **Solution:** The 5th percentile corresponds to the 25th observation out of 500 (500 * 0.05 = 25). ``` VaR_95% = -(-2.8%) = 2.8% ``` Interpretation: On 95% of days, the loss is expected not to exceed 2.8% based on the historical distribution. ## Common Pitfalls - **Not annualizing volatility correctly:** Volatility scales with the square root of time (multiply by sqrt(N)), not linearly. Multiplying daily vol by 252 instead of sqrt(252) produces wildly inflated numbers. - **Using calendar days vs trading days inconsistently:** Use 252 trading days (not 365 calendar days) for equity markets when annualizing. Bond markets and some international markets may differ. - **Survivorship bias in historical data:** Data sets that exclude delisted or failed securities understate realized risk. - **Lookback period sensitivity:** A 1-year lookback captures different risk regimes than a 5-year lookback. Always state the lookback window and consider whether it spans relevant market conditions. - **Confusing VaR confidence level direction:** 95% VaR corresponds to the 5th percentile of returns (the loss tail). The "95%" refers to the confidence level, not the percentile of gains. - **Log returns vs simple returns:** For volatility estimation, log returns are preferred because they are additive across time. For reporting cumulative performance, simple returns are more intuitive. ## Cross-References - **performance-metrics** (wealth-management plugin): Uses volatility, downside deviation, tracking error, and max drawdown as denominators in risk-adjusted ratios (Sharpe, Sortino, Information Ratio, Calmar). - **forward-risk** (wealth-management plugin): Historical VaR and historical volatility serve as inputs to forward-looking VaR models and stress tests. - **volatility-modeling** (wealth-management plugin): EWMA and GARCH models extend the simple historical volatility estimators covered here into forecasting frameworks. - **retirement-decumulation** (wealth-management plugin): drawdown and volatility history drives sequence-of-returns risk in withdrawal-phase portfolios ## Running the script Run with `uv run scripts/historical_risk.py` (the PEP 723 header resolves numpy automatically) or with `python3 scripts/historical_risk.py` after `pip install numpy scipy`. A bare run prints a full risk analysis (annualized and Parkinson volatility, maximum drawdown with timing, 95%/99% historical VaR, downside deviation, semi-variance, tracking error, rolling volatility) on seeded synthetic data with an injected drawdown event. Use `--verify` to assert outputs match this skill's worked examples and the demo's expected values (exit code 0 on PASS) and `--help` for an overview of the class. The file is primarily meant to be imported as a module (e.g., `from historical_risk import HistoricalRiskAnalyzer`). [View on SkillFed](https://skillfed.io/JoelLewis/finance_skills/historical-risk) · [View on GitHub](https://github.com/JoelLewis/finance_skills)