Przeszukiwanie siatki

Uczenie maszynowe z PySpark

Andrew Collier

Data Scientist, Fathom Data

Wybór optymalnej wartości parametru

Samochody ponownie

cars.select('mass', 'cyl', 'consumption').show(5)

+------+---+-----------+
|  mass|cyl|consumption|
+------+---+-----------+
|1451.0|  6|       9.05|
|1129.0|  4|       6.53|
|1399.0|  4|       7.84|
|1147.0|  4|       7.84|
|1111.0|  4|       9.05|
+------+---+-----------+

Zużycie paliwa z wyrazem wolnym

Regresja liniowa z wyrazem wolnym. Dopasowanie do danych treningowych.

regression = LinearRegression(labelCol='consumption', fitIntercept=True)
regression = regression.fit(cars_train)

Oblicz RMSE na danych testowych.

evaluator.evaluate(regression.transform(cars_test))

# RMSE for model with an intercept
0.745974203928479

Zużycie paliwa bez wyrazu wolnego

Regresja liniowa bez wyrazu wolnego. Dopasowanie do danych treningowych.

regression = LinearRegression(labelCol='consumption', fitIntercept=False)
regression = regression.fit(cars_train)

Oblicz RMSE na danych testowych.

# RMSE for model without an intercept (second model)
0.852819012439

# RMSE for model with an intercept    (first model)
0.745974203928

Siatka parametrów

from pyspark.ml.tuning import ParamGridBuilder

# Create a parameter grid builder
params = ParamGridBuilder()

# Add grid points
params = params.addGrid(regression.fitIntercept, [True, False])

# Construct the grid
params = params.build()


# How many models?
print('Number of models to be tested: ', len(params))

Number of models to be tested:  2

Przeszukiwanie siatki z walidacją krzyżową

Utwórz walidator krzyżowy i dopasuj do danych treningowych.

cv = CrossValidator(estimator=regression,
                    estimatorParamMaps=params,
                    evaluator=evaluator)
cv = cv.setNumFolds(10).setSeed(13).fit(cars_train)

Jaka jest wartość RMSE po walidacji krzyżowej dla każdego modelu?

cv.avgMetrics

[0.800663722151, 0.907977823182]

Najlepszy model i parametry

# Access the best model
cv.bestModel

Można też użyć obiektu walidatora krzyżowego.

predictions = cv.transform(cars_test)

Pobierz najlepszy parametr.

cv.bestModel.explainParam('fitIntercept')

'fitIntercept: whether to fit an intercept term (default: True, current: True)'

Bardziej rozbudowana siatka

params = ParamGridBuilder() \
            .addGrid(regression.fitIntercept, [True, False]) \

            .addGrid(regression.regParam, [0.001, 0.01, 0.1, 1, 10]) \

            .addGrid(regression.elasticNetParam, [0, 0.25, 0.5, 0.75, 1]) \
            .build()

Ile modeli teraz?

print ('Number of models to be tested: ', len(params))

Number of models to be tested:  50

Znajdź najlepsze parametry!

Uczenie maszynowe z PySpark