BLOG

Behind the numbers.

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

How to Get Historical Market Cap Data in Python
How Much of the S&P 500 Survives 20 Years? Index Turnover Analysis in Python
Does Joining the S&P 500 Bring New Institutional Owners? 13F Event Study in Python
Is Volatility Seasonal? Calendar Month Analysis of Realized Volatility in Python
Does Fast Revenue Growth Force Companies to Borrow? Cash Funding Analysis in Python
How Far Back Does SEC EDGAR Data Go?
Are One-Time Charges Really One-Time? Charge Frequency Analysis in Python
Does Buying the Dip Work? Short-Term Reversal by Volatility Regime in Python
Alpha Vantage vs Massive vs xfinlink for Fundamentals
How Long Does a Stock Take to Recover From a 50% Fall? Drawdown Analysis in Python
Do Companies That Shrink Their Share Count Outperform? Net Buyback Yield in Python
How Much Does the Dow's Price Weighting Distort It? Index Weighting Analysis in Python
How to Get SEC Form 4 Insider Trading Data in Python
Can Anything Predict Next Month's Stock Returns? Out-of-Sample R-Squared Testing in Python
How Much of a Stock's Return Comes From Its Sector? Variance Decomposition in Python
Altman Z-Score: Where To Get It in Python
Do Value Screens Agree on Which Stocks Are Cheap? Multiple Overlap Analysis in Python
Annual vs Quarterly Financial Data: Which to Use
Do High Returns on Capital Persist? ROIC Fade Analysis in Python
Does Past Beta Predict Future Beta? Beta Stability Testing in Python
Do Defensive Sectors Actually Defend? Up and Down Capture in Python
How to Choose a Financial Data API
Do Small Caps Actually Beat Large Caps? Size Premium Test in Python
What If You Miss the Market's Best Days? Extreme-Day Analysis in Python
Does Rebalancing Add Return? Fixed-Weight vs Drift Portfolios in Python
Which S&P 500 Companies Are Closest to Default? Merton Distance-to-Default in Python
What Happens to Stocks Removed From the S&P 500? Replacement Pair Analysis in Python
Does the Golden Cross Work? 50/200 Moving Average Crossover Backtest in Python
Financial Data for Academic Finance Research
Does Skipping the Most Recent Month Improve Momentum? S&P 500 Decile Sorts in Python
Does Cointegration Survive Out of Sample? Pairs Trading Validation in Python
Which Dividends Are Not Covered by Cash? Free-Cash-Flow Coverage Screening in Python
GICS vs SIC vs NAICS: Which Industry Classification to Use
How Many Stocks Does It Take to Diversify? Random Portfolio Simulation in Python
How Many Days of Data Does a Volatility Estimate Need? Range-Based Estimators in Python
How Much of S&P 500 Cash Flow Is Stock Compensation? Cross-Sectional Analysis in Python
What Is a 13F Filing? Institutional Holdings Explained
Does Revenue Growth Explain Profit Growth? Cross-Sectional Decomposition in Python
How Much of the Nasdaq 100 Is Already in the S&P 500? Index Overlap Analysis in Python
How Often Does a 99% Value-at-Risk Limit Actually Break? VaR Backtesting in Python
Real-Time vs End-of-Day Market Data: Which Do You Need?
How Concentrated Are S&P 500 Earnings? Point-in-Time Index Analysis in Python
Does Volatility Scale With the Square Root of Time? Variance Ratio Test in Python
Does Goodwill Distort the Price-to-Book Screen? Goodwill-Adjusted Valuation in Python
How Are Shares Outstanding Reported (and Why They Disagree)
Do Low-Volatility Stocks Deliver Better Risk-Adjusted Returns? S&P 500 Quintile Sorts in Python
Does Trend Following Beat Buy and Hold? Time-Series Momentum in Python
Has the Stock-Bond Correlation Flipped? 60/40 Portfolio Risk in Python
What API to Use for a Stock Screener
Which Assets Hedge Inflation Shocks? Macro Factor Betas in Python
Does Covariance Shrinkage Beat the Sample Covariance? Minimum-Variance Portfolios in Python
Do Faster Inventory Turns Mean Thinner Margins? Gross Margin Return on Inventory in Python
SEC EDGAR API vs Fundamentals API: Which to Use
Does Fast Earnings Growth Persist? Rank Correlation Analysis in Python
Split Adjustment Explained: Adjusted Close vs Close
Which Trading Day of the Month Pays Best? Turn-of-the-Month Analysis in Python
Can You Use Yahoo Finance Data Commercially?
How Many Independent Bets Does a Nine-Sector Portfolio Give You? Eigenvalue Analysis in Python
Which Volatility Forecast Wins One Month Ahead? HAR vs EWMA in Python
How Concentrated Are Institutional Equity Portfolios? Form 13F Concentration Analysis in Python
Do Stocks Earn Their Returns Overnight or Intraday? Return Decomposition in Python
When Do Corporate Insiders Actually Trade? Form 4 Timing Analysis in Python
Data Requirements for Backtesting a Trading Strategy
What Is Survivorship Bias in Backtesting?
Do High Dividend Yields Come From Bigger Payouts or Falling Prices? Yield Decomposition in Python
Do Stocks Fall Harder Than They Rise? Downside Beta vs Upside Beta in Python
Does Volatility Targeting Improve Sharpe Ratios? Seven-Asset Backtest in Python
Free Stock Market Data APIs: What You Actually Get
How to Give an LLM Financial Data With an MCP Server
Does Post-Earnings Announcement Drift Survive Real Filing Dates? PEAD Event Study in Python
Do Insider Buying Clusters Predict Returns? Signal Testing in Python
Does the S&P 500 Index Effect Still Exist? Event Study in Python
Are Companies Leaving the S&P 500 Faster Than They Used To? Index Survival Analysis in Python
Does Gross Profitability Predict Stock Returns? Quintile Factor Test in Python
What Growth Rate Is the Market Pricing In? Reverse DCF in Python
Does a Strong Balance Sheet Cushion Drawdowns? Leverage and Downside Risk in Python
Does Ticker Recycling Corrupt a Mean-Reversion Backtest? Entity-Resolved Z-Scores in Python
How Much of a Growth Screen's Backtested Edge Is Survivorship Bias? Point-in-Time Index Testing in Python
Why Do Leveraged ETFs Decay? Measuring Volatility Drag in Python
Are Consumer Staples Margins Shrinking Under Inflation? Gross Margin Trend Analysis in Python
Do Weak Jobs Reports Predict Market Drawdowns? NFP Surprise Event Study in Python
Is the Rotation From Tech to Industrials Backed by Earnings? Relative EPS Growth Analysis in Python
Is the Semiconductor Rally Broadening Beyond NVIDIA? Return Dispersion Analysis in Python
Which Stocks Benefit Most When Oil Prices Fall? Oil Beta Screening in Python
Do Bond Returns Predict Stock Returns? Granger Causality Test in Python
Which Stocks Actually Drive Portfolio Returns? Shapley Value Attribution in Python
Does "Sell in May" Still Work? Calendar Anomaly Backtest in Python
How to Build Complete Price History Through Ticker Changes? Entity Resolution in Python
Are KO and PEP Cointegrated? Pairs Trading Signal Construction in Python
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 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

