using System; using L=Science.Biology.PopulationGenetics; namespace ScienceTest.BiologyTest.PopulationGeneticsTest { /// /// One gene is located in chromosome X. /// public class ColorBlind { private string result; public string Result { get{return result;} } public ColorBlind() { } public void Compute() { L.Chromosome[] Ac = new L.Chromosome[22]; Ac[0] = new L.Chromosome(); Ac[1] = new L.Chromosome(); Ac[2] = new L.Chromosome(); Ac[3] = new L.Chromosome(); Ac[4] = new L.Chromosome(); Ac[5] = new L.Chromosome(); Ac[6] = new L.Chromosome(); Ac[7] = new L.Chromosome(); Ac[8] = new L.Chromosome(); Ac[9] = new L.Chromosome(); Ac[10] = new L.Chromosome(); Ac[11] = new L.Chromosome(); Ac[12] = new L.Chromosome(); Ac[13] = new L.Chromosome(); Ac[14] = new L.Chromosome(); Ac[15] = new L.Chromosome(); Ac[16] = new L.Chromosome(); Ac[17] = new L.Chromosome(); Ac[18] = new L.Chromosome(); Ac[19] = new L.Chromosome(); Ac[20] = new L.Chromosome(); Ac[21] = new L.Chromosome(); // give one gene L.Gene g0 = new L.Gene(); g0.NumberOfAllele = 2; g0.StartBase = 1200; g0.EndBase = 1530; g0.FindMutation(); g0.From = 1; g0.To = 0; g0.ModifiedValue = 0.0; g0.SetMutation(); L.Chromosome X = new L.Chromosome(g0); // end of giving gene L.Chromosome Y = new L.Chromosome(); L.Mitochondria mito = new L.Mitochondria(); L.SocialGene[] sg = new L.SocialGene[1]; sg[0] = new L.SocialGene(); sg[0].NumberOfAllele = 1; L.ChromosomePair[] Cp = new L.ChromosomePair[22]; for(int i = 0; i < 22; i++) Cp[i] = new L.ChromosomePair(Ac[i]); L.ChromosomePair XX = new L.ChromosomePair(X); L.GenotypeFemale fg = new L.GenotypeFemale(Cp,XX,mito,sg); L.GenotypeMale mg = new L.GenotypeMale(Cp,X,Y,mito,sg); L.Mating mating = new L.Mating(); L.Mutation mut = new L.Mutation(fg,mg); L.Reproduction reprod = new L.Reproduction(fg,mg); L.Selection select = new L.Selection(fg,mg); select.WhichFemale = 2; select.ModifiedValue = 0.0012; select.SetDisadvantageXX(); select.WhichMale = 1; select.ModifiedValue = 0.0012; select.SetDisadvantageXY(); L.PopulationEquation eq = new L.PopulationEquation(mating, mut, reprod, select); eq.HowManyTimes = 1000000; eq.Compute(); foreach(double k in eq.AdultFemalePopulation) result += Convert.ToString(k)+"\r\n"; result += "------\r\n"; foreach(double k in eq.AdultMalePopulation) result += Convert.ToString(k)+"\r\n"; result += "------\r\n"; result += eq.InfantMortality(); } } }