private List<SimulationObject> GeneratePersons(SpatialZone currZone, int numPerson,
Person initAgent, bool warmUpStatus,
OutputFileWriter currWriter)
{
int seltdDim = 0;
List<ConditionalDistribution> condList = currZone.GetPersonDataCollectionsList();
var generatedAgents = new List<SimulationObject>();
Person prevAgent = initAgent;
ImportanceSampler currImpSampler = new ImportanceSampler();
int iter = 0;
if(warmUpStatus == true)
{
iter = Constants.WARMUP_ITERATIONS;
}
else
{
iter = Constants.SKIP_ITERATIONS * numPerson;
}
Person newAgent = new Person();
StringBuilder builder = new StringBuilder();
for(int i = 0; i < iter; i++)
{
seltdDim = randGen.NextInRange(0, condList.Count - 1);
DiscreteCondDistribution currDist =
(DiscreteCondDistribution)condList[seltdDim];
/*if (currDist.GetDimensionName() == "HouseholdSize2")
{
newAgent = (Person) currImpSampler.GetNextAgent(
currZone.myHhldSize2Marginal,
currDist, currDist.GetDimensionName(),
(SimulationObject) prevAgent, currZone);
}
else if (currDist.GetDimensionName() == "Age")
{
newAgent = (Person)currImpSampler.GetNextAgent(
currZone.myAgeMarginal,
currDist, currDist.GetDimensionName(),
(SimulationObject)prevAgent, currZone);
}
else*/
if(currDist.GetDimensionName() == "Sex")
{
newAgent = (Person)currImpSampler.GetNextAgent(
currZone.mySexMarginal,
currDist, currDist.GetDimensionName(),
(SimulationObject)prevAgent, currZone);
}
/*else if (currDist.GetDimensionName() == "EducationLevel")
{
newAgent = (Person)currImpSampler.GetNextAgent(
currZone.myEducationMarginal,
currDist, currDist.GetDimensionName(),
(SimulationObject)prevAgent, currZone);
}*/
else
{
List<KeyValPair> currComm = currDist.GetCommulativeValue(
prevAgent.GetNewJointKey(currDist.GetDimensionName())
, currZone);
newAgent = (Person)GenerateNextAgent(currComm,
(SimulationObject)prevAgent,
currDist.GetDimensionName());
}
prevAgent = newAgent;
if(warmUpStatus == false && (i % Constants.SKIP_ITERATIONS == 0))
{
//generatedAgents.Add(newAgent);
builder.Append((int)newAgent.GetAge()); builder.Append(',');
builder.Append(newAgent.GetZoneID()); builder.Append(',');
builder.Append((int)newAgent.GetSex()); builder.Append(',');
builder.Append((int)newAgent.GetHhldSize()); builder.Append(',');
builder.Append((int)newAgent.GetEducationLevel());
currWriter.WriteToFile(builder.ToString());
builder.Clear();
}
}
return generatedAgents;
}
private List<double> GetUtilityValuesForSex(Person person)
{
var currValues = new List<double>(2);
currValues.Add(1);
int twoAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren);
int twoPlusAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdults);
int twoPlusAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdultsChildren);
int threePlusCars = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.ThreeOrMore);
int lNo = Convert.ToInt32(person.GetDrivingLicense() == DrivingLicense.No);
int homeFT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.FullTimeHome);
int homePT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTimeHome);
int PT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTime);
int town = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Townhouse);
int manConst = Convert.ToInt32(person.GetOccupation() == Occupation.Construction);
int profMan = Convert.ToInt32(person.GetOccupation() == Occupation.TradesProfessional);
int retail = Convert.ToInt32(person.GetOccupation() == Occupation.TechnicalRetailSalesServiceNotEmployedUnknown);
int oNE = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.Unemployed);
int eleven = Convert.ToInt32(person.GetAge() == Age.ElevenToThirteen);
int sixteen = Convert.ToInt32(person.GetAge() == Age.SixteenToSeventeen);
int eighteen = Convert.ToInt32(person.GetAge() == Age.EighteenToTwentyFive);
int twentysix = Convert.ToInt32(person.GetAge() == Age.TwentySixToThirty);
int fortyone = Convert.ToInt32(person.GetAge() == Age.FortyOneToFifty);
int fiftyone = Convert.ToInt32(person.GetAge() == Age.FiftyOneToFiftyFour);
int fiftyfive = Convert.ToInt32(person.GetAge() == Age.FiftyfiveToSixtyFour);
int sixtyfivePlus = Convert.ToInt32(person.GetAge() == Age.MoreThanSixtyFive);
currValues.Add(Math.Exp(0.436 * 1.00 + -0.261 * twoAdultsChildren + -0.174 * twoPlusAdults
+ -0.222 * twoPlusAdultsChildren + 0.0472 * threePlusCars + 0.443 * lNo + -0.267 * homeFT
+ 0.530 * homePT + 0.900 * PT + 0.0508 * town + -2.33 * manConst + -0.733 * profMan
+ -0.824 * retail + -0.329 * oNE + -0.451 * eleven + -0.354 * sixteen + -0.154 * eighteen
+ 0.371 * twentysix + 0.428 * fortyone + 0.326 * fiftyone + 0.216 * fiftyfive + -0.0818 * sixtyfivePlus));
return currValues;
}
private List<double> GetUtilityValuesForOccupation(Person person, SpatialZone curZ)
{
var currValues = new List<double>(5);
currValues.Add(1);
int twoAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren);
int twoPlusAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdultsChildren);
int homeFT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.FullTimeHome);
int homePT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTimeHome);
int PT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTime);
int apart = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Apartment);
int eighteen = Convert.ToInt32(person.GetAge() == Age.EighteenToTwentyFive);
int twentysix = Convert.ToInt32(person.GetAge() == Age.TwentySixToThirty);
int fortyone = Convert.ToInt32(person.GetAge() == Age.FortyOneToFifty);
int fiftyone = Convert.ToInt32(person.GetAge() == Age.FiftyOneToFiftyFour);
int fiftyfive = Convert.ToInt32(person.GetAge() == Age.FiftyfiveToSixtyFour);
int thirtyone = Convert.ToInt32(person.GetAge() == Age.ThirtyOneToForty);
int threePlusCars = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.ThreeOrMore);
double z_manConst = curZ.GetPersonOccuMarginal().GetValue("1");
int oneAdultChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.OneAdultOneChild);
int twoAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults);
int twoPlusAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdults
|| person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults);
int oneCar = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.OneCar);
int twoCars = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.TwoCars);
currValues.Add(Math.Exp(-1.95 * 1.00 + -0.746 * oneAdultChildren + 0.280 * twoAdults + 0.594 * twoAdultsChildren
+ 0.556 * twoPlusAdults + 0.827 * twoPlusAdultsChildren + 0.247 * apart + 0.0751 * oneCar + 0.0867 * threePlusCars
+ 2.13 * homeFT + 0.237 * homePT + 0.841 * PT + -0.106 * eighteen + -0.123 * twentysix + -0.137 * thirtyone
+ 0.145 * fiftyone + 0.140 * fiftyfive + 16.1 * z_manConst));
int town = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Townhouse);
int fourteen = Convert.ToInt32(person.GetAge() == Age.FourteenToFifteen);
int sixteen = Convert.ToInt32(person.GetAge() == Age.SixteenToSeventeen);
int sixtyfivePlus = Convert.ToInt32(person.GetAge() == Age.MoreThanSixtyFive);
double z_profMan = curZ.GetPersonOccuMarginal().GetValue("2");
currValues.Add(Math.Exp(-0.648 * 1.00 + -0.0753 * twoAdults + 0.121 * twoAdultsChildren + -0.243 * twoPlusAdults
+ -0.110 * twoPlusAdultsChildren + -0.265 * apart + -0.147 * town + 0.271 * oneCar + 0.337 * twoCars + 0.303 * threePlusCars
+ 2.31 * homeFT + -1.02 * fourteen + -0.866 * sixteen + -0.389 * eighteen + 0.0738 * thirtyone + 0.229 * sixtyfivePlus
+ 7.22 * z_profMan));
int eleven = Convert.ToInt32(person.GetAge() == Age.ElevenToThirteen);
double z_retail = curZ.GetPersonOccuMarginal().GetValue("3");
currValues.Add(Math.Exp(-1.05 * 1.00 + 0.308 * oneAdultChildren + 0.205 * twoAdults + 0.479 * twoAdultsChildren + 0.379 * twoPlusAdults
+ 0.781 * twoPlusAdultsChildren + 0.220 * apart + 0.108 * town + -0.152 * oneCar + -0.245 * twoCars + -0.27 * threePlusCars
+ 2.85 * homeFT + 1.48 * homePT + 2.30 * PT + 3.90 * eleven + 2.20 * fourteen + 2.22 * sixteen + 0.507 * eighteen
+ 0.160 * twentysix + 0.0563 * fiftyone + 0.129 * fiftyfive + 0.398 * sixtyfivePlus + 6.32 * z_retail));
double z_oNE = curZ.GetPersonOccuMarginal().GetValue("4");
currValues.Add(Math.Exp(0.743 * 1.00 + 2.37 * oneAdultChildren + 0.298 * twoAdults + 2.47 * twoAdultsChildren + 0.437 * twoPlusAdults
+ 2.12 * twoPlusAdultsChildren + 0.114 * apart + -0.0672 * town + -0.63 * oneCar + -1.08 * twoCars + -1.33 * threePlusCars
+ 8.71 * eleven + 5.34 * fourteen + 3.87 * sixteen + 0.969 * eighteen + -0.842 * twentysix + -1.70 * thirtyone
+ -0.528 * fiftyone + 0.517 * fiftyfive + 3.08 * sixtyfivePlus + 0.854 * z_oNE));
return currValues;
}
private List<double> GetUtilityValuesForEmploymentStatus(Person person)
{
var currValues = new List<double>(5);
currValues.Add(1);
int twoAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren);
int twoPlusAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdults);
int threePlusCars = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.ThreeOrMore);
int lNo = Convert.ToInt32(person.GetDrivingLicense() == DrivingLicense.No);
int apart = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Apartment);
int town = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Townhouse);
int profMan = Convert.ToInt32(person.GetOccupation() == Occupation.TradesProfessional);
int retail = Convert.ToInt32(person.GetOccupation() == Occupation.TechnicalRetailSalesServiceNotEmployedUnknown);
int eighteen = Convert.ToInt32(person.GetAge() == Age.EighteenToTwentyFive);
int twentysix = Convert.ToInt32(person.GetAge() == Age.TwentySixToThirty);
int thirtyone = Convert.ToInt32(person.GetAge() == Age.ThirtyOneToForty);
int female = Convert.ToInt32(person.GetSex() == Sex.Female);
int twoAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults);
int oneCar = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.OneCar);
int twoPlusAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdultsChildren);
currValues.Add(Math.Exp(-2.34 * 1.00 - 0.122 * twoAdults + -0.335 * twoAdultsChildren + -0.453 * twoPlusAdults + -0.248 * twoPlusAdultsChildren
+ -0.236 * apart + -0.323 * town + 0.104 * oneCar + 0.126 * threePlusCars + 0.522 * profMan + 0.802 * retail
+ -1.32 * eighteen + -0.711 * twentysix + -0.325 * thirtyone + -0.252 * female + 0.247 * lNo));
int manConst = Convert.ToInt32(person.GetOccupation() == Occupation.Construction);
int eleven = Convert.ToInt32(person.GetAge() == Age.ElevenToThirteen);
int fiftyone = Convert.ToInt32(person.GetAge() == Age.FiftyOneToFiftyFour);
int sixteen = Convert.ToInt32(person.GetAge() == Age.SixteenToSeventeen);
currValues.Add(Math.Exp(-3.25 * 1.00 + -0.430 * twoAdultsChildren + -0.573 * twoPlusAdults + -0.456 * twoPlusAdultsChildren
+ -0.249 * town + 0.0881 * oneCar + -0.757 * manConst + 0.204 * profMan + 0.351 * retail + 8.78 * eleven + 3.50 * sixteen
+ -0.367 * eighteen + -1.06 * twentysix + -0.865 * thirtyone + -0.508 * fiftyone + 0.542 * female));
int fourteen = Convert.ToInt32(person.GetAge() == Age.FourteenToFifteen);
currValues.Add(Math.Exp(2.01 * 1.00 + -0.689 * twoPlusAdults + -0.0915 * twoPlusAdultsChildren + 0.0922 * apart
+ -0.240 * town + 0.398 * oneCar + -0.126 * threePlusCars + -17.7 * manConst + -17.9 * profMan + 14.6 * eleven
+ 5.15 * fourteen + 6.33 * sixteen + 0.646 * eighteen + -1.15 * twentysix + -1.45 * thirtyone + -1.31 * fiftyone
+ -0.280 * female + 2.15 * lNo + -19.6 * retail));
currValues.Add(Math.Exp(-2.69 * 1.00 + 0.116 * twoAdults + -0.232 * twoAdultsChildren + 0.0514 * twoPlusAdults + 0.188 * apart
+ -0.103 * town + 0.0953 * oneCar + -0.164 * threePlusCars + -0.194 * manConst + 0.0903 * profMan + 0.914 * retail
+ 3.93 * fourteen + 5.99 * sixteen + 1.74 * eighteen + -0.317 * thirtyone + -0.392 * fiftyone + 0.853 * female
+ 0.462 * lNo + 7.15 * eleven));
return currValues;
}
private List<double> GetUtilityValuesForDriverLicense(Person person)
{
var currValues = new List<double>(2);
currValues.Add(1);
int twoAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren);
int twoPlusAdultsChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.TwoAdultsChildren
|| person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdultsChildren);
int homeFT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.FullTimeHome);
int homePT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTimeHome);
int PT = Convert.ToInt32(person.GetEmploymentStatus() == EmploymentStatus.PartTime);
int apart = Convert.ToInt32(person.GetHhld().GetDwellingType() == DwellingType.Apartment);
int manConst = Convert.ToInt32(person.GetOccupation() == Occupation.Construction);
int retail = Convert.ToInt32(person.GetOccupation() == Occupation.TechnicalRetailSalesServiceNotEmployedUnknown);
int sixteen = Convert.ToInt32(person.GetAge() == Age.SixteenToSeventeen);
int eighteen = Convert.ToInt32(person.GetAge() == Age.EighteenToTwentyFive);
int twentysix = Convert.ToInt32(person.GetAge() == Age.TwentySixToThirty);
int fortyone = Convert.ToInt32(person.GetAge() == Age.FortyOneToFifty);
int oneAdultChildren = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.OneAdultOneChild);
int twoAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults);
int twoPlusAdults = Convert.ToInt32(person.GetHhld().GetHhldSize() == HouseholdSize.ThreeOrMoreAdults
|| person.GetHhld().GetHhldSize() == HouseholdSize.Twoadults);
int oneCar = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.OneCar);
int twoCars = Convert.ToInt32(person.GetHhld().GetNumOfCars() == NumOfCars.TwoCars);
currValues.Add(Math.Exp(-0.833 * 1.00 + 2.17 * oneAdultChildren + -0.103 * twoAdults + 1.85 * twoAdultsChildren
+ 1.02 * twoPlusAdults + 2.17 * twoPlusAdultsChildren + 0.693 * apart + -0.392 * oneCar + -0.743 * twoCars
+ -1.35 * homeFT + -1.26 * homePT + -0.220 * PT + -1.95 * manConst + -1.22 * retail + 0.811 * sixteen
+ -0.621 * eighteen + -1.52 * twentysix + -1.80 * fortyone));
return currValues;
}
private void updateCounts(Person person)
{
{
switch (person.GetAge())
{
case (Age.LessThanEleven):
{
count_less11++;
break;
}
case (Age.ElevenToThirteen):
{
count_11++;
break;
}
case (Age.FourteenToFifteen):
{
count_14++;
break;
}
case (Age.SixteenToSeventeen):
{
count_16++;
break;
}
case (Age.EighteenToTwentyFive):
{
count_18++;
break;
}
case (Age.TwentySixToThirty):
{
count_26++;
break;
}
case (Age.ThirtyOneToForty):
{
count_31++;
break;
}
case (Age.FortyOneToFifty):
{
count_41++;
break;
}
case (Age.FiftyOneToFiftyFour):
{
count_51++;
break;
}
case (Age.FiftyfiveToSixtyFour):
{
count_55++;
break;
}
case (Age.MoreThanSixtyFive):
{
count_65++;
break;
}
default:
break;
}
switch (person.GetEmploymentStatus())
{
case (EmploymentStatus.PartTime):
count_part_time++;
break;
case (EmploymentStatus.FullTimeHome):
count_full_time_home++;
break;
case (EmploymentStatus.FullTime):
count_full_time++;
break;
case (EmploymentStatus.PartTimeHome):
count_part_time_home++;
break;
case (EmploymentStatus.Unemployed):
count_unemployed++;
break;
default:
break;
}
switch (person.GetSex())
{
case (Sex.Male):
count_male++;
break;
case (Sex.Female):
count_female++;
break;
default:
break;
}
switch (person.GetDrivingLicense())
{
case (DrivingLicense.Yes):
count_driving_license_yes++;
break;
case (DrivingLicense.No):
count_driving_license_no++;
break;
default:
break;
}
switch (person.GetOccupation())
{
case (Occupation.GeneralOffice):
{
count_general_office++;
break;
}
case (Occupation.ClericalManufacturing):
{
count_clerical_manufacturing++;
break;
}
case (Occupation.TradesProfessional):
{
count_professional_man++;
break;
}
case (Occupation.TechnicalRetailSalesServiceNotEmployedUnknown):
{
count_retail++;
break;
}
case (Occupation.Construction):
{
count_costruction_man++;
break;
}
}
}
}
