Score and Predict Large Datasets¶
Sometimes you’ll train on a smaller dataset that fits in memory, but need to predict or score for a much larger (possibly larger than memory) dataset. Perhaps your learning curve has leveled off, or you only have labels for a subset of the data.
In this situation, you can use ParallelPostFit to parallelize and distribute the scoring or prediction steps.
from dask.distributed import Client, progress # Scale up: connect to your own cluster with bmore resources # see http://dask.pydata.org/en/latest/setup.html client = Client(processes=False, threads_per_worker=4, n_workers=1, memory_limit='2GB') client
import numpy as np import dask.array as da from sklearn.datasets import make_classification
We’ll generate a small random dataset with scikit-learn.
X_train, y_train = make_classification( n_features=2, n_redundant=0, n_informative=2, random_state=1, n_clusters_per_class=1, n_samples=1000) X_train[:5]
array([[ 1.53682958, -1.39869399], [ 1.36917601, -0.63734411], [ 0.50231787, -0.45910529], [ 1.83319262, -1.29808229], [ 1.04235568, 1.12152929]])
And we’ll clone that dataset many times with
y_large represent our larger than memory dataset.
# Scale up: increase N, the number of times we replicate the data. N = 100 X_large = da.concatenate([da.from_array(X_train, chunks=X_train.shape) for _ in range(N)]) y_large = da.concatenate([da.from_array(y_train, chunks=y_train.shape) for _ in range(N)]) X_large
dask.array<concatenate, shape=(100000, 2), dtype=float64, chunksize=(1000, 2)>
Since our training dataset fits in memory, we can use a scikit-learn
estimator as the actual estimator fit during traning. But we know that
we’ll want to predict for a large dataset, so we’ll wrap the
scikit-learn estimator with
from sklearn.linear_model import LogisticRegressionCV from dask_ml.wrappers import ParallelPostFit
clf = ParallelPostFit(LogisticRegressionCV(cv=3))
Now we’ll call
clf.fit. Dask-ML does nothing here, so this step can
only use datasets that fit in memory.
ParallelPostFit(estimator=LogisticRegressionCV(Cs=10, class_weight=None, cv=3, dual=False, fit_intercept=True, intercept_scaling=1.0, max_iter=100, multi_class='warn', n_jobs=None, penalty='l2', random_state=None, refit=True, scoring=None, solver='lbfgs', tol=0.0001, verbose=0), scoring=None)
Now that training is done, we’ll turn to predicting for the full (larger than memory) dataset.
y_pred = clf.predict(X_large) y_pred
dask.array<predict, shape=(100000,), dtype=int64, chunksize=(1000,)>
y_pred is Dask arary. Workers can write the predicted values to a shared file system, without ever having to collect the data on a single machine.
Or we can check the models score on the entire large dataset. The computation will be done in parallel, and no single machine will have to hold all the data.