Do Bond Returns Predict Stock Returns? Granger Causality Test in Python

What's the question?

A longstanding question in macro finance is whether bond market movements contain information about future equity returns. The logic is intuitive: bonds reflect expectations about interest rates, inflation, and credit conditions. If those expectations shift before equity markets react, then lagged bond returns would carry predictive power for stock returns.

The Granger causality test provides a formal statistical framework for this question. Developed by Clive Granger in 1969, the test asks whether the past values of one time series improve the forecast of another, beyond what the second series' own history can explain. Importantly, "Granger causality" does not imply true causation. It measures temporal precedence and incremental predictive content. If lagged bond returns significantly reduce the forecast error for stock returns, bonds are said to Granger-cause stocks.

The test works by estimating two models. The restricted model regresses stock returns on their own lags only. The unrestricted model adds lagged bond returns as additional regressors. An F-test then compares the residual sum of squares. A low p-value (below 0.05) indicates that the bond return lags contribute statistically significant explanatory power.

The approach

The test uses monthly returns for five ETFs spanning the fixed-income spectrum: TLT (20+ year Treasury bonds), IEF (7-10 year Treasury bonds), SHY (1-3 year Treasury bonds), and AGG (broad aggregate bond index), with SPY representing the equity market.

  1. Pull 7 years of monthly closing prices for all five ETFs
  2. Compute monthly returns from closing prices
  3. Run the Granger causality test in both directions: bond returns predicting stock returns, and stock returns predicting bond returns
  4. Test lag orders 1 through 4 (1 to 4 months of history) to check whether predictive content emerges at longer horizons

Duration diversity matters. If the yield curve transmits information to equity markets, the effect may appear at different lag structures for long-duration bonds (TLT) versus short-duration instruments (SHY).

Code

import xfinlink as xfl
import pandas as pd
from statsmodels.tsa.stattools import grangercausalitytests

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

tickers = ["SPY", "TLT", "IEF", "SHY", "AGG"]
df = xfl.prices(tickers, period="7y", interval="1mo", fields=["close"])

returns = {}
for t in tickers:
    sub = df[df["ticker"] == t].sort_values("date").copy()
    sub["return"] = sub["close"].pct_change()
    sub = sub.dropna(subset=["return"])
    returns[t] = sub[["date", "return"]].set_index("date")
    returns[t].columns = [t]

merged = returns["SPY"]
for t in tickers[1:]:
    merged = merged.join(returns[t], how="inner")

