3,293 answered market questions

WITH
closes AS
(
    SELECT
        date,
        toFloat64(any(close)) AS px
    FROM global_markets.stocks_daily_aggs
    WHERE ticker = 'SPY'
      AND date >= '2021-06-01'
    GROUP BY date
),
daily_returns AS
(
    SELECT
        date,
        pow(log(px / lagInFrame(px) OVER (ORDER BY date ROWS BETWEEN 1 PRECEDING AND CURRENT ROW)), 2) AS r2
    FROM closes
),
forward_rv AS
(
    SELECT
        date,
        100 * sqrt(252 * avg(r2) OVER (ORDER BY date ROWS BETWEEN 1 FOLLOWING AND 21 FOLLOWING)) AS rv_pct,
        count(r2) OVER (ORDER BY date ROWS BETWEEN 1 FOLLOWING AND 21 FOLLOWING)                 AS sessions_ahead
    FROM daily_returns
    WHERE isFinite(r2)
),
implied AS
(
    SELECT
        date,
        100 * avg(toFloat64(implied_volatility)) AS iv_pct
    FROM global_markets.options_greeks
    WHERE underlying_symbol = 'SPY'
      AND iv_converged = 1
      AND volume > 0
      AND days_to_expiry BETWEEN 20 AND 45
      AND abs(toFloat64(strike_price) / toFloat64(underlying_close) - 1) < 0.05
      AND date >= '2021-08-01'
    GROUP BY date
    HAVING count() >= 20
),
paired AS
(
    SELECT
        i.date              AS date,
        i.iv_pct            AS iv_pct,
        f.rv_pct            AS rv_pct,
        i.iv_pct - f.rv_pct AS premium_pct
    FROM implied AS i
    INNER JOIN forward_rv AS f ON f.date = i.date
    WHERE f.sessions_ahead = 21
)
SELECT
    formatDateTime(argMin(date, premium_pct), '%b %e, %Y') AS label,
    round(argMin(iv_pct, premium_pct), 2)                  AS implied_vol_pct,
    round(argMin(rv_pct, premium_pct), 2)                  AS realized_vol_pct,
    round(-1 * min(premium_pct), 2)                        AS shortfall_pct
FROM paired
GROUP BY toStartOfMonth(date)
ORDER BY min(premium_pct) ASC
LIMIT 10