[옵션] 11번의 과거 시장폭락 사태를 분석해보자

스샷 첨부하고 싶었는데, 캡처도 못하게 틀어 막아놨다. 계좌정보가 보이는 화면도 아니고 호가창인데 대체 왜? 애초에 네이키드 매도는 기대조차도 하지도 않았다. 들어가보니 역시나 콜/풋 스프레드 매매나 아이언 콘도르, 버터플라이 같은 손실이 제한적인 멀티레그 전략들만 도입해놨다. 좋아, 그거라도 해준게 어디냐? 거기까진 이해하겠다 근데, 스프레드 간격을 지들 맘대로 4달러가 최대로 막아놨음 만약 스프레드 간격 자유롭게 설정하도록 해놨다가 손실구간 넓으면 자기들 돈 떼일까봐 걱정하는 건 알겠는데 스프레드 간격 4달러로 뭘 어떤 전략을 짜라는 거냐 도대체? 아니 어차피 다 손실 제한되는 스프레드나 박스 매매들 뿐인데 그럼 증거금을 최대 손실금액만큼 납입하게 하고 스프레드 간격은 내가 원하는 만큼 알아서 정하게 하면 되잖아? 그걸 왜 니네가 정하는데 ㅋㅋㅋㅋ 어차피 돈 떼이는게 걱정되는 거면, 스프레드 넓게 잡게 해줬다가 떼이나, 계약수 무지막지하게 올렸다가 떼이나 똑같은데 왜 니들맘대로 스프레드 간격을 제한을 두는거냐고 ㅋㅋㅋㅋ HTS 지원도 안됨 핸드폰 MTS로 딸깍딸깍 대충 빨리 매매하고 빨리 수수료나 많이 내라 이거 같은데? 가격 차트도 안보여줌 살다살다 매매하는데 가격정보도 안보여주면 뭘 기준으로 어떻게 매매하라는 거냐 도대체? 더러워서 그냥 안하고 만다. 지들 맘대로 팔다리 다짤라놓고 우리가 정해주는 대로만 매매하라는 심보가 어이가 없어서 앱 삭제해버렸다. 해외옵션 불모지에 좀 변화가 생기나 했더니, 차트도 그릭스도 제대로 안보여주고 HTS 지원도 안되고 핸드폰으로만 대충 깜깜이로 회전율만 높여서 개인투자자들 수수료 뽑으려고 구조 짜는게 너무 대놓고라 욕지거리가 올라온다. 유진투자선물 계좌가 이미 있는데, 멀티레그 전용 계좌 따로 개설하라고 할 때부터 알아봤다. 에혀

휴대폰으로 문자가 왔는데, 스캠인지 진짜인지 헷갈려서 링크 클릭을 못하겠다. UX/UI는 진짜 구려도 개별주식 옵션이 그나마 (거의) 전종목 지원돼서 그나마 쓰고 있는 증권사이긴 한데... 개별종목 가치평가하고 나서도 적정 주가 올 때까지 Cash secured put 못해서 한이었는데... 이 한을 이제 푸나요?? 근데 다시 읽어보니 네이키드 매도는 또 안되고 합성포지션만 된다는 듯. 아마 네이키드 매도 포지션 열어줬다가 갑작스런 이벤트로 개인이 드러누우면 돈 떼일 우려가 있어서 손익이 제한되는 합성만 열어주는 듯 12가지 전략이 뭔지는 아직 모르겠지만, 같은 맥락으로 아마 스트래들/스트랭글 매도같은 상하방 손실이 무한히 열려있는 건 제외일 듯 그래도 해외 옵션불모지 한국에서 그게 어디냐 싶기도하다. 그와중에 또 신규계좌 별도개설 필요... 제발 통합좀 해라. 주식계좌, 선물 옵션 계좌, 거기에 멀티레그 계좌까지 세개를 돌리라고?