bond_etfs = ["TLT", "IEF", "SHY", "AGG"]
for bond in bond_etfs:
    test_data = merged[["SPY", bond]].copy()
    gc = grangercausalitytests(test_data, maxlag=4, verbose=False)
    for lag in range(1, 5):
        f_stat = gc[lag][0]["ssr_ftest"][0]
        p_val = gc[lag][0]["ssr_ftest"][1]
        sig = "YES" if p_val < 0.05 else "no"
        print(f"{bond:<6} lag={lag}  F={f_stat:.3f}  p={p_val:.4f}  {sig}")

Full script with formatting and visualisation: bond-returns-granger-causality-stock-python.py

Output

Bar chart showing Granger causality p-values for bond ETFs predicting stock returns at lags 1 through 4, with a 5 percent significance threshold line
Observation period: 2019-08 to 2026-06
Monthly observations: 83

Bond ETF   Lag       F-stat    p-value   Significant?
-------------------------------------------------------
TLT        1          2.416     0.1241             no
TLT        2          1.898     0.1569             no
TLT        3          1.539     0.2117             no
TLT        4          1.136     0.3467             no
IEF        1          2.364     0.1282             no
IEF        2          1.846     0.1649             no
IEF        3          2.466     0.0689             no
IEF        4          1.708     0.1579             no
SHY        1          0.974     0.3268             no
SHY        2          1.455     0.2399             no
SHY        3          4.564     0.0055            YES
SHY        4          3.307     0.0154            YES
AGG        1          3.198     0.0775             no
AGG        2          2.788     0.0679             no
AGG        3          2.586     0.0596             no
AGG        4          1.837     0.1314             no

Reverse test: Do stock returns Granger-cause bond returns?
Bond ETF   Lag       F-stat    p-value   Significant?
-------------------------------------------------------
TLT        1          0.072     0.7898             no
TLT        2          0.129     0.8791             no
TLT        3          0.188     0.9043             no
TLT        4          0.151     0.9618             no
IEF        1          0.458     0.5004             no
IEF        2          0.231     0.7942             no
IEF        3          0.242     0.8666             no
IEF        4          0.170     0.9528             no
SHY        1          1.778     0.1863             no
SHY        2          0.970     0.3836             no
SHY        3          0.541     0.6559             no
SHY        4          0.379     0.8228             no
AGG        1          0.056     0.8143             no
AGG        2          0.165     0.8482             no
AGG        3          0.359     0.7827             no
AGG        4          0.295     0.8806             no

What this tells us

Long-duration Treasuries (TLT, IEF) and the broad bond aggregate (AGG) do not Granger-cause equity returns at any lag tested. P-values range from 0.07 to 0.35 — none crossing the 0.05 threshold. The bond market's most volatile and widely watched segment, the long end of the yield curve, does not contain statistically significant predictive information for monthly stock returns over this sample.

The exception is SHY, the 1-3 year Treasury ETF, which achieves significance at lag 3 (p = 0.006) and lag 4 (p = 0.015). Short-term Treasuries are the segment of the curve most directly influenced by Federal Reserve policy expectations. A shift in SHY returns reflects changing expectations for near-term fed funds rate decisions. The 3-4 month lag aligns with the typical horizon over which monetary policy expectations crystallize and propagate into risk asset repricing.

The reverse test is uniformly insignificant. Equity returns do not Granger-cause bond returns at any lag or for any bond ETF. The relationship, to the extent it exists, is one-directional: short-term rate expectations may lead equities, but equities do not lead bonds.

AGG approaches significance at lags 1 and 2 (p-values of 0.08 and 0.07) without crossing the threshold. AGG blends duration exposures, diluting whatever signal SHY carries at longer lags while picking up partial sensitivity to rate expectations through its short-duration allocation.

So what?

For portfolio managers and systematic traders, the results suggest that monitoring short-term Treasury returns — not the long bond — may offer a marginal timing edge for equity allocation decisions at the quarterly horizon. SHY captures the market's real-time pricing of near-term monetary policy, and that pricing appears to lead stock market moves by 3 to 4 months.

This does not mean building a trading strategy on the SHY signal alone. With only 83 observations and two significant lags out of 16 tests, the finding is suggestive rather than definitive, and multiple testing inflates the chance of a false positive. The practical application is as a supplementary input: if short-term bond returns have been consistently positive (rates falling, rate cuts expected) over the past quarter, it is a tailwind for equities. If short-term bonds have sold off (rates rising, tightening expected), the headwind may take another quarter to arrive in stock prices.

The absence of long-bond predictive power is itself informative. Duration risk in TLT and IEF is driven by term premium and inflation expectations — factors that the equity market prices concurrently rather than with a lag. The short end, by contrast, carries cleaner monetary policy information that takes time to filter through credit conditions, corporate earnings expectations, and investor positioning.

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

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