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
|
function score = cond_mutual_info_score(i,si,j,sj,c,sc,data)
% G = cond_mutual_info_score(i,si,j,sj,c,sc,data)
% Only for tabular node which values are 1,2,...,size .
% si is size of node i, sj is size of node j, sc is the size of node c.
% data(i,m) is the node i in the case m.
%
%
% pphilippe.leray@univ-nantes.fr, francois.olivier.c.h@gmail.com
[n N]=size(data);
Pc=hist(data(c,:),1:sc)/N;
score=0;
for cvalue=1:sc,
ind=find(data(c,:)==cvalue);
Nj=hist(data(j,ind),1:sj);
Ni=hist(data(i,ind),1:si);
NiNj=Ni'*Nj;
for k=1:si
ind2=find(data(i,ind)==k) ;
Nij(k,:) = hist(data(j,ind(ind2)),1:sj);
end
% sommons les valeurs non-infinies:
ind=find(NiNj~=0 & Nij~=0);
score=score+Pc(cvalue)*sum(sum(Nij(ind).*log(N*Nij(ind)./NiNj(ind))/N));
end
|