0. 들어가며 먼저, CBOE에서 발간한 따끈따끈한 보고서를 소개해주신 Valley의 옵션장인 asymmetry님께 감사드립니다. Cboe - VIX Index Decomposition 내용이 궁금해서, 제가 보려는 목적으로 번역본도 올려보았는데, 조금 복잡하긴 하지만 너무 흥미로워서 제가 이해한 것이 맞는지도 의견을 나눠볼 겸, 나름대로의 해석을 덧붙여 볼까 합니다. 저는 현업도, 옵션을 본격적으로 매매하는 전문 트레이더도 아니고, 그저 옵션을 공부하면서 정리하기 위해 쓴 글이니, 제가 틀리게 알고 있는 부분이 있으면 가감없는 댓글 부탁드립니다. 번역본을 먼저 한번 훑어보시고, 대략 이해가 가시면 이하 내용은 동어반복이라 굳이 읽어보실 필요가 없을 것 같습니다. 1. VIX 상승의 구성요소 분해 우리는 VIX를 단순히 공포지수 혹은 내재변동성 지수라고만 알고 있고, S&P500이 하락하면 VIX는 상승한다 정도로만 알고 있습니다. 조금 더 나아간다면, S&P500 옵션시장 참여자들이 생각하는 향후 30일간의 S&P500의 기대변동성의 연율화된 값이라고 설명할 수 있겠습니다. 그리고 실제로 S&P500과 VIX의 상관관계는 거의 -1에 육박한다고 알려져 있습니다. 따라서, 단순히 "시장이 떨어져서 VIX가 올랐다"라는 해석 정도로는 시장에 대한 특별한 인사이트를 읽을 수가 없습니다. 저 역시도 VIX에 대한 이해도가 부족해서 딱 그 정도 레벨에 머물러 있었고, 기껏해야 SPX와 VIX의 다이버전스를 통해 "추세 반전이 임박했다" 정도로만 해석하는 데에 그쳤을 뿐 그 이상의 인사이트를 뽑아내기에는 어려움이 있었습니다. "당신은 몇차원 시장참여자인가?"라는 질문과 파생상품 시장 그런데 이번에 CBOE에서 8월 1일자로 발간한 VIX Index Decomposition, 즉 VIX 지수의 분해라는 보고서를 통해 그 실마리를 조금 잡을 것도 같다는 생각을 하게 되었습니다. 보고서의 요지는 VIX지수가 상승했을 때, 그 상승의 기여도를 여러가지 하위요소로 분해해서 바라보면 단순히 VIX 수치 이상의 인사이트를 얻을 수 있다. 라는 것입니다. 보고서에서 CBOE는 VIX 변화의 구성요소를 크게 아래와 같이 분해합니다. 행사가 변화에 의한 VIX 변화량 (ATM 구간) 스큐 곡선의 평행이동에 의한 VIX 변화량 시장의 심리에 의한 변화량 (지수와 VIX의 동반 움직임) 외가 콜옵션 및 풋옵션 기울기 변화에 의한 VIX 변화량 (15Δ ~ 45Δ 구간) 극외가 콜옵션 및 풋옵션 곡률 변화에 의한 VIX 변화량 (1Δ ~ 15Δ 구간) 즉, 내재변동성 커브(IV Curve)를 행사 가능성인 델타를 기준으로 나누어서 해석해보는 것입니다. 하나씩 살펴보겠습니다. 1-1. 행사가 변화에 의한 VIX 변화량 (ATM 구간) VIX 변화의 첫번째 구성요소는 행사가 자체가 변화함에 따라 내재변동성(이하 IV) 역시 정해진 IV 커브(IV Curve) 위에서의 움직임이는 데에서 오는 변화입니다. 아래 차트를 보면, X축은 옵션의 행사가이고 Y축은 해당 행사가에서의 내재변동성(%)를 나타냅니다. 연한색은 앤캐리 청산 사태 폭락 전인 8월 2일의 IV커브, 진한색은 폭락 후인 8월 5일의 IV커브입니다. 여기서 곡선 자체가 위로 이동한것은 제외하고, 연한 곡선위에서 8월 2일자에서 8월 5일자로 이동한 것을 확인해보겠습니다.(빨간색 화살표) X축이 옵션의 행사가라고 했는데, 폭락전인 8월 2일의 등가격은 5,346이었습니다. 그리고 8월 5일 폭락 후 등가격은 5,186이었습니다. 그럼 순수하게 S&P500이 5,346에서 5,186까지 떨어지면서, 즉 X축이 움직이면서 내재변동성이 상승한 변화량을 측정하려면, 이동전의 Y값과 이동후의 Y값의 차이를 구하면 됩니다. 따라서, 연한색 곡선위에서 X축 행사가 5,346에 해당하는 Y값(IV)는 18.99 연한색 곡선위에서 X축 행사가 5,186에 해당하는 Y값(IV)는 21.55 즉, IV곡선의 모양과 위치가 아무것도 변하지 않았다고 가정했을 때 순수하게 행사가의 변화 때문에 발생한 IV의 변화량은 21.55-18.99 = 2.57이라는 의미입니다. 이것이 고정 행사가(Sticky strike)당 예상 움직임이라고 보고서에서는 표현하고 있습니다. IV 커브를 고정시킨 상태에서 순수 행사가의 움직임으로부터 기인한 IV변화인이 고정행사가 당 예상 움직임은 변동성이 낮은 체제(Regime)에서는 꽤나 예측가능성 있게 움직입니다. 따라서 IV 커브가 상하좌우로 요동치지 않는 안정된 변동성 체제하에서는 SPX가 하락하는 만큼 비례해서 정직하게 움직이기 때문에 예측가능한 VIX 체제를 만듭니다. 아래 표를 봐도, ① VIX레벨이 낮을수록, ② VIX 움직임의 폭이 작을 수록 고정행사가당 예상움직임의 설명력이 높은 것을 확인할 수 있습니다. 1-2. ...

