Monte Carlo Risk Simulation — Pagmomodelo ng Mga Resulta ng Portfolio
Ang Monte Carlo simulations ay nagpapahintulot sa iyo na mag-model ng 100,000+ potensyal na mga landas sa hinaharap ng iyong portfolio, na nagbibigay-katuturan sa mga sukat ng panganib tulad ng Value at Risk (VaR) at inaasahang pagkukulang. Sa halip na mag-assume ng isang resulta lamang, makikita mo ang distribusyon ng mga posibilidad.
Halimbawa: Magpatakbo ng 100k simulations ng iyong trading strategy. Makikita mo: 68% ng mga resulta ay kumikita ng $10-$50, 25% ay lugi ng $10-$30, 5% ay lugi ng $50-$200, at 2% ay kumikita ng $100+. Ngayon ay naiintindihan mo na ang iyong risk profile.
Basic Geometric Brownian Motion (GBM)
Modelo ang mga presyo sa hinaharap bilang mga random walk na may drift at volatility:
import numpy as np
def simulate_price_paths(S0, mu, sigma, T, dt, n_sims=10000):
steps = int(T / dt)
paths = np.zeros((n_sims, steps))
paths[:, 0] = S0
for i in range(1, steps):
Z = np.random.normal(0, 1, n_sims)
dS = mu * paths[:, i-1] * dt + sigma * paths[:, i-1] * np.sqrt(dt) * Z
paths[:, i] = paths[:, i-1] + dS
return paths
Portfolio Risk Metrics
Value at Risk (VaR)
Ang pagkawala na lalampas mo lamang sa 5% ng mga senaryo (95% confidence):
def calculate_var(simulations, confidence=0.95):
pnl = simulations[:, -1] - simulations[:, 0]
var = np.percentile(pnl, (1 - confidence) * 100)
return var
Expected Shortfall (CVaR)
Average na pagkawala sa pinakamasamang 5% ng mga kaso:
def calculate_cvar(simulations, confidence=0.95):
pnl = simulations[:, -1] - simulations[:, 0]
var_threshold = np.percentile(pnl, (1 - confidence) * 100)
worst_5_pct = pnl[pnl <= var_threshold]
cvar = np.mean(worst_5_pct)
return cvar
Strategy-Specific Simulation
Simulate ang iyong aktwal na trading rules, hindi lamang mga price paths:
def simulate_strategy(capital, win_rate, avg_win, avg_loss, trades_per_day=10, days=30, n_sims=10000):
results = np.zeros(n_sims)
for sim in range(n_sims):
equity = capital
n_trades = trades_per_day * days
for _ in range(n_trades):
if np.random.rand() < win_rate:
pnl = np.random.normal(avg_win, avg_win * 0.2)
else:
pnl = -np.random.normal(avg_loss, avg_loss * 0.2)
equity += pnl
if equity < capital * 0.5:
break
results[sim] = equity
return results
Stress Testing with Smart Money Signals
Modelo ang pinakamasamang mga senaryo kapag bumagsak ang Smart Money confidence:
def stress_test_scenarios(base_win_rate):
scenario_normal = simulate_strategy(10000, win_rate=base_win_rate)
scenario_medium = simulate_strategy(10000, win_rate=base_win_rate * 0.95)
scenario_veto = simulate_strategy(10000, win_rate=0.45)
return {
'normal': {'var': np.percentile(scenario_normal, 5), 'mean': np.mean(scenario_normal)},
'medium': {'var': np.percentile(scenario_medium, 5), 'mean': np.mean(scenario_medium)},
'veto': {'var': np.percentile(scenario_veto, 5), 'mean': np.mean(scenario_veto)}
}
Interpretation
Mula sa 10k simulations, matututunan mo:
- Median outcome (50th percentile)
- Best case (95th percentile)
- Worst case (5th percentile)
- Probability of ruin (equity < 0)
- Expected drawdown length
Stress-test strategies with Smart Money signals
Unawain ang portfolio risk sa ilalim ng normal, deteriorating, at VETO Smart Money conditions. Ang aming API ay nagbibigay sa iyo ng mga signal; ang Monte Carlo ay nagmo-modelo ng mga resulta.
Risk Model Today →