Changeset 9338 for branches/HeuristicLab.Problems.GaussianProcessTuning/GaussianProcessDemo/Form1.cs

Timestamp:

03/31/13 13:17:49 (11 years ago)

Author:

gkronber

Message:

#1967: minor adaptations necessary for the EuroCAST presentation

Location:

branches/HeuristicLab.Problems.GaussianProcessTuning

Files:

• . (modified) (1 prop)
• GaussianProcessDemo/Form1.cs (modified) (7 diffs)

branches/HeuristicLab.Problems.GaussianProcessTuning

@@ -4 +4 @@
 *.testsettings
 *.user
+bin

@@ -4 +4 @@
 *.testsettings
 *.user
+bin

branches/HeuristicLab.Problems.GaussianProcessTuning/GaussianProcessDemo/Form1.cs

@@ -10 +10 @@
 using HeuristicLab.Algorithms.DataAnalysis;
 using HeuristicLab.Core;
+using HeuristicLab.Data;
 using HeuristicLab.Problems.DataAnalysis;
 using HeuristicLab.Problems.Instances.DataAnalysis;
@@ -27 +28 @@
 
       var sum = new CovarianceSum();
-      sum.Terms.Add(new CovariancePeriodic());
+      var t = new CovarianceSquaredExponentialIso();
+      t.InverseLengthParameter.Value = new DoubleValue(1.0 / Math.Exp(-2));
+      sum.Terms.Add(t);
       sum.Terms.Add(new CovarianceNoise());
       this.covFunction = sum;
@@ -38 +41 @@
     private void UpdateSliders() {
       flowLayoutPanel1.Controls.Clear();
+      flowLayoutPanel1.Controls.Add(dataButton);
       for (int i = 0; i < covFunction.GetNumberOfParameters(1); i++) {
         var sliderControl = new TrackBar();
@@ -49 +53 @@
 
     private void InitData() {
-      int n = 20;
+      int n = 200;
       data = new List<List<double>>();
       data.Add(ValueGenerator.GenerateSteps(0, 1, 1.0 / n).ToList());
@@ -69 +73 @@
       }
 
+      var trainingData = new List<List<double>>();
+      var trainingIndices = RandomEnumerable.SampleRandomWithoutRepetition(Enumerable.Range(0, y.Count), random, 10);
+      var trainingY = trainingIndices.Select(i => y[i]).ToList();
+      var trainingX = trainingIndices.Select(i => data[0][i]).ToList();
+      trainingData.Add(trainingY);
+      trainingData.Add(trainingX);
+      
+      //chart1.Series[2].Points.Clear();
+      //foreach (var p in trainingY.Zip(trainingX, (t, x) => new { t, x })) {
+      //  chart1.Series[2].Points.AddXY(p.x, p.t);
+      //}
+
       var allData = new List<List<double>>();
       allData.Add(y);
       allData.Add(data[0]);
       var variableNames = new string[] { "y", "x" };
-      var ds = new Dataset(variableNames, allData);
-      var rows = Enumerable.Range(0, data[0].Count);
-      var correctModel = new GaussianProcessModel(ds, variableNames.First(), variableNames.Skip(1), rows, hyp, new MeanZero(),
+      var fullDataSet = new Dataset(variableNames, allData);
+      var trainingDataSet = new Dataset(variableNames, trainingData);
+      var trainingRows = Enumerable.Range(0, trainingIndices.Count());
+      var fullRows = Enumerable.Range(0, data[0].Count);
+      var correctModel = new GaussianProcessModel(fullDataSet, variableNames.First(), variableNames.Skip(1), fullRows, hyp, new MeanZero(),
                                                 (ICovarianceFunction)covFunction.Clone());
-      var yPred = correctModel.GetEstimatedValues(ds, rows);
+      var yPred = correctModel.GetEstimatedValues(fullDataSet, fullRows);
       chart1.Series[1].Points.Clear();
       foreach (var p in yPred.Zip(data[0], (t, x) => new { t, x })) {
@@ -86 +104 @@
     private double[] GetSliderValues() {
       var hyp = new List<double>();
-      foreach (var slider in flowLayoutPanel1.Controls.Cast<TrackBar>()) {
+      foreach (var slider in flowLayoutPanel1.Controls.OfType<TrackBar>()) {
         Console.Write(slider.Value / 10.0 + " ");
         hyp.Add(slider.Value / 10.0);
@@ -94 +112 @@
       return hyp.ToArray();
     }
+
+    private void dataButton_Click(object sender, EventArgs e) {
+      var dataForm = new Form();
+      var dataTextField = new TextBox();
+      dataTextField.Multiline = true;
+      dataTextField.Text = DataToText();
+      dataTextField.Dock = DockStyle.Fill;
+      dataForm.Controls.Add(dataTextField);
+      dataForm.ShowDialog();
+    }
+
+    private string DataToText() {
+      var str = new StringBuilder();
+      foreach (var p in chart1.Series[1].Points) {
+        str.AppendLine(p.XValue + "\t" + p.YValues.First());
+      }
+      return str.ToString();
+    }
   }
 }

@@ -4 +4 @@
 *.testsettings
 *.user
+bin