Monte Carlo Simulation–Based Portfolio Sustainability Analysis
Comprehensive methodology and detailed explanation of Monte Carlo simulation for portfolio analysis
1. Introduction
Long-term portfolio planning under uncertainty is a central problem in financial economics, particularly in the context of retirement planning, endowment management, and sustainable withdrawal strategies. Deterministic forecasting techniques are insufficient to capture the stochastic nature of financial markets, inflation dynamics, and sequence-of-returns risk. Consequently, probabilistic modeling approaches—most notably Monte Carlo simulation—are widely employed to assess the distribution of potential future portfolio outcomes.
This chapter presents a comprehensive Monte Carlo simulation–based analysis of portfolio performance and sustainability under systematic withdrawals. The methodology, assumptions, metrics, and interpretations are documented in full detail to ensure transparency, reproducibility, and academic rigor.
2. Monte Carlo Simulation Framework
2.1 Conceptual Overview
A Monte Carlo simulation is a stochastic numerical technique that estimates the probability distribution of outcomes by repeatedly sampling from random variables governed by specified statistical properties. In the context of portfolio analysis, each simulation represents a plausible future realization of asset returns, inflation, and portfolio cashflows.
Rather than producing a single forecasted outcome, the simulation generates a large ensemble of potential return paths, enabling probabilistic inference regarding portfolio growth, risk, and failure rates.
3. Return Generation Model
3.1 Historical Returns–Based Simulation
The simulation employs a historical bootstrap model to generate future asset returns. Annual portfolio returns are randomly sampled (with replacement) from observed historical returns spanning January 1972 to December 2024.
Rationale
- Preserves empirical distribution characteristics (fat tails, skewness)
- Retains real-world extreme events
- Avoids parametric assumptions about return normality
Assumption
Future market behavior is statistically consistent with historical observations.
4. Portfolio Configuration
4.1 Initial Conditions
- Initial portfolio value: $1,000,000
- Investment horizon: 30 years
- Number of simulations: 10,000
4.2 Asset Allocation
The portfolio consists of a single asset class:
| Asset Class | Allocation |
|---|---|
| US Stock Market | 100% |
This allocation represents a high-risk, high-return strategy with no diversification benefits.
5. Inflation Modeling
5.1 Historical Inflation Process
Inflation is modeled using historical CPI-U data for the United States over the same period as asset returns.
- Mean inflation rate: 3.93%
- Standard deviation: 1.31%
Inflation realizations are sampled concurrently with asset returns and incorporate historical correlations between equity returns and inflation.
6. Withdrawal Policy
6.1 Fixed Real Withdrawal Strategy
The portfolio is subjected to a fixed annual withdrawal of $45,000, adjusted for inflation and withdrawn at the end of each year.
This withdrawal policy is representative of retirement income strategies and introduces sequence-of-returns risk, whereby early adverse returns disproportionately increase the probability of portfolio failure.
7. Rebalancing Assumptions
The portfolio is rebalanced annually to maintain the target asset allocation. Rebalancing constrains portfolio risk drift and ensures consistency with the specified investment policy.
8. Definition of Portfolio Survival
A simulation is classified as successful if the portfolio balance remains strictly positive after all scheduled withdrawals. The probability of success is defined as the proportion of simulations satisfying this condition.
In the present analysis, 84.10% of simulated portfolios survived the full 30-year horizon.
9. Performance Metrics
9.1 Return Measures
9.1.1 Time-Weighted Rate of Return (TWRR)
TWRR measures the compound growth rate of the portfolio independent of external cashflows and reflects pure investment performance.
9.1.2 Nominal and Real Returns
Nominal returns are measured in unadjusted monetary units, whereas real returns are adjusted for inflation to reflect changes in purchasing power.
9.2 Risk Measures
9.2.1 Volatility
Volatility is computed as the annualized standard deviation of monthly returns and serves as a primary measure of return uncertainty.
9.2.2 Maximum Drawdown
Maximum drawdown quantifies the largest observed peak-to-trough decline in portfolio value. Metrics are computed both including and excluding the effect of withdrawals.
9.3 Risk-Adjusted Performance
9.3.1 Sharpe Ratio
The Sharpe ratio measures excess return per unit of total risk relative to a risk-free benchmark.
9.3.2 Sortino Ratio
The Sortino ratio refines the Sharpe ratio by penalizing only downside volatility, making it particularly suitable for asymmetric return distributions.
10. Withdrawal Sustainability Metrics
10.1 Safe Withdrawal Rate (SWR)
The safe withdrawal rate is defined as the maximum inflation-adjusted withdrawal (expressed as a percentage of the initial portfolio value) that avoids portfolio depletion over the simulation horizon.
10.2 Perpetual Withdrawal Rate (PWR)
The perpetual withdrawal rate represents the maximum annual withdrawal percentage that preserves the inflation-adjusted portfolio principal indefinitely.
11. Distributional Results and Probabilities
11.1 Percentile Analysis
Simulation outcomes are reported across percentiles (10th, 25th, 50th, 75th, and 90th), enabling interpretation of worst-case, median, and best-case scenarios.
11.2 Return and Loss Probabilities
The probability tables report the likelihood of achieving returns above specified thresholds and experiencing losses of given magnitudes, both with and without cashflows.
12. Interpretation of Results
The results indicate substantial dispersion in long-term outcomes, highlighting the inherent uncertainty of equity-dominated portfolios under systematic withdrawals. While median outcomes suggest strong real wealth accumulation, downside scenarios reveal nontrivial failure risk driven primarily by early negative return sequences.
13. Limitations
- Results are hypothetical and not predictive
- Structural market changes are not modeled
- Taxes, transaction costs, and behavioral responses are excluded
- Assumes stable statistical relationships over time
14. Conclusion
Monte Carlo simulation provides a robust framework for evaluating portfolio sustainability under uncertainty. The analysis demonstrates the trade-off between growth potential and withdrawal risk, emphasizing the importance of diversification, withdrawal discipline, and inflation awareness in long-term financial planning.
Mathematical Appendices
Appendix A: Return Compounding
A.1 Portfolio Value Evolution
Vt+1 = (Vt - Wt) × (1 + Rt)where Vt = portfolio value at time t, Wt = withdrawal at time t, Rt = portfolio return at time t
Appendix B: Inflation Adjustment
Wt = W0 × ∏i=1t (1 + Ii)where Ii is the inflation rate in year i.
Appendix C: Time-Weighted Rate of Return
TWRR = (∏i=1n (1 + ri))1/n - 1Appendix D: Real Returns
rreal = (1 + rnominal) / (1 + i) - 1Appendix E: Volatility
σannual = σmonthly × √12Appendix F: Maximum Drawdown
Max Drawdown = maxt ((Pt - Vt) / Pt)where Pt is the historical peak value.
Appendix G: Sharpe Ratio
Sharpe = (E[Rp] - Rf) / σpAppendix H: Sortino Ratio
Sortino = (E[Rp] - Rf) / σdownsideAppendix I: Safe Withdrawal Rate
SWR = W / V0Appendix J: Monte Carlo Estimation of Success Probability
P(Success) = (1/N) × Σi=1N 1(Vi,T > 0)