미국 ETF,
이름 끝의 (H)

# -*- coding: utf-8 -*-
r"""미국 ETF — 「이름 끝의 H, 붙은 걸 살까 그냥 살까」 편 재현 코드 (AlgoLab, 2026-09-18)

영상에서 말한 숫자를 **이 파일 하나로** 그대로 다시 냅니다.

  python run_fx.py

필요한 것: 파이썬 표준 라이브러리만. (numpy·pandas 불필요)
데이터: 옆 폴더 `data_fx/` 에 그날 시세를 얼려 두었습니다. 그 폴더가 있으면 그것을 쓰고,
없으면 야후 파이낸스에서 새로 받습니다. **새로 받으면 분배금 보정이 갱신돼 끝자리가
조금 달라질 수 있습니다** — 영상과 똑같은 숫자를 보려면 얼린 시세를 그대로 두십시오.

무엇을 견줬나
  · 같은 지수를 담는 국내상장 ETF 두 개 — 이름 끝에 (H) 가 붙은 것과 안 붙은 것
  · 짝비교 5번 (상품 8개)
      S&P500    360750(그냥) ↔ 143850(H·선물)  ·  379800(그냥) ↔ 219480(H·선물)
      S&P500    360750(그냥) ↔ 448290(H·현물)  ← 선물 구조를 걷어낸 짝
      나스닥100  133690(그냥) ↔ 304940(H·선물)  ·  379810(그냥) ↔ 304940(H·선물)
  · ★ (H) 상품 넷 중 셋이 선물로 지수를 따라가는 상품입니다(상품명 조회로 확인).
    그래서 선물이 아닌 현물 환헤지 상품(448290)으로 한 짝을 더 만들었습니다.
  · 헤드라인은 **구간이 가장 긴 짝** — 고르는 규칙을 먼저 정하고 결과를 봤습니다
    (나스닥100 · 133690 ↔ 304940 · 2018-08-30 ~ 2026-09-17 · 8.05년)

정직하게 잰 방법
  · 분배금 보정(총수익) 종가 · 사고팔 때 편도 0.1%
  · 두 종목이 모두 거래된 날만 씁니다(같은 날짜·같은 구간)
  · 적립식 — 매달 첫 거래일 종가에 100만원. 거치식 결과도 같이 냅니다
  · (H) 여부는 **상품 이름이 아니라 데이터로** 갈랐습니다. 하루치 수익률을
    「미국 지수 수익률」과 「원/달러 수익률」에 최소제곱으로 회귀해, 원/달러 계수가
    0.12 아래면 (H), 위면 그냥형입니다. 실제로는 -0.03~-0.01 과 +0.32~+0.38 로
    중간값 없이 두 무리로 갈립니다.

안 넣은 것: 호가 스프레드·괴리율·세금(양도세·배당소득세)·연금계좌 과세이연.
매수만 하고 하락에 거는 기능은 없습니다. 여덟 종목 다 지금 상장돼 있습니다(생존편향).
환율이 위기에 오른다는 것은 지난 8년의 성질입니다. 투자 권유가 아닙니다.
"""
import csv
import datetime
import io
import json
import os
import sys

# 한국어 윈도 기본 콘솔(cp949)에서 대시·화살표가 UnicodeEncodeError 를 낸다.
# 시청자가 받아 가는 파일이라 첫 줄에서 죽으면 안 된다.
try:
    sys.stdout.reconfigure(encoding="utf-8", errors="replace")
except Exception:
    pass

ASOF = "2026-09-17"          # 영상을 만든 날의 직전 거래일. None 으로 바꾸면 오늘까지
COST = 0.001                 # 편도 0.1%
DEP = 1000000.0              # 적립식 매달 납입
HERE = os.path.dirname(os.path.abspath(__file__))
FROZEN = os.path.join(HERE, "data_fx")

PAIRS = [("S&P500",   "360750.KS", "143850.KS", "선물"),
         ("S&P500",   "379800.KS", "219480.KS", "선물"),
         ("S&P500",   "360750.KS", "448290.KS", "현물"),
         ("나스닥100", "133690.KS", "304940.KS", "선물"),
         ("나스닥100", "379810.KS", "304940.KS", "선물")]
# 국내 **정식 종목명**이다. 해외 데이터 제공처의 영문명은 143850 에서 (H) 를 빠뜨린다.
NAMES = {"133690": "TIGER 미국나스닥100 (현물·환노출)",
         "379810": "KODEX 미국나스닥100 (현물·환노출)",
         "360750": "TIGER 미국S&P500 (현물·환노출)",
         "379800": "KODEX 미국S&P500 (현물·환노출)",
         "304940": "KODEX 미국나스닥100선물(H) (선물·환헤지)",
         "143850": "TIGER 미국S&P500선물(H) (선물·환헤지)",
         "219480": "KODEX 미국S&P500선물(H) (선물·환헤지)",
         "448290": "TIGER 미국S&P500(H) (현물·환헤지)"}
IDX = {"S&P500": "^GSPC", "나스닥100": "^NDX"}   # ^IXIC 는 나스닥 종합 — 100 이 아니다
FX = "KRW=X"


# ── 시세 ──────────────────────────────────────────────────────────────
def _path(sym):
    safe = sym.replace(".", "_").replace("^", "_").replace("=", "_")
    return os.path.join(FROZEN, safe + "_1d.csv")


def load(sym):
    if not os.path.exists(_path(sym)):
        fetch(sym)
    out = []
    with io.open(_path(sym), encoding="utf-8") as f:
        for row in csv.reader(f):
            if not row or row[0] == "date":
                continue
            d = row[0][:10]
            if ASOF and d > ASOF:
                continue
            try:
                out.append((d, float(row[4])))
            except (ValueError, IndexError):
                pass
    if not out:
        raise SystemExit("시세가 비었습니다: " + sym)
    return out


def fetch(sym):
    import urllib.parse
    import urllib.request
    os.makedirs(FROZEN, exist_ok=True)
    url = ("https://query1.finance.yahoo.com/v8/finance/chart/%s"
           "?range=30y&interval=1d&events=div%%2Csplit" % urllib.parse.quote(sym))
    req = urllib.request.Request(url, headers={"User-Agent": "Mozilla/5.0"})
    j = json.load(urllib.request.urlopen(req, timeout=60))["chart"]["result"][0]
    ts, q = j["timestamp"], j["indicators"]["quote"][0]
    adj = (j["indicators"].get("adjclose") or [{}])[0].get("adjclose") or q["close"]
    with io.open(_path(sym), "w", encoding="utf-8", newline="") as f:
        w = csv.writer(f)
        w.writerow(["date", "open", "high", "low", "close", "volume"])
        for i, t in enumerate(ts):
            if q["close"][i] is None or adj[i] is None:
                continue
            k = adj[i] / q["close"][i]
            d = datetime.datetime.utcfromtimestamp(t).strftime("%Y-%m-%d")
            w.writerow([d,
                        q["open"][i] * k if q["open"][i] else adj[i],
                        q["high"][i] * k if q["high"][i] else adj[i],
                        q["low"][i] * k if q["low"][i] else adj[i],
                        adj[i], q["volume"][i] or 0])
    print("  받았습니다 %-12s" % sym)


# ── 계산 ──────────────────────────────────────────────────────────────
def 보유(dates, px):
    g = [1.0 - COST]
    for i in range(1, len(px)):
        g.append(g[-1] * px[i] / px[i - 1])
    d0 = datetime.date(*map(int, dates[0].split("-")))
    d1 = datetime.date(*map(int, dates[-1].split("-")))
    yrs = (d1 - d0).days / 365.25
    cagr = g[-1] ** (1 / yrs) - 1
    peak, mdd = g[0], 0.0
    for v in g:
        peak = max(peak, v)
        mdd = min(mdd, v / peak - 1)
    return {"CAGR": cagr, "MDD": mdd, "Calmar": cagr / abs(mdd), "배수": g[-1], "년": yrs}


def 적립(dates, px):
    first, seen = set(), set()
    for i, d in enumerate(dates):
        if d[:7] not in seen:
            seen.add(d[:7])
            first.add(i)
    units, bal = 0.0, []
    for i in range(len(dates)):
        if i in first:
            units += DEP * (1 - COST) / px[i]
        bal.append(units * px[i])
    peak, mdd = 0.0, 0.0
    for v in bal:
        peak = max(peak, v)
        if peak > 0:
            mdd = min(mdd, v / peak - 1)
    cf = [-DEP if i in first else 0.0 for i in range(len(dates))]
    cf[-1] += bal[-1]
    lo, hi = -0.99, 3.0
    for _ in range(200):
        mid = (lo + hi) / 2
        v = sum(x / ((1 + mid) ** (i / 252.0)) for i, x in enumerate(cf))
        lo, hi = (mid, hi) if v > 0 else (lo, mid)
    return {"최종잔고": bal[-1], "총납입": DEP * len(first), "배수": bal[-1] / (DEP * len(first)),
            "IRR": (lo + hi) / 2, "MDD": mdd, "납입횟수": len(first)}


def 환율베타(etf, idx, fx):
    """일간수익률 r_etf = a + b1*r_idx + b2*r_fx 를 최소제곱으로 푼다(3x3 정규방정식)."""
    em = {d: i for i, (d, _) in enumerate(etf)}
    im = {d: i for i, (d, _) in enumerate(idx)}
    fm = {d: i for i, (d, _) in enumerate(fx)}
    days = sorted(set(em) & set(im) & set(fm))
    X, y = [], []
    for k in range(1, len(days)):
        p, q = days[k - 1], days[k]
        if em[q] - em[p] != 1 or im[q] < 1 or im[p] < 1:
            continue
        X.append([1.0,
                  idx[im[q] - 1][1] / idx[im[p] - 1][1] - 1,
                  fx[fm[q]][1] / fx[fm[p]][1] - 1])
        y.append(etf[em[q]][1] / etf[em[p]][1] - 1)
    n = 3
    A = [[sum(X[r][i] * X[r][j] for r in range(len(X))) for j in range(n)] for i in range(n)]
    b = [sum(X[r][i] * y[r] for r in range(len(X))) for i in range(n)]
    for i in range(n):                                   # 가우스 소거
        p = max(range(i, n), key=lambda r: abs(A[r][i]))
        A[i], A[p] = A[p], A[i]
        b[i], b[p] = b[p], b[i]
        for r in range(i + 1, n):
            f = A[r][i] / A[i][i]
            for c in range(i, n):
                A[r][c] -= f * A[i][c]
            b[r] -= f * b[i]
    x = [0.0] * n
    for i in range(n - 1, -1, -1):
        x[i] = (b[i] - sum(A[i][c] * x[c] for c in range(i + 1, n))) / A[i][i]
    return x[2], len(X)


def won(v, 올림=False):
    # 영상·썸네일은 **화면에 뜬 두 금액의 차이**로 말한다(2억 5,387만 − 1억 9,775만 = 5,612만원).
    # 차이만 그 규약을 따르고, 잔고는 그대로 버린다.
    man = int(round(v / 10000)) if 올림 else int(v // 10000)
    e, m = divmod(man, 10000)
    return ("%d억 %s만원" % (e, format(m, ","))) if e else ("%s만원" % format(man, ","))


def main():
    syms = sorted({s for _, a, b, _c in PAIRS for s in (a, b)} | set(IDX.values()) | {FX})
    if not os.path.isdir(FROZEN):
        print("얼린 시세가 없어 야후에서 받습니다 (%d종목)" % len(syms))
    series = {s: load(s) for s in syms}
    fx = series[FX]

    print("\n── (H) 여부를 이름이 아니라 데이터로 가른다 ──────────────────")
    print("   일간수익률의 원/달러 계수 — 0 에 붙으면 (H), 0.3 을 넘으면 그냥형")
    for name, idxsym in IDX.items():
        for sym in sorted({s for n, a, b, _c in PAIRS if n == name for s in (a, b)}):
            beta, cnt = 환율베타(series[sym], series[idxsym], fx)
            print("   %-12s 원/달러 계수 %+.3f  (표본 %5d일)  → %-9s %s"
                  % (sym, beta, cnt, "(H) 환헤지" if beta < 0.12 else "그냥형",
                     NAMES.get(sym.split(".")[0], "")))

    print("\n── 짝비교 %d번 ───────────────────────────────────────────" % len(PAIRS))
    승 = {"보유 해마다": 0, "보유 낙폭": 0, "적립 남은 돈": 0, "적립 낙폭": 0}
    rows = []
    for name, un, hd, 구조 in PAIRS:
        du, dh = dict(series[un]), dict(series[hd])
        days = sorted(set(du) & set(dh))
        pu, ph = [du[d] for d in days], [dh[d] for d in days]
        hu, hh = 보유(days, pu), 보유(days, ph)
        au, ah = 적립(days, pu), 적립(days, ph)
        승["보유 해마다"] += hu["CAGR"] > hh["CAGR"]
        승["보유 낙폭"] += hu["MDD"] > hh["MDD"]
        승["적립 남은 돈"] += au["최종잔고"] > ah["최종잔고"]
        승["적립 낙폭"] += au["MDD"] > ah["MDD"]
        rows.append((name, un, hd, days, hu, hh, au, ah, 구조))
        print("\n  %s  %s(그냥) ↔ %s(H·%s)   %s ~ %s  %.2f년"
              % (name, un[:6], hd[:6], 구조, days[0], days[-1], hu["년"]))
        print("    거치식  그냥 해마다 %5.2f%% · 가장 많이 빠졌을 때 %6.2f%%"
              "   |  (H) 해마다 %5.2f%% · %6.2f%%"
              % (hu["CAGR"] * 100, hu["MDD"] * 100, hh["CAGR"] * 100, hh["MDD"] * 100))
        print("    적립식  그냥 %14s · 가장 많이 빠졌을 때 %6.2f%%"
              "   |  (H) %14s · %6.2f%%"
              % (won(au["최종잔고"]), au["MDD"] * 100, won(ah["최종잔고"]), ah["MDD"] * 100))
        print("            매달 100만원 %d번 · 넣은 돈 %s" % (au["납입횟수"], won(au["총납입"])))

    n = len(PAIRS)
    print("\n── 집계 — 그냥형(H 없는 쪽)이 이긴 횟수 ──────────────────────")
    for k, v in 승.items():
        print("   %-12s %d/%d" % (k, v, n))

    head = max(rows, key=lambda r: r[4]["년"])          # 규칙: 구간이 가장 긴 짝
    name, un, hd, days, hu, hh, au, ah, 구조 = head
    print("\n── 영상이 말한 숫자 (헤드라인 = 구간이 가장 긴 짝) ───────────")
    print("   %s · %s(그냥) ↔ %s(H·%s) · %s ~ %s · %.2f년"
          % (name, un[:6], hd[:6], 구조, days[0], days[-1], hu["년"]))
    print("   ※ 이 짝이 다섯 짝 중 가장 크게 벌어진 짝입니다. 나머지 넷은 그보다 작습니다:")
    for nm, u, h, dd, _hu, _hh, _au, _ah, st in rows:
        if (nm, u, h) != (name, un, hd):
            print("      %-8s %s↔%s(%s)  %s" % (nm, u[:6], h[:6], st, won(_au["최종잔고"] - _ah["최종잔고"], 올림=True)))
    print("   매달 100만원씩 %d번 · 넣은 돈 %s" % (au["납입횟수"], won(au["총납입"])))
    print("   H 없는 쪽 남은 돈   %s" % won(au["최종잔고"]))
    print("   H 붙은 쪽 남은 돈   %s" % won(ah["최종잔고"]))
    print("   차이               %s   (영상·썸네일 표기)" % won(au["최종잔고"] - ah["최종잔고"], 올림=True))
    print("   가장 많이 빠졌을 때  H 없는 쪽 %.2f%% · H 붙은 쪽 %.2f%%"
          % (au["MDD"] * 100, ah["MDD"] * 100))
    print("\n   ※ 영상은 만원 미만을 버리고 낙폭은 소수를 버려 말합니다"
          " — 둘 다 우리 주장에 불리한 쪽입니다.")


if __name__ == "__main__":
    main()