Changeset 3905 for trunk/sources/HeuristicLab.Problems.DataAnalysis.Regression/3.3/Symbolic/Analyzers/BestSymbolicRegressionSolutionAnalyzer.cs

Timestamp:

06/08/10 17:47:15 (14 years ago)

Author:

gkronber

Message:

Implemented calculation of variable impacts as integral over relative variable frequencies. #1036

File:

• trunk/sources/HeuristicLab.Problems.DataAnalysis.Regression/3.3/Symbolic/Analyzers/BestSymbolicRegressionSolutionAnalyzer.cs (modified) (6 diffs)

trunk/sources/HeuristicLab.Problems.DataAnalysis.Regression/3.3/Symbolic/Analyzers/BestSymbolicRegressionSolutionAnalyzer.cs

@@ -35 +35 @@
 using HeuristicLab.Problems.DataAnalysis;
 using HeuristicLab.Problems.DataAnalysis.Evaluators;
+using HeuristicLab.Analysis;
 
 namespace HeuristicLab.Problems.DataAnalysis.Regression.Symbolic.Analyzers {
@@ -43 +44 @@
     private const string SymbolicExpressionTreeInterpreterParameterName = "SymbolicExpressionTreeInterpreter";
     private const string BestSolutionInputvariableCountResultName = "Variables used by best solution";
+    private const string VariableFrequenciesParameterName = "VariableFrequencies";
+    private const string VariableImpactsResultName = "Integrated variable frequencies";
     private const string BestSolutionParameterName = "BestSolution";
 
@@ -55 +58 @@
       get { return (ILookupParameter<SymbolicRegressionSolution>)Parameters[BestSolutionParameterName]; }
     }
+    public ILookupParameter<DataTable> VariableFrequenciesParameter {
+      get { return (ILookupParameter<DataTable>)Parameters[VariableFrequenciesParameterName]; }
+    }
     #endregion
     #region properties
@@ -63 +69 @@
       get { return SymbolicExpressionTreeParameter.ActualValue; }
     }
+    public DataTable VariableFrequencies {
+      get { return VariableFrequenciesParameter.ActualValue; }
+    }
     #endregion
 
@@ -69 +78 @@
       Parameters.Add(new ScopeTreeLookupParameter<SymbolicExpressionTree>(SymbolicExpressionTreeParameterName, "The symbolic expression trees to analyze."));
       Parameters.Add(new ValueLookupParameter<ISymbolicExpressionTreeInterpreter>(SymbolicExpressionTreeInterpreterParameterName, "The interpreter that should be used for the analysis of symbolic expression trees."));
+      Parameters.Add(new LookupParameter<DataTable>(VariableFrequenciesParameterName, "The variable frequencies table to use for the calculation of variable impacts"));
       Parameters.Add(new LookupParameter<SymbolicRegressionSolution>(BestSolutionParameterName, "The best symbolic regression solution."));
+    }
+
+    [StorableHook(HookType.AfterDeserialization)]
+    private void Initialize() {
+      if (!Parameters.ContainsKey(VariableFrequenciesParameterName)) {
+        Parameters.Add(new LookupParameter<DataTable>(VariableFrequenciesParameterName, "The variable frequencies table to use for the calculation of variable impacts"));
+      }
     }
 
@@ -89 +106 @@
         if (Results.ContainsKey(BestSolutionInputvariableCountResultName)) {
           Results[BestSolutionInputvariableCountResultName].Value = new IntValue(model.InputVariables.Count());
+          Results[VariableImpactsResultName].Value = CalculateVariableImpacts();
         } else {
           Results.Add(new Result(BestSolutionInputvariableCountResultName, new IntValue(model.InputVariables.Count())));
+          Results.Add(new Result(VariableImpactsResultName, CalculateVariableImpacts()));
         }
       }
       return BestSolutionParameter.ActualValue;
+    }
+
+    private DoubleMatrix CalculateVariableImpacts() {
+      if (VariableFrequencies != null) {
+        var impacts = new DoubleMatrix(VariableFrequencies.Rows.Count, 1, new string[] { "Impact" }, VariableFrequencies.Rows.Select(x => x.Name));
+        int rowIndex = 0;
+        foreach (var dataRow in VariableFrequencies.Rows) {
+          string variableName = dataRow.Name;
+          double integral = 0;
+          if (dataRow.Values.Count > 1) {
+            double baseline = dataRow.Values.First();
+            integral = (from value in dataRow.Values
+                        select value - baseline)
+                              .Sum();
+          }
+          impacts[rowIndex++, 0] = integral;
+        }
+        return impacts;
+      } else return new DoubleMatrix(1, 1);
     }