Skip to content

Commit a0c7a8a

Browse files
committed
Added both symmetric and asymmetric version to figure
1 parent fa9eefe commit a0c7a8a

File tree

1 file changed

+21
-2
lines changed

1 file changed

+21
-2
lines changed

demo_successive_norm.m

Lines changed: 21 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -22,16 +22,22 @@
2222
results{2}.method = 'Row-Column Standardize';
2323
results{2}.output = standard_correlation_sn(X)
2424

25+
results{3}.method = 'Row-Column Sym Standardize';
26+
results{3}.output = standard_correlation_sym_sn(X)
2527

2628
figure;
27-
subplot(1,2,1);
29+
subplot(1,3,1);
2830
imagesc(results{1}.output.corr); axis equal image;
2931
colormap('winter');
3032
title(results{1}.method);
31-
subplot(1,2,2);
33+
subplot(1,3,2);
3234
imagesc(results{2}.output.corr); axis equal image;
3335
colormap('winter')
3436
title(results{2}.method);
37+
subplot(1,3,3);
38+
imagesc(results{3}.output.corr); axis equal image;
39+
colormap('winter')
40+
title(results{3}.method);
3541

3642
fname = ['tmp' filesep datestr(now,'dd-mmm-yyyy-HHMMSS')];
3743
if(~exist(fname,'dir'))
@@ -58,6 +64,18 @@
5864

5965
output = struct();
6066

67+
X = standardize.successive_normalize(X);
68+
69+
ggmobj = GGM(X);
70+
ggmobj.MLECovEstimate();
71+
output.corr = ggmobj.Sigma;
72+
73+
end
74+
75+
function output = standard_correlation_sym_sn(X)
76+
77+
output = struct();
78+
6179
X = standardize.successive_normalize(X, ...
6280
struct('method','sym'));
6381

@@ -66,3 +84,4 @@
6684
output.corr = ggmobj.Sigma;
6785

6886
end
87+

0 commit comments

Comments
 (0)