Free cookie consent management tool by TermsFeed Policy Generator

source: branches/FitnessLandscapeAnalysis/HeuristicLab.Analysis.FitnessLandscape/Analysis/UpDownWalkAnalyzer.cs @ 7128

Last change on this file since 7128 was 7128, checked in by epitzer, 12 years ago

#1696 Integrate fitness landscape analysis plugins from Heureka! repository.

File size: 10.3 KB
Line 
1#region License Information
2/* HeuristicLab
3 * Copyright (C) 2002-2010 Heuristic and Evolutionary Algorithms Laboratory (HEAL)
4 *
5 * This file is part of HeuristicLab.
6 *
7 * HeuristicLab is free software: you can redistribute it and/or modify
8 * it under the terms of the GNU General Public License as published by
9 * the Free Software Foundation, either version 3 of the License, or
10 * (at your option) any later version.
11 *
12 * HeuristicLab is distributed in the hope that it will be useful,
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15 * GNU General Public License for more details.
16 *
17 * You should have received a copy of the GNU General Public License
18 * along with HeuristicLab. If not, see <http://www.gnu.org/licenses/>.
19 */
20#endregion
21
22using System.Collections.Generic;
23using System.Drawing;
24using System.Linq;
25using HeuristicLab.Collections;
26using HeuristicLab.Common;
27using HeuristicLab.Core;
28using HeuristicLab.Data;
29using HeuristicLab.Operators;
30using HeuristicLab.Optimization.Operators;
31using HeuristicLab.Parameters;
32using HeuristicLab.Persistence.Default.CompositeSerializers.Storable;
33
34namespace HeuristicLab.Analysis.FitnessLandscape.Analysis {
35
36  [StorableClass]
37  public class UpDownWalkAnalyzer : AlgorithmOperator, IQualityTrailAnalyzer {
38
39    #region Parameters
40    public LookupParameter<DataTable> QualityTrailParameter {
41      get { return (LookupParameter<DataTable>)Parameters["Quality Trail"]; }
42    }
43    public LookupParameter<VariableCollection> ResultsParameter {
44      get { return (LookupParameter<VariableCollection>)Parameters["Results"]; }
45    }
46    public LookupParameter<DataTable> UpDownStepsParameter {
47      get { return (LookupParameter<DataTable>)Parameters["UpDownSteps"]; }
48    }
49    public LookupParameter<DataTable> UpDownLevelsParameter {
50      get { return (LookupParameter<DataTable>)Parameters["UpDownLevels"]; }
51    }
52    public LookupParameter<DoubleValue> UpWalkLengthParameter {
53      get { return (LookupParameter<DoubleValue>)Parameters["UpWalkLength"]; }
54    }
55    public LookupParameter<DoubleValue> DownWalkLengthParameter {
56      get { return (LookupParameter<DoubleValue>)Parameters["DownWalkLength"]; }
57    }
58    public LookupParameter<DoubleValue> UpWalkLenVarParameter {
59      get { return (LookupParameter<DoubleValue>)Parameters["UpWalkLenVar"]; }
60    }
61    public LookupParameter<DoubleValue> DownWalkLenVarParameter {
62      get { return (LookupParameter<DoubleValue>)Parameters["DownWalkLenVar"]; }
63    }
64    public LookupParameter<DoubleValue> UpperLevelParameter {
65      get { return (LookupParameter<DoubleValue>)Parameters["UpperLevel"]; }
66    }
67    public LookupParameter<DoubleValue> LowerLevelParameter {
68      get { return (LookupParameter<DoubleValue>)Parameters["LowerLevel"]; }
69    }
70    public LookupParameter<DoubleValue> UpperVarianceParameter {
71      get { return (LookupParameter<DoubleValue>)Parameters["UpperVariance"]; }
72    }
73    public LookupParameter<DoubleValue> LowerVarianceParameter {
74      get { return (LookupParameter<DoubleValue>)Parameters["LowerVariance"]; }
75    }
76    #endregion
77
78    [StorableConstructor]
79    protected UpDownWalkAnalyzer(bool deserializing) : base(deserializing) { }
80    protected UpDownWalkAnalyzer(UpDownWalkAnalyzer original, Cloner cloner) : base(original, cloner) { }
81
82    public UpDownWalkAnalyzer() {
83      Parameters.Add(new LookupParameter<DataTable>("Quality Trail", "The qualities of the solutions"));
84      Parameters.Add(new LookupParameter<VariableCollection>("Results", "The collection of all results of this algorithm"));
85      Parameters.Add(new LookupParameter<DataTable>("UpDownSteps", "Distribution of upward and downward steps between extremes."));
86      Parameters.Add(new LookupParameter<DataTable>("UpDownLevels", "Distribution of upper and lower levels in up-down walks."));
87      Parameters.Add(new LookupParameter<DoubleValue>("DownWalkLength", "Average downward walk length."));
88      Parameters.Add(new LookupParameter<DoubleValue>("UpWalkLength", "Average updward walk length."));
89      Parameters.Add(new LookupParameter<DoubleValue>("UpWalkLenVar", "Upward walk length variance."));
90      Parameters.Add(new LookupParameter<DoubleValue>("DownWalkLenVar", "Downward walk length variance."));
91      Parameters.Add(new LookupParameter<DoubleValue>("UpperLevel", "Average maximum fitness value."));
92      Parameters.Add(new LookupParameter<DoubleValue>("LowerLevel", "Average minimum fitness value."));
93      Parameters.Add(new LookupParameter<DoubleValue>("LowerVariance", "Lower level variance."));
94      Parameters.Add(new LookupParameter<DoubleValue>("UpperVariance", "Upper level variance."));
95   
96
97      var resultsCollector = new ResultsCollector();
98      resultsCollector.CollectedValues.Add(new LookupParameter<DataTable>(UpDownStepsParameter.Name));
99      resultsCollector.CollectedValues.Add(new LookupParameter<DataTable>(UpDownLevelsParameter.Name));
100      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(DownWalkLengthParameter.Name));
101      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(UpWalkLengthParameter.Name));
102      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(UpWalkLenVarParameter.Name));
103      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(DownWalkLenVarParameter.Name));
104      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(UpperLevelParameter.Name));
105      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(LowerLevelParameter.Name));
106      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(UpperVarianceParameter.Name));
107      resultsCollector.CollectedValues.Add(new LookupParameter<DoubleValue>(LowerVarianceParameter.Name));
108
109      OperatorGraph.InitialOperator = resultsCollector;
110      resultsCollector.Successor = null;
111    }
112
113    public override IDeepCloneable Clone(Cloner cloner) {
114      return new UpDownWalkAnalyzer(this, cloner);
115    }
116
117    public override IOperation Apply() {
118      var stepsTable = GetOrCreateTable(UpDownStepsParameter, "Up", "Down");
119      var levelsTable = GetOrCreateTable(UpDownLevelsParameter, "Top", "Bottom");
120      DataTable qualityTrail = QualityTrailParameter.ActualValue;
121      if (qualityTrail != null && qualityTrail.Rows.Count > 0) {
122        var qualities = qualityTrail.Rows.First().Values.ToList();
123        if (qualities.Count > 2) {
124          var extremes = qualities
125            .Delta((a, b) => new {a, b, diff = b - a})
126            .Select((p, i) => new {p.a, p.b, p.diff, i = i+1})
127            .Delta((s1, s2) => new {
128              s1.i,
129              value = s2.a,
130              top = s1.diff >= 0 && s2.diff < 0,
131              bottom = s1.diff <= 0 && s2.diff > 0
132            })
133            .Where(x => x.top || x.bottom)
134            .GroupConsecutive(x => x.bottom)
135            .Select(g => g.Count() == 1
136                           ? g.First()
137                           : (g.First().bottom
138                                ? g.OrderBy(x => x.value).First()
139                                : g.OrderByDescending(x => x.value).First())).ToList();
140          var maxima = extremes.Where(x => x.top).ToList();
141          var minima = extremes.Where(x => x.bottom).ToList();
142          var tops = Enumerable.Repeat(new {length = 0, value = 0.0}, 0).ToList();
143          var bottoms = tops;
144          if (maxima.Count > 0 && minima.Count > 0) {
145            if (maxima.First().i < minima.First().i) {
146              bottoms = maxima.Zip(minima, (t, b) => new {length = b.i - t.i, b.value}).ToList();
147              minima.Insert(0, new {i = -1, value = 0.0, top = false, bottom = false});
148              tops = maxima.Zip(minima, (t, b) => new {length = t.i - b.i, t.value}).ToList();
149            } else {
150              tops = maxima.Zip(minima, (t, b) => new {length = t.i - b.i, t.value}).ToList();
151              maxima.Insert(0, new {i = -1, value = 0.0, top = false, bottom = false});
152              bottoms = maxima.Zip(minima, (t, b) => new {length = b.i - t.i, b.value}).ToList();
153            }
154            if (tops.Count > 0) {
155              var topLengths = tops.Select(t => (double) t.length).ToList();
156              var topVals = tops.Select(t => t.value).ToList();
157              ReplaceHistogram(stepsTable.Rows["Up"], topLengths);
158              ReplaceHistogram(levelsTable.Rows["Top"], topVals);
159              UpperLevelParameter.ActualValue = new DoubleValue(topVals.Average());
160              UpperVarianceParameter.ActualValue = new DoubleValue(topVals.Variance());
161              UpWalkLengthParameter.ActualValue = new DoubleValue(topLengths.Average());
162              UpWalkLenVarParameter.ActualValue = new DoubleValue(topLengths.Variance());
163            }
164            if (bottoms.Count > 0) {
165              var bottomLengths = bottoms.Select(b => (double) b.length).ToList();
166              var bottomVals = bottoms.Select(b => b.value).ToList();
167              ReplaceHistogram(stepsTable.Rows["Down"], bottomLengths);
168              ReplaceHistogram(levelsTable.Rows["Bottom"], bottomVals);
169              LowerLevelParameter.ActualValue = new DoubleValue(bottomVals.Average());
170              LowerVarianceParameter.ActualValue = new DoubleValue(bottomVals.Variance());
171              DownWalkLengthParameter.ActualValue = new DoubleValue(bottomLengths.Average());
172              DownWalkLenVarParameter.ActualValue = new DoubleValue(bottomLengths.Variance());
173            }
174          }
175        }
176      }
177      return base.Apply();
178    }
179
180    private void ReplaceHistogram(DataRow row, IEnumerable<double> values) {
181      row.VisualProperties.ChartType = DataRowVisualProperties.DataRowChartType.Histogram;
182      row.VisualProperties.Bins = 25;
183      row.Values.Clear();
184      row.Values.AddRange(values);
185    }
186
187    private DataTable GetOrCreateTable(LookupParameter<DataTable> tableParameter, params string[] rowNames) {
188      DataTable table = tableParameter.ActualValue;
189      if (table == null) {
190        table = new DataTable(tableParameter.Name, tableParameter.Description);
191        tableParameter.ActualValue = table;
192        foreach (var name in rowNames) {
193          table.Rows.Add(new DataRow(name));
194        }
195      }
196      return table;
197    }
198  }
199}
Note: See TracBrowser for help on using the repository browser.