public double ParentalSetProbability(int doubleDominant, int dominantHetero, int dominantRecessive, int doubleHetero, int heteroRecessive, int doubleRecessive)
{
var model = new ExpectedOffspringModel(doubleDominant, doubleHetero, dominantRecessive, doubleHetero, heteroRecessive, doubleRecessive);
var parentalProbabilitySum = model.Parents.Values.Sum();
return(Math.Round(parentalProbabilitySum, Constants.RoundingPrecision));
}
public double ExpectedOffspring_Tests(int doubleDominant, int dominantHetero, int dominantRecessive, int doubleHetero, int heteroRecessive,
int doubleRecessive, Genotype genotype, float offspringPerPair)
{
var model = new ExpectedOffspringModel(doubleDominant, dominantHetero, dominantRecessive, doubleHetero, heteroRecessive, doubleRecessive);
return(model.ExpectedOffspring(genotype, offspringPerPair));
}
static void Main()
{
const string parent = @"C:\Users\rianjs\Desktop";
const string filename = "overlap-graphs3.txt";
//var model = new ExpectedOffspringModel(19843, 16233, 18989, 19312, 16213, 17310); //155311
var model = new ExpectedOffspringModel(1, 0, 0, 1, 0, 1); //3.5
var pSum = model.Parents.Values.Sum();
Console.WriteLine("Probability sum: {0}", pSum); //OK
var hetero = model.ExpectedOffspring(Genotype.Heterozygous, 2);
var dominant = model.ExpectedOffspring(Genotype.Dominant, 2);
var recessive = model.ExpectedOffspring(Genotype.Recessive, 2);
var answer = hetero + dominant;
Console.WriteLine("Direct answer: {0:N}", answer);
Console.ReadLine();
}
public int Constructor_Tests(int doubleDominant, int dominantHetero, int dominantRecessive, int doubleHetero, int heteroRecessive, int doubleRecessive)
{
var model = new ExpectedOffspringModel(doubleDominant, dominantHetero, dominantRecessive, doubleHetero, heteroRecessive, doubleRecessive);
return(0);
}
