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 ClassAllocation
US Stock Market100%

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 - 1

Appendix D: Real Returns

rreal = (1 + rnominal) / (1 + i) - 1

Appendix E: Volatility

σannual = σmonthly × √12

Appendix F: Maximum Drawdown

Max Drawdown = maxt ((Pt - Vt) / Pt)

where Pt is the historical peak value.

Appendix G: Sharpe Ratio

Sharpe = (E[Rp] - Rf) / σp

Appendix H: Sortino Ratio

Sortino = (E[Rp] - Rf) / σdownside

Appendix I: Safe Withdrawal Rate

SWR = W / V0

Appendix J: Monte Carlo Estimation of Success Probability

P(Success) = (1/N) × Σi=1N 1(Vi,T > 0)

Implementation Verification

A detailed comparison of the implementation with this documentation is available, including verification of formulas, methodology alignment, and any differences.