comparison core/magnatagatune/sim_from_comparison_fair_components.m @ 0:e9a9cd732c1e tip

first hg version after svn

+function [partBinTrn, partBinTst, partBinNoTrn] = sim_from_comparison_fair_components(comparison, comparison_ids, k, trainpart, filename)
+%
+% [partBinTrn, partBinTst, partBinNoTrn] =
+% sim_from_comparison_fair_components(comparison, comparison_ids, k, trainpart, [filename])
+
+% creates a cross-validation partitioning of the 
+% similarity data in "multiG", PRESERVING the 
+% connected components in it during partitioning
+
+% ---
+% get the similarity multigraph and remove cycles
+% ---
+cprint(2, 'creating graph')
+% Gm = ClipSimGraphMulti(comparison, comparison_ids);
+% Gm.remove_cycles_length2; 
+cprint(2, 'loading Multigraph for Similarity Constraints')
+load('comp_SimGraphMulti.mat', 'G');
+
+% ---
+% Note: we get the connected components in the graph
+% and filter out those who have only one node
+% ---
+cprint(2, 'extracting connected components')
+[Gs, s, id] = connected_components(G);
+
+valid = find(s > 1);
+Gsv = Gs(valid);
+
+% ---
+% We randomise the graph triplet order,
+% as well as the in-component
+% constraint succession be randomised here.
+% ---
+datPermu = randperm(numel(Gsv));
+Gsv = Gsv(datPermu);
+
+conPermu = zeros(numel(Gsv),3);
+for i = 1:numel(Gsv)
+    conPermu(i,:) = randperm(3);
+end
+
+% ---
+% NOTE: we try the easy route: partition the graphs
+% and look at which constraints  balance we end up with
+% ---
+P = cvpartition(numel(Gsv), 'k', k);
+
+% ---
+% here we export the graphs similarity test sets
+% ---
+cprint(2, 'export test similarity')
+partBinTst = {};
+for i = 1:P.NumTestSets  % test runs
+    partBinTst{i} = zeros(0, 3);
+    
+    tmp_idx = find(P.test(i));
+    for j = 1:numel(tmp_idx); % componens
+        
+        % ---
+        % get the graphs which are associated 
+        % to this set and save them into a new bin.
+        % ---
+         [weights, a, b, c] = Gsv(tmp_idx(j)).similarities();
+         partBinTst{i} = [partBinTst{i}; [a' b' c' weights]];
+    end
+end
+
+
+% ---
+% Note: This uses a "truly" increasing training set
+% to do the partial training partition
+% ---
+cprint(2, 'export train similarity')
+for m = 1:numel(trainpart)
+
+    Ptrain(m) = cvpartition_trunctrain_incsubsets(P, trainpart(m));
+end
+
+% ---
+% here we export the graph's similarity training sets
+% ---
+partBinTrn = {};
+for i = 1:P.NumTestSets % train runs
+
+    for m = 1:numel(trainpart) % increasing training sets
+        partBinTrn{i,m} = zeros(0, 3);
+    
+        tmp_idx = find(Ptrain(m).training(i));
+        for j = 1:numel(tmp_idx); % components 
+       
+            % ---
+            % get the graphs which are associated 
+            % to this set and save them into a new bin.
+            % ---
+             [weights, a, b, c] = Gsv(tmp_idx(j)).similarities();
+             
+             % ---
+             % NOTE: WE apply the inner-triplet permutation,
+             % and truncate it where necessary
+             % ---
+             tmp_permu = conPermu(tmp_idx(j),:);
+             if numel(a) < 3
+                 tmp_permu = tmp_permu(tmp_permu <= numel(a));
+             end
+             
+             a = a(tmp_permu);
+             b = b(tmp_permu);
+             c = c(tmp_permu);
+             weights = weights(tmp_permu);
+             
+             % save the clips
+             partBinTrn{i,m} = [partBinTrn{i,m}; [a' b' c' weights]];
+        end
+    end
+end
+
+partBinNoTrn = {};
+for i = 1:P.NumTestSets % train runs
+
+    for m = 1:numel(trainpart) % increasing training sets
+        partBinNoTrn{i,m} = zeros(0, 3);
+    
+        tmp_idx = find(~Ptrain(m).training(i) & ~Ptrain(m).test(i));
+        for j = 1:numel(tmp_idx); % components 
+       
+            % ---
+            % get the graphs which are associated 
+            % to this set and save them into a new bin.
+            % ---
+             [weights, a, b, c] = Gsv(tmp_idx(j)).similarities();
+             
+             % ---
+             % NOTE: WE apply the inner-triplet permutation,
+             % and truncate it where necessary
+             % ---
+             tmp_permu = conPermu(tmp_idx(j),:);
+             if numel(a) < 3
+                 tmp_permu = tmp_permu(tmp_permu <= numel(a));
+             end
+             
+             a = a(tmp_permu);
+             b = b(tmp_permu);
+             c = c(tmp_permu);
+             weights = weights(tmp_permu);
+
+             % save the clips
+             partBinNoTrn{i,m} = [partBinNoTrn{i,m}; [a' b' c' weights]];
+        end
+    end
+end
+
+if nargin == 5
+    save(filename, 'partBinTrn', 'partBinTst', 'partBinNoTrn')
+end
+
+end