-
Notifications
You must be signed in to change notification settings - Fork 2
/
calculateNMI_semiPI.m
69 lines (66 loc) · 1.36 KB
/
calculateNMI_semiPI.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
66
67
68
function [nmi] = calculateNMI_semiPI(cluster,true_cluster,seeds)
n = size(cluster,1);
K = length(unique(cluster));
I = 0;
for i=1:K %cluster
for j=1:K %label
wkcj = 0;
wk = 0;
cj = 0;
N = 0;
for p=1:n
if ~ismember(p,seeds)
N = N + 1;
if (cluster(p) == i) && (true_cluster(p) == j)
wkcj = wkcj + 1;
end;
if cluster(p) == i
wk = wk + 1;
end;
if true_cluster(p) == j
cj = cj + 1;
end;
end;
end;
if wkcj > 0
I = I + (wkcj/N)*log((N*wkcj)/(wk*cj) + realmin);
end;
end;
end;
H1 = 0;
for i=1:K
wk = 0;
N = 0;
for p=1:n
if ~ismember(p,seeds)
N = N + 1;
if cluster(p) == i
wk = wk + 1;
end;
end;
end;
if wk > 0
H1 = H1 + (wk/N)*log(wk/N + realmin);
end;
end;
H1 = H1 * (-1);
H2 = 0;
for j=1:K
cj = 0;
N = 0;
for p=1:n
if ~ismember(p,seeds)
N = N + 1;
if true_cluster(p) == j
cj = cj + 1;
end;
end;
end;
if cj > 0
H2 = H2 + (cj/N) * log(cj/N + realmin);
end;
end;
H2 = H2 * (-1);
nmi = 2 * I / (H1 + H2);
% fprintf('%f %f %f', I, H1, H2);
end