-
Notifications
You must be signed in to change notification settings - Fork 25
/
shap_example.py
70 lines (55 loc) · 2.42 KB
/
shap_example.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
# Copyright 2023 Neal Lathia
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
import shap
from libraries.util.datasets import load_regression_dataset
from libraries.util.domains import DIABETES_DOMAIN
from sklearn.ensemble import GradientBoostingRegressor
from sklearn.metrics import mean_squared_error
from sklearn.pipeline import Pipeline
from modelstore.model_store import ModelStore
EXPLAINER_DOMAIN = f"{DIABETES_DOMAIN}-explainer"
def _train_example_model() -> Pipeline:
X_train, X_test, y_train, y_test = load_regression_dataset()
# Train a model using an sklearn pipeline
params = {
"n_estimators": 250,
"max_depth": 4,
"min_samples_split": 5,
"learning_rate": 0.01,
"loss": "ls",
}
model = GradientBoostingRegressor(**params)
model.fit(X_train, y_train)
results = mean_squared_error(y_test, model.predict(X_test))
print(f"🔍 Trained model MSE={results}.")
explainer = shap.TreeExplainer(model)
# Example only
shap_values = explainer.shap_values(X_test)[0]
print(f"🔍 Shap values={shap_values[:10]}.")
return explainer
def train_and_upload(modelstore: ModelStore) -> dict:
# Train a model and return the explainer
explainer = _train_example_model()
# Upload the explainer to the model store
print(f'⤴️ Uploading the explainer to the "{EXPLAINER_DOMAIN}" domain.')
meta_data = modelstore.upload(EXPLAINER_DOMAIN, explainer=explainer)
return meta_data
def load_and_test(modelstore: ModelStore, model_domain: str, model_id: str):
# Load the explainer back into memory!
print(f'⤵️ Loading the explainer "{model_domain}" domain model={model_id}')
explainer = modelstore.load(model_domain, model_id)
# Run some example predictions
_, X_test, _, _ = load_regression_dataset()
shap_values = explainer.shap_values(X_test)[0]
print(f"🔍 Shap values={shap_values[:10]}.")