function D = ID3(train_features, train_targets, params, region)
% Classify using Quinlan's ID3 algorithm
% Inputs:
% features - Train features
% targets - Train targets
% params - [Number of bins for the data, Percentage of incorrectly assigned samples at a node]
% region - Decision region vector: [-x x -y y number_of_points]
%
% Outputs
% D - Decision sufrace
[Ni, M] =
size(train_features); %·μ??DDêyNioíáDêyM
%Get parameters
[Nbins, inc_node] = process_params(params);
inc_node = inc_node*M/100;
%For the decision region
N = region(5);
mx = ones(N,1) * linspace
(region(1),region(2),N); %linspace(?eê??μ£????1?μ£??a????êy)
my = linspace (region(3),region(4),N)' * ones(1,N);
flatxy = [mx(:), my(:)]';
%Preprocessing
[f, t, UW, m] = PCA(train_features,
train_targets, Ni, region);
train_features = UW * (train_features -
m*ones(1,M));
flatxy = UW * (flatxy - m*ones(1,N^2));