static void Main(string[] args)
{
DataTable table = LoadData();
TGPConfig config = new TGPConfig("TGPConfig.xml");
TGPPop <TGPSolution> pop = new TGPPop <TGPSolution>(config);
pop.OperatorSet.AddOperator(new TGPOperator_Plus());
pop.OperatorSet.AddOperator(new TGPOperator_Minus());
pop.OperatorSet.AddOperator(new TGPOperator_Division());
pop.OperatorSet.AddOperator(new TGPOperator_Multiplication());
pop.OperatorSet.AddOperator(new TGPOperator_Power());
pop.OperatorSet.AddIfltOperator();
for (int i = 1; i < 10; ++i)
{
pop.ConstantSet.AddConstant(string.Format("C{0}", i), i);
}
pop.VariableSet.AddVariable("X1");
pop.VariableSet.AddVariable("X2");
pop.CreateFitnessCase += (index) =>
{
MexicanHatFitnessCase fitness_case = new MexicanHatFitnessCase();
fitness_case.X1 = double.Parse(table.Rows[index]["X1"].ToString());
fitness_case.X2 = double.Parse(table.Rows[index]["X2"].ToString());
fitness_case.Y = double.Parse(table.Rows[index]["Y"].ToString());
return(fitness_case);
};
pop.GetFitnessCaseCount += () =>
{
return(table.Rows.Count);
};
pop.EvaluateObjectiveForSolution += (fitness_cases, solution, objective_index) =>
{
double fitness = 0;
for (int i = 0; i < fitness_cases.Count; i++)
{
MexicanHatFitnessCase fitness_case = (MexicanHatFitnessCase)fitness_cases[i];
double correct_y = fitness_case.Y;
double computed_y = fitness_case.ComputedY;
fitness += (correct_y - computed_y) * (correct_y - computed_y);
}
return(fitness);
};
pop.BreedInitialPopulation();
while (!pop.IsTerminated)
{
pop.Evolve();
Console.WriteLine("Mexican Hat Symbolic Regression Generation: {0}", pop.CurrentGeneration);
Console.WriteLine("Global Fitness: {0}\tCurrent Fitness: {1}", pop.GlobalBestProgram.Fitness.ToString("0.000"), pop.FindFittestProgramInCurrentGeneration().Fitness.ToString("0.000"));
Console.WriteLine("Global Best Solution:\n{0}", pop.GlobalBestProgram);
}
Console.WriteLine(pop.GlobalBestProgram.ToString());
}
static void Main(string[] args)
{
DataTable table = LoadData("dataset.txt");
TGPConfig config = new TGPConfig("TGPConfig.xml");
//config.CrossoverRate = 0.35;
//config.MicroMutationRate = 0.15;
//config.MacroMutationRate = 0.45;
//config.ReproductionRate = 0.05;
config.ElitismRatio = 0.1;
TGPPop <TGPSolution> pop = new TGPPop <TGPSolution>(config);
pop.PopulationReplacement = TGPPop <TGPSolution> .PopulationReplacementMode.TinyGP;
pop.OperatorSet.AddOperator(new TGPOperator_Plus());
pop.OperatorSet.AddOperator(new TGPOperator_Minus());
pop.OperatorSet.AddOperator(new TGPOperator_Division());
pop.OperatorSet.AddOperator(new TGPOperator_Multiplication());
pop.OperatorSet.AddOperator(new TGPOperator_Sin());
pop.OperatorSet.AddOperator(new TGPOperator_Cos());
pop.OperatorSet.AddIfgtOperator();
for (int i = 1; i < 10; ++i)
{
pop.ConstantSet.AddConstant(string.Format("C{0}", i), i);
}
pop.VariableSet.AddVariable("x");
pop.VariableSet.AddVariable("y");
pop.CreateFitnessCase += (index) =>
{
SpiralFitnessCase fitness_case = new SpiralFitnessCase();
fitness_case.X = double.Parse(table.Rows[index]["X"].ToString());
fitness_case.Y = double.Parse(table.Rows[index]["Y"].ToString());
fitness_case.Label = int.Parse(table.Rows[index]["Label"].ToString());
return(fitness_case);
};
pop.GetFitnessCaseCount += () =>
{
return(table.Rows.Count);
};
pop.EvaluateObjectiveForSolution += (fitness_cases, solution, objective_index) =>
{
double fitness = 0;
for (int i = 0; i < fitness_cases.Count; i++)
{
SpiralFitnessCase fitness_case = (SpiralFitnessCase)fitness_cases[i];
int correct_y = fitness_case.Label;
int computed_y = fitness_case.ComputedLabel;
fitness += (correct_y == computed_y) ? 0 : 1;
}
return(fitness);
};
pop.BreedInitialPopulation();
while (!pop.IsTerminated)
{
pop.Evolve();
Console.WriteLine("Spiral Classification Generation: {0}", pop.CurrentGeneration);
Console.WriteLine("Global Fitness: {0}\tCurrent Fitness: {1}", pop.GlobalBestProgram.Fitness, pop.FindFittestProgramInCurrentGeneration().Fitness);
Console.WriteLine("Global Best Solution:\n{0}", pop.GlobalBestProgram);
}
Console.WriteLine(pop.GlobalBestProgram.ToString());
}