Community Scripts

# ╔══════════════════════════════════════════════════════════════╗
# ║                  STRATEGY REQUEST LOG                       ║
# ╚══════════════════════════════════════════════════════════════╝
# Generated        : 2026-05-08 02:08:02
# Model            : XGBoost
# Feature Eng.     : Auto-add features: ON
# Signal / Entry   : —
# Optimization     : —
# Risk Mgmt        : —
# Risk Filter      : —
# ══════════════════════════════════════════════════════════════

# ============================================================
# SECTION 0 — IMPORTS & CONSTANTS
import numpy as np
import pandas as pd

DATA_PATH = "/root/Desktop/QuantifyMe/data/ohlc/USDJPY_1min.parquet"
START_DATE = "2026-05-04 00:00:00"
END_DATE   = "2026-05-07 00:00:00"
VALIDATION_DATE = ""
TRAIN_SPLIT = 0.6993736951983298

# SECTION 1 — FEATURE ENGINEERING
def feature_engineering(df, close, open_, high, low):
    # --- Returns over multiple horizons ---
    for n in [1, 3, 5, 10, 20]:
        df[f"ret_{n}"] = close.pct_change(n)

    # --- RSI 14 ---
    delta = close.diff()
    gain = delta.clip(lower=0)
    loss = -delta.clip(upper=0)
    avg_gain = gain.ewm(com=13, min_periods=14).mean()
    avg_loss = loss.ewm(com=13, min_periods=14).mean()
    rs = avg_gain / (avg_loss + 1e-10)
    df["rsi_14"] = 100.0 - (100.0 / (1.0 + rs))

    # --- RSI derived features ---
    df["rsi_14_zscore"] = (df["rsi_14"] - df["rsi_14"].rolling(50).mean()) / (df["rsi_14"].rolling(50).std() + 1e-10)
    df["rsi_overbought"] = np.where(df["rsi_14"] > 70, 1, 0)
    df["rsi_oversold"]   = np.where(df["rsi_14"] < 30, 1, 0)

    # --- Bollinger Bands 20, 2 ---
    bb_mid   = close.rolling(20).mean()
    bb_std   = close.rolling(20).std()
    bb_upper = bb_mid + 2.0 * bb_std
    bb_lower = bb_mid - 2.0 * bb_std
    df["bb_mid"]   = bb_mid
    df["bb_upper"] = bb_upper
    df["bb_lower"] = bb_lower
    df["bb_width"]  = (bb_upper - bb_lower) / (bb_mid + 1e-10)
    df["bb_pct_b"]  = (close - bb_lower) / (bb_upper - bb_lower + 1e-10)
    df["bb_above"]  = np.where(close > bb_upper, 1, 0)
    df["bb_below"]  = np.where(close < bb_lower, 1, 0)

    # --- Moving Averages ---
    for w in [50, 100, 200]:
        df[f"sma_{w}"] = close.rolling(w).mean()
        df[f"price_vs_sma_{w}"] = (close - df[f"sma_{w}"]) / (df[f"sma_{w}"] + 1e-10)

    # --- MA crossover signals ---
    df["sma50_vs_sma100"] = np.where(df["sma_50"] > df["sma_100"], 1, -1)
    df["sma50_vs_sma200"] = np.where(df["sma_50"] > df["sma_200"], 1, -1)
    df["sma100_vs_sma200"] = np.where(df["sma_100"] > df["sma_200"], 1, -1)

    # --- ATR 14 ---
    tr = pd.concat([
        high - low,
        (high - close.shift(1)).abs(),
        (low  - close.shift(1)).abs()
    ], axis=1).max(axis=1)
    df["atr_14"] = tr.ewm(com=13, min_periods=14).mean()
    df["natr_14"] = df["atr_14"] / (close + 1e-10)

    # --- Momentum / rate of change ---
    for n in [5, 10, 20]:
        df[f"mom_{n}"] = close - close.shift(n)
        df[f"roc_{n}"] = (close - close.shift(n)) / (close.shift(n) + 1e-10)

    # --- Volume features (if volume exists) ---
    if "volume" in df.columns:
        vol = df["volume"].replace(0, np.nan)
        df["vol_sma_20"]   = vol.rolling(20).mean()
        df["vol_ratio_20"] = vol / (df["vol_sma_20"] + 1e-10)
    else:
        df["vol_ratio_20"] = 1.0

    # --- Price spread & body features ---
    df["hl_spread"]  = (high - low) / (close + 1e-10)
    df["body_ratio"] = (close - open_).abs() / (high - low + 1e-10)
    df["upper_wick"] = (high - pd.concat([close, open_], axis=1).max(axis=1)) / (high - low + 1e-10)
    df["lower_wick"] = (pd.concat([close, open_], axis=1).min(axis=1) - low) / (high - low + 1e-10)
    df["bull_candle"] = np.where(close > open_, 1, 0)

    # --- Lagged returns for autocorrelation signal ---
    for lag in [1, 2, 3, 5]:
        df[f"ret1_lag{lag}"] = df["ret_1"].shift(lag)

    # --- Rolling volatility ---
    df["vol_10"]  = df["ret_1"].rolling(10).std()
    df["vol_20"]  = df["ret_1"].rolling(20).std()
    df["vol_50"]  = df["ret_1"].rolling(50).std()

    # --- Z-score of close over 20 and 50 bars ---
    df["zscore_20"] = (close - close.rolling(20).mean()) / (close.rolling(20).std() + 1e-10)
    df["zscore_50"] = (close - close.rolling(50).mean()) / (close.rolling(50).std() + 1e-10)

    # --- Relative distance of price from BB bands ---
    df["dist_upper"] = (bb_upper - close) / (close + 1e-10)
    df["dist_lower"] = (close - bb_lower) / (close + 1e-10)

    # --- EMA 9 and 21 for short-term momentum ---
    df["ema_9"]  = close.ewm(span=9, min_periods=9).mean()
    df["ema_21"] = close.ewm(span=21, min_periods=21).mean()
    df["ema9_vs_ema21"] = np.where(df["ema_9"] > df["ema_21"], 1, -1)
    df["price_vs_ema9"]  = (close - df["ema_9"]) / (df["ema_9"] + 1e-10)
    df["price_vs_ema21"] = (close - df["ema_21"]) / (df["ema_21"] + 1e-10)

    # --- MACD-like signal ---
    ema_12 = close.ewm(span=12, min_periods=12).mean()
    ema_26 = close.ewm(span=26, min_periods=26).mean()
    macd_line   = ema_12 - ema_26
    signal_line = macd_line.ewm(span=9, min_periods=9).mean()
    df["macd"]         = macd_line
    df["macd_signal"]  = signal_line
    df["macd_hist"]    = macd_line - signal_line
    df["macd_cross"]   = np.where(macd_line > signal_line, 1, -1)

    # --- Fill NaN from warm-up periods ---
    df = df.bfill().ffill()

    return df

# SECTION 2 — STRATEGY CONFIG
def strategy_config():
    return {
        "title":            "USD/JPY Multi-MA + RSI/BB XGBoost Sharpe",
        "model_type":       "XGBClassifier",
        "model_params":     {
            "n_estimators":     400,
            "max_depth":        4,
            "learning_rate":    0.04,
            "subsample":        0.75,
            "colsample_bytree": 0.75,
            "min_child_weight": 5,
            "gamma":            0.1,
            "reg_alpha":        0.1,
            "reg_lambda":       1.5,
            "objective":        "binary:logistic",
            "tree_method":      "hist",
            "random_state":     42,
        },
        "signal_threshold": 0.55,
        "direction":        "both",
        "stop_loss":        0.0008,
        "take_profit":      0.0016,
        "cooldown":         0,
        "max_positions":    1,
        "on_opposite":      "close_only",
        "session_filter":   None,
        "min_atr":          None,
        "trend_filter":     None,
        "target_horizon":   5,
        "objective":        (
            "Maximize Sharpe ratio on USD/JPY 1-min data using XGBoost with "
            "returns, RSI, Bollinger Bands, multiple MAs (50/100/200), MACD, "
            "ATR, and candle-body features. Stop-loss and take-profit set at "
            "a 1:2 risk/reward to filter noise and improve Sharpe. n_estimators "
            "and moderate depth balance bias-variance. Regularization (alpha/lambda) "
            "reduces overfitting on short date range."
        ),
        "notes":            (
            "Target horizon of 5 bars (5 minutes) is chosen to capture short-term "
            "directional moves on 1-min data without excessive label noise. "
            "colsample_bytree and subsample add stochasticity to reduce variance. "
            "close_only on opposite signal avoids whipsaw from rapid reversals. "
            "No session filter applied since USD/JPY has liquidity around the clock."
        ),
    }