-
Notifications
You must be signed in to change notification settings - Fork 2
/
solution_evaluation.m
65 lines (60 loc) · 1.26 KB
/
solution_evaluation.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
function [NC,Sil,Silmin,NCfix] = solution_evaluation(data,M,labels,NC,NCfix,simatrix,nrow,type,cut)
if type == 1
type = 'euclidean';
else
type = 'correlation';
end
dim = 0;
if simatrix
% M = simatrix_make(data,type,nrow);
Ms = NaN*ones(nrow,nrow);
for i=1:nrow
Ms(i,i) = 0;
end
for i=1:size(M,1)
ni=M(i,1);
nj=M(i,2);
Y = -M(i,3);
Ms(ni,nj)= Y;
if Y < dim
dim = Y;
end
end
end
if dim < 0
Ms = Ms - dim + 1;
for i=1:nrow
Ms(i,i) = 0;
end
end
dim = length(NC);
Sil =[];
Sildelete = [];
for i = 1:dim
Y = labels(:,i);
if simatrix
Smax = silhouette2(data, Y, Ms);
else
Smax = silhouette(data, Y, type);
end
dn = isfinite(Smax);
Sil(i) = mean(Smax(dn));
[C, Y, dmax] = ind2cluster(Y);
dmax = min(dmax);
Sildelete(i) = dmax < cut;
Q =[];
for j = 1:length(C)
R = C{j};
R = Smax(R);
dn = isfinite(R);
Q(j) = mean(R(dn));
end
Silmin(i) = min(Q);
end
Sildelete = find(Sildelete);
if length(Sildelete) < dim
Sil(Sildelete) = [];
Silmin(Sildelete) = [];
NC(Sildelete) = [];
NCfix(Sildelete) = [];
end