To check out this repository please hg clone the following URL, or open the URL using EasyMercurial or your preferred Mercurial client.
root / _FullBNT / BNT / graph / moralize.m @ 8:b5b38998ef3b
History | View | Annotate | Download (322 Bytes)
| 1 |
function [M, moral_edges] = moralize(G) |
|---|---|
| 2 |
% MORALIZE Ensure that for every child, all its parents are married, and drop directionality of edges. |
| 3 |
% [M, moral_edges] = moralize(G) |
| 4 |
|
| 5 |
M = G; |
| 6 |
n = length(M); |
| 7 |
for i=1:n |
| 8 |
fam = family(G,i); |
| 9 |
M(fam,fam)=1; |
| 10 |
end |
| 11 |
M = setdiag(M,0); |
| 12 |
moral_edges = sparse(triu(max(0,M-G),1)); |