mirror of
https://github.com/webclinic017/drift.git
synced 2026-07-28 03:08:01 +00:00
Mark Aron Szulyovszky
5c4a5b0cf1
* refactor(Project): removed regression method (we can still use regression models, but we'll need map them to classes later) * fix(Training): removed mistakenly left in `method` parameter
255 lines
9.7 KiB
Python
255 lines
9.7 KiB
Python
import pandas as pd
|
|
import numpy as np
|
|
from utils.types import DataSource, FeatureExtractor
|
|
from utils.helpers import deduplicate_indexes, drop_columns_if_exist
|
|
from data_loader.collections import DataCollection
|
|
from typing import Literal
|
|
import ray
|
|
import os
|
|
from config.hashing import hash_data_config
|
|
from diskcache import Cache
|
|
cache = Cache(".cachedir/data")
|
|
|
|
def load_data(**kwargs) -> tuple[pd.DataFrame, pd.Series, pd.Series]:
|
|
hashed = hash_data_config(kwargs)
|
|
if hashed in cache:
|
|
return cache.get(hashed)
|
|
else:
|
|
return_value = __load_data(**kwargs)
|
|
cache[hashed] = return_value
|
|
return return_value
|
|
|
|
def __load_data(assets: DataCollection,
|
|
other_assets: DataCollection,
|
|
exogenous_data: DataCollection,
|
|
target_asset: DataSource,
|
|
load_non_target_asset: bool,
|
|
log_returns: bool,
|
|
forecasting_horizon: int,
|
|
own_features: list[tuple[str, FeatureExtractor, list[int]]],
|
|
other_features: list[tuple[str, FeatureExtractor, list[int]]],
|
|
exogenous_features: list[tuple[str, FeatureExtractor, list[int]]],
|
|
index_column: Literal['date', 'int'],
|
|
no_of_classes: Literal['two', 'three-balanced', 'three-imbalanced'],
|
|
narrow_format: bool = False
|
|
) -> tuple[pd.DataFrame, pd.Series, pd.Series]:
|
|
"""
|
|
Loads asset data from the specified path.
|
|
Returns:
|
|
- DataFrame `X` with all the training data
|
|
- Series `y` with the target asset returns shifted by 1 day OR if it's a classification problem, the target class)
|
|
- Series `forward_returns` with the target asset returns shifted by 1 day
|
|
"""
|
|
|
|
target_file = [f for f in assets if f[1].startswith(target_asset[1])]
|
|
assert len(target_file) == 1, "There should be exactly one target file"
|
|
other_files = [f for f in assets if load_non_target_asset == True and f[1].startswith(target_asset[1]) == False]
|
|
files = other_files + other_assets
|
|
|
|
target_asset_future = [__load_df.remote(
|
|
data_source=data_source,
|
|
prefix=data_source[1],
|
|
returns='log_returns' if log_returns else 'returns',
|
|
feature_extractors=own_features,
|
|
narrow_format=narrow_format,
|
|
) for data_source in target_file]
|
|
target_asset_df = ray.get(target_asset_future)
|
|
|
|
target_asset_only_returns_future = __load_df.remote(
|
|
data_source=target_file[0],
|
|
prefix=target_file[0][1],
|
|
returns='returns',
|
|
feature_extractors=[],
|
|
narrow_format=narrow_format,
|
|
)
|
|
df_target_asset_only_returns = ray.get(target_asset_only_returns_future)
|
|
|
|
asset_futures = [__load_df.remote(
|
|
data_source=data_source,
|
|
prefix=data_source[1],
|
|
returns='log_returns' if log_returns else 'returns',
|
|
feature_extractors=other_features,
|
|
narrow_format=narrow_format,
|
|
) for data_source in files]
|
|
asset_dfs = ray.get(asset_futures)
|
|
|
|
exogenous_futures = [__load_df.remote(
|
|
data_source=data_source,
|
|
prefix=data_source[1],
|
|
returns='none',
|
|
feature_extractors=exogenous_features,
|
|
narrow_format=narrow_format,
|
|
) for data_source in exogenous_data]
|
|
exogenous_dfs = ray.get(exogenous_futures)
|
|
|
|
dfs = target_asset_df + asset_dfs + exogenous_dfs
|
|
dfs = [deduplicate_indexes(df) for df in dfs]
|
|
target_df = dfs[0]
|
|
if narrow_format:
|
|
dfs = pd.concat([df.sort_index().reindex(target_df.index) for df in dfs], axis=0).fillna(0.)
|
|
else:
|
|
dfs = pd.concat([df.sort_index().reindex(target_df.index) for df in dfs], axis=1).fillna(0.)
|
|
|
|
dfs.index = pd.DatetimeIndex(dfs.index)
|
|
|
|
if index_column == 'int':
|
|
dfs.reset_index(drop=True, inplace=True)
|
|
df_target_asset_only_returns.reset_index(drop=True, inplace=True)
|
|
|
|
if narrow_format:
|
|
dfs = dfs.drop(index=dfs.index[0], axis=0)
|
|
df_target_asset_only_returns = df_target_asset_only_returns.drop(index=dfs.index[0], axis=0)
|
|
|
|
## Create target
|
|
target_col = 'target'
|
|
returns_col = target_asset[1] + '_returns'
|
|
forward_returns = __create_target_cum_forward_returns(df_target_asset_only_returns, returns_col, forecasting_horizon)
|
|
dfs[target_col] = __create_target_classes(dfs, returns_col, forecasting_horizon, no_of_classes)
|
|
# we need to drop the last row, because we forward-shift the target (see what happens if you call .shift[-1] on a pd.Series)
|
|
dfs = dfs.iloc[:-forecasting_horizon]
|
|
forward_returns = forward_returns.iloc[:-forecasting_horizon]
|
|
|
|
X = dfs.drop(columns=[target_col])
|
|
y = dfs[target_col]
|
|
|
|
return X, y, forward_returns
|
|
|
|
|
|
@ray.remote
|
|
def __load_df(data_source: DataSource,
|
|
prefix: str,
|
|
returns: Literal['none', 'price', 'returns', 'log_returns'],
|
|
feature_extractors: list[tuple[str, FeatureExtractor, list[int]]],
|
|
narrow_format: bool = False) -> pd.DataFrame:
|
|
df = pd.read_csv(os.path.join(data_source[0], data_source[1] + '.csv'), header=0, index_col=0).fillna(0)
|
|
|
|
if returns == 'log_returns':
|
|
df['returns'] = np.log(df['close']).diff(1)
|
|
elif returns == 'price':
|
|
df['returns'] = df['close']
|
|
elif returns == 'returns':
|
|
df['returns'] = df['close'].pct_change()
|
|
|
|
df = __apply_feature_extractors(df, log_returns=True if returns == 'log_returns' else False, feature_extractors = feature_extractors)
|
|
|
|
df = df.replace([np.inf, -np.inf], 0.)
|
|
df = drop_columns_if_exist(df, ['open', 'high', 'low', 'close', 'volume'])
|
|
|
|
if narrow_format:
|
|
df["ticker"] = np.repeat(prefix, df.shape[0])
|
|
else:
|
|
df.columns = [prefix + "_" + c if 'date' not in c else c for c in df.columns]
|
|
return df
|
|
|
|
|
|
def __apply_feature_extractors(df: pd.DataFrame,
|
|
log_returns: bool,
|
|
feature_extractors: list[tuple[str, FeatureExtractor, list[int]]]) -> pd.DataFrame:
|
|
|
|
for name, extractor, periods in feature_extractors:
|
|
for period in periods:
|
|
features = extractor(df, period, log_returns)
|
|
if type(features) == pd.DataFrame:
|
|
df = pd.concat([df, features], axis=1)
|
|
elif type(features) == pd.Series:
|
|
df[name + '_' + str(period)] = features
|
|
else:
|
|
assert False, "Feature extractor must return a pd.DataFrame or pd.Series"
|
|
return df
|
|
|
|
|
|
def __create_target_cum_forward_returns(df: pd.DataFrame, source_column: str, period: int) -> pd.Series:
|
|
assert period > 0
|
|
return df[source_column].shift(-period)
|
|
|
|
|
|
def __create_target_classes(df: pd.DataFrame, source_column: str, period: int, no_of_classes: Literal["two", "three"]) -> pd.Series:
|
|
assert period > 0
|
|
|
|
def get_class_binary(x: float) -> int:
|
|
return -1 if x <= 0.0 else 1
|
|
|
|
def get_class_threeway_balanced(series: pd.Series) -> pd.Series:
|
|
|
|
def get_bins_threeway(x):
|
|
bins = pd.qcut(df[source_column], 3, retbins=True, duplicates = 'drop')[1]
|
|
|
|
if len(bins) != 4:
|
|
# if we don't have enough data for the quantiles, we'll need to add hard-coded values
|
|
lower_bound = bins[0]
|
|
upper_bound = bins[-1]
|
|
bins = [lower_bound] + [-0.02, 0.02] + [upper_bound]
|
|
return bins
|
|
bins = get_bins_threeway(series)
|
|
|
|
def map_class_threeway(current_value):
|
|
lower_threshold = bins[1]
|
|
upper_threshold = bins[2]
|
|
if current_value <= lower_threshold:
|
|
return -1
|
|
elif current_value > lower_threshold and current_value < upper_threshold:
|
|
return 0
|
|
else:
|
|
return 1
|
|
return series.map(map_class_threeway)
|
|
|
|
def get_class_threeway_imbalanced(series: pd.Series) -> pd.Series:
|
|
|
|
def get_bins_threeway(x):
|
|
bins = pd.qcut(df[source_column], 4, retbins=True, duplicates = 'drop')[1]
|
|
|
|
if len(bins) != 5:
|
|
# if we don't have enough data for the quantiles, we'll need to add hard-coded values
|
|
lower_bound = bins[0]
|
|
upper_bound = bins[-1]
|
|
bins = [lower_bound] + [-0.02, 0.0, 0.02] + [upper_bound]
|
|
return bins
|
|
bins = get_bins_threeway(series)
|
|
|
|
def map_class_threeway(current_value):
|
|
lower_threshold = bins[1]
|
|
upper_threshold = bins[3]
|
|
if current_value <= lower_threshold:
|
|
return -1
|
|
elif current_value > lower_threshold and current_value < upper_threshold:
|
|
return 0
|
|
else:
|
|
return 1
|
|
return series.map(map_class_threeway)
|
|
|
|
target_column = df[source_column].shift(-period)
|
|
|
|
if no_of_classes == "three-balanced":
|
|
return get_class_threeway_balanced(target_column)
|
|
elif no_of_classes == "three-imbalanced":
|
|
return get_class_threeway_imbalanced(target_column)
|
|
else:
|
|
return target_column.map(get_class_binary)
|
|
|
|
|
|
def datasource_to_file(data_source: DataSource) -> str:
|
|
return data_source[0] + '/' + data_source[1] + '.csv'
|
|
|
|
|
|
|
|
def load_only_returns(assets: DataCollection, index_column: Literal['date', 'int'], returns: Literal['price', 'returns']) -> pd.DataFrame:
|
|
|
|
assets_future = [__load_df.remote(
|
|
data_source=data_source,
|
|
prefix=data_source[1],
|
|
returns=returns,
|
|
feature_extractors=[],
|
|
narrow_format=False,
|
|
) for data_source in assets]
|
|
target_asset_df = ray.get(assets_future)
|
|
|
|
dfs = [deduplicate_indexes(df) for df in target_asset_df]
|
|
dfs = pd.concat(dfs, axis=1)
|
|
# dfs = dfs.applymap(lambda x: np.nan if x == 0 else x)
|
|
dfs.index = pd.DatetimeIndex(dfs.index)
|
|
|
|
if index_column == 'int':
|
|
dfs.reset_index(drop=True, inplace=True)
|
|
|
|
return dfs
|