BLOG

Behind the numbers.

Code examples, market analysis, and data quality deep-dives.

Which Commodities Have the Strongest Momentum? Rotation Backtest in Python
Which Commodity ETFs Have the Worst Tail Risk? Expected Shortfall in Python
Are Gold Miners Leveraged Gold Bets? Rolling Beta Analysis in Python
Does the Base-Metals-to-Gold Ratio Lead Cyclical Stocks? Signal Test in Python
Can Risk Parity Tame Commodity Volatility? Portfolio Optimization in Python
Are Power Stocks Becoming an AI Infrastructure Trade? Momentum Screening in Python
Which AI Chip Stocks Have Margin Momentum? Profitability Trend Analysis in Python
Which AI Stocks Are Cheapest Relative to Growth? Growth-Adjusted Valuation in Python
Does AI Stock Leadership Persist? Momentum Backtest in Python
Which AI Stocks Have the Cleanest Balance Sheets? Net Cash Screening in Python
Can Risk Parity Reduce Mega-Cap Drawdowns? Portfolio Optimization in Python
Which Growth Stocks Are Self-Funding? Cash-Flow Quality Screening in Python
Which Sectors Struggle When the Dollar Rallies? Sector Rotation Analysis in Python
Do Cheap Stocks Hold Up When Bonds Sell Off? Valuation Rotation in Python
Does the Nasdaq 100 Have Better Growth Quality Than the Dow? Index Constituent Analysis in Python
Do Healthcare Cash-Flow Margins Predict Returns? Signal Evaluation in Python
Which Dividend Stocks Survive a Cash-Flow Stress Test? Dividend Screening in Python
Does Heavy Insider Selling Predict Weak Returns? Insider Flow Test in Python
Can Quality Screens Reduce Small-Cap Balance-Sheet Risk? Russell 2000 Test in Python
Which Retailers Have Positive Operating Leverage? Margin Screening in Python
Is MSTR a Leveraged Bitcoin Proxy? Rolling Beta Analysis in Python
Is Micron's Memory Cycle Recovering? Inventory and Margin Forecasting in Python
Which Sectors Work When Bonds Rally? Rate-Sensitive Rotation in Python
Do One-Month Price Extremes Reverse? Signal Evaluation in Python
Do Low-Volatility S&P 500 Stocks Reduce Drawdowns? Factor Test in Python
Is AI Capex Paying Back Fast Enough? Revenue Hurdle Forecasting in Python
Could Shorter AI Asset Lives Hit Earnings? Depreciation Stress Test in Python
How Much AI Capex Risk Can a Portfolio Remove? Constrained Optimization in Python
Is the AI Capex Trade Crowded? Rolling Volatility and Sector Rotation in Python
Did the AI Boom Come From Existing S&P 500 Members? Point-in-Time Momentum Test in Python
Is AI Revenue Circular? Customer-Vendor Capex Loop Analysis in Python
Is the AI Trade Connected to Private Credit? Rolling Correlation Network in Python
Is Apollo More Balance-Sheet Sensitive Than Peers? Leverage Screen in Python
Are AI Earnings Supported by Cash Flow? Accrual and Capex Screen in Python
Can Defensive Stocks Hedge AI Drawdowns? Basket Regime Test in Python
How Fast Does the Market Price In Fed Decisions? FOMC Event Study in Python
How Much Are Options Sellers Overpaid? The Variance Risk Premium in Python
Which Companies Have the Worst Earnings Quality? Sloan Accrual Screen with Geographic Revenue Data in Python
Does the Oil-to-Gold Ratio Signal Recessions? XLE/GLD Backtest in Python
Is AI Spending Crowding Out Free Cash Flow? Capex Sustainability Across the Mag 7 in Python
Does a Long Energy / Short Bonds Portfolio Capture Inflation Surprises? Factor Construction in Python
Can a Hidden Markov Model Detect Oil Market Regimes? HMM Analysis in Python
Do Grain Prices Predict Food Inflation? Granger Causality Test in Python
Does the Corporate Credit Spread Predict Stock Market Crashes? BAA-AAA Spread Analysis in Python
Do Oil Stocks Hedge Inflation? Rolling Beta Analysis in Python
Which Stocks Are Most Rate-Sensitive? Equity Duration via Bond Beta in Python
Which Companies Have the Highest Accrual Ratios? Earnings Quality Screening in Python
Is Alpha Persistent or Decaying? Rolling Sharpe Ratio Analysis in Python
Are Markets Trending or Mean-Reverting? Hurst Exponent Analysis in Python
Is Consumer Discretionary vs Staples a Leading Indicator? XLY/XLP Ratio Analysis in Python
Does Heavy Capex Predict Future Stock Returns? Capital Expenditure Analysis in Python
How to Estimate Cost of Equity Using CAPM in Python
Is Volatility Predictable? Testing for Volatility Clustering in Python
Which Industrials Are Overleveraged? Net Debt to EBITDA Screening in Python
GM Before and After Bankruptcy: Why Entity Resolution Matters for Financial Data
What Is Adjusted Beta? Merrill Lynch Beta Shrinkage in Python
How Good Is a Stock Pick? Information Ratio and Tracking Error in Python
Do Stock Returns Follow a Normal Distribution? Testing for Fat Tails in Python
Which Large Caps Have the Highest Free Cash Flow Yield? FCF Screening in Python
Which Sectors Won Over 5 Years? Sector Rotation Analysis in Python
How to Forecast Stock Volatility with GARCH Models in Python
Are Stock Prices Mean-Reverting? Augmented Dickey-Fuller Test in Python
How to Calculate CAPM Alpha and Beta with Regression in Python
How to Compare Sector Sharpe Ratios and Sortino Ratios in Python
DELL: Why Stitching Historical Price Data Together Is Wrong
How to Analyze Drawdown and Recovery for Bank Stocks in Python
How to Screen SaaS Stocks by Revenue Growth and Cash Flow in Python
How to Screen REITs by Dividend Yield and Valuation in Python
How Correlated Are the Magnificent 7? Intra-Group Correlation in Python
AAPL vs XOM: Do Individual Stocks Have Seasonal Patterns?
How to Rank Large-Cap Stocks by Momentum in Python
How to Build a Multi-Endpoint Financial Dashboard in Python
How to Compare Volatility Across Energy Stocks in Python
How to Screen Healthcare Stocks by Valuation in Python
How to Build a Sector Correlation Matrix for Portfolio Diversification in Python
How to Find Oversold and Overbought Stocks Using Z-Scores in Python
How to Measure Earnings Quality: Cash Flow vs Net Income in Python
How to Build a Multi-Factor Stock Screen in Python (Value + Momentum + Quality)
How to Build a Simple DCF Model for Any Stock in Python
How to Screen Tech Stocks by Revenue Growth in Python
How to Screen Stocks by Balance Sheet Health in Python
Is "Sell in May" Real? SPY Monthly Seasonality Over 10 Years
How to Compare Sector Performance YTD Using Python
How to Track S&P 500 Additions and Removals Over Time in Python
How to Screen Dividend Stocks by Yield and Quality in Python
How to Calculate Max Drawdown and Recovery Time for Any Stock in Python
How to Compare Profitability Across Mega-Cap Tech Stocks in Python
Why Ticker Symbols Are Unreliable: The Recycling Problem Every Quant Should Know
How to Calculate and Compare Stock Volatility in Python
How to Screen Blue-Chip Stocks by P/E Ratio in Python
How to Track Companies Through Ticker Changes, Bankruptcies, and Renames in Python
S&P 500 Turnover: How Much the Index Has Changed Since 2010
How to Calculate Stock Beta and Correlation in Python
← All articles

Which Commodity ETFs Have the Worst Tail Risk? Expected Shortfall in Python

What's the question?

Volatility is an average measure of uncertainty. It treats upside and downside movement symmetrically. Commodity investors often need something more specific: how bad are the worst days?

Expected shortfall answers that question. It measures the average return inside the worst part of the return distribution. A 5% expected shortfall is the average daily return among the worst 5% of days. It is stricter than value at risk, which only marks the cutoff point. Expected shortfall asks what happens after that cutoff is breached.

The question is which commodity ETFs carry the most severe downside tails, and whether the ranking differs from ordinary volatility.

The approach

The universe is GLD, SLV, USO, UNG, DBA, DBB, and DBC. Built from SEC EDGAR public filings and market data, the screen uses five years of daily returns.

  1. Pull five years of daily returns
  2. Compute annualized volatility
  3. Compute the 5% value-at-risk cutoff
  4. Compute expected shortfall as the average return below that cutoff
  5. Compare maximum drawdown, worst day, skewness, and kurtosis

Skewness measures whether extreme moves are tilted positive or negative. Kurtosis measures whether the distribution has unusually large outliers relative to a normal distribution.

Code

import xfinlink as xfl
import pandas as pd
import numpy as np

xfl.set_api_key("YOUR_API_KEY")  # free at https://xfinlink.com/signup

tickers = ["GLD", "SLV", "USO", "UNG", "DBA", "DBB", "DBC"]
prices = xfl.prices(tickers, period="5y", fields=["return_daily"])
returns = prices.pivot_table(index="date", columns="ticker", values="return_daily").dropna()

rows = []
for ticker in tickers:
    r = returns[ticker]
    cutoff = r.quantile(0.05)
    tail = r[r <= cutoff]
    wealth = (1 + r).cumprod()
    rows.append({
        "ticker": ticker,
        "volatility": r.std() * np.sqrt(252),
        "var_5": cutoff,
        "expected_shortfall_5": tail.mean(),
        "max_drawdown": (wealth / wealth.cummax() - 1).min(),
        "kurtosis": r.kurt(),
    })

print(pd.DataFrame(rows).sort_values("expected_shortfall_5"))

Full script with formatting and visualisation: commodity-tail-risk-expected-shortfall-python.py

Output

Expected shortfall ranking for commodity ETFs across the worst five percent of daily returns
=== Commodity ETF Tail Risk Screen ===
Universe: GLD, SLV, USO, UNG, DBA, DBB, DBC
Sample: 2021-06-21 to 2026-06-17 (1254 trading days)
Tail metric: expected shortfall across the worst 5% of daily returns

Tail-risk ranking:
UNG   vol=64.1%  VaR_5= -6.5%  ES_5= -8.6%  max_drawdown=-92.5%  worst_day=2026-02-02 -24.9%  skew=-0.10  kurtosis=+1.82
SLV   vol=36.5%  VaR_5= -3.2%  ES_5= -5.2%  max_drawdown=-45.4%  worst_day=2026-01-30 -28.5%  skew=-1.73  kurtosis=+22.24
USO   vol=36.3%  VaR_5= -3.5%  ES_5= -5.2%  max_drawdown=-36.2%  worst_day=2022-03-09 -11.7%  skew=-0.17  kurtosis=+3.07
DBC   vol=19.2%  VaR_5= -1.9%  ES_5= -2.9%  max_drawdown=-27.4%  worst_day=2022-03-09  -7.9%  skew=-0.59  kurtosis=+3.51
DBB   vol=20.2%  VaR_5= -2.0%  ES_5= -2.8%  max_drawdown=-35.0%  worst_day=2022-03-09  -8.4%  skew=-0.18  kurtosis=+2.43
GLD   vol=18.2%  VaR_5= -1.7%  ES_5= -2.7%  max_drawdown=-24.5%  worst_day=2026-01-30 -10.3%  skew=-0.72  kurtosis=+7.59
DBA   vol=14.0%  VaR_5= -1.4%  ES_5= -2.0%  max_drawdown=-15.9%  worst_day=2024-05-13  -4.9%  skew=-0.23  kurtosis=+1.59

Highest tail risk: UNG with ES_5= -8.6%
Lowest tail risk:  DBA with ES_5= -2.0%

What this tells us

UNG has the highest tail risk by a large margin. Its annualized volatility is 64.1%, expected shortfall is -8.6%, and maximum drawdown is -92.5%. Natural gas exposure is not simply volatile; it has repeated severe downside events.

SLV and USO have similar expected shortfall, both around -5.2%, but the shape of the risk differs. SLV shows much higher kurtosis at +22.24, meaning its tail risk is more concentrated in extreme outliers. USO has a lower kurtosis reading but still carries large daily downside moves.

DBA is the lowest-risk instrument in this sample. Its expected shortfall is -2.0%, volatility is 14.0%, and maximum drawdown is -15.9%. Broad commodity exposure through DBC is safer than energy or silver, but it still has a -27.4% maximum drawdown.

So what?

Commodity allocation should not be sized from headline exposure alone. A 10% position in UNG does not carry the same risk as a 10% position in DBA or GLD. Expected shortfall makes that difference explicit.

For portfolio construction, the useful rule is to size commodity positions by tail loss rather than by capital dollars. Instruments with high expected shortfall need smaller weights, tighter drawdown limits, or explicit hedges. The worst days, not the average days, determine whether the commodity sleeve survives stress.

Built with xfinlink — free financial data API for Python. pip install xfinlink

Built with xfinlink — free financial data API for Python. pip install xfinlink
← All articles