/// <summary>
///
/// </summary>
/// <param name="args">Array of four mortality parameters (n, l, m, h)</param>
/// <param name="logToFile">If false, program only outputs to Console. If true, program outputs to file rather than Console</param>
private static void batchMode(string[] args)
{
// Get mortality parameters
MortalityParam mortMaternal, mortSiblingComp, mortSiblingHelp, mortChildInfluence, mortMotherInfluence;
StreamWriter sw = null;
// run backwards and forwards using passed mortality params (no user input required)
// population name
string population = args[0];
//
mortMaternal = parseArgs(args[1]);
mortSiblingComp = parseArgs(args[2]);
mortSiblingHelp = parseArgs(args[3]);
// Siler params
double silerA1 = double.Parse(args[4]);
double silerB1 = double.Parse(args[5]);
double silerA2 = double.Parse(args[6]);
double silerA3 = double.Parse(args[7]);
double silerB3 = double.Parse(args[8]);
string births = "";
string popMatrix = "";
string rFile = "";
if (args.Length == 10)
{
births = args[9];
}
else if (args.Length >= 11)
{
births = args[9];
popMatrix = args[10];
rFile = args[11];
}
// Return a string of mortality parameter values (e.g. medium, low, medium, high returns "MLMH")
string sp1 = Enum.GetName(typeof(MortalityParam), mortMaternal).Substring(0, 1);
string sp2 = Enum.GetName(typeof(MortalityParam), mortSiblingComp).Substring(0, 1);
string sp3 = Enum.GetName(typeof(MortalityParam), mortSiblingHelp).Substring(0, 1);
string p = population + "-" + sp1 + sp2 + sp3; // +sp4 + sp5;
// set up log file
string logFilename = "Log-" + p + ".log";
Console.WriteLine("(Running in batch mode)");
Console.WriteLine();
// create and initialise model
Console.WriteLine("-- initialising model");
BirthIntervals sdp = new BirthIntervals();
Console.WriteLine("-- model initialised");
// set parameters of model
sdp.population = population;
sdp.childbirthMortality = mortMaternal;
sdp.siblingCompetition = mortSiblingComp;
sdp.siblingHelp = mortSiblingHelp;
// set Siler params
sdp.silerA1 = silerA1;
sdp.silerB1 = silerB1;
sdp.silerA2 = silerA2;
sdp.silerA3 = silerA3;
sdp.silerB3 = silerB3;
Console.WriteLine(PrintParams(mortMaternal, mortSiblingComp, mortSiblingHelp));
if (args.Length < 10)
{
// calculate optimal birth intervals
Console.WriteLine("-- beginning backward iteration");
sdp.runBackward();
Console.WriteLine("-- backward iteration finished");
// output birth decisions and fitness here
Console.WriteLine("-- writing birth decisions to disk");
sdp.outputBirthDecisions();
sdp.outputGrowthRate();
Console.WriteLine("-- birth decisions written to disk");
Console.WriteLine("-- beginning forward iteration");
sdp.runForward();
Console.WriteLine("-- forward iteration finished");
Console.WriteLine("-- writing population data to disk");
sdp.outputGrowthRate();
Console.WriteLine("-- population data written to disk");
}
if (args.Length == 10)
{
Console.WriteLine("-- importing births from {0}", births);
sdp.importBirthDecisionsFromFile(births);
Console.WriteLine("-- birth decisions imported");
Console.WriteLine("-- beginning forward iteration");
sdp.runForward();
Console.WriteLine("-- forward iteration finished");
Console.WriteLine("-- writing population data to disk");
sdp.outputGrowthRate();
Console.WriteLine("-- population data written to disk");
}
else if (args.Length >= 11)
{
Console.WriteLine("-- importing births from {0}", births);
sdp.importBirthDecisionsFromFile(births);
Console.WriteLine("-- birth decisions imported");
Console.WriteLine("-- importing population matrix from {0}", popMatrix);
sdp.importPopulationFromFile(popMatrix);
Console.WriteLine("-- population matrix imported");
Console.WriteLine("-- importing population growth rate from {0}", rFile);
sdp.importPopulationGrowth(rFile);
Console.WriteLine("-- population growth rate imported");
Console.WriteLine("-- importing population matrix from {0}", popMatrix);
sdp.importPopulationFromFile(popMatrix);
Console.WriteLine("-- population matrix imported");
}
Console.WriteLine("-- calculating stats");
sdp.calcStats();
Console.WriteLine("-- calculated stats");
Console.WriteLine("-- writing stats data to disk");
sdp.outputStats();
Console.WriteLine("-- stats data written to disk");
}
private static void interactiveMode()
{
Console.WriteLine("(Running in interactive mode)");
Console.WriteLine();
Console.WriteLine("Please enter mortality parameters (h=high, m=medium, l=low, n=none):");
Console.Write("Risk of mortality during childbirth: ");
MortalityParam p1 = parseArgs(Console.ReadLine());
Console.Write("Intensity of sibling competition: ");
MortalityParam p2 = parseArgs(Console.ReadLine());
Console.Write("Level of juvenile help: ");
MortalityParam p3 = parseArgs(Console.ReadLine());
Console.WriteLine(PrintParams(p1, p2, p3));
Console.Write("Enter Siler parameter a1: ");
double silerA1 = double.Parse(Console.ReadLine());
Console.Write("Enter Siler parameter b1: ");
double silerB1 = double.Parse(Console.ReadLine());
Console.Write("Enter Siler parameter a2: ");
double silerA2 = double.Parse(Console.ReadLine());
Console.Write("Enter Siler parameter a3: ");
double silerA3 = double.Parse(Console.ReadLine());
Console.Write("Enter Siler parameter b3: ");
double silerB3 = double.Parse(Console.ReadLine());
bool done = false;
// create and initialise model
sdpInteractive = new BirthIntervals();
// set parameters of model
sdpInteractive.childbirthMortality = p1;
sdpInteractive.siblingCompetition = p2;
sdpInteractive.siblingHelp = p3;
// set Siler params
sdpInteractive.silerA1 = silerA1;
sdpInteractive.silerB1 = silerB1;
sdpInteractive.silerA2 = silerA2;
sdpInteractive.silerA3 = silerA3;
sdpInteractive.silerB3 = silerB3;
do
{
// Main menu
Console.WriteLine();
Console.WriteLine("1 = calculate optimal interbirth intervals (backward iteration)");
Console.WriteLine("2 = calculate population distributions (forward iteration)");
Console.WriteLine("3 = calculate stats");
Console.WriteLine("4 = full model run");
Console.Write("Please select an option (0 = quit): ");
string input = Console.ReadLine();
int option;
try
{
option = int.Parse(input);
}
catch (FormatException e)
{
Console.WriteLine("Invalid option. Try again");
continue;
}
if (option == 0)
{
done = true;
continue;
}
else if (option == 1)
{
back();
}
else if (option == 2)
{
forward();
}
else if (option == 3)
{
stats();
}
else if (option == 4)
{
back();
forward();
stats();
}
else
{
Console.WriteLine("{0} is an invalid option.", option);
}
} while (!done);
}
