深入浅出特征工程|深入浅出特征工程 -- 基于 OpenMLDB 的实践指南(下)
1. 准备知识
在上期系列文章中(深入浅出特征工程 -- 基于 OpenMLDB 的实践指南(上)),我们介绍了特征工程的基础概念、实践工具,以及基本的基于单表的特征脚本开发。在本篇文章中,我们将基于主表和副表,去展开详细介绍更加复杂和强大的基于多表的特征脚本开发。同时,我们依然依托 OpenMLDB 所提供的 SQL 语法进行特征工程脚本示例,关于 OpenMLDB 的更多信息可以访问我们的 GitHub repo,以及文档网站。
OpenMLDB GitHub repo
OpenMLDB 文档 · 看云
如果你想运行本篇教程所举例的 SQL,请按照以下两个步骤做准备:
- 推荐使用 OpenMLDB docker 镜像在单机版下运行本教程,运行方式参考 OpenMLDB 快速上手。如果使用集群版,请使用离线模式(SET @@execute_mode='offline' )。集群版普通线上模式仅支持简单的数据预览功能,因此无法运行教程中大部分的 SQL。
- 本教程相关的所有数据以及导入操作脚本可以在这里下载。
表一:主表,用户交易表 t1
文章图片
副表:表二,商户流水表 t2
文章图片
在传统关系数据库中,为了取得多表的信息,最常用的方式是使用 join 进行拼接。但是对于特征工程的需求来说,数据库的 join 并不能非常高效的满足需求。最主要的原因是我们的主表样本表有一个用于模型训练的 label 列,其每一个值只能对应一行数据记录。所以实际中我们希望在 join 以后,结果表格的行数需要和主表的行数保持一致。
2. 副表单行特征 2.1 LAST JOIN OpenMLDB 目前支持LAST JOIN来进行类似数据库的 join 操作。LAST JOIN 可以看作一种特殊的 LEFT JOIN。在满足 JOIN 条件的前提下,左表的每一行拼取一条符合条件的最后一行。LAST JOIN分为无序拼接,和有序拼接。我们用更简单的表为例,假设表 s1,s2 的 schema 均为
(id int, col1 string, std_ts timestamp)
那么,我们可以做这样的join操作:
-- des c: 基于 ORDER BY 的有序 LAST JOIN 拼接
SELECT * FROM s1 LAST JOIN s2 ORDER BY s2.std_ts ON s1.col1 = s2.col1;
如下图,
LAST JOIN
时配置 ORDER BY
,则右表按 std_ts
排序,拼接最后一条命中的数据行。以左表第二行为例,符合条件的右表有 2 条,按 std_ts
排序后,选择最后一条 3, b, 2020-05-20 10:11:13
。文章图片
? LAST JOIN 过程示例
文章图片
? LAST JOIN 拼接结果
3. 副表多行聚合特征 【深入浅出特征工程|深入浅出特征工程 -- 基于 OpenMLDB 的实践指南(下)】OpenMLDB 针对副表拼接场景,扩展了标准的 WINDOW 语法,新增了 WINDOW UNION 的特性,支持从副表拼接多条数据形成副表窗口。在副表拼接窗口的基础上,可以方便构建副表多行聚合特征。同样地,构造副表多行聚合特征也需要完成两个步骤:
- 步骤一:定义副表拼接窗口。
- 步骤二:在副表拼接窗口上构造副表多行聚合特征。
副表拼接窗口的语法定义为:
window window_name as (UNION other_table PARTITION BY key_col ORDER BY order_col ROWS_RANGE|ROWS BETWEEN StartFrameBound AND EndFrameBound)
其中,最基本的不可或缺的语法元素包括:
UNION other_table
:other_table
是指进行 WINDOW UNION 的副表。
主表和副表需要保持schema一致。大部分情况下,主表和副表的schema都是不同的。因此,我们可以通过对主表和副表进行列筛选和默认列配置来保证参与窗口计算的主表和副表schema一致。列筛选还可以去掉无用列,只在关键列上做
WINDOW UNION 和聚合。PARTITION BY key_col
: 表示按列key_col
从副表拼接匹配数据。ORDER BY order_col
:表示副表拼接数据集按照order_col
列进行排序- `ROWS_RANGE BETWEEN StartFrameBound AND
EndFrameBound`: 表示副表拼接窗口的时间区间 StartFrameBound
表示该窗口的上界
UNBOUNDED PRECEDING
: 无上界。time_expression PRECEDING
: 如果是时间区间,可以定义时间偏移,如30d preceding
表示窗口上界为当前行的时间-30天。
EndFrameBound
表示该时间窗口的下界。
CURRENT ROW
: 当前行 time_expressionPRECEDING
: 如果是时间区间,可以定义时间偏移,如1d PRECEDING。这表示窗口下界为当前行的时间-1天。
ROWS BETWEEN StartFrameBound AND EndFrameBound
: 表示副表拼接窗口的时间区间StartFrameBound
表示该窗口的上界。
UNBOUNDED PRECEDING
: 无上界。number PRECEDING
: 如果是条数区间,可以定义时间条数。如,100 PRECEDING
表示窗口上界为的当前行的前100行。
EndFrameBound
表示该时间窗口的下界。
CURRENT ROW
: 当前行number PRECEDING
: 如果是条数窗口,可以定义时间条数。如,`1
- 配置窗口区间界时,请注意:
- OpenMLDB 目前无法支持当前行以后的时间作为上界和下界。如`1d
FOLLOWING`。换言之,我们只能处理历史时间窗口。这也基本满足大部分的特征工程的应用场景。 - OpenMLDB的下界时间必须>=上界时间
- OpenMLDB 的下界的条数必须<=上界条数
- OpenMLDB 目前无法支持当前行以后的时间作为上界和下界。如`1d
INSTANCE_NOT_IN_WINDOW
:标记为副表拼接窗口。主表除了当前行以外,其他数据不进入窗口。
以下通过具体例子来描述 WINDOW UNION 的拼接窗口定义操作。对于前面所述为用户交易表 t1,我们需要定义商户流水表 t2 的副表上拼接窗口,该拼接是基于
mid
进行。由于 t1 和 t2 的schema不同,所以我们首先分别从 t1 和 t2 抽取相同的列,对于不存在的列,可以配置缺省值。其中, mid
列用于两个表的拼接,所以是必须的;其次,时间戳的列(t1 中的 trans_time
,t2 中的 purchase_time
)包含时序信息,在定义时间窗口时候也是必须的;其余列按照聚合函数需要,进行必要的筛选保留。 以下 SQL 和示意图为从 t1 抽取必要列,生成 t11。
(select id, mid, trans_time as purchase_time, 0.0 as purchase_amt, "" as purchase_type from t1) as t11
文章图片
? 从 t1 抽取必要列,生成 t11
以下 SQL 和示意图为从 t2 抽取必要列,生成 t22。
(select 0L as id, mid, purchase_time, purchase_amt, purchase_type from t2) as t22
文章图片
? 从 t2 抽取必要列,生成 t22
可以看到,分别完成抽取以后生成的表格 t11 和 t22,已经具有了相同的 schema,两者可以进行逻辑上的 UNION 操作。但是在 OpenMLDB 中,WINDOW UNION 并不是真的为了进行传统数据库中的 UNION 操作,而是为了对于 t11 中的每一个样本行,去构建副表 t22 上的时间窗口。我们按照商户ID
mid
,对 t11 中的每一行数据,从 t22 中获取对应的拼接数据,然后按消费时间(purchase_time
) 排序,构造副表拼接窗口。比如我们定义一个 w_t2_10d
的窗口:不包含主表除了当前行以外的数据行,加上副表通过 mid 拼接上的十天以内的数据,示意图如下所示。可以看到,黄色和蓝色阴影部分,分别定义了样本 6 和样本 9 的副表拼接窗口。?
文章图片
? 副表拼接窗口示例
该窗口定义过程的 SQL 脚本如下所示(注意,这还不是一个完整的 SQL):
(SELECT id, mid, trans_time as purchase_time, 0.0 as purchase_amt, "" as purchage_type FROM t1) as t11
window w_t2_10d as (
UNION (SELECT 0L as id, mid, purchase_time, purchase_amt, purchase_type FROM t2) as t22
PARTITION BY mid ORDER BY purchase_time
ROWS_RANGE BETWEEN 10d PRECEDING AND 1 PRECEDING INSTANCE_NOT_IN_WINDOW)
3.2 步骤二:构建副表多行聚合特征 对于副表拼接窗口进行多行聚合函数加工,构造多行副表聚合特征,使得最后生成的行数和主表相同。以简单聚合函数为例,我们可以构造样本的副表拼接特征:商户的最近10天的零售总额
w10d_merchant_purchase_amt_sum
,商户的最近10天消费总次数 w10d_merchant_purchase_count
。以下 SQL 基于上面 3.1 中所定义的副表拼接窗口,构建多行聚合特征。SELECT
id,
-- 样本商户最近10天零售总额
sum(purchase_amt) over w_t2_10d as w10d_merchant_purchase_amt_sum,
-- 样本商户最近10天交易次数
count(purchase_amt) over w_t2_10d as w10d_merchant_purchase_count
FROM
(SELECT id, mid, trans_time as purchase_time, 0.0 as purchase_amt, "" as purchase_type FROM t1) as t11
window w_t2_10d as (
UNION (SELECT 0L as id, mid, purchase_time, purchase_amt, purchase_type FROM t2) as t22
PARTITION BY mid ORDER BY purchase_time
ROWS_RANGE BETWEEN 10d PRECEDING AND 1 PRECEDING INSTANCE_NOT_IN_WINDOW)
4. 特征组构建 一般而言,一个完整特征抽取脚本将抽取几十、上百,甚至几百个特征。我们可以根据特征类型、特征关联的表和窗口将这些特征分成若干组,然后将每一组特征放置到不同的SQL子查询里;最后将这些子查询按主表ID拼接在一起。本节,我们将承接前面的例子,演示如果将各种特征拼接在一起形成一个特征大宽表。
首先,我们将特征分成3组:
文章图片
然后,我们使用OpenMLDB SQL将同一组特征构建在同一个子查询中:
- 特征组 1
-- 主表单行特征
SELECT
-- 列直取
id, uid, trans_type,
-- 单行时间特征:day of week
dayofweek(trans_time) as f_trans_day_of_week,
-- 单行时间特征:交易日
day(trans_time) as f_trans_day,
-- 单行时间特征:交易小时
hour(trans_time) as f_trans_hour,
-- 单行时间特征:交易分钟
minute(trans_time),
-- 单行数学特征:交易额向上取整后取对数
log(ceiling(trans_amt)) as f_trans_amt_log,
-- 单行字符串特征:卡号前四位
substr(cardno, 1, 4),
-- 副表特征
t2.purchase_time as f_purchase_time,
t2.purchase_amt as f_purchase_amt,
t2.purchase_type as f_purchase_typeFROM t1 LAST JOIN t2 ORDER BY t2.purchase_time ON t1.mid = t2.mid
- 特征组 2
-- 主表窗口特征
SELECT
id as out2id,
-- 最近30天POS交易总金额
sum_where(trans_amt, trans_type = "POS") over w30d as w30d_sum_pos_trans_amt,
-- 最近30天的最大POS交易金额
max_where(trans_amt, trans_type = "POS") over w30d as w30d_max_pos_trans_amt,
-- 最近30天的单次平均POS交易金额
avg_where(trans_amt, trans_type = "POS") over w30d as w30d_avg_pos_trans_amt,
-- 最近30天的POS交易总次数
count_where(trans_amt, trans_type = "POS") over w30d as w30d_count_pos_trans_amt,-- 最近一周总交易金额总额
sum(trans_amt) over w7d as w7d_sum_trans_amt,
-- 最近一周的总交易次数
count(trans_amt) over w7d as w7d_count_trans_amt,
from t1
-- 最近30天的窗口
window w30d as (PARTITION BY uid ORDER BY trans_time ROWS_RANGE BETWEEN 30d PRECEDING AND CURRENT ROW),
-- 最近一周的窗口
w7d as (PARTITION BY uid ORDER BY trans_time ROWS_RANGE BETWEEN 7d PRECEDING AND CURRENT ROW)
- 特征组 3
-- 副表聚合特征
SELECT
id as out3id,
-- 样本商户最近一周零售总额
sum(purchase_amt) over w7d_merchant as w7d_merchant_purchase_amt_sum,
-- 样本商户最近一周交易次数
count(purchase_amt) over w7d_merchant as w7d_merchant_purchase_count,
-- 主表商户最近一周的流水
FROM (select id, mid, cardno as card, trans_time as purchase_time, 0.0 as purchase_amt, "" as purchase_type from t1) as t11
window w7d_merchant as (UNION (select 0L as id, mid, card, purchase_time, purchase_amt, purchase_type from t2) as t22 PARTITION BY mid ORDER BY purchase_time ROWS_RANGE BETWEEN 30d PRECEDING AND 1 PRECEDING INSTANCE_NOT_IN_WINDOW)
最后,将这三组特征,按主表ID拼接在一起:
SELECT * FROM
-- 特征组1
(
-- 主表单行特征
SELECT
-- 列直取
id, uid, trans_type,
-- 单行时间特征:day of week
dayofweek(trans_time) as f_trans_day_of_week,
-- 单行时间特征:交易日
day(trans_time) as f_trans_day,
-- 单行时间特征:交易小时
hour(trans_time) as f_trans_hour,
-- 单行时间特征:交易分钟
minute(trans_time),
-- 单行数学特征:交易额向上取整后取对数
log(ceiling(trans_amt)) as f_trans_amt_log,
-- 单行字符串特征:卡号前四位
substr(cardno, 1, 4),-- 副表特征
t2.purchase_time as f_purchase_time,
t2.purchase_amt as f_purchase_amt,
t2.purchase_type as f_purchase_typeFROM t1 LAST JOIN t2 ORDER BY t2.purchase_time ON t1.mid = t2.mid
) as out1 LAST JOIN
-- 特征组2
(
-- 主表窗口特征
SELECT
id as out2id,
-- 最近30天POS交易总金额
sum_where(trans_amt, trans_type = "POS") over w30d as w30d_sum_pos_trans_amt,
-- 最近30天的最大POS交易金额
max_where(trans_amt, trans_type = "POS") over w30d as w30d_max_pos_trans_amt,
-- 最近30天的单次平均POS交易金额
avg_where(trans_amt, trans_type = "POS") over w30d as w30d_avg_pos_trans_amt,
-- 最近30天的POS交易总次数
count_where(trans_amt, trans_type = "POS") over w30d as w30d_count_pos_trans_amt,-- 最近一周总交易金额总额
sum(trans_amt) over w7d as w7d_sum_trans_amt,
-- 最近一周的总交易次数
count(trans_amt) over w7d as w7d_count_trans_amt
from t1
-- 最近30天的窗口
window w30d as (PARTITION BY uid ORDER BY trans_time ROWS_RANGE BETWEEN 30d PRECEDING AND CURRENT ROW),
-- 最近一周的窗口
w7d as (PARTITION BY uid ORDER BY trans_time ROWS_RANGE BETWEEN 7d PRECEDING AND CURRENT ROW)
) as out2 ON out1.id = out2.out2id LAST JOIN
-- 特征组3
(
-- 副表聚合特征
SELECT
id as out3id,
-- 样本商户最近一周零售总额
sum(purchase_amt) over w7d_merchant as w7d_merchant_purchase_amt_sum,
-- 样本商户最近一周交易次数
count(purchase_amt) over w7d_merchant as w7d_merchant_purchase_count
-- 主表商户最近一周的流水
FROM (select id, mid, cardno as card, trans_time as purchase_time, 0.0 as purchase_amt, "" as purchase_type from t1) as t11
window w7d_merchant as (UNION (select 0L as id, mid, card, purchase_time, purchase_amt, purchase_type from t2) as t22 PARTITION BY mid ORDER BY purchase_time ROWS_RANGE BETWEEN 30d PRECEDING AND 1 PRECEDING INSTANCE_NOT_IN_WINDOW)
) as out3 ON out1.id = out3.out3id;
5. 扩展阅读 开源机器学习数据库OpenMLDB v0.4.0产品介绍
想快速试用 OpenMLDB 来开始写特征计算脚本?快来看一下 OpenMLDB 快速上手
完整的 SQL 语法参考文档:中国镜像站点,国际站点
本系列教程上篇:4PD开发者社区:深入浅出特征工程 – 基于 OpenMLDB 的实践指南(上)
推荐阅读
- 【Java分享客栈】我有一个朋友,和前端工程师联调接口被狠狠鄙视了一番。
- java|温故知新-什么是软件工程
- 深入浅出归并排序(一)
- 图像特征描述|边缘检测之Sobel检测算子
- 5年工作经验的Android工程师面试经验总结,阿里+腾讯+字节跳动面试真题汇总
- 计算机四级网络工程——操作系统——第一章操作系统概述
- Unity3D|Unity Addressables 深入浅出(一)
- Unity3D|【游戏开发高阶】从零到一教你Unity使用ToLua实现热更新(含Demo工程 | LuaFramework | 增量 | HotUpdate)
- Unity3D|【游戏开发实战】Unity UGUI实现循环复用列表,显示巨量列表信息,含Demo工程源码
- 芯片工程师太贵(贵你妹啊!)