imp = xlsread('Academia.xlsx')
eva = evalclusters(imp, 'kmeans', 'CalinskiHarabasz','Klist',[1:15])
eva2 = evalclusters(imp, 'kmeans', 'Silhouette','Klist',[1:15])
eva3 = evalclusters(imp, 'kmeans', 'DaviesBouldin','Klist',[1:15])
subplot(2,2,1)
plot(eva)
title('CalinskiHarabasz')
subplot(2,2,2)
plot(eva2)
title('Silhouette')
subplot(2,2,3)
plot(eva3)
title('DaviesBouldin')
[idx,C]=kmeans(imp,4)
[numInst,numDims] = size(imp);
clr = lines(4);
figure, hold on
scatter3(imp(:,1), imp(:,2), imp(:,3), 36, clr(idx,:), 'Marker','.')
scatter3(C(:,1), C(:,2), C(:,3), 100, clr, 'Marker','o', 'LineWidth',3)
hold off
view(3), axis vis3d, box on, rotate3d on
xlabel('flexibilidade'), ylabel('velocidade'), zlabel('forca')