Changeset 9384 for branches/HeuristicLab.Analysis.AlgorithmBehavior

Timestamp:

04/19/13 13:47:47 (12 years ago)

Author:

ascheibe

Message:

#1886 added calculation of selection intensity

File:

• branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/SelectionPressureAnalyzer.cs (modified) (9 diffs)

branches/HeuristicLab.Analysis.AlgorithmBehavior/HeuristicLab.Analysis.AlgorithmBehavior.Analyzers/3.3/SelectionPressureAnalyzer.cs

@@ -22 +22 @@
 
 using System;
+using System.Collections.Generic;
 using System.Linq;
 using HeuristicLab.Common;
@@ -65 +66 @@
     }
     [Storable]
-    private ScatterPlotHelper selPressurePlot;
+    private ScatterPlotHelper selPressurePlot, selIntensityPlot;
     [Storable]
     private int cnt = 0;
@@ -79 +80 @@
       lastGeneration = original.lastGeneration;
       selPressurePlot = (ScatterPlotHelper)original.selPressurePlot.Clone(cloner);
+      selIntensityPlot = (ScatterPlotHelper)original.selIntensityPlot.Clone(cloner);
       SolutionQualityName = original.SolutionQualityName;
     }
@@ -93 +95 @@
 
       selPressurePlot = new ScatterPlotHelper(false, true, true, true);
+      selIntensityPlot = new ScatterPlotHelper(false, true, true, true);
     }
 
@@ -101 +104 @@
     protected override void InitializeAction() {
       selPressurePlot.InitializePlot(Results, "Selection Pressure", "Solution Index", "Quality Difference");
+      selIntensityPlot.InitializePlot(Results, "Selection Intensity", "Solution Index", "Scaled Quality Difference");
       ParentsQualityParameter.ActualName = SolutionQualityName;
 
@@ -121 +125 @@
 
         double oldPopQuality = 0.0;
+        List<double> qValues = new List<double>();
         int popSize = 0;
         foreach (IScope oldSolScope in oldPop.SubScopes) {
           double curQuality = ((DoubleValue)oldSolScope.Variables[SolutionQualityName].Value).Value;
           oldPopQuality += curQuality;
+          qValues.Add(curQuality);
           popSize++;
         }
@@ -138 +144 @@
               selPressurePlot.Max = bkQuality - wkQuality;
               selPressurePlot.Min = 0;
+              selIntensityPlot.Max = bkQuality - wkQuality;
+              selIntensityPlot.Min = 0;
             }
           } else {
@@ -143 +151 @@
               selPressurePlot.Max = wkQuality - bkQuality;
               selPressurePlot.Min = 0;
+              selIntensityPlot.Max = wkQuality - bkQuality;
+              selIntensityPlot.Min = 0;
             }
           }
         }
 
-        Point2D<double> popQualityPoint;
+        Point2D<double> popQualityPoint, selectionIntensityPoint;
         if (MaximizationParameter.ActualValue.Value) {
           popQualityPoint = new Point2D<double>(cnt, quality - (oldPopQuality / popSize));
+          selectionIntensityPoint = new Point2D<double>(cnt, (quality - (oldPopQuality / popSize) / qValues.StandardDeviation()));
         } else {
           popQualityPoint = new Point2D<double>(cnt, (oldPopQuality / popSize) - quality);
+          selectionIntensityPoint = new Point2D<double>(cnt, ((oldPopQuality / popSize) - quality) / qValues.StandardDeviation());
         }
 
         selPressurePlot.AddPoint(curGenStr, popQualityPoint);
+        selIntensityPlot.AddPoint(curGenStr, selectionIntensityPoint);
       }
 
@@ -167 +180 @@
     public override void ClearState() {
       selPressurePlot.CleanUp();
+      selIntensityPlot.CleanUp();
     }