1. VIX [S&P500과 VIX(파란색, 상하반전)] 주식시장은 전고점을 넘겨 상승했지만, 옵션시장(VIX, 상하반전)은 전고점을 회복하지 못함. 즉, 상하반전된 차트이므로 주가가 상승하는 것에 비해 내재변동성은 더이상 하락하지 못하고 바닥이 Sticky한 상태라고 해석할 수 있음 7월 9일 상호관세 유예 마감일을 의식한 것일까? 이런 경우는 당연한 소리일 수도 있지만 둘 중 한가지로 귀결됨. 옵션시장 참여자들이 틀렸다는 것이 밝혀져서 불확실성 해소에 따라 VIX가 급격히 하락하여 주가와 상관관계를 다시 회복하거나 주식시장 참여자들이 틀렸다는 것이 밝혀져서 VIX의 방향으로 주가가 수렴하여 하락하거나 2. VVIX [S&P500과 VIX(파란색, 상하반전), ...

Asymmetry님의 [Repository] #002 옵션내재확률 구하기글을 참조하여, 만기별 옵션 내재확률분포를 시각화하고, 인터랙티브 차트로 구성하는 코드를 작성해봄 만기를 세개까지 동시에 비교할 수 있도록 하고, 현물 가격을 세로선으로 표시하도록 수정 굳이 인터랙티브 차트를 추가 구성해 본 이유는 마우스를 갖다대면 내재확률(PDF)와 누적확률(CDF)가 나오기 때문에 만기에 특정 행사가 내에서 마감할 확률의 구간 추정에 용이하기 때문일 것 같아서임. spx_quotedata 받아서 어떻게든 혼자 해보려고 Cursor 고문하다가 포기하려던 무렵, Asymmetry님이 옛날에 올려주셨던 코드를 발견해서 마개조함. (감사합니다 비대칭좌...) 이제 이를 어떻게 해석해서 실전적인 마켓뷰로 반영할 수 있는지 연구해보자. 코드전문 import pandas as pd import numpy as np from scipy.interpolate import interp1d, UnivariateSpline from scipy.ndimage import gaussian_filter1d from scipy.stats import norm, skew, kurtosis import matplotlib.pyplot as plt from matplotlib.ticker import FuncFormatter from scipy.optimize import brentq from scipy.integrate import cumulative_trapezoid import datetime from datetime import datetime import plotly.graph_objects as go from plotly.subplots import make_subplots def is_third_friday(date): return date.weekday() == 4 and 15 <= date.day <= 21 # Time difference calculation: seconds -> minutes -> years def calculate_time_to_expiry(expiration_datetime): current_datetime = datetime.now() time_difference = (expiration_datetime - current_datetime).total_seconds() / 60 # Convert seconds to minutes time_in_years = time_difference / 525600 # Convert minutes to years return time_in_years # Modified function to handle expiration filtering and time delta calculation def process_expiration_data(option_data, target_dates): results = {} spotPrice = get_spot_price(file_path) for target_date in target_dates: print(f"\nProcessing data for expiration date: {target_date}") # Filter for expiration date and identify if it's the third Friday data_date = option_data[option_data['Expiration Date'] == target_date].copy() if is_third_friday(target_date): # Print a message confirming that we are excluding SPXW print(f"Third Friday identified for {target_date}, only taking SPX (excluding SPXW).") data_date = data_date[~data_date['Calls'].str.contains('SPXW')] # Exclude SPXW if third Friday else: print(f"Not the third Friday. Both SPX and SPXW are included for {target_date}.") # Apply expiration datetime and calculate time to expiry in years data_date['Expiration DateTime'] = data_date.apply(assign_expiration_datetime, axis=1) expiration_datetime = data_date['Expiration DateTime'].iloc[0] T = calculate_time_to_expiry(expiration_datetime) # Time to expiration in years print(f"Time to Expiration: {T} years") # Clean and filter valid data data_date = data_date.dropna(subset=['Bid', 'Ask', 'Bid.1', 'Ask.1']) data_date = data_date[(data_date['Bid'] > 0) & (data_date['Ask'] > 0) & (data_date['Bid.1'] > 0) & (data_date['Ask.1'] > 0)] # Continue with further steps such as calculating mid prices, implied volatility, etc. # ... # Store results for plotting and comparison later results[target_date] = { 'spot_price': spotPrice, 'time_to_expiry': T, # Other results you compute... } return results def assign_expiration_datetime(row): exp_date = row['Expiration Date'] if 'SPXW' in row['Calls']: # PM-settled SPXW (weekly options expire at 4:00 PM) return datetime(exp_date.year, exp_date.month, exp_date.day, 16, 0) else: # AM-settled SPX (standard options expire at 9:30 AM) return datetime(exp_date.year, exp_date.month, exp_date.day, 9, 30) def implied_volatility(option_price, S, K, T, r, option_type='call'): def bs_price(sigma): if sigma <= 0: return np.nan d1 = (np.log(S / K) + (r + 0.5 * sigma**2) * T) / (sigma * np.sqrt(T)) d2 = d1 - sigma * np.sqrt(T) if option_type == 'call': price = S * norm.cdf(d1) - K * np.exp(-r * T) * norm.cdf(d2) else: price = K * np.exp(-r * T) * norm.cdf(-d2) - S * norm.cdf(-d1) return price def objective_function(sigma): return bs_price(sigma) - option_price try: implied_vol = brentq(objective_function, 1e-6, 5, maxiter=1000) except (ValueError, RuntimeError): implied_vol = np.nan return implied_vol def bs_price(S, K, T, r, sigma, option_type='call'): sigma = max(sigma, 1e-6) d1 = (np.log(S / K) + (r + 0.5 * sigma**2) * T) / (sigma * np.sqrt(T)) d2 = d1 - sigma * np.sqrt(T) if option_type == 'call': price = S * norm.cdf(d1) - K * np.exp(-r * T) * norm.cdf(d2) else: price = K * np.exp(-r * T) * norm.cdf(-d2) - S * norm.cdf(-d1) return price def calculate_implied_pdf(prices, strikes): interp_func = interp1d(strikes, prices, kind='cubic', fill_value='extrapolate') fine_strikes = np.linspace(min(strikes), max(strikes), 500) prices_fine = interp_func(fine_strikes) price_d2 = np.gradient(np.gradient(prices_fine, fine_strikes), fine_strikes) return price_d2, fine_strikes def normalize_pdf(pdf, strikes): pdf = np.maximum(pdf, 0) total_prob_before = np.trapezoid(pdf, strikes) if total_prob_before == 0: total_prob_before = 1e-6 pdf_normalized = pdf / total_prob_before cdf = cumulative_trapezoid(pdf_normalized, strikes, initial=0) cdf_normalized = cdf / cdf[-1] return pdf_normalized, cdf_normalized def normalize_pdf_within_zoomed_range(pdf, strikes, total_prob_zoomed): pdf_normalized = pdf / total_prob_zoomed cdf_zoomed = cumulative_trapezoid(pdf_normalized, strikes, initial=0) cdf_zoomed_normalized = cdf_zoomed / cdf_zoomed[-1] return pdf_normalized, cdf_zoomed_normalized def get_spot_price(filename): with open(filename, 'r', encoding='utf-8') as optionsFile: optionsFileData = optionsFile.readlines() spotLine = optionsFileData[1] spotPrice_str = spotLine.split('Last:')[1].split('Change:')[0].strip() spotPrice = float(spotPrice_str.replace(',', '')) return spotPrice def find_percentile(strikes, cdf, percentile): idx = np.abs(cdf - percentile / 100).argmin() return strikes[idx] # Load the data file_path = r'spx_quotedata.csv' option_data = pd.read_csv(file_path, skiprows=3) # Parse and filter by expiration date option_data['Expiration Date'] = pd.to_datetime(option_data['Expiration Date']).dt.date # Define the three target dates target_dates = [pd.Timestamp('2025-06-27').date(), pd.Timestamp('2025-07-18').date(), pd.Timestamp('2025-08-15').date()] results = process_expiration_data(option_data, target_dates) # Get spot ...