forked from animatedread/Warrior_EA
6192393511
Hypotheses were registered in FX_PLAN.md before each round. Trend, breakout, cross reversion, hour seasonality, month-end USD, carry-cross dip-buy, metals dip-buy and flight-to-safety all fail the bar. The weekend-gap fade looked like the best result of the project on bar data (OOS t 20, 28/28 pairs) and loses on real ticks (EURCHF PF 0.52, AUDNZD PF 0.53): the Sunday-open spread is as wide as the gap. WarriorGapFade is kept as the research artifact that proved it and is flagged DO NOT TRADE. Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
65 lines
2.9 KiB
Python
65 lines
2.9 KiB
Python
"""
|
|
Is the START of the month special for metals, or is it just drift?
|
|
|
|
Every trading-day-of-month position gets the same treatment: the 2-day window
|
|
entered at the close of TDOM k-1 and exited at the close of TDOM k+1. If the
|
|
first-2-days window is merely riding the 2016-26 metals rally, every position
|
|
earns about the drift and TDOM 1-2 is not special. If it is a genuine turn-of-
|
|
month effect, it should rank at or near the top of ~20 positions in BOTH IS and
|
|
OOS, and beat the drift baseline by a margin a permutation test respects.
|
|
"""
|
|
|
|
from __future__ import annotations
|
|
|
|
import sys
|
|
import numpy as np
|
|
|
|
sys.path.insert(0, __file__.rsplit("\\", 1)[0] if "\\" in __file__ else ".")
|
|
import fx_screen as fs # noqa: E402
|
|
|
|
SPLIT = np.datetime64("2016-01-01")
|
|
|
|
|
|
def tdom(ts):
|
|
"""1-based trading-day-of-month for each bar."""
|
|
mon = ts.astype("datetime64[D]").astype("datetime64[M]")
|
|
out = np.ones(len(ts), int)
|
|
for i in range(1, len(ts)):
|
|
out[i] = 1 if mon[i] != mon[i - 1] else out[i - 1] + 1
|
|
return out
|
|
|
|
|
|
def profile(sym, width=2):
|
|
d = fs.load(sym, "D1")
|
|
ts, c = d["ts"], d["c"]
|
|
k = tdom(ts)
|
|
r = np.full(len(c), np.nan)
|
|
r[width:] = (c[width:] - c[:-width]) / c[:-width] # window ENDING at bar i
|
|
end_k = k # TDOM of the exit bar
|
|
res = {}
|
|
for pos in range(width, 21): # exit on TDOM pos -> window TDOM pos-1..pos
|
|
m = (end_k == pos) & np.isfinite(r)
|
|
res[pos] = (r[m & (ts < SPLIT)], r[m & (ts >= SPLIT)])
|
|
return res, np.nanmean(r[ts < SPLIT]), np.nanmean(r[ts >= SPLIT])
|
|
|
|
|
|
if __name__ == "__main__":
|
|
rng = np.random.default_rng(0)
|
|
for sym in sys.argv[1:] or ["XAUUSD", "XAGUSD", "XPTUSD", "XPDUSD", "XAUEUR"]:
|
|
res, dis, doos = profile(sym)
|
|
is_m = {p: v[0].mean() * 1e4 for p, v in res.items()}
|
|
oo_m = {p: v[1].mean() * 1e4 for p, v in res.items()}
|
|
rank_is = sorted(is_m, key=is_m.get, reverse=True).index(2) + 1
|
|
rank_oo = sorted(oo_m, key=oo_m.get, reverse=True).index(2) + 1
|
|
#--- permutation: how often does a RANDOM set of month-windows of the same
|
|
#--- size beat the TDOM-1..2 excess over drift, OOS
|
|
oos_all = np.concatenate([v[1] for v in res.values()])
|
|
tom = res[2][1]
|
|
exc = tom.mean() - oos_all.mean()
|
|
perm = np.mean([rng.choice(oos_all, len(tom), replace=False).mean() - oos_all.mean() >= exc
|
|
for _ in range(5000)])
|
|
print(f"{sym}: TDOM1-2 window IS {is_m[2]:+6.1f} bp (rank {rank_is}/19, all-window mean {dis * 1e4:+.1f})"
|
|
f" OOS {oo_m[2]:+6.1f} bp (rank {rank_oo}/19, all-window mean {doos * 1e4:+.1f})"
|
|
f" OOS excess {exc * 1e4:+.1f} bp, permutation p = {perm:.4f}")
|
|
top = sorted(oo_m.items(), key=lambda x: -x[1])[:4]
|
|
print(f" OOS top windows (exit TDOM: bp): {[(p, round(v, 1)) for p, v in top]}")
|