https://github.com/chen0040/cs-grammatical-evolution
Grammatical evolution implemented using C#
https://github.com/chen0040/cs-grammatical-evolution
grammatical-evolution symbolic-regression
Last synced: 10 months ago
JSON representation
Grammatical evolution implemented using C#
- Host: GitHub
- URL: https://github.com/chen0040/cs-grammatical-evolution
- Owner: chen0040
- License: mit
- Created: 2017-10-31T01:11:12.000Z (over 8 years ago)
- Default Branch: master
- Last Pushed: 2017-12-03T01:44:56.000Z (over 8 years ago)
- Last Synced: 2025-08-27T19:28:44.244Z (11 months ago)
- Topics: grammatical-evolution, symbolic-regression
- Language: C#
- Size: 641 KB
- Stars: 4
- Watchers: 1
- Forks: 1
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE
Awesome Lists containing this project
README
# cs-grammatical-evolution
Grammatical evolution implemented using C#
# Install
```bash
Install-Package cs-grammatical-evolution -Version 1.0.4
```
# Usage
The sample code belows show how to use the Grammatical Evolution to solve the spiral classification problem:
```cs
class Program
{
static DataTable LoadData(string filename)
{
DataTable table = new DataTable();
table.Columns.Add("X");
table.Columns.Add("Y");
table.Columns.Add("Label");
int line_count = 0;
using (StreamReader reader = new StreamReader(filename))
{
string line = reader.ReadLine();
int.TryParse(line, out line_count);
while ((line = reader.ReadLine()) != null)
{
string[] elements = line.Split(new char[] { '\t' });
double x, y;
int label;
double.TryParse(elements[0].Trim(), out x);
double.TryParse(elements[1].Trim(), out y);
int.TryParse(elements[2].Trim(), out label);
table.Rows.Add(x, y, label);
}
}
return table;
}
static void CreateGrammar(GEPop pop)
{
pop.GERuleSet.RemoveAllRules();
//pop.GERuleSet.AddRule("", "", "", "", "", "");
pop.GERuleSet.AddRule("", "", "", "", "", "");
pop.GERuleSet.AddRule("", "", "", "", "");
//pop.GERuleSet.AddRule("", "", "", "", "", "", "");
pop.GERuleSet.AddRule("", "");
pop.GERuleSet.AddRule("", "+");
pop.GERuleSet.AddRule("", "-");
pop.GERuleSet.AddRule("", "/");
pop.GERuleSet.AddRule("", "*");
//pop.GERuleSet.AddRule("", ",");
//pop.GERuleSet.AddRule("", "Pow");
pop.GERuleSet.AddRule("", "Sin");
//pop.GERuleSet.AddRule("", "Tan");
pop.GERuleSet.AddRule("", "Cos");
pop.GERuleSet.AddRule("", "(");
pop.GERuleSet.AddRule("", ")");
pop.GERuleSet.AddRule("", "X");
pop.GERuleSet.AddRule("", "Y");
pop.GERuleSet.AddRule("", "C1");
}
static void CreateTerminalSet(GEPop pop)
{
pop.OperatorSet.AddOperator("+");
pop.OperatorSet.AddOperator("-");
pop.OperatorSet.AddOperator("/");
pop.OperatorSet.AddOperator("*");
pop.OperatorSet.AddOperator("Sin");
pop.OperatorSet.AddOperator("Cos");
pop.OperatorSet.AddOperator("(");
pop.OperatorSet.AddOperator(")");
//pop.OperatorSet.AddOperator("Tan");
//pop.OperatorSet.AddOperator("Pow");
pop.ConstantSet.AddConstant("C1", 8);
pop.VariableSet.AddVariable("X");
pop.VariableSet.AddVariable("Y");
pop.BuildTerminalSet();
}
static void Main(string[] args)
{
DataTable table = LoadData("dataset.txt");
GEConfig config = new GEConfig("GEConfig.xml");
GEPop pop = new GEPop(config);
CreateTerminalSet(pop);
//User can load BNF grammar from file or create programmatic ally
pop.GERuleSet.LoadGrammar("grammar.bnf");
Console.WriteLine(pop.GERuleSet.ToString());
//CreateGrammar(pop);
//Xianshun: Say if you implement a operator named "Move(X, Y)",
//Which returns a value 1 if X > Y otherwise 0
//use the following code to enable this use-defined operator in the BNF
/*
pop.EvaluateFunction += (function_name, function_args) =>
{
if (function_name == "Move")
{
double X = Convert.ToDouble(function_args.Parameters[0].Evaluate());
double Y = Convert.ToDouble(function_args.Parameters[1].Evaluate());
if (X > Y)
{
function_args.Result = 1;
}
else
{
function_args.Result = 0;
}
}
};
*/
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, s, 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("Current Best Solution:\n{0}", pop.FindFittestProgramInCurrentGeneration());
}
Console.WriteLine(pop.GlobalBestProgram.ToString());
}
}
```
The GEConfig.xml and its child configuration files will be automatically generated if they do not exist, otherwise the configuration will be loaded from the existing GEConfig.xml and its child configuration files.