Scribd
Upload
97 views

Machine Learning in MATLAB: Roland Memisevic

This document provides an introduction to machine learning concepts and how to implement them in MATLAB. It discusses representing data as matrices, linear models using vector multiplication (wTx), and learning algorithms like gradient descent and perceptron learning. Gradient descent is used to learn the weights for linear regression by minimizing the error between predictions and targets. Perceptron learning is applied to learn a hyperplane that can perform linear classification. Visualization tools in MATLAB are demonstrated for gaining insights into machine learning problems and solutions.

Uploaded by

Rupesh Sushir
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
97 views

Machine Learning in MATLAB: Roland Memisevic

This document provides an introduction to machine learning concepts and how to implement them in MATLAB. It discusses representing data as matrices, linear models using vector multiplication (wTx), and learning algorithms like gradient descent and perceptron learning. Gradient descent is used to learn the weights for linear regression by minimizing the error between predictions and targets. Perceptron learning is applied to learn a hyperplane that can perform linear classification. Visualization tools in MATLAB are demonstrated for gaining insights into machine learning problems and solutions.

Uploaded by

Rupesh Sushir
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 16

Machine Learning in MATLAB

Roland Memisevic

January 25, 2007


Last time

I Working with MATLAB.


I Vectors, matrices, operations.
I Control structures.
I Scripts and functions.
I Plotting.
I Slicing/Logical indexing.
A few random things
I Get the size of an object with size. Takes an optional
argument to specify the dimension (without, it returns an
array with the sizes of all dimensions). Example:
A = [1 2 3
4 5 6]
size(A,1) = 2
size(A,2) = 3
size(A) = [2 3]

I Use ... to let commands span several lines.


I Clear memory by deleting some unused variables using clear
variable-name
I who and whos show the variables in the workspace.
I Use drawnow to draw now (flush the event queue) ...
I Use save to save the current workspace. load to get it back.
repmat

I Use repmat to build new matrices from smaller ones.


Example:
x = [1;2]
repmat(x,1,3) = [1 1 1
2 2 2]
repmat(x,2,2) = [1 1
2 2
1 1
2 2]

I Used mostly to vectorize things (more in a moment...)


Representing data

I We are often dealing with lots of data, and normally that data
is given in the form of real-vectors. If not, we can often
massage it accordingly... (performing feature extraction).
I It is convenient to use matrices to store data, for example by
stacking the vectors column-wise:

X = (x1 x2 . . . xN )

I To display up to 3-d data given in this form, you could use


scatter(X(1,:),X(2,:)) % display 2-d dataset
scatter3(X(1,:),X(2,:),X(3,:)) % display 3-d dataset
High dimensions

I We can index arrays linearly (regardless of the number of


dimensions). Example: If A=[1 2;3 4], then A(3) = 3.
I Summarizing operations, such as sum, mean, etc. can be
applied to arrays of any dimension (not just vectors). Use an
additional argument to specify over which dimension to
perform the operation. Example: If A = [1 2 ;3 4], then
sum(A,1) = [4 6]
sum(A,2) = [3
7]

I Many operations on vectors and matrices extend to objects


with more than 1 or 2 dimensions. Example: Construct a
random 5 5 5-array with randn(5,5,5).
Linear models

I Many machine learning methods are based on simple


neurons computing:

wTx ( =
X
wk xk )
k
I Perceptrons, back-prop networks, support vector machines,
logistic regression, PCA, (nearest neighbors), ...
Vectorizing
I Applying the linear function to data-points stacked
column-wise in a matrix X is simply Y = wT X .
In MATLAB:
Y = w*X;
I In MATLAB writing stuff in matrix form can be faster than
using loops. Referred to as vectorization.
I Another example. Suppose you want to mean center a set of
vectors stored in X . Instead of
m = mean(X,2);
for i = 1 : size(X,2)
X(:,ii) = X(:,ii) - m;
end

we could write:
X = X - repmat(mean(X,2),1,size(X,2));
Linear regression

I Problem: We are given input/output data (real values). Want


to learn the underlying function.
I Assumption for now f (x) = wT x and y = f (x) + noise
I We can use MATLAB to get a quick idea of whats going on:
X = randn(2,1000); % make some input data
w = [1;1]; % define a linear function
n = randn(1,1000)*0.1; % produce some noise
Y = w*X + n; % produce outputs
scatter3(X(1,:),X(2,:),Y);
Learning

I Unfortunately, in reality we dont know w...


I Given only a set of examples yi and xi , what would be a
reasonable guess for w?
I Standard approach: Minimize the sum of squared errors
X
E (w) = (wT xi yi )2
i

I There are many ways of finding the w that minimizes E (w).


I One very good way is this:
I Start with a random guess winit .
E
I Follow the negative gradient w , until the error stops
changing a lot.
Learning
E P Tx
I The gradient is w =2 i (w i yi )xi .
I So we just have to iterate:
X
w w 2 (wT xi yi )xi ,
i

where  is a small learning rate, without which we will


overshoot the minimum.
I With vectorization, learning takes about 5 lines in MATLAB:
for iteration = 1 : 5000 %in practice: until stopping
%criterion satisfied
grad = 2*sum(repmat(w*X-Y,size(X,1),1).*X,2);
w = w - epsilon * grad;
err = sum((Y - w*X).^2) %just to check
end
Visualizing the solution

I To see what the learned model does on the training data:


hold on; %keep old plot
Y_hat = w_learned*X;
scatter3(X(1,:),X(2,:),Y_hat, r);
Linear classification

I Classification: We are given some data with associated


class-labels, and want to learn the underlying function, that
maps inputs to labels.
I Examples: Spam filtering, face recognition, intrusion
detection, ... and many many many more.
I For now we consider only binary problems (2 classes). Most
ideas can be quite easily extended to more than 2.
I Again, we will use a linear model. In particular:

f (x) = sgn(wT x)

I This means that we try to separate classes using a


hyper-plane.
Linear classification
I Again we can use MATLAB to visualize whats going on:
X = [randn(3,200)-ones(3,200)*1.8 ...
randn(3,200)+ones(3,200)*1.8]; %produce some inputs
Y = [zeros(1,200), ones(1,200)]; %produce some labels
scatter3(X(1,:),X(2,:),X(3,:), 80, Y, filled);

I How well does some random w do on the training set?


w = randn(3,1);
Y_random = sign(w*X);
scatter3(X(1,:),X(2,:),X(3,:),80,Y_random, filled);
hold on;
plot3([0 w(1)], [0 w(2)], [0 w(3)], k); %show w
hold off;
sum(Y_random~=Y)/200 %error rate
Learning
I How can we learn a good hyper-plane?
I There are many possible ways.
I An extremely powerful one is perceptron learning:
I Start with some initial guess for w. Then iterate, picking
training examples (in any order):
I if sgn(wT xi ) = yi , then do nothing
I otherwise, set w := w xi , depending on what the output was
I In MATLAB:
for iteration = 1 : 100
for ii = 1 : size(X,2)
if sign(w*X(:,ii)) ~= Y(ii)
w = w + X(:,ii) * Y(ii);
end
end
end
Conclusions

I A lot of machine learning is based on the simple neuron:


wT x
I We have looked at basic regression and classification.
I Usually a few lines in MATLAB.
I A couple of things were oversimplified here. For example, in
practice we would adapt the learning rate in gradient descent,
add an extra input-dimension for the bias, etc.
I Can be easily applied to real data: E.g. Spam, ...

You might also like