++++Notebook converted from Jupyter for blog publishing.
02-KNN-Exercise-Solutions
KNN Project Exercise - Solutions
Due to the simplicity of KNN for Classification, let's focus on using a PipeLine and a GridSearchCV tool, since these skills can be generalized for any model.
The Sonar Data
Detecting a Rock or a Mine
Sonar (sound navigation ranging) is a technique that uses sound propagation (usually underwater, as in submarine navigation) to navigate, communicate with or detect objects on or under the surface of the water, such as other vessels.
The data set contains the response metrics for 60 separate sonar frequencies sent out against a known mine field (and known rocks). These frequencies are then labeled with the known object they were beaming the sound at (either a rock or a mine).
Our main goal is to create a machine learning model capable of detecting the difference between a rock or a mine based on the response of the 60 separate sonar frequencies.
Data Source: https://archive.ics.uci.edu/ml/datasets/Connectionist+Bench+(Sonar,+Mines+vs.+Rocks) (opens in a new tab)
Complete the Tasks in bold
TASK: Run the cells below to load the data.
import numpy as np
import pandas as pd
import seaborn as sns
import matplotlib.pyplot as pltdf = pd.read_csv('../DATA/sonar.all-data.csv')df.head()Freq_1
Freq_2
Freq_3
Freq_4
Freq_5Data Exploration
TASK: Create a heatmap of the correlation between the difference frequency responses.
# CODE HEREplt.figure(figsize=(8,6))
sns.heatmap(df.corr(),cmap='coolwarm')<AxesSubplot:>
TASK: What are the top 5 correlated frequencies with the target\label?
Note: You many need to map the label to 0s and 1s.
Additional Note: We're looking for absolute correlation values.
#CODE HEREdf['Target'] = df['Label'].map({'R':0,'M':1})np.abs(df.corr()['Target']).sort_values().tail(6)Freq_45 0.339406
Freq_10 0.341142
Freq_49 0.351312
Freq_12 0.392245
Freq_11 0.432855Train | Test Split
Our approach here will be one of using Cross Validation on 90% of the dataset, and then judging our results on a final test set of 10% to evaluate our model.
TASK: Split the data into features and labels, and then split into a training set and test set, with 90% for Cross-Validation training, and 10% for a final test set.
# CODE HEREfrom sklearn.model_selection import train_test_splitX = df.drop(['Target','Label'],axis=1)
y = df['Label']X_cv, X_test, y_cv, y_test = train_test_split(X, y, test_size=0.1, random_state=42)TASK: Create a PipeLine that contains both a StandardScaler and a KNN model
from sklearn.preprocessing import StandardScaler
from sklearn.neighbors import KNeighborsClassifierscaler = StandardScaler()knn = KNeighborsClassifier()operations = [('scaler',scaler),('knn',knn)]from sklearn.pipeline import Pipelinepipe = Pipeline(operations)TASK: Perform a grid-search with the pipeline to test various values of k and report back the best performing parameters.
from sklearn.model_selection import GridSearchCVk_values = list(range(1,30))param_grid = {'knn__n_neighbors': k_values}full_cv_classifier = GridSearchCV(pipe,param_grid,cv=5,scoring='accuracy')full_cv_classifier.fit(X_cv,y_cv)GridSearchCV(cv=5,
estimator=Pipeline(steps=[('scaler', StandardScaler()),
('knn', KNeighborsClassifier())]),
param_grid={'knn__n_neighbors': [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,
12, 13, 14, 15, 16, 17, 18, 19,full_cv_classifier.best_estimator_.get_params(){'memory': None,
'steps': [('scaler', StandardScaler()),
('knn', KNeighborsClassifier(n_neighbors=1))],
'verbose': False,
'scaler': StandardScaler(),(HARD) TASK: Using the .cv_results_ dictionary, see if you can create a plot of the mean test scores per K value.
#CODE HEREfull_cv_classifier.cv_results_['mean_test_score']array([0.84537696, 0.78065434, 0.77524893, 0.75917496, 0.75931721,
0.74822191, 0.75945946, 0.71664296, 0.7113798 , 0.68421053,
0.70042674, 0.68435277, 0.68449502, 0.67908962, 0.69530583,
0.68990043, 0.7113798 , 0.70042674, 0.72204836, 0.67908962,
0.70071124, 0.69530583, 0.69530583, 0.68463727, 0.68477952,scores = full_cv_classifier.cv_results_['mean_test_score']
plt.plot(k_values,scores,'o-')
plt.xlabel("K")
plt.ylabel("Accuracy")Text(0, 0.5, 'Accuracy')
Final Model Evaluation
TASK: Using the grid classifier object from the previous step, get a final performance classification report and confusion matrix.
#Code Herepred = full_cv_classifier.predict(X_test)from sklearn.metrics import classification_report,confusion_matrix,accuracy_scoreconfusion_matrix(y_test,pred)array([[12, 1],
[ 1, 7]], dtype=int64)print(classification_report(y_test,pred)) precision recall f1-score support
M 0.92 0.92 0.92 13
R 0.88 0.88 0.88 8