clc
close all
clear all

F={'randtop'; 'gridtop';'hextop'};
topologyFcn=F{1};
G={'linkdist';'mandist';'boxdist'};
distanceFcn=G{1};

%z generatora randn
d1 = load('d1.mat');
d2 = load('d2.mat');
d3 = load('d3.mat');
d1 = d1.d1;
d2 = d2.d2;
d3 = d3.d3;

% d1=randn(2,100);
% d2=randn(2,100)+2; %<- to zmieni wykres, powinny trochę nacgodzić
% d3=randn(2,100);%<- to zmieni wykres, powinny trochę nacgodzić
d3(2,:)=d3(2,:)+3;%<- to zmieni wykres, powinny trochę nacgodzić
d=[d1,d2,d3];

D=[ones(1,100) 2*ones(1,100) 3*ones(1,100)];
D=D';

dimensions = [1 12];
epoki = 100;
[wynik, cmdx] = sieci(d, D, dimensions, epoki, topologyFcn, distanceFcn);
wynik

figure;
gscatter(d(1,:), d(2,:),cmdx);

Data = load('Data_microarray.mat');

d1 = Data.Data(:,1);
d2 = Data.Data(:,2);
d = [d1,d2]';
D = Data.Data(:,3)';

[wynik, cmdx] = sieci(d,D,dimensions, epoki, topologyFcn, distanceFcn);
wynik

figure;
gscatter(d(1,:), d(2,:),cmdx);