辅导CS4006 Adaboost Classifiers调试Python

Adaboost。

Requirement

In this assignment, your task is to train an AdaBoost classifier on synthetic data. For reference, you are provided with the posterior P(y = 1 | x), with x regularly sampled over the domain X = [0, 1] × [0, 1] , so that you can see, in the end, how the output of the AdaBoost classifier better approximates the posterior at each round.

Please read the assignment entirely before you start coding, in order to get a sense of how it is organized. In particular, note that the AdaBoost algorithm is only run at the very last cell of the “Train the classifier” section. Before that, a number of functions are defined, one of which you need to complete.

Question 1

Fill in the missing parts to implement the Adaboost algorithm described in class (slide 64 of the course). This involves iterating over the following steps:

a. Find the best weak learner h at each round.
b. Using the weak learner’s weighted error e, compute t.
c. Update the weight distribution D of the training samples.

Question 2

Modify your loop to compute the loss at each round. Then, plot E and make sure that it is monotonically decreasing with time. Verify that E provides an upper bound for the number of errors.

Question 3

First show the approximate posterior of your strong learner side-by-side with the original posterior. Then, show the approximate posteriors for each step at which the learner’s response has been saved. Make sure that they look increasingly similar to the original posterior.

A word on notation

The response of a weak learner h for the sample x is h(x) ∈ {−1, 1}.
At each round we find the best weak learner. The overall response of the strong learner at round t for the sample x. In order to be coherent with the weak learner’s expression, we can also define H(x) = sign(f(x)) ∈ {1, 1}, which can also be called the overall response. However, in this assignment, we are only interested in f .

Code

Imports

1
2
3

import matplotlib.pyplot as plt
import numpy as np
from construct_data import construct_data

Visualize training data and posterior

features, labels, posterior = construct_data(500, 'train', 'nonlinear', plusminus=True)

# Extract features for both classes
features_pos = features[labels == 1]
features_neg = features[labels != 1]

# Display data
fig = plt.figure(figsize=plt.figaspect(0.3))
ax = fig.add_subplot(1, 2, 1)
ax.scatter(features_pos[:, 0], features_pos[:, 1], c="red", label="Positive class")
ax.scatter(features_neg[:, 0], features_neg[:, 1], c="blue", label="Negative class")

ax.set_title("Training data")
ax.set_xlabel("Feature 1")
ax.set_ylabel("Feature 2")
ax.legend()

ax = fig.add_subplot(1, 2, 2)
ax.imshow(posterior, extent=[0, 1, 0, 1], origin='lower')
ax.set_title("Posterior of the positive class $P(y=1 \mid x)$")
ax.set_xlabel("Feature 1")
ax.set_ylabel("Feature 2")

plt.show()

Train the classifier

Weak learner evaluation

The weak learner we use for this classification problem is a decision stump (see slide 63 of the course), whose response is defined as h(x) = s(2[xd ≥ θ] − 1), where

d is the the dimension along which the decision is taken,
[·] is 1 if is true and 0 otherwise,
θ is the threshold applied along dimension d and
s ∈ {−1, 1} is the polarity of the decision stump (this is a multiplicative factor, not a func- tion!).

For example, if s = 1, the decision stump will consider that all samples whose d-th feature is greater than θ are in the positive class (h(x) = +1), and all samples with a feature strictly lower than are in the negative class (h(x) = -1).

def evaluate_stump(features, coordinate_wl, polarity_wl, theta_wl):
    """Evaluate the stump's response for each point."""
    feature_slice = features[:, coordinate_wl]
    weak_learner_output = polarity_wl * (2*(feature_slice gt;= theta\_wl) - 1)

    return weak_learner_output

def evaluate_stump_on_grid(x_rng, y_rng, coordinate_wl, polarity_wl, theta_wl):
    """Evaluate the stump's response for each point on a rectangular grid."""
    feature_slice = np.meshgrid(x_rng, y_rng)[coordinate_wl]
    weak_learner_on_grid = polarity_wl * (2*(feature_slice gt;= theta_wl) - 1)

    return weak_learner_on_grid

Finding the best weak learner

At each round of AdaBoost, the samples are reweighted, thus producing a new classification prob- lem, where the samples with a larger weigth count more in the classification error. The first step of a new round is to find the weak learner with the best performance for this new problem, that is, with the smallest classification error: