Changeset 13495 for branches/HeuristicLab.EvolutionTracking/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Tracking/Analyzers/SymbolicDataAnalysisFragmentLengthAnalyzer.cs

Timestamp:

01/11/16 21:26:00 (8 years ago)

Author:

bburlacu

Message:

#1772: Improve the average operator improvement and average fragment length analyzers. Update genealogy analyzer to always evaluate the intermediate children.

File:

• branches/HeuristicLab.EvolutionTracking/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Tracking/Analyzers/SymbolicDataAnalysisFragmentLengthAnalyzer.cs (modified) (2 diffs)

branches/HeuristicLab.EvolutionTracking/HeuristicLab.Problems.DataAnalysis.Symbolic/3.4/Tracking/Analyzers/SymbolicDataAnalysisFragmentLengthAnalyzer.cs

@@ -61 +61 @@
       updateCounter.Value = 1;
 
-      if (PopulationGraph == null)
+      if (PopulationGraph == null || Generation.Value == 0)
         return base.Apply();
 
-      var children = PopulationGraph.GetByRank(Generation.Value).Where(x => x.InDegree > 0).ToList();
-      double avgFragmentLength = 0;
-      if (children.Any())
-        avgFragmentLength = children.Select(x => (IFragment<ISymbolicExpressionTreeNode>)x.InArcs.Last().Data).Where(f => f != null).Average(x => x.Root.GetLength());
+      // consider all fragments (including those of the intermediate crossover children)
+      double averageCrossoverFragmentLength = 0;
+      double averageMutationFragmentLength = 0;
+      int crossoverChildren = 0;
+      int mutationChildren = 0;
+
+      var vertices = PopulationGraph.Vertices.Where(x => x.Rank > Generation.Value - 1);
+      foreach (var v in vertices) {
+        if (!v.InArcs.Any() || v.InArcs.Last().Data == null)
+          continue;
+        var fragment = (IFragment<ISymbolicExpressionTreeNode>)v.InArcs.Last().Data;
+        if (v.InDegree == 2) {
+          averageCrossoverFragmentLength += fragment.Root.GetLength();
+          crossoverChildren++;
+        } else {
+          averageMutationFragmentLength += fragment.Root.GetLength();
+          mutationChildren++;
+        }
+      }
+      var averageFragmentLength = (averageCrossoverFragmentLength + averageMutationFragmentLength) / (crossoverChildren + mutationChildren);
+      averageCrossoverFragmentLength /= crossoverChildren;
+      averageMutationFragmentLength /= mutationChildren;
 
       DataTable table;
@@ -73 +91 @@
         table = new DataTable("Average fragment length");
         var row = new DataRow("Average fragment length") { VisualProperties = { StartIndexZero = true } };
-        row.Values.Add(avgFragmentLength);
+        row.Values.Add(averageFragmentLength);
+        table.Rows.Add(row);
+        row = new DataRow("Average crossover fragment length") { VisualProperties = { StartIndexZero = true } };
+        row.Values.Add(averageCrossoverFragmentLength);
+        table.Rows.Add(row);
+        row = new DataRow("Average mutation fragment length") { VisualProperties = { StartIndexZero = true } };
+        row.Values.Add(averageMutationFragmentLength);
         table.Rows.Add(row);
         Results.Add(new Result("AverageFragmentLength", table));
       } else {
         table = (DataTable)Results["AverageFragmentLength"].Value;
-        table.Rows["Average fragment length"].Values.Add(avgFragmentLength);
+        table.Rows["Average fragment length"].Values.Add(averageFragmentLength);
+        table.Rows["Average crossover fragment length"].Values.Add(averageCrossoverFragmentLength);
+        table.Rows["Average mutation fragment length"].Values.Add(averageMutationFragmentLength);
       }
       return base.Apply();