と。

統計学は趣味、マーケティングは義務。

pipelineライブラリKedroを使ってみた話

02_機械学習

「パイプラインって何？」

仕事でも機械学習の案件がちょっと増えてきたというのと、
kaggleもベースラインくらいは自動的にsubmitできるところまで持っていきたいって思ったので、
pipelineを作ろうと言うことになりました。

ただ、私はエンジニアリング畑ではないので、ゼロから作れる自信がありません。
困ったなー困ったなーと思っていたところに、こんなQiitaを見かけました。

qiita.com

なるほど、いろいろあるんだな、となりました。
この中から今回はKedroを導入しとりあえず触ってみたのでレポします。

【2020年2月9日追記】

この記事を書いたら下のQiitaを書いたMinamiさんがこんなことをおっしゃいまして

「そんな褒め過ぎですよ」とか返そうとか思ったら、
KedroのProduct ManagerのYetundeさんがこんなことをおっしゃいました。

「リソースガイドとしてKedroのFAQに載せるよ」って書いてる……！！！？？？！？！？！！？

私は技術力がなく、この記事も個人的なメモ扱いだったのですが、
開発に携わっている人たちにこうして受け入れてもらえているのは素直に光栄に思います。

私も私で、kedroはいろいろトライアンドエラーしようと思っています。
ということで以下にもREADMEとして記事を転載しておきます。 多分フォルダが増えたりREADMEが加筆修正されたりすると思います。

github.com

pipelineとは

機械学習のタスクをまとめ上げて、自動化するやつです(ざっくり)

データの前処理やパラメータのチューニングなど、
機械学習の運用はいろんな煩雑化要素が現れて悲しくなります。
このあたりを含めて整理し、結果の再現やメンテナンス性を
確保するための試み、という感じです。

Kedroとは

もうgithubに飛んだほうが早い github.com

Quantumblack社が手がけるオープンソースのpipelineライブラリです。

なんでKedroなの？

他のライブラリが悪い、というわけではなく、理由は主に3点です。

  • ドキュメントが丁寧

  • 開発に携わってる人が優しい

  • ロゴ・バナーがかっこいい(重要)

ドキュメントが丁寧

get startedはどのライブラリにもありますが、ステップバイステップで解説してくれており、
ファイルの構造やいつも触る機会のないyamlファイルなどの対応関係がとてもわかりやすかったです。
pipelineはおろか、Pythonも基本対話型で使っちゃうような私にとっては、この点はすごい重要でした。

開発者やコミュニティがニュービーに優しい

Twitterで「このやり方あとで調べよう」つったら
メインテインしてる人がリプくれたり、
開発してる人から反応もらえたり、
悩んでいたコーディングの実装例を示してくれたりして、
びっくりした反面めちゃくちゃ嬉しかったというのがあります。
多分他のライブラリを使ってても、誰かしら何かしら教えてくれるとは思うんですが、
こうしていろいろ教えてくれるんだと思うだけで結構心強かったりします。

ロゴ・バナーがかっこいい

このツイートが全てです*1。

twitter.com

バナーみたとき「あ、かっこいいし使おう」ってなっちゃいました。中二病なので。†でHN囲むタイプの。
†きぬいと†。
つらくなってきた。

使ってみた

ドキュメントがあるのでぶっちゃけここから先は蛇足でしかないんですが、使ってみたレポートです。
実装githubは以下。
github.com

Kedroはpipでインストールできます。

pip install kedro

使ったデータ

前ブログで使ったやつを改造したデータを使ってます。
乱数で作ってるので、学習済みモデルも公開して問題ないっちゃ問題ない。何にも使えないので。

プロジェクト立ち上げはkedro new

Rを使っていれば割と自然なプロジェクトという概念、Kedroにもあります*2。
コマンドラインにkedro newと入力すると対話的に「プロジェクト名」などを入力できます。
これが済むと、作業ディレクトリ上にkedroのプロジェクトフォルダができます。しゅごい。

実行準備はkedro install

プロジェクトを立ち上げたらkedro installと入力するとなんかわにゃわにゃ動きます。
動き終わると「データを置く場所」「前処理をする場所」「モデルを構築する場所」などが
大体作られます。すごい。

フォルダ構成

きぬいとのgitではこんな感じになっています。
__init__とか__pycache__とかもありますが省略しています。
フォルダ構成の時点できれいなのも好きです。
ただ、フォルダ量は膨大なので、各々の説明は公式ドキュメント等を参考にしてください。

kedro_classification
├── README.md
├── conf
│   ├── README.md
│   ├── base
│   │   ├── catalog.yml
│   │   ├── credentials.yml
│   │   ├── logging.yml
│   │   └── parameters.yml
│   └── local
├── data
│   ├── 01_raw
│   ├── 02_intermediate
│   ├── 03_primary
│   ├── 04_features
│   ├── 05_model_input
│   ├── 06_models
│   ├── 07_model_output
│   └── 08_reporting
├── docs
│   └── source
│       ├── conf.py
│       └── index.rst
├── kedro_cli.py
├── logs
│   ├── errors.log
│   ├── info.log
│   └── journals
├── notebooks
├── references
├── results
├── setup.cfg
└── src
    ├── kedro_classification
    │   ├── nodes
    │   ├── pipeline.py
    │   ├── pipelines
    │   │   ├── data_engineering
    │   │   │   ├── nodes.py
    │   │   │   └── pipeline.py
    │   │   └── data_science
    │   │       ├── nodes.py
    │   │       └── pipeline.py
    │   └── run.py
    ├── requirements.txt
    ├── setup.py
    └── tests
        └── test_run.py

pipeline構築

上記のフォルダ構成ですが、第一にdata/01_rawに使うデータを突っ込むというところをやります。
今回はダミーデータがぶち込まれています。 それ以外では、confにあるyamlファイルと、src/[プロジェクト名]/pipelines/data_engineering、
同じくdata_scienceのpythonスクリプト(nodes.py、pipeline.py)を編集していきます。 全部編集したら1階層上にもpipeline.pyがあるので、これも編集します。
この辺の名前はおそらく参照さえしっかりしていればわかりやすく名前を付けても良いと思います。

conf内のyamlファイルでは、中間テーブルの出力やモデルの出力をcatalog.yml、
モデルのパラメータなどをparameters.ymlに格納します。
これで、nodes.pyで引数として指定すれば、Kedro側でよしなに引っ張ってきてくれます。 nodes.pyには、具体的なデータの前処理やモデル関数の定義を書き、
pipeline.pyには、それらの入出力を定義します。

data_engineering

このフォルダでは主に前処理を行います。
今回はどちらかというとモデリングの部分のトライアルがメインだったので、
適当な前処理になっちゃってます。

import pandas as pd
import numpy as np


def preprocessing(usedata: pd.DataFrame) -> pd.DataFrame:
    for i in range(1,40):
        var_name = 'Var.' + str(i)
        usedata[var_name] = 1 - usedata[var_name]
    return usedata

もちろん複数関数を定義して、適宜適用することもできます。

そしてこれを、格納している元データに適用し、
適用結果を中間テーブルとして吐き出すためのスクリプトがpipeline.py

from kedro.pipeline import node, Pipeline
from kedro_classification.pipelines.data_engineering.nodes import preprocessing

def create_pipeline(**kwargs):
    print('loading create_pipeline in pipeline.py....')
    return Pipeline(
        [
            node(
                func=preprocessing,
                inputs='usedata',
                outputs='preprocessed_Data',
                name='preprocessed_Data',
            ),
        ]
    )

実装もnodes.pyの関数を持ってきて、
出力はpreprocessed_Dataとして出す、というシンプルなものです*3。
pipeline.pyはdata_scienceでも同様の記法で書きます。なので省略です。

data_science

公式ドキュメントでは線型回帰の実装があったので、背伸びしてLightGBMの実装をやってみました。
コード(一部)は以下。

def LightGBM_model(
    data: pd.DataFrame,
    parameters: Dict
    ) -> lgb.LGBMRegressor:
    
    ### define classes
    regressor = lgb.LGBMRegressor()
    y = data['y']
    X = data.drop(['y', 'ID'], axis=1)
    ### hyperparameters from parameters.yml
    lgb_params = {
            'num_iterations'       : parameters['n_estimators'],
            'boosting_type'        : parameters['boosting_type'],
            'objective'            : parameters['objective'],
            'metric'               : parameters['metric'],
            'num_leaves'           : parameters['num_leaves'],
            'learning_rate'        : parameters['learning_rate'],
            'max_depth'            : parameters['max_depth'],
            'verbosity'            : parameters['verbose'],
            'early_stopping_round' : parameters['early_stopping_rounds'],
            'seed'                 : parameters['seed']
            }

    fold = KFold(n_splits=parameters['folds'], random_state=parameters['random_state'])

    oof_pred = np.zeros(len(X))
    ### run model with kfold
    for k,   (train_index, valid_index) in enumerate(fold.split(X, y)):
        #print(train_index)
        X_train, X_valid = X.iloc[train_index], X.iloc[valid_index]
        y_train, y_valid = y.iloc[train_index], y.iloc[valid_index]
    
        lgb_train = lgb.Dataset(X_train, y_train)
        lgb_valid = lgb.Dataset(X_valid, y_valid)

        regressor = lgb.train(lgb_params, lgb_train, valid_sets=lgb_valid, verbose_eval=False)
        
        y_train_pred = regressor.predict(X_train, num_iteration=regressor.best_iteration)
        y_valid_pred = regressor.predict(X_valid, num_iteration=regressor.best_iteration)

        auc_train = roc_auc_score(y_train, y_train_pred)
        auc_valid = roc_auc_score(y_valid, y_valid_pred)
        print('Early stopping round is: {iter}'.format(iter=regressor.current_iteration()))
        print('Fold {n_folds}: train AUC is {train: .3f} valid AUC is {valid: .3f}'.format(n_folds=k+1, train=auc_train, valid=auc_valid))
    
    return regressor


def evaluate_LightGBM_model(regressor: lgb.basic.Booster, X_test: np.ndarray, y_test: np.ndarray): 
    y_pred = regressor.predict(X_test, num_iteration=regressor.best_iteration)
    print('y predicted!')
    print(type(y_pred)) 
    #y_pred = np.argmax(y_pred, axis=1)
    #roc_curve = r
    score  = roc_auc_score(y_test, y_pred)
    logger = logging.getLogger(__name__)
    logger.info('AUC is %.3f.', score)

汚い実装で恥ずかしいですが、こんな感じで書ければ大丈夫です。
parametersを入れておけば、対応するハイパーパラメータをKedroがよしなに持ってきてくれます。
pipeline.pyでこれらの関数と入出力を指定すればOK。

# coding: utf-8

from kedro.pipeline import node, Pipeline
from typing import Dict, Any, List
from kedro_classification.pipelines.data_science.nodes import (
    split_data,
    Linear_Regression_model,
    LightGBM_model,
    evaluate_LightGBM_model
)

def create_pipeline(**kwargs):
    return Pipeline(
        [
            node(
                func=split_data,
                inputs=['preprocessed_Data', 'parameters'],
                outputs=['X_train', 'X_test', 'y_train', 'y_test'],
            ),            
            node(
                func=LightGBM_model,
                inputs=['preprocessed_Data', 'parameters'],
                outputs='regressor',
                name='regressor',
            ),
            node(
                func=evaluate_LightGBM_model,
                inputs=['regressor', 'X_test', 'y_test'],
                outputs=None,
            ),
        ]
    )

最後に、src直下のpipeline.pyを編集します。
ここで前処理からモデル実行・出力に至るまでの流れを指定します。

from typing import Dict

from kedro.pipeline import Pipeline
from kedro_classification.pipelines.data_engineering import pipeline as de
from kedro_classification.pipelines.data_science import pipeline as ds

def create_pipelines(**kwargs) -> Dict[str, Pipeline]:
    de_pipeline = de.create_pipeline()
    ds_pipeline = ds.create_pipeline()

    return {
        'de' : de_pipeline,
        '__default__': de_pipeline + ds_pipeline
    }

kedro test

ここまでやったら(あるいはここまでに至るまでのどこかで)テストを実行できます。
コマンドラインでkedro testを実行するだけでどこで詰むか詰まないかがわかります。 すごい。

kedro run

テストが上手く行ったら走らせます。kedro runで走ります。
すごい。

おわりに

長くなりました。
ここまでの実装は、ドキュメントを見ながら試行錯誤し、3日位で無理なく進められました。
今回はローカルでの実装ですが、もちろんKedroはAWSやGCPでも動きますし、
様々なDBからデータを獲得できます。
データは分析が主で、パイプライン構築以前にこういうエンジニアリングの経験が浅い私にも
どうにか形にできるところまで道案内してくれるので、とても良いライブラリでした。

一回ベースラインができたとはいえ、まだまだ足りないことも多いです。
現状はLighGBMのモデル関数の中でCVの手法をねじ込んでいる部分の改善や、
複数の前処理関数を組合せたpipeline構築など、
脱Tutorialな課題はいっぱいあります。
ただ、これはKedro側の問題というよりはnodes.py内でどう実装するかの問題なので、
Kedroは割と柔軟にここを受け入れてくれると信じています(？)

*1:このあと「でもさっきバナー変えようってPR送っちゃった」って開発主担当の人が言ってくれたんだけどどのみちかっこいいからいい

*2:多分大体のライブラリにある

*3:中間テーブルはcatalog.ymlで管理されます