clear
load count.dat
X = count;
typy_odleglosci = {'sqEuclidean'; 'cityblock'; 'cosine'; 'correlation'; 'Hamming'};
rodzaj = {'Euclidean'; 'sqEuclidean'; 'cityblock' ; 'cosine'; 'correlation';'Hamming'};
wynik = zeros(8,8);
sup = 0;
supI = 0;
supRodzaj = char(rodzaj(1));
for i = 1:6
    for j = 1:4
        K = kmeans(X,i,'distance',char(typy_odleglosci(j)));
        [silh,h] = silhouette(X,K,char(typy_odleglosci(j)));
        S = mean(silh);
        if(S > sup)
            sup = S;
            supI= K;
            supRodzaj = char(rodzaj(j));
        end
        wynik(i,j) = S;
    end
end
[silh,h] = silhouette(X,supI,supRodzaj)