function test

load 2moons;

Y = produce_labelY(y);

%{

for alpha = 0 : 0.1 : 1.0

    for sigma = 1.5 : 0.1 : 2.5;

       LGC(x, Y, sigma, alpha);

    end

end

%}

%{

for sigma = 0.05 : 0.01 : 0.20

LGC(x, Y, sigma, 0.9);

end

%}

%for alpha = 0.0 : 0.1 : 0.9

%{

for alpha = 0.9 : 0.01 : 0.99

LGC(x, Y, 0.11, alpha);

end

%}

%{

for sigma = 0.05 : 0.005 : 0.09

LGC(x, Y, sigma, 0.99);

end

%}

%{

min_error = 200;

best_s = 0;

best_a = 0;

for sigma = 0 : 0.005 : 0.1

    for alpha = 0.9 : 0.01 : 0.99

        error_number = LGC(x, Y, -y, sigma, alpha, 0);

        if error_number < min_error

            min_error = error_number;

            best_s = sigma;

            best_a = alpha;

        end

    end

end

error_number = LGC(x, Y, -y, best_s, best_a, 1);

%}

% The last proper parameter tuning

min_error = 200;

best_s = 0;

best_a = 0;

for sigma = 0.07 : 0.001 : 0.08

    for alpha = 0.95 : 0.001 : 0.999

        %error_number = LGC(x, Y, -y, sigma, alpha, 0, 1);

        error_number = LGC(x, Y, -y, sigma, alpha, 0, 0);

        if error_number < min_error

            min_error = error_number;

            best_s = sigma;

            best_a = alpha;

        end

    end

end

error_number = LGC(x, Y, -y, best_s, best_a, 1, 0)

%error_rate = LGC(x, Y, -y, 0.07, 0.996, 1, 0);

function Y = produce_labelY(y)

n = size(y, 1);

Y = zeros(n, 2);

%{

label_number = 40;

index = find(y == 1);

Y( index(find(index <= label_number)), 1) = 1;

index = find(y == -1);

Y( index(find(index <= label_number)), 2) = 1;

%}

Y(53, 1) = 1;      %Positive example

Y(143, 2) = 1;     %Active exampl

————————————————