📈 TinyCTA
A Lightweight Python Package for Commodity Trading Advisor Strategies.
Quick Links: 📚 Repository • 📦 PyPI • 🐛 Issues • 💬 Discussions
TinyCTA provides essential tools for quantitative finance and algorithmic trading, particularly for trend-following strategies. The package includes:
- Polars-based signal processing: oscillators, moving-average crossovers, and volatility-adjusted returns
- Robust volatility estimation via rolling median absolute deviation
- Linear algebra utilities that handle matrices with missing values
- Matrix shrinkage techniques commonly used in portfolio optimization
This package is designed to be the foundation for implementing CTA strategies in just a few lines of code, hence the name "TinyCTA".
📖 New here? Follow the end-to-end CTA tutorial
to go from raw prices through signals and the Engine to cash positions.
pip install tinyctaThe core install keeps a minimal dependency footprint (numpy, polars, pydantic, cvx-linalg).
The optional Optuna-based hyperparameter-optimisation layer (tinycta.hyper) is installed via
the hyper extra:
pip install "tinycta[hyper]"Clone the repository and install using the provided Makefile:
git clone https://github.com/tschm/tinycta.git
cd tinycta
make installThis will install uv (a fast Python package installer) and create a virtual environment with all dependencies.
import polars as pl
from tinycta.osc import osc
prices = pl.DataFrame({"A": [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]})
result = prices.with_columns(osc(pl.col("A"), fast=2, slow=6).alias("osc_A"))import polars as pl
from tinycta.ewma import ma_cross
prices = pl.DataFrame({"A": [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]})
result = prices.with_columns(
ma_cross(pl.col("A"), fast=2, slow=6).alias("sig_A")
)import polars as pl
from tinycta.util import vol_adj, adj_log_prices
prices = pl.DataFrame({"A": [100, 101, 99, 102, 98, 103]})
result = prices.with_columns(
vol_adj(pl.col("A"), vola=3, clip=4.2).alias("vol_adj_A"),
adj_log_prices(pl.col("A"), vola=3, clip=4.2).alias("adj_log_A"),
)moving_absolute_deviation is an outlier-resistant rolling volatility estimate (a scaled
rolling median absolute deviation of log returns), and shrink2id blends a square matrix
towards the identity — a common regulariser for noisy covariance/correlation matrices.
import numpy as np
import polars as pl
from tinycta.signal import moving_absolute_deviation, shrink2id
prices = pl.DataFrame({"A": [100.0, 101.0, 99.0, 102.0, 98.0, 103.0, 97.0, 104.0]})
vol = prices.with_columns(
moving_absolute_deviation(pl.col("A"), com=2).alias("vol_A")
)
print(sorted(vol.columns))
print(vol.shape)
cov = np.array([[4.0, 1.0], [1.0, 9.0]])
print(shrink2id(cov, lamb=0.5)) # blend halfway towards the identity['A', 'vol_A']
(8, 2)
[[2.5 0.5]
[0.5 5. ]]
import numpy as np
from tinycta.linalg import solve
matrix = np.array([[1.0, 0.5], [0.5, 1.0]])
rhs = np.array([1.0, 2.0])
solution = solve(matrix, rhs)
print(np.round(solution, 10) + 0)[0. 2.]
The Engine turns aligned price and expected-return (mu) frames into
correlation-shrinkage-optimized cash positions. It is configured by a validated Config.
import polars as pl
from tinycta.config import Config
from tinycta.engine import Engine
dates = list(range(1, 9))
prices = pl.DataFrame({"date": dates, "A": [100.0, 101.0, 102.0, 101.5, 103.0, 104.0, 103.5, 105.0]})
mu = pl.DataFrame({"date": dates, "A": [0.0, 0.1, 0.2, 0.1, 0.2, 0.3, 0.2, 0.3]})
cfg = Config(vola=3, corr=3, clip=4.2, shrink=0.5)
positions = Engine(prices=prices, mu=mu, cfg=cfg).cash_position
# .cash_position mirrors the input frame: one row per date, one column per asset.
print(sorted(positions.columns))
print(positions.shape)['A', 'date']
(8, 2)
Config is a frozen Pydantic model: vola, corr (must be >= vola) and clip must be
positive, and shrink must lie in [0, 1].
tinycta.hyper.optimize runs an Optuna study over a function that builds a portfolio from a
trial and scores it by Sharpe ratio, returning a frozen Study.
from tinycta.hyper import optimize
def suggest_portfolio(trial):
fast = trial.suggest_int("fast", 2, 20)
slow = trial.suggest_int("slow", fast + 1, 100)
# ... build and return a jquantstats Portfolio from the suggested params ...
return build_portfolio(fast, slow)
study = optimize(suggest_portfolio, n_trials=100, seed=42)
print(study.best_params, study.best_value)tinycta.hyper.get_config bundles the config sections and a configured logger for a notebook
experiment. It reads a shared config.yml and, for an experiment named name, an optional
experiment-specific config/{name}.yml; each data / params / optuna section is taken from
config.yml when present, otherwise from the sibling file. All three sections are optional.
# config.yml — every section is optional
data:
output_path: output # output dir, relative to the notebooks directory (default: "output")
params:
fast: 12 # arbitrary experiment parameters
optuna:
n_trials: 100 # arbitrary Optuna settingsPaths resolve relative to the notebooks directory — the parent of config.yml, or its
grandparent when config.yml lives inside a config/ subdirectory. Outputs are written to
{notebooks}/{output_path}/{name}/, and a config-supplied output_path is confined to the
notebooks directory (a traversing or absolute path raises ValueError). Set the
NOTEBOOK_OUTPUT_FOLDER environment variable to override the output directory entirely — this
explicit operator override is trusted and not confined.
from tinycta.hyper import get_config
cfg = get_config("my_experiment") # reads ./config.yml (+ ./config/my_experiment.yml)
cfg.logger.info("run starting") # loguru logger, also writing to output.log
fast = cfg.params["fast"] # config sections as plain dictsosc(x, fast, slow, min_samples=1)— analytically scaled EWMA-difference oscillator (Polars)ma_cross(prices, fast, slow, min_samples=1)— sign of fast-vs-slow EWM crossover: -1, 0, or +1 (Polars)vol_adj(x, vola, clip, min_samples=1)— clipped, volatility-adjusted log returns (Polars)adj_log_prices(x, vola, clip, min_samples=1)— cumulative sum of volatility-adjusted log returns (Polars)
moving_absolute_deviation(price, com=32)— robust rolling volatility estimate via median absolute deviation (Polars)shrink2id(matrix, lamb=1.0)— shrink a matrix towards the identity matrix
valid(matrix)— extract the finite subset of a matrix by filtering NaN rows/columnsa_norm(vector, matrix=None)— matrix-norm of a vectorinv_a_norm(vector, matrix=None)— inverse matrix-norm of a vectorsolve(matrix, rhs)— solve a linear system, handling matrices with NaN values
Config(vola, corr, clip, shrink)— frozen Pydantic config;corr >= vola,vola/corr/clip > 0,shrink ∈ [0, 1]Engine(prices, mu, cfg)— correlation-aware position optimizer;.cash_positionreturns per-asset cash positions.assets,.ret_adj,.vola,.cor— intermediate per-asset/per-timestamp quantities
optimize(suggest_portfolio_fn, n_trials=100, seed=42)— run an Optuna study scored by Sharpe; returns aStudyStudy— frozen result wrapper exposingbest_params,best_value,n_completed,n_trials, and.plot(output_dir)get_config(name, config_path=None)— load mergeddata/params/optunasections and a configured logger; returns anExperimentConfigExperimentConfig—NamedTuplebundlingname,logger, and the optionalparams,optunaanddatasections
make installmake testmake fmtmake cleanTinyCTA is licensed under the MIT License. See the LICENSE file for details.
Contributions are welcome! Please feel free to submit a Pull Request.