public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
}
private bool Equals(ISampleItem other) {
if (ReferenceEquals(null, other)) {
return false;
}
if (ReferenceEquals(this, other)) {
return true;
}
return other.ParcelId == ParcelId;
}
Task IDataService.SaveItem(ISampleItem item) => MethodRunner.RunAsync(() =>
{
var dto = _sampleMapper.MapToSampleItemDto(item);
if (item.IsNew)
{
_sampleProvider.CreateItem(dto);
}
else
{
_sampleProvider.UpdateItem(dto);
}
});
private bool Equals(ISampleItem other)
{
if (ReferenceEquals(null, other))
{
return(false);
}
if (ReferenceEquals(this, other))
{
return(true);
}
return(other.ParcelId == ParcelId);
}
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
ChoiceProbabilityCalculator.Alternative alternative = sampleItem.Alternative;
if (!alternative.Available)
{
return;
}
IHouseholdWrapper household = _tour.Household;
IPersonWrapper person = _tour.Person;
//var personDay = _tour.PersonDay;
bool householdHasChildren = household.HasChildren;
bool householdHasNoChildren = householdHasChildren ? false : true;
IParcelWrapper destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
int jointTourFlag = (_tour.JointTourSequence > 0).ToFlag();
double fastestTravelTime =
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, _tour.OriginParcel, destinationParcel).Variable +
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, destinationParcel, _tour.OriginParcel).Variable;
if (fastestTravelTime >= _maxAvailableMinutes)
{
alternative.Available = false;
return;
}
alternative.Choice = destinationParcel;
//double tourLogsum = 0;
double jointTourLogsum = 0;
double workTourLogsum = 0;
double otherTourLogsum = 0;
double subtourLogsum = 0;
if (_tour.IsHomeBasedTour)
{
if (_tour.DestinationPurpose == Global.Settings.Purposes.Work)
{
//var destinationArrivalTime = ChoiceModelUtility.GetDestinationArrivalTime(Global.Settings.Models.WorkTourModeModel);
//var destinationDepartureTime = ChoiceModelUtility.GetDestinationDepartureTime(Global.Settings.Models.WorkTourModeModel);
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeTimeModel>().RunNested(_tour, destinationParcel, _tour.Household.VehiclesAvailable, _tour.Person.GetTransitFareDiscountFraction());
workTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
// else if (_tour.DestinationPurpose == Global.Settings.Purposes.Escort) {
// var nestedAlternative = (Global.ChoiceModelDictionary.Get("HEscortTourModeModel") as HEscortTourModeModel).RunNested(_tour, destinationParcel);
// tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
// }
else
{
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <OtherHomeBasedTourModeTimeModel>().RunNested(_tour, destinationParcel, _tour.Household.VehiclesAvailable, _tour.Person.GetTransitFareDiscountFraction());
if (_tour.JointTourSequence > 0)
{
jointTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
else
{
otherTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
}
}
else
{
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <WorkBasedSubtourModeTimeModel>().RunNested(_tour, destinationParcel, _tour.Household.VehiclesAvailable, _tour.Person.GetTransitFareDiscountFraction());
subtourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
//var purpose = _tour.TourPurposeSegment;
//var carOwnership = person.CarOwnershipSegment;
//var votSegment = _tour.VotALSegment;
//var transitAccess = destinationParcel.TransitAccessSegment();
//var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][purpose][carOwnership][votSegment][transitAccess];
//var aggregateLogsumHomeBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][transitAccess];
//var aggregateLogsumWorkBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.Work_BASED][carOwnership][votSegment][transitAccess];
double distanceFromOrigin = _tour.OriginParcel.DistanceFromOrigin(destinationParcel, _tour.DestinationArrivalTime);
double piecewiseDistanceFrom5To10Miles = Math.Max(0, Math.Min(distanceFromOrigin - .5, 1 - .5));
double piecewiseDistanceFrom10MilesToInfinity = Math.Max(0, distanceFromOrigin - 1);
double piecewiseDistanceFrom0To1Mile = Math.Min(distanceFromOrigin, .10);
double piecewiseDistanceFrom1To5Miles = Math.Max(0, Math.Min(distanceFromOrigin - .1, .5 - .1));
double piecewiseDistanceFrom1To3AndAHalfMiles = Math.Max(0, Math.Min(distanceFromOrigin - .1, .35 - .1));
double piecewiseDistanceFrom3AndAHalfTo10Miles = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
double distanceFromOriginLog = Math.Log(1 + distanceFromOrigin);
//var distanceFromWorkLog = person.UsualWorkParcel.DistanceFromWorkLog(destinationParcel, 1);
double distanceFromSchoolLog = person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
double timePressure = Math.Log(1 - fastestTravelTime / _maxAvailableMinutes);
// log transforms of buffers for Neighborhood effects
double logOfOnePlusEmploymentEducationBuffer1 = Math.Log(destinationParcel.EmploymentEducationBuffer1 + 1.0);
double logOfOnePlusEmploymentFoodBuffer1 = Math.Log(destinationParcel.EmploymentFoodBuffer1 + 1.0);
double logOfOnePlusEmploymentGovernmentBuffer1 = Math.Log(destinationParcel.EmploymentGovernmentBuffer1 + 1.0);
double logOfOnePlusEmploymentOfficeBuffer1 = Math.Log(destinationParcel.EmploymentOfficeBuffer1 + 1.0);
double logOfOnePlusEmploymentRetailBuffer1 = Math.Log(destinationParcel.EmploymentRetailBuffer1 + 1.0);
double logOfOnePlusEmploymentServiceBuffer1 = Math.Log(destinationParcel.EmploymentServiceBuffer1 + 1.0);
double logOfOnePlusEmploymentMedicalBuffer1 = Math.Log(destinationParcel.EmploymentMedicalBuffer1 + 1.0);
double logOfOnePlusEmploymentIndustrial_Ag_ConstructionBuffer1 = Math.Log(destinationParcel.EmploymentIndustrialBuffer1 + destinationParcel.EmploymentAgricultureConstructionBuffer1 + 1.0);
double logOfOnePlusEmploymentTotalBuffer1 = Math.Log(destinationParcel.EmploymentTotalBuffer1 + 1.0);
double logOfOnePlusHouseholdsBuffer1 = Math.Log(destinationParcel.HouseholdsBuffer1 + 1.0);
double logOfOnePlusStudentsK12Buffer1 = Math.Log(destinationParcel.StudentsK8Buffer1 + destinationParcel.StudentsHighSchoolBuffer1 + 1.0);
double logOfOnePlusStudentsUniversityBuffer1 = Math.Log(destinationParcel.StudentsUniversityBuffer1 + 1.0);
// var EMPHOU_B = Math.Log(destinationParcel.EmploymentTotalBuffer1 + destinationParcel.HouseholdsBuffer1 + 1.0);
//var logOfOnePlusParkingOffStreetDailySpacesBuffer1 = Math.Log(1 + destinationParcel.ParkingOffStreetPaidDailySpacesBuffer1);
// connectivity attributes
//var c34Ratio = destinationParcel.C34RatioBuffer1();
//var carCompetitionFlag = FlagUtility.GetCarCompetitionFlag(carOwnership); // exludes no cars
//var noCarCompetitionFlag = FlagUtility.GetNoCarCompetitionFlag(carOwnership);
//var noCarsFlag = FlagUtility.GetNoCarsFlag(carOwnership);
alternative.AddUtilityTerm(2, sampleItem.AdjustmentFactor);
alternative.AddUtilityTerm(3, workTourLogsum);
alternative.AddUtilityTerm(4, otherTourLogsum);
alternative.AddUtilityTerm(5, jointTourLogsum);
alternative.AddUtilityTerm(6, subtourLogsum);
//subpopulation-specific terms
//alternative.AddUtilityTerm(260, _secondaryFlag * _workOrSchoolPatternFlag * piecewiseDistanceFrom0To1Mile);
//alternative.AddUtilityTerm(261, _secondaryFlag * _workOrSchoolPatternFlag * piecewiseDistanceFrom1To5Miles);
//alternative.AddUtilityTerm(262, _secondaryFlag * _workOrSchoolPatternFlag * piecewiseDistanceFrom5To10Miles);
//alternative.AddUtilityTerm(263, _secondaryFlag * _workOrSchoolPatternFlag * piecewiseDistanceFrom10MilesToInfinity);
alternative.AddUtilityTerm(260, _secondaryFlag * _workOrSchoolPatternFlag * distanceFromOriginLog);
//alternative.AddUtilityTerm(264, _secondaryFlag * _otherPatternFlag * piecewiseDistanceFrom0To1Mile);
//alternative.AddUtilityTerm(265, _secondaryFlag * _otherPatternFlag * piecewiseDistanceFrom1To5Miles);
//alternative.AddUtilityTerm(266, _secondaryFlag * _otherPatternFlag * piecewiseDistanceFrom5To10Miles);
//alternative.AddUtilityTerm(267, _secondaryFlag * _otherPatternFlag * piecewiseDistanceFrom10MilesToInfinity);
alternative.AddUtilityTerm(264, _secondaryFlag * _otherPatternFlag * distanceFromOriginLog);
alternative.AddUtilityTerm(268, (!_tour.IsHomeBasedTour).ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(269, household.Has0To15KIncome.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(270, household.HasMissingIncome.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(271, person.IsRetiredAdult.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(272, person.IsUniversityStudent.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(273, person.IsChildAge5Through15.ToFlag() * distanceFromOriginLog);
//alternative.AddUtilityTerm(274, person.IsChildUnder5.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(275, (_tour.IsHomeBasedTour).ToFlag() * timePressure);
alternative.AddUtilityTerm(276, (_tour.IsHomeBasedTour).ToFlag() * distanceFromSchoolLog);
//alternative.AddUtilityTerm(14, distanceFromWorkLog);
//alternative.AddUtilityTerm(277, (carCompetitionFlag + noCarCompetitionFlag) * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
//alternative.AddUtilityTerm(278, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
//alternative.AddUtilityTerm(279, carCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
//alternative.AddUtilityTerm(280, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
//alternative.AddUtilityTerm(281, noCarsFlag * c34Ratio);
//alternative.AddUtilityTerm(282, noCarCompetitionFlag * c34Ratio);
//alternative.AddUtilityTerm(283, (carCompetitionFlag + noCarCompetitionFlag) * logOfOnePlusParkingOffStreetDailySpacesBuffer1);
alternative.AddUtilityTerm(284, jointTourFlag * piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(285, jointTourFlag * piecewiseDistanceFrom1To5Miles);
alternative.AddUtilityTerm(286, jointTourFlag * piecewiseDistanceFrom5To10Miles);
alternative.AddUtilityTerm(287, jointTourFlag * piecewiseDistanceFrom10MilesToInfinity);
if (_tour.DestinationPurpose == Global.Settings.Purposes.Work)
{
alternative.AddUtilityTerm(10, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(11, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(12, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(13, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(14, aggregateLogsumWorkBased);
// Neighborhood
alternative.AddUtilityTerm(20, logOfOnePlusEmploymentEducationBuffer1);
//alternative.AddUtilityTerm(21, logOfOnePlusEmploymentGovernmentBuffer1);
alternative.AddUtilityTerm(22, logOfOnePlusEmploymentOfficeBuffer1);
//alternative.AddUtilityTerm(23, logOfOnePlusEmploymentServiceBuffer1);
alternative.AddUtilityTerm(24, logOfOnePlusEmploymentMedicalBuffer1);
alternative.AddUtilityTerm(25, logOfOnePlusHouseholdsBuffer1);
alternative.AddUtilityTerm(26, logOfOnePlusStudentsUniversityBuffer1);
alternative.AddUtilityTerm(27, logOfOnePlusStudentsK12Buffer1);
alternative.AddUtilityTerm(28, logOfOnePlusEmploymentIndustrial_Ag_ConstructionBuffer1);
// Size terms
alternative.AddUtilityTerm(30, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(31, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(32, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(33, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(34, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(35, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(36, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(37, destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(38, destinationParcel.Households);
alternative.AddUtilityTerm(39, destinationParcel.GetStudentsK12());
//(for application) alternative.AddUtilityTerm(40, destinationParcel.StudentsUniversity);
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Escort)
{
alternative.AddUtilityTerm(50, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(51, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(52, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(53, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(54, aggregateLogsumHomeBased);
// Neighborhood
alternative.AddUtilityTerm(60, householdHasNoChildren.ToFlag() * logOfOnePlusEmploymentGovernmentBuffer1);
alternative.AddUtilityTerm(61, householdHasNoChildren.ToFlag() * logOfOnePlusHouseholdsBuffer1);
//alternative.AddUtilityTerm(62, householdHasChildren.ToFlag() * logOfOnePlusHouseholdsBuffer1);
alternative.AddUtilityTerm(63, householdHasChildren.ToFlag() * logOfOnePlusStudentsK12Buffer1);
alternative.AddUtilityTerm(64, logOfOnePlusEmploymentTotalBuffer1);
// Size terms
alternative.AddUtilityTerm(70, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(71, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(72, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(73, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(74, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(75, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(76, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(77, (!householdHasChildren).ToFlag() * (destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction));
alternative.AddUtilityTerm(78, (!householdHasChildren).ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(79, (!householdHasChildren).ToFlag() * destinationParcel.GetStudentsK12());
alternative.AddUtilityTerm(80, householdHasChildren.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(81, householdHasChildren.ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(82, householdHasChildren.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(83, householdHasChildren.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(84, householdHasChildren.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(85, householdHasChildren.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(86, householdHasChildren.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(87, householdHasChildren.ToFlag() * (destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction));
alternative.AddUtilityTerm(88, householdHasChildren.ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(89, householdHasChildren.ToFlag() * destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.PersonalBusiness)
{
alternative.AddUtilityTerm(90, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(91, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(92, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(93, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(94, aggregateLogsumHomeBased);
// Neighborhood
alternative.AddUtilityTerm(100, logOfOnePlusEmploymentEducationBuffer1);
alternative.AddUtilityTerm(101, logOfOnePlusEmploymentOfficeBuffer1);
//alternative.AddUtilityTerm(102, logOfOnePlusEmploymentServiceBuffer1);
//alternative.AddUtilityTerm(103, logOfOnePlusEmploymentMedicalBuffer1);
alternative.AddUtilityTerm(104, logOfOnePlusHouseholdsBuffer1);
//alternative.AddUtilityTerm(105, logOfOnePlusStudentsUniversityBuffer1);
alternative.AddUtilityTerm(106, logOfOnePlusEmploymentGovernmentBuffer1);
alternative.AddUtilityTerm(107, logOfOnePlusEmploymentRetailBuffer1);
// Size terms
alternative.AddUtilityTerm(110, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(111, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(112, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(113, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(114, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(115, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(116, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(117, destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(118, destinationParcel.Households);
//(for application) alternative.AddUtilityTerm(119, destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Shopping)
{
alternative.AddUtilityTerm(120, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(121, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(122, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(123, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(124, aggregateLogsumHomeBased);
// Neighborhood
alternative.AddUtilityTerm(130, logOfOnePlusEmploymentEducationBuffer1);
alternative.AddUtilityTerm(131, logOfOnePlusEmploymentRetailBuffer1);
// Size terms
alternative.AddUtilityTerm(140, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(141, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(142, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(143, destinationParcel.EmploymentService);
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Meal)
{
alternative.AddUtilityTerm(150, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(151, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(152, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(153, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(154, aggregateLogsumHomeBased);
// Neighborhood
alternative.AddUtilityTerm(156, logOfOnePlusEmploymentFoodBuffer1);
//alternative.AddUtilityTerm(157, logOfOnePlusEmploymentRetailBuffer1);
alternative.AddUtilityTerm(158, logOfOnePlusEmploymentServiceBuffer1);
// Size terms
alternative.AddUtilityTerm(160, destinationParcel.EmploymentFood);
//(for application) alternative.AddUtilityTerm(161, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(162, destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(163, destinationParcel.Households);
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Social)
{
alternative.AddUtilityTerm(170, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(171, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(172, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(173, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(174, aggregateLogsumHomeBased);
// Neighborhood
//alternative.AddUtilityTerm(180, logOfOnePlusEmploymentOfficeBuffer1);
//alternative.AddUtilityTerm(181, logOfOnePlusEmploymentServiceBuffer1);
//alternative.AddUtilityTerm(182, logOfOnePlusHouseholdsBuffer1);
//alternative.AddUtilityTerm(183, logOfOnePlusStudentsK12Buffer1);
//alternative.AddUtilityTerm(184, logOfOnePlusStudentsUniversityBuffer1);
//alternative.AddUtilityTerm(185, logOfOnePlusEmploymentTotalBuffer1);
// Size terms
alternative.AddUtilityTerm(190, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(191, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(192, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(193, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(194, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(195, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(196, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(197, destinationParcel.Households);
alternative.AddUtilityTerm(198, destinationParcel.StudentsUniversity);
//(for application) alternative.AddUtilityTerm(199, destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Recreation)
{
alternative.AddUtilityTerm(200, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(201, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(202, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(203, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(204, aggregateLogsumHomeBased);
// Neighborhood
alternative.AddUtilityTerm(210, logOfOnePlusEmploymentOfficeBuffer1);
alternative.AddUtilityTerm(211, logOfOnePlusEmploymentServiceBuffer1);
alternative.AddUtilityTerm(212, logOfOnePlusHouseholdsBuffer1);
//alternative.AddUtilityTerm(213, logOfOnePlusStudentsK12Buffer1);
alternative.AddUtilityTerm(214, logOfOnePlusStudentsUniversityBuffer1);
//alternative.AddUtilityTerm(215, logOfOnePlusEmploymentTotalBuffer1);
// Size terms
alternative.AddUtilityTerm(220, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(221, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(222, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(223, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(224, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(225, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(226, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(227, destinationParcel.Households);
alternative.AddUtilityTerm(228, destinationParcel.StudentsUniversity);
alternative.AddUtilityTerm(229, destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Medical)
{
alternative.AddUtilityTerm(230, piecewiseDistanceFrom0To1Mile);
alternative.AddUtilityTerm(231, piecewiseDistanceFrom1To3AndAHalfMiles);
alternative.AddUtilityTerm(232, piecewiseDistanceFrom3AndAHalfTo10Miles);
alternative.AddUtilityTerm(233, piecewiseDistanceFrom10MilesToInfinity);
//alternative.AddUtilityTerm(234, aggregateLogsumHomeBased);
// Neighborhood
//alternative.AddUtilityTerm(240, logOfOnePlusEmploymentEducationBuffer1);
//alternative.AddUtilityTerm(241, logOfOnePlusEmploymentOfficeBuffer1);
//alternative.AddUtilityTerm(242, logOfOnePlusEmploymentServiceBuffer1);
//alternative.AddUtilityTerm(243, logOfOnePlusEmploymentMedicalBuffer1);
alternative.AddUtilityTerm(244, logOfOnePlusHouseholdsBuffer1);
//alternative.AddUtilityTerm(245, logOfOnePlusStudentsUniversityBuffer1);
//alternative.AddUtilityTerm(246, logOfOnePlusEmploymentGovernmentBuffer1);
alternative.AddUtilityTerm(247, logOfOnePlusEmploymentRetailBuffer1);
// Size terms
//(for application) alternative.AddUtilityTerm(250, destinationParcel.EmploymentEducation);
//(for application) alternative.AddUtilityTerm(251, destinationParcel.EmploymentFood);
//(for application) alternative.AddUtilityTerm(252, destinationParcel.EmploymentGovernment);
//(for application) alternative.AddUtilityTerm(253, destinationParcel.EmploymentOffice);
//(for application) alternative.AddUtilityTerm(254, destinationParcel.EmploymentRetail);
//(for application) alternative.AddUtilityTerm(255, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(256, destinationParcel.EmploymentMedical);
//(for application) alternative.AddUtilityTerm(257, destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
//(for application) alternative.AddUtilityTerm(258, destinationParcel.Households);
//(for application) alternative.AddUtilityTerm(259, destinationParcel.GetStudentsK12());
}
}
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
ChoiceProbabilityCalculator.Alternative alternative = sampleItem.Alternative;
//remove nesting for estimation of conditional MNL
if (!Global.Configuration.IsInEstimationMode || ESTIMATE_NESTED_MODEL)
{
alternative.AddNestedAlternative(_sampleSize + 2, 0, THETA_PARAMETER);
}
if (!alternative.Available)
{
return;
}
IParcelWrapper destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
// var destinationZoneTotals = ChoiceModelRunner.ZoneTotals[destinationParcel.ZoneId];
alternative.Choice = destinationParcel;
ChoiceProbabilityCalculator.Alternative nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeModel>().RunNested(_person, _person.Household.ResidenceParcel, destinationParcel, _destinationArrivalTime, _destinationDepartureTime, _person.Household.HouseholdTotals.DrivingAgeMembers);
double workTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
int votSegment = _person.Household.GetVotALSegment();
int taSegment = destinationParcel.TransitAccessSegment();
double aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][Global.Settings.CarOwnerships.OneOrMoreCarsPerAdult][votSegment][taSegment];
double distanceFromOrigin = _person.Household.ResidenceParcel.DistanceFromOrigin(destinationParcel, 1);
double distance1 = Math.Min(distanceFromOrigin, .35);
double distance2 = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
double distance3 = Math.Max(0, distanceFromOrigin - 1);
double distance20 = Math.Max(0, distanceFromOrigin - 2);
double distance30 = Math.Max(0, distanceFromOrigin - 3);
double distance40 = Math.Max(0, distanceFromOrigin - 4);
double distanceLog = Math.Log(1 + distanceFromOrigin);
double distanceFromSchool = _person.IsFullOrPartTimeWorker ? 0 : _person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
// parcel buffers
double educationBuffer = Math.Log(destinationParcel.EmploymentEducationBuffer2 + 1);
double governmentBuffer = Math.Log(destinationParcel.EmploymentGovernmentBuffer2 + 1);
double officeBuffer = Math.Log(destinationParcel.EmploymentOfficeBuffer2 + 1);
double serviceBuffer = Math.Log(destinationParcel.EmploymentServiceBuffer2 + 1);
double householdsBuffer = Math.Log(destinationParcel.HouseholdsBuffer2 + 1);
// var retailBuffer = Math.Log(destinationParcel.EmploymentRetailBuffer2 + 1);
double industrialAgricultureConstructionBuffer = Math.Log(destinationParcel.EmploymentIndustrialBuffer2 + destinationParcel.EmploymentAgricultureConstructionBuffer2 + 1);
double foodBuffer = Math.Log(destinationParcel.EmploymentFoodBuffer2 + 1);
double medicalBuffer = Math.Log(destinationParcel.EmploymentMedicalBuffer2 + 1);
double employmentTotalBuffer = Math.Log(destinationParcel.EmploymentTotalBuffer2 + 1);
double studentsUniversityBuffer = Math.Log(destinationParcel.StudentsUniversityBuffer2 + 1);
double studentsK12Buffer = Math.Log(destinationParcel.StudentsK8Buffer2 + destinationParcel.StudentsHighSchoolBuffer2 + 1);
// var mixedUse4Index = destinationParcel.MixedUse4Index2();
//size attributes (derived)
double employmentIndustrialAgricultureConstruction = destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction;
// parking attributes
double parcelParkingDensity = destinationParcel.ParcelParkingPerTotalEmployment();
// connectivity attributes
double c34Ratio = destinationParcel.C34RatioBuffer1();
alternative.AddUtilityTerm(1, sampleItem.AdjustmentFactor); //constrain coeff to 1
alternative.AddUtilityTerm(2, _person.IsFulltimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(3, _person.IsPartTimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(4, _person.IsNotFullOrPartTimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(5, distanceLog); // for distance calibration
alternative.AddUtilityTerm(6, _person.IsFulltimeWorker.ToFlag() * distance1);
alternative.AddUtilityTerm(7, _person.IsFulltimeWorker.ToFlag() * distance2);
alternative.AddUtilityTerm(8, _person.IsFulltimeWorker.ToFlag() * distance3);
alternative.AddUtilityTerm(9, _person.IsPartTimeWorker.ToFlag() * distanceLog);
alternative.AddUtilityTerm(10, _person.IsNotFullOrPartTimeWorker.ToFlag() * distanceLog);
alternative.AddUtilityTerm(11, _person.Household.Has0To15KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(12, _person.Household.Has50To75KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(13, _person.Household.Has75To100KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(14, _person.IsFemale.ToFlag() * distanceLog);
alternative.AddUtilityTerm(15, _person.IsStudentAge.ToFlag() * distanceFromSchool);
alternative.AddUtilityTerm(16, _person.IsFulltimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(17, _person.IsPartTimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(18, _person.IsNotFullOrPartTimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(19, parcelParkingDensity);
alternative.AddUtilityTerm(20, c34Ratio);
//extra additive distance terms for calibrating longer distances
alternative.AddUtilityTerm(46, distance20);
alternative.AddUtilityTerm(47, distance30);
alternative.AddUtilityTerm(48, distance40);
//Neighborhood
alternative.AddUtilityTerm(21, _person.Household.HasValidIncome.ToFlag() * serviceBuffer);
alternative.AddUtilityTerm(22, _person.Household.HasValidIncome.ToFlag() * educationBuffer);
alternative.AddUtilityTerm(23, _person.Household.HasValidIncome.ToFlag() * foodBuffer);
alternative.AddUtilityTerm(24, _person.Household.HasValidIncome.ToFlag() * governmentBuffer);
alternative.AddUtilityTerm(25, _person.Household.HasValidIncome.ToFlag() * officeBuffer);
alternative.AddUtilityTerm(26, _person.Household.HasValidIncome.ToFlag() * medicalBuffer);
alternative.AddUtilityTerm(27, _person.Household.HasValidIncome.ToFlag() * householdsBuffer);
alternative.AddUtilityTerm(28, _person.Household.HasValidIncome.ToFlag() * studentsUniversityBuffer);
alternative.AddUtilityTerm(29, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * studentsK12Buffer);
alternative.AddUtilityTerm(30, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * studentsUniversityBuffer);
alternative.AddUtilityTerm(31, _person.Household.HasValidIncome.ToFlag() * _person.IsPartTimeWorker.ToFlag() * industrialAgricultureConstructionBuffer);
alternative.AddUtilityTerm(32, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * foodBuffer);
alternative.AddUtilityTerm(33, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * medicalBuffer);
alternative.AddUtilityTerm(34, _person.IsFulltimeWorker.ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * employmentTotalBuffer);
alternative.AddUtilityTerm(35, _person.IsNotFullOrPartTimeWorker.ToFlag() * _person.Household.HasIncomeUnder50K.ToFlag() * governmentBuffer);
alternative.AddUtilityTerm(36, _person.IsNotFullOrPartTimeWorker.ToFlag() * _person.Household.HasIncomeUnder50K.ToFlag() * employmentTotalBuffer);
//Size
alternative.AddUtilityTerm(51, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(52, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(53, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(54, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(55, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(56, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(57, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(58, _person.Household.HasValidIncome.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(59, _person.Household.HasValidIncome.ToFlag() * destinationParcel.StudentsUniversity);
alternative.AddUtilityTerm(60, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(61, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(62, _person.Household.HasValidIncome.ToFlag() * _person.IsPartTimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(63, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(64, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(65, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(66, _person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(67, _person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(68, _person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(69, _person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(70, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(71, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(72, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(73, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(74, _person.IsFulltimeWorker.ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(75, _person.IsFulltimeWorker.ToFlag() * (!_person.Household.Has75KPlusIncome).ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(76, _person.IsPartTimeWorker.ToFlag() * (!_person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(77, (!_person.IsFulltimeWorker).ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(78, _person.IsNotFullOrPartTimeWorker.ToFlag() * (!_person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(79, _person.Household.HasMissingIncome.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(80, _person.Household.HasMissingIncome.ToFlag() * destinationParcel.StudentsUniversity);
//add any region-specific new terms in region-specific class, using coefficient numbers 91-100 or other unused variable #
_parentClass.RegionSpecificCustomizations(alternative, _person, destinationParcel);
// OD shadow pricing
if (Global.Configuration.ShouldUseODShadowPricing)
{
int res = _person.Household.ResidenceParcel.District;
int des = destinationParcel.District;
//var first = res <= des? res : des;
//var second = res <= des? des : res;
double shadowPriceConfigurationParameter = res == des ? Global.Configuration.WorkLocationOOShadowPriceCoefficient : Global.Configuration.WorkLocationODShadowPriceCoefficient;
int odShadowPriceF12Value = MAX_REGULAR_PARAMETER + Global.Configuration.NumberOfODShadowPricingDistricts * (res - 1) + des;
alternative.AddUtilityTerm(odShadowPriceF12Value, shadowPriceConfigurationParameter);
}
// set shadow price depending on persontype and add it to utility
// we are using the sampling adjustment factor assuming that it is 1
if (Global.Configuration.ShouldUseShadowPricing)
{
alternative.AddUtilityTerm(1, destinationParcel.ShadowPriceForEmployment);
}
}
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
var alternative = sampleItem.Alternative;
if (!alternative.Available)
{
return;
}
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
// var destinationZoneTotals = ChoiceModelRunner.ZoneTotals[destinationParcel.ZoneId];
alternative.Choice = destinationParcel;
var schoolTourLogsum = 0D;
var nestedAlternative = Global.ChoiceModelSession.Get <SchoolTourModeTimeModel>().RunNested(_person, _person.Household.ResidenceParcel, destinationParcel, _destinationArrivalTime, _destinationDepartureTime, _person.Household.HouseholdTotals.DrivingAgeMembers, 0.0);
schoolTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
var votSegment = _person.Household.GetVotALSegment();
var taSegment = destinationParcel.TransitAccessSegment();
var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][Global.Settings.CarOwnerships.OneOrMoreCarsPerAdult][votSegment][taSegment];
var distanceFromOrigin = _person.Household.ResidenceParcel.DistanceFromOrigin(destinationParcel, 1);
var distance1 = Math.Min(distanceFromOrigin, .1);
var distance2 = Math.Max(0, Math.Min(distanceFromOrigin - .1, .5 - .1));
var distance3 = Math.Max(0, distanceFromOrigin - .5);
var distanceLog = Math.Log(1 + distanceFromOrigin);
var distanceFromWork = _person.IsFullOrPartTimeWorker ? _person.UsualWorkParcel.DistanceFromWorkLog(destinationParcel, 1) : 0;
// var millionsSquareFeet = destinationZoneTotals.MillionsSquareFeet();
// zone densities
// var eduDensity = destinationZoneTotals.GetEmploymentEducationDensity(millionsSquareFeet);
// var govDensity = destinationZoneTotals.GetEmploymentGovernmentDensity(millionsSquareFeet);
// var offDensity = destinationZoneTotals.GetEmploymentOfficeDensity(millionsSquareFeet);
// var serDensity = destinationZoneTotals.GetEmploymentServiceDensity(millionsSquareFeet);
// var houDensity = destinationZoneTotals.GetHouseholdsDensity(millionsSquareFeet);
// parcel buffers
var educationBuffer1 = Math.Log(destinationParcel.EmploymentEducationBuffer1 + 1);
//var governmentBuffer1 = Math.Log(destinationParcel.EmploymentGovernmentBuffer1 + 1);
//var officeBuffer1 = Math.Log(destinationParcel.EmploymentOfficeBuffer1 + 1);
//var serviceBuffer1 = Math.Log(destinationParcel.EmploymentServiceBuffer1 + 1);
//var householdsBuffer1 = Math.Log(destinationParcel.HouseholdsBuffer1 + 1);
//var retailBuffer1 = Math.Log(destinationParcel.EmploymentRetailBuffer1 + 1);
//var industrialAgricultureConstructionBuffer1 = Math.Log(destinationParcel.EmploymentIndustrialBuffer1 + destinationParcel.EmploymentAgricultureConstructionBuffer1 + 1);
//var foodBuffer1 = Math.Log(destinationParcel.EmploymentFoodBuffer1 + 1);
//var medicalBuffer1 = Math.Log(destinationParcel.EmploymentMedicalBuffer1 + 1);
//var employmentTotalBuffer1 = Math.Log(destinationParcel.EmploymentTotalBuffer1 + 1);
var studentsUniversityBuffer1 = Math.Log(destinationParcel.StudentsUniversityBuffer1 + 1);
var studentsK8Buffer1 = Math.Log(destinationParcel.StudentsK8Buffer1 + 1);
var studentsHighSchoolBuffer1 = Math.Log(destinationParcel.StudentsHighSchoolBuffer1 + 1);
//var educationBuffer2 = Math.Log(destinationParcel.EmploymentEducationBuffer2 + 1);
//var governmentBuffer2 = Math.Log(destinationParcel.EmploymentGovernmentBuffer2 + 1);
//var officeBuffer2 = Math.Log(destinationParcel.EmploymentOfficeBuffer2 + 1);
//var serviceBuffer2 = Math.Log(destinationParcel.EmploymentServiceBuffer2 + 1);
var householdsBuffer2 = Math.Log(destinationParcel.HouseholdsBuffer2 + 1);
//var retailBuffer2 = Math.Log(destinationParcel.EmploymentRetailBuffer2 + 1);
//var industrialAgricultureConstructionBuffer2 = Math.Log(destinationParcel.EmploymentIndustrialBuffer2 + destinationParcel.EmploymentAgricultureConstructionBuffer2 + 1);
//var foodBuffer2 = Math.Log(destinationParcel.EmploymentFoodBuffer2 + 1);
//var medicalBuffer2 = Math.Log(destinationParcel.EmploymentMedicalBuffer2 + 1);
var employmentTotalBuffer2 = Math.Log(destinationParcel.EmploymentTotalBuffer2 + 1);
var studentsUniversityBuffer2 = Math.Log(destinationParcel.StudentsUniversityBuffer2 + 1);
//var studentsK8Buffer2 = Math.Log(destinationParcel.StudentsK8Buffer2 + 1);
//var studentsHighSchoolBuffer2 = Math.Log(destinationParcel.StudentsHighSchoolBuffer2 + 1);
// var educationBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentEducationBuffer1 - destinationParcel.EmploymentEducation) + 1);
// var governmentBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentGovernmentBuffer1 - destinationParcel.EmploymentGovernment) + 1);
// var officeBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentOfficeBuffer1 - destinationParcel.EmploymentOffice) + 1);
// var serviceBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentServiceBuffer1 - destinationParcel.EmploymentService) + 1);
// var householdsBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.HouseholdsBuffer1 - destinationParcel.Households) + 1);
// var retailBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentRetailBuffer1 - destinationParcel.EmploymentRetail) + 1);
// var industrialAgricultureConstructionBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentIndustrialBuffer1 + destinationParcel.EmploymentAgricultureConstructionBuffer1
// - destinationParcel.EmploymentIndustrial - destinationParcel.EmploymentAgricultureConstruction) + 1);
// var foodBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentFoodBuffer1 - destinationParcel.EmploymentFood) + 1);
// var medicalBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentMedicalBuffer1 - destinationParcel.EmploymentMedical) + 1);
// var employmentTotalBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentTotalBuffer1 - destinationParcel.EmploymentTotal) + 1);
// var studentsUniversityBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.StudentsUniversityBuffer1 - destinationParcel.StudentsUniversity) + 1);
// var studentsK8Buffer1 = Math.Log(Math.Max(0.0, destinationParcel.StudentsK8Buffer1 - destinationParcel.StudentsK8) + 1);
// var studentsHighSchoolBuffer1 = Math.Log(Math.Max(0.0, destinationParcel.StudentsHighSchoolBuffer1 - destinationParcel.StudentsHighSchool) + 1);
//
// var educationBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentEducationBuffer2 - destinationParcel.EmploymentEducation) + 1);
// var governmentBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentGovernmentBuffer2 - destinationParcel.EmploymentGovernment) + 1);
// var officeBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentOfficeBuffer2 - destinationParcel.EmploymentOffice) + 1);
// var serviceBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentServiceBuffer2 - destinationParcel.EmploymentService) + 1);
// var householdsBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.HouseholdsBuffer2 - destinationParcel.Households) + 1);
// var retailBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentRetailBuffer2 - destinationParcel.EmploymentRetail) + 1);
// var industrialAgricultureConstructionBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentIndustrialBuffer2 + destinationParcel.EmploymentAgricultureConstructionBuffer2
// - destinationParcel.EmploymentIndustrial - destinationParcel.EmploymentAgricultureConstruction) + 1);
// var foodBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentFoodBuffer2 - destinationParcel.EmploymentFood) + 1);
// var medicalBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentMedicalBuffer2 - destinationParcel.EmploymentMedical) + 1);
// var employmentTotalBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.EmploymentTotalBuffer2 - destinationParcel.EmploymentTotal) + 1);
// var studentsUniversityBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.StudentsUniversityBuffer2 - destinationParcel.StudentsUniversity) + 1);
// var studentsK8Buffer2 = Math.Log(Math.Max(0.0, destinationParcel.StudentsK8Buffer2 - destinationParcel.StudentsK8) + 1);
// var studentsHighSchoolBuffer2 = Math.Log(Math.Max(0.0, destinationParcel.StudentsHighSchoolBuffer2 - destinationParcel.StudentsHighSchool) + 1);
alternative.AddUtilityTerm(1, sampleItem.AdjustmentFactor);
alternative.AddUtilityTerm(2, _person.IsChildUnder5.ToFlag() * schoolTourLogsum);
alternative.AddUtilityTerm(3, _person.IsChildAge5Through15.ToFlag() * schoolTourLogsum);
alternative.AddUtilityTerm(4, _person.IsDrivingAgeStudent.ToFlag() * schoolTourLogsum);
alternative.AddUtilityTerm(5, _person.IsUniversityStudent.ToFlag() * schoolTourLogsum);
alternative.AddUtilityTerm(6, (!_person.IsStudentAge).ToFlag() * schoolTourLogsum);
alternative.AddUtilityTerm(7, _person.IsChildUnder5.ToFlag() * distance1);
alternative.AddUtilityTerm(8, _person.IsChildUnder5.ToFlag() * distance2);
alternative.AddUtilityTerm(9, _person.IsChildUnder5.ToFlag() * distance3);
alternative.AddUtilityTerm(10, _person.IsChildAge5Through15.ToFlag() * distance1);
alternative.AddUtilityTerm(11, _person.IsChildAge5Through15.ToFlag() * distance2);
alternative.AddUtilityTerm(12, _person.IsChildAge5Through15.ToFlag() * distance3);
alternative.AddUtilityTerm(13, _person.IsDrivingAgeStudent.ToFlag() * distanceLog);
alternative.AddUtilityTerm(14, _person.IsUniversityStudent.ToFlag() * distanceLog);
alternative.AddUtilityTerm(15, (!_person.IsStudentAge).ToFlag() * distanceLog);
alternative.AddUtilityTerm(16, (!_person.IsStudentAge).ToFlag() * distanceFromWork);
alternative.AddUtilityTerm(17, _person.IsChildUnder5.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(18, _person.IsChildAge5Through15.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(19, _person.IsDrivingAgeStudent.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(20, _person.IsUniversityStudent.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(21, (!_person.IsStudentAge).ToFlag() * aggregateLogsum);
//Neighborhood
alternative.AddUtilityTerm(30, _person.IsChildUnder5.ToFlag() * householdsBuffer2);
alternative.AddUtilityTerm(31, _person.IsChildUnder5.ToFlag() * studentsHighSchoolBuffer1);
alternative.AddUtilityTerm(32, _person.IsChildUnder5.ToFlag() * employmentTotalBuffer2);
alternative.AddUtilityTerm(33, _person.IsChildAge5Through15.ToFlag() * studentsK8Buffer1);
alternative.AddUtilityTerm(34, _person.IsDrivingAgeStudent.ToFlag() * studentsHighSchoolBuffer1);
alternative.AddUtilityTerm(35, _person.IsUniversityStudent.ToFlag() * educationBuffer1);
alternative.AddUtilityTerm(36, _person.IsAdult.ToFlag() * studentsUniversityBuffer1);
alternative.AddUtilityTerm(37, _person.IsAdult.ToFlag() * studentsUniversityBuffer2);
alternative.AddUtilityTerm(38, _person.IsAdult.ToFlag() * studentsK8Buffer1);
//Size
alternative.AddUtilityTerm(61, _person.IsChildUnder5.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(62, _person.IsChildUnder5.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(63, _person.IsChildUnder5.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(64, _person.IsChildUnder5.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(65, _person.IsChildUnder5.ToFlag() * 10.0 * destinationParcel.Households);
alternative.AddUtilityTerm(66, _person.IsChildUnder5.ToFlag() * destinationParcel.StudentsK8);
alternative.AddUtilityTerm(67, _person.IsChildAge5Through15.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(68, _person.IsChildAge5Through15.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(69, _person.IsChildAge5Through15.ToFlag() * destinationParcel.StudentsHighSchool);
alternative.AddUtilityTerm(70, _person.IsChildAge5Through15.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(71, _person.IsChildAge5Through15.ToFlag() * 10.0 * destinationParcel.Households);
alternative.AddUtilityTerm(72, _person.IsChildAge5Through15.ToFlag() * destinationParcel.StudentsK8);
alternative.AddUtilityTerm(73, _person.IsDrivingAgeStudent.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(74, _person.IsDrivingAgeStudent.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(75, _person.IsDrivingAgeStudent.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(76, _person.IsDrivingAgeStudent.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(77, _person.IsDrivingAgeStudent.ToFlag() * 10.0 * destinationParcel.Households);
alternative.AddUtilityTerm(78, _person.IsDrivingAgeStudent.ToFlag() * destinationParcel.StudentsHighSchool);
alternative.AddUtilityTerm(79, _person.IsAdult.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(80, _person.IsAdult.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(81, _person.IsAdult.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(82, _person.IsAdult.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(83, _person.IsAdult.ToFlag() * destinationParcel.StudentsUniversity);
alternative.AddUtilityTerm(84, _person.IsAdult.ToFlag() * destinationParcel.StudentsHighSchool);
// set shadow price depending on persontype and add it to utility
// we are using the sampling adjustment factor assuming that it is 1
alternative.AddUtilityTerm(1, _person.IsAdult ? destinationParcel.ShadowPriceForStudentsUniversity : destinationParcel.ShadowPriceForStudentsK12);
//remove nesting for estimation of conditional MNL
alternative.AddNestedAlternative(_sampleSize + 2, 0, 99);
}
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
var alternative = sampleItem.Alternative;
if (!alternative.Available)
{
return;
}
var household = _tour.Household;
var person = _tour.Person;
// var personDay = _tour.PersonDay;
var householdHasChildren = household.HasChildren;
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
var fastestTravelTime =
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, _tour.OriginParcel, destinationParcel).Variable +
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, destinationParcel, _tour.OriginParcel).Variable;
if (fastestTravelTime >= _maxAvailableMinutes)
{
alternative.Available = false;
return;
}
alternative.Choice = destinationParcel;
double tourLogsum;
if (_tour.IsHomeBasedTour)
{
if (_tour.DestinationPurpose == Global.Settings.Purposes.Escort)
{
var nestedAlternative = Global.ChoiceModelSession.Get <EscortTourModeModel>().RunNested(_tour, destinationParcel);
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
else
{
var nestedAlternative = Global.ChoiceModelSession.Get <OtherHomeBasedTourModeModel>().RunNested(_tour, destinationParcel);
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
}
else
{
var nestedAlternative = Global.ChoiceModelSession.Get <WorkBasedSubtourModeModel>().RunNested(_tour, destinationParcel);
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
//var purpose = _tour.TourPurposeSegment;
var carOwnership = person.GetCarOwnershipSegment();
//var votSegment = _tour.VotALSegment;
//var transitAccess = destinationParcel.TransitAccessSegment();
//var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][purpose][carOwnership][votSegment][transitAccess];
//var aggregateLogsumHomeBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][transitAccess];
//var aggregateLogsumWorkBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.Work_BASED][carOwnership][votSegment][transitAccess];
var distanceFromOrigin = _tour.OriginParcel.DistanceFromOrigin(destinationParcel, _tour.DestinationArrivalTime);
var distanceFromOrigin0 = Math.Max(0, Math.Min(distanceFromOrigin - .5, 1 - .5));
var distanceFromOrigin3 = Math.Max(0, distanceFromOrigin - 1);
var distanceFromOrigin4 = Math.Min(distanceFromOrigin, .10);
var distanceFromOrigin5 = Math.Max(0, Math.Min(distanceFromOrigin - .1, .5 - .1));
var distanceFromOrigin8 = Math.Max(0, Math.Min(distanceFromOrigin - .1, .35 - .1));
var distanceFromOrigin9 = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
var distanceFromOriginLog = Math.Log(1 + distanceFromOrigin);
var distanceFromWorkLog = person.UsualWorkParcel.DistanceFromWorkLog(destinationParcel, 1);
var distanceFromSchoolLog = person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
var timePressure = Math.Log(1 - fastestTravelTime / _maxAvailableMinutes);
// log transforms of buffers for Neighborhood effects
var empEduBuffer = Math.Log(destinationParcel.EmploymentEducationBuffer1 + 1.0);
var empFooBuffer = Math.Log(destinationParcel.EmploymentFoodBuffer1 + 1.0);
// var EMPGOV_B = Math.Log(destinationParcel.EmploymentGovernmentBuffer1 + 1.0);
var empOfcBuffer = Math.Log(destinationParcel.EmploymentOfficeBuffer1 + 1.0);
var empRetBuffer = Math.Log(destinationParcel.EmploymentRetailBuffer1 + 1.0);
var empSvcBuffer = Math.Log(destinationParcel.EmploymentServiceBuffer1 + 1.0);
var empMedBuffer = Math.Log(destinationParcel.EmploymentMedicalBuffer1 + 1.0);
// var EMPIND_B = Math.Log(destinationParcel.EmploymentIndustrialBuffer1 + destinationParcel.EmploymentAgricultureConstructionBuffer1 + 1.0);
var empTotBuffer = Math.Log(destinationParcel.EmploymentTotalBuffer1 + 1.0);
var housesBuffer = Math.Log(destinationParcel.HouseholdsBuffer1 + 1.0);
var studK12Buffer = Math.Log(destinationParcel.StudentsK8Buffer1 + destinationParcel.StudentsHighSchoolBuffer1 + 1.0);
var studUniBuffer = Math.Log(destinationParcel.StudentsUniversityBuffer1 + 1.0);
// var EMPHOU_B = Math.Log(destinationParcel.EmploymentTotalBuffer1 + destinationParcel.HouseholdsBuffer1 + 1.0);
// connectivity attributes
var c34Ratio = destinationParcel.C34RatioBuffer1();
var carCompetitionFlag = FlagUtility.GetCarCompetitionFlag(carOwnership); // exludes no cars
var noCarCompetitionFlag = FlagUtility.GetNoCarCompetitionFlag(carOwnership);
var noCarsFlag = FlagUtility.GetNoCarsFlag(carOwnership);
alternative.AddUtilityTerm(1, sampleItem.AdjustmentFactor);
alternative.AddUtilityTerm(2, (_tour.IsHomeBasedTour).ToFlag() * timePressure);
alternative.AddUtilityTerm(3, _secondaryFlag * _workOrSchoolPatternFlag * distanceFromOrigin0);
alternative.AddUtilityTerm(4, _secondaryFlag * _otherPatternFlag * distanceFromOrigin5);
alternative.AddUtilityTerm(5, _secondaryFlag * _otherPatternFlag * distanceFromOrigin0);
alternative.AddUtilityTerm(6, _secondaryFlag * _otherPatternFlag * distanceFromOrigin3);
alternative.AddUtilityTerm(7, (!_tour.IsHomeBasedTour).ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(8, household.Has0To15KIncome.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(9, household.HasMissingIncome.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(10, person.IsRetiredAdult.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(11, person.IsChildAge5Through15.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(12, person.IsChildUnder5.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(13, (_tour.IsHomeBasedTour).ToFlag() * distanceFromSchoolLog);
alternative.AddUtilityTerm(14, distanceFromWorkLog);
alternative.AddUtilityTerm(15, carCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
alternative.AddUtilityTerm(16, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixInParcel());
alternative.AddUtilityTerm(17, carCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
alternative.AddUtilityTerm(18, noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
alternative.AddUtilityTerm(19, noCarsFlag * c34Ratio);
if (_tour.DestinationPurpose == Global.Settings.Purposes.Escort)
{
alternative.AddUtilityTerm(20, tourLogsum);
alternative.AddUtilityTerm(21, distanceFromOrigin4);
alternative.AddUtilityTerm(22, distanceFromOrigin8);
alternative.AddUtilityTerm(23, distanceFromOrigin9);
// Neighborhood
alternative.AddUtilityTerm(24, householdHasChildren.ToFlag() * studK12Buffer);
alternative.AddUtilityTerm(25, empTotBuffer);
// Size terms
alternative.AddUtilityTerm(71, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(72, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(73, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(74, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(75, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(76, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(77, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(78, (!householdHasChildren).ToFlag() * destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(79, (!householdHasChildren).ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(80, householdHasChildren.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(81, householdHasChildren.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(82, householdHasChildren.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(83, householdHasChildren.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(84, householdHasChildren.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(85, householdHasChildren.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(86, householdHasChildren.ToFlag() * destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(87, householdHasChildren.ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(88, householdHasChildren.ToFlag() * destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.PersonalBusiness || _tour.DestinationPurpose == Global.Settings.Purposes.Medical)
{
alternative.AddUtilityTerm(26, tourLogsum);
alternative.AddUtilityTerm(27, distanceFromOrigin4);
alternative.AddUtilityTerm(28, distanceFromOrigin8);
alternative.AddUtilityTerm(29, distanceFromOrigin9);
alternative.AddUtilityTerm(30, distanceFromOrigin3);
// Neighborhood
alternative.AddUtilityTerm(31, empEduBuffer);
alternative.AddUtilityTerm(32, empSvcBuffer);
alternative.AddUtilityTerm(33, empMedBuffer);
alternative.AddUtilityTerm(34, housesBuffer); // also psrc
alternative.AddUtilityTerm(35, studUniBuffer);
// Size terms
alternative.AddUtilityTerm(89, destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(90, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(91, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(92, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(93, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(94, destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(95, destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(96, destinationParcel.Households);
alternative.AddUtilityTerm(97, destinationParcel.GetStudentsK12());
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Shopping)
{
alternative.AddUtilityTerm(36, tourLogsum);
alternative.AddUtilityTerm(37, distanceFromOrigin4);
alternative.AddUtilityTerm(38, distanceFromOrigin8);
alternative.AddUtilityTerm(39, distanceFromOrigin9);
alternative.AddUtilityTerm(40, distanceFromOrigin3);
// Neighborhood
alternative.AddUtilityTerm(41, empEduBuffer); // also psrc
alternative.AddUtilityTerm(42, empRetBuffer); // also psrc
// Size terms
alternative.AddUtilityTerm(98, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(99, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(100, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(101, destinationParcel.EmploymentService);
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Meal)
{
alternative.AddUtilityTerm(43, tourLogsum);
alternative.AddUtilityTerm(44, distanceFromOrigin4);
alternative.AddUtilityTerm(45, distanceFromOrigin8);
alternative.AddUtilityTerm(46, distanceFromOrigin9);
alternative.AddUtilityTerm(47, distanceFromOrigin3);
// Neighborhood
alternative.AddUtilityTerm(48, empFooBuffer); // psrc
// Size terms
alternative.AddUtilityTerm(102, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(103, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(104, destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(105, destinationParcel.Households);
}
else if (_tour.DestinationPurpose == Global.Settings.Purposes.Social || _tour.DestinationPurpose == Global.Settings.Purposes.Recreation)
{
alternative.AddUtilityTerm(49, tourLogsum);
alternative.AddUtilityTerm(50, distanceFromOrigin4);
alternative.AddUtilityTerm(51, distanceFromOrigin8);
alternative.AddUtilityTerm(52, distanceFromOrigin9);
alternative.AddUtilityTerm(53, distanceFromOrigin3);
// Neighborhood
alternative.AddUtilityTerm(54, empOfcBuffer); // also psrc
alternative.AddUtilityTerm(55, empSvcBuffer); // also psrc
alternative.AddUtilityTerm(56, housesBuffer); // also psrc
alternative.AddUtilityTerm(57, studUniBuffer); // psrc
// Size terms
alternative.AddUtilityTerm(106, destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(107, destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(108, destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(109, destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(110, destinationParcel.EmploymentService);
alternative.AddUtilityTerm(111, destinationParcel.Households);
alternative.AddUtilityTerm(112, destinationParcel.StudentsUniversity);
alternative.AddUtilityTerm(113, destinationParcel.GetStudentsK12());
}
//add any region-specific new terms in region-specific class, using coefficient numbers 114-120, or other unused variable #
_parentClass.RegionSpecificOtherTourDistrictCoefficients(alternative, _tour, destinationParcel);
// OD shadow pricing
if (Global.Configuration.ShouldUseODShadowPricing)
{
var ori = _tour.OriginParcel.District;
var des = destinationParcel.District;
//var first = res <= des? res : des;
//var second = res <= des? des : res;
var shadowPriceConfigurationParameter = ori == des ? Global.Configuration.OtherTourDestinationOOShadowPriceCoefficient : Global.Configuration.OtherTourDestinationODShadowPriceCoefficient;
var odShadowPriceF12Value = MAX_REGULAR_PARAMETER + Global.Configuration.NumberOfODShadowPricingDistricts * (ori - 1) + des;
alternative.AddUtilityTerm(odShadowPriceF12Value, shadowPriceConfigurationParameter);
}
}
Task IDataService.DeleteItem(ISampleItem item) => MethodRunner.RunAsync(() =>
{
_sampleProvider.DeleteItem(item.Id);
_items.Remove(item);
});
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
var alternative = sampleItem.Alternative;
var household = _tour.Household;
var person = _tour.Person;
var personDay = _tour.PersonDay;
// var householdHasChildren = household.HasChildren;
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
var excludedParcel = person.UsualWorkParcel == null || person.UsualWorkParcelId == household.ResidenceParcelId || _tour.DestinationPurpose != Global.Settings.Purposes.Work || _tour.GetTourCategory() == Global.Settings.TourCategories.WorkBased ? null : person.UsualWorkParcel;
var usualWorkParcel = (excludedParcel != null && excludedParcel.Id == destinationParcel.Id); // only 1 for oddball alternative on tours with oddball alternative
var usualWorkParcelFlag = usualWorkParcel.ToFlag();
// Comment out these nesting calls when estimating the conditional flat model
// model is NL with oddball alternative
if (usualWorkParcelFlag == 0)
{
// this alternative is in the non-oddball nest
alternative.AddNestedAlternative(_sampleSize + 2, 0, 60); // associates alternative with non-oddball nest
}
else
{
// this is the oddball alternative
alternative.AddNestedAlternative(_sampleSize + 3, 1, 60); // associates alternative with oddball nest
}
if (!alternative.Available)
{
return;
}
// use this block of code to eliminate the oddball alternative for estimation of the conditional model
//if (usualWorkParcelFlag == 1) {
// alternative.Available = false;
//
// return;
//}
var fastestTravelTime =
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, _tour.OriginParcel, destinationParcel).Variable +
ImpedanceRoster.GetValue("ivtime", Global.Settings.Modes.Hov3, Global.Settings.PathTypes.FullNetwork, Global.Settings.ValueOfTimes.DefaultVot, _fastestAvailableTimeOfDay, destinationParcel, _tour.OriginParcel).Variable;
if (fastestTravelTime >= _maxAvailableMinutes)
{
alternative.Available = false;
return;
}
alternative.Choice = destinationParcel;
double tourLogsum;
if (_tour.IsHomeBasedTour)
{
var nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeModel>().RunNested(_tour, destinationParcel, household.VehiclesAvailable, person.GetTransitFareDiscountFraction());
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
else
{
var nestedAlternative = Global.ChoiceModelSession.Get <WorkBasedSubtourModeModel>().RunNested(_tour, destinationParcel);
tourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
}
//var purpose = _tour.TourPurposeSegment;
//var nonWorkPurpose = _tour.DestinationPurpose != Global.Settings.Purposes.Work;
var carOwnership = person.GetCarOwnershipSegment();
//var votSegment = _tour.VotALSegment;
//var transitAccess = destinationParcel.TransitAccessSegment();
//var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][purpose][carOwnership][votSegment][transitAccess];
//var aggregateLogsumHomeBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][carOwnership][votSegment][transitAccess];
//var aggregateLogsumWorkBased = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.Work_BASED][carOwnership][votSegment][transitAccess];
var distanceFromOrigin = _tour.OriginParcel.DistanceFromOrigin(destinationParcel, _tour.DestinationArrivalTime);
// var distanceFromOrigin0 = Math.Max(0, Math.Min(distanceFromOrigin - .5, 1 - .5));
// var distanceFromOrigin3 = Math.Max(0, distanceFromOrigin - 1);
// var distanceFromOrigin4 = Math.Min(distanceFromOrigin, .10);
// var distanceFromOrigin5 = Math.Max(0, Math.Min(distanceFromOrigin - .1, .5 - .1));
// var distanceFromOrigin8 = Math.Max(0, Math.Min(distanceFromOrigin - .1, .35 - .1));
// var distanceFromOrigin9 = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
var distanceFromOriginLog = Math.Log(1 + distanceFromOrigin);
var distanceFromWorkLog = person.UsualWorkParcel.DistanceFromWorkLog(destinationParcel, 1);
var distanceFromSchoolLog = person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
// var millionsSquareFeet = destinationZoneTotals.MillionsSquareFeet();
var timePressure = Math.Log(1 - fastestTravelTime / _maxAvailableMinutes);
// log transforms of buffers for Neighborhood effects
var empEduBuffer = Math.Log(destinationParcel.EmploymentEducationBuffer1 + 1.0);
// var EMPFOO_B = Math.Log(destinationParcel.EmploymentFoodBuffer1 + 1.0);
// var EMPGOV_B = Math.Log(destinationParcel.EmploymentGovernmentBuffer1 + 1.0);
var empOfcBuffer = Math.Log(destinationParcel.EmploymentOfficeBuffer1 + 1.0);
// var EMPRET_B = Math.Log(destinationParcel.EmploymentRetailBuffer1 + 1.0);
// var EMPSVC_B = Math.Log(destinationParcel.EmploymentServiceBuffer1 + 1.0);
// var EMPMED_B = Math.Log(destinationParcel.EmploymentMedicalBuffer1 + 1.0);
var empIndBuffer = Math.Log(destinationParcel.EmploymentIndustrialBuffer1 + destinationParcel.EmploymentAgricultureConstructionBuffer1 + 1.0);
// var EMPTOT_B = Math.Log(destinationParcel.EmploymentTotalBuffer1 + 1.0);
var housesBuffer = Math.Log(destinationParcel.HouseholdsBuffer1 + 1.0);
// var STUDK12B = Math.Log(destinationParcel.StudentsK8Buffer1 + destinationParcel.StudentsHighSchoolBuffer1 + 1.0);
// var STUDUNIB = Math.Log(destinationParcel.StudentsUniversityBuffer1 + 1.0);
// var EMPHOU_B = Math.Log(destinationParcel.EmploymentTotalBuffer1 + destinationParcel.HouseholdsBuffer1 + 1.0);
// parking attributes
// var parcelParkingDensity = destinationParcel.ParcelParkingPerTotalEmployment();
// var zoneParkingDensity = destinationParcel.ZoneParkingPerTotalEmploymentAndK12UniversityStudents(destinationZoneTotals, millionsSquareFeet);
// var ParkingPaidDailyLogBuffer1 = Math.Log(1 + destinationParcel.ParkingOffStreetPaidDailySpacesBuffer1);
// connectivity attributes
// var c34Ratio = destinationParcel.C34RatioBuffer1();
// var carDeficitFlag = FlagUtility.GetCarDeficitFlag(carOwnership); // includes no cars
// var carCompetitionFlag = FlagUtility.GetCarCompetitionFlag(carOwnership); // exludes no cars
var noCarCompetitionFlag = FlagUtility.GetNoCarCompetitionFlag(carOwnership);
// var noCarsFlag = FlagUtility.GetNoCarsFlag(carOwnership);
// Usual location attributes
alternative.AddUtilityTerm(1, usualWorkParcelFlag);
alternative.AddUtilityTerm(2, person.IsPartTimeWorker.ToFlag() * usualWorkParcelFlag);
alternative.AddUtilityTerm(3, person.IsStudentAge.ToFlag() * usualWorkParcelFlag);
alternative.AddUtilityTerm(4, _primaryFlag * personDay.TwoOrMoreWorkToursExist().ToFlag() * usualWorkParcelFlag);
alternative.AddUtilityTerm(5, personDay.WorkStopsExist().ToFlag() * usualWorkParcelFlag);
alternative.AddUtilityTerm(6, _secondaryFlag * usualWorkParcelFlag);
// non-usual location attributes
alternative.AddUtilityTerm(11, (!usualWorkParcel).ToFlag() * sampleItem.AdjustmentFactor);
alternative.AddUtilityTerm(12, _tour.IsHomeBasedTour.ToFlag() * (!usualWorkParcel).ToFlag() * timePressure);
alternative.AddUtilityTerm(13, (!usualWorkParcel).ToFlag() * person.IsFulltimeWorker.ToFlag() * tourLogsum);
alternative.AddUtilityTerm(14, (!usualWorkParcel).ToFlag() * person.IsPartTimeWorker.ToFlag() * tourLogsum);
alternative.AddUtilityTerm(15, (!usualWorkParcel).ToFlag() * person.IsNotFullOrPartTimeWorker.ToFlag() * tourLogsum);
alternative.AddUtilityTerm(16, (!usualWorkParcel).ToFlag() * person.IsRetiredAdult.ToFlag() * distanceFromOriginLog);
alternative.AddUtilityTerm(17, (!usualWorkParcel).ToFlag() * distanceFromWorkLog);
alternative.AddUtilityTerm(18, (!usualWorkParcel).ToFlag() * person.IsStudentAge.ToFlag() * distanceFromSchoolLog);
alternative.AddUtilityTerm(19, (!usualWorkParcel).ToFlag() * noCarCompetitionFlag * destinationParcel.ParkingHourlyEmploymentCommercialMixBuffer1());
// non-usual location Neighborhood attributes
alternative.AddUtilityTerm(31, (!usualWorkParcel).ToFlag() * empEduBuffer);
alternative.AddUtilityTerm(32, (!usualWorkParcel).ToFlag() * empOfcBuffer);
alternative.AddUtilityTerm(33, (!usualWorkParcel).ToFlag() * housesBuffer);
alternative.AddUtilityTerm(34, (!usualWorkParcel).ToFlag() * empIndBuffer);
// non-usual location Size terms (consider conditioning these by fulltime, parttime, notFTPT, an income (see original sacog spec)
alternative.AddUtilityTerm(40, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(41, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(42, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(43, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(44, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(45, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(46, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(47, (!usualWorkParcel).ToFlag() * destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction);
alternative.AddUtilityTerm(48, (!usualWorkParcel).ToFlag() * destinationParcel.Households);
alternative.AddUtilityTerm(49, (!usualWorkParcel).ToFlag() * destinationParcel.StudentsUniversity);
// // usual location size term
// alternative.AddUtilityTerm(50, usualWorkParcelFlag * 1);
// // Comment out these nesting calls when estimating the conditional flat model
// // model is NL with oddball alternative
// if (usualWorkParcelFlag == 0) {
// // this alternative is in the non-oddball nest
// alternative.AddNestedAlternative(_sampleSize + 2, 0, 60); // associates alternative with non-oddball nest
// }
// else {
// // this is the oddball alternative
// alternative.AddNestedAlternative(_sampleSize + 3, 1, 60); // associates alternative with oddball nest
// }
}
public void SetUtilities(ISampleItem sampleItem, int sampleFrequency)
{
if (sampleItem == null)
{
throw new ArgumentNullException("sampleItem");
}
var alternative = sampleItem.Alternative;
if (!alternative.Available)
{
return;
}
var destinationParcel = ChoiceModelFactory.Parcels[sampleItem.ParcelId];
// var destinationZoneTotals = ChoiceModelRunner.ZoneTotals[destinationParcel.ZoneId];
alternative.Choice = destinationParcel;
var workTourLogsum = 0D;
var nestedAlternative = Global.ChoiceModelSession.Get <WorkTourModeTimeModel>().RunNested(_person, _person.Household.ResidenceParcel, destinationParcel, _destinationArrivalTime, _destinationDepartureTime, _person.Household.HouseholdTotals.DrivingAgeMembers, 0.0);
workTourLogsum = nestedAlternative == null ? 0 : nestedAlternative.ComputeLogsum();
var votSegment = _person.Household.GetVotALSegment();
var taSegment = destinationParcel.TransitAccessSegment();
var aggregateLogsum = Global.AggregateLogsums[destinationParcel.ZoneId][Global.Settings.Purposes.HomeBasedComposite][Global.Settings.CarOwnerships.OneOrMoreCarsPerAdult][votSegment][taSegment];
var distanceFromOrigin = _person.Household.ResidenceParcel.DistanceFromOrigin(destinationParcel, 1);
var distance1 = Math.Min(distanceFromOrigin, .35);
var distance2 = Math.Max(0, Math.Min(distanceFromOrigin - .35, 1 - .35));
var distance3 = Math.Max(0, distanceFromOrigin - 1);
var distanceLog = Math.Log(1 + distanceFromOrigin);
var distanceFromSchool = _person.IsFullOrPartTimeWorker ? 0 : _person.UsualSchoolParcel.DistanceFromSchoolLog(destinationParcel, 1);
// parcel buffers
var educationBuffer = Math.Log(destinationParcel.EmploymentEducationBuffer2 + 1);
var governmentBuffer = Math.Log(destinationParcel.EmploymentGovernmentBuffer2 + 1);
var officeBuffer = Math.Log(destinationParcel.EmploymentOfficeBuffer2 + 1);
var serviceBuffer = Math.Log(destinationParcel.EmploymentServiceBuffer2 + 1);
var householdsBuffer = Math.Log(destinationParcel.HouseholdsBuffer2 + 1);
// var retailBuffer = Math.Log(destinationParcel.EmploymentRetailBuffer2 + 1);
var industrialAgricultureConstructionBuffer = Math.Log(destinationParcel.EmploymentIndustrialBuffer2 + destinationParcel.EmploymentAgricultureConstructionBuffer2 + 1);
var foodBuffer = Math.Log(destinationParcel.EmploymentFoodBuffer2 + 1);
var medicalBuffer = Math.Log(destinationParcel.EmploymentMedicalBuffer2 + 1);
var employmentTotalBuffer = Math.Log(destinationParcel.EmploymentTotalBuffer2 + 1);
var studentsUniversityBuffer = Math.Log(destinationParcel.StudentsUniversityBuffer2 + 1);
var studentsK12Buffer = Math.Log(destinationParcel.StudentsK8Buffer2 + destinationParcel.StudentsHighSchoolBuffer2 + 1);
// var mixedUse4Index = destinationParcel.MixedUse4Index2();
//size attributes (derived)
var employmentIndustrialAgricultureConstruction = destinationParcel.EmploymentIndustrial + destinationParcel.EmploymentAgricultureConstruction;
// parking attributes
var parcelParkingDensity = destinationParcel.ParcelParkingPerTotalEmployment();
// connectivity attributes
var c34Ratio = destinationParcel.C34RatioBuffer1();
alternative.AddUtilityTerm(1, sampleItem.AdjustmentFactor);
alternative.AddUtilityTerm(2, _person.IsFulltimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(3, _person.IsPartTimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(4, _person.IsNotFullOrPartTimeWorker.ToFlag() * workTourLogsum);
alternative.AddUtilityTerm(5, distanceLog); // for distance calibration
alternative.AddUtilityTerm(6, _person.IsFulltimeWorker.ToFlag() * distance1);
alternative.AddUtilityTerm(7, _person.IsFulltimeWorker.ToFlag() * distance2);
alternative.AddUtilityTerm(8, _person.IsFulltimeWorker.ToFlag() * distance3);
alternative.AddUtilityTerm(9, _person.IsPartTimeWorker.ToFlag() * distanceLog);
alternative.AddUtilityTerm(10, _person.IsNotFullOrPartTimeWorker.ToFlag() * distanceLog);
alternative.AddUtilityTerm(11, _person.Household.Has0To15KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(12, _person.Household.Has50To75KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(13, _person.Household.Has75To100KIncome.ToFlag() * distanceLog);
alternative.AddUtilityTerm(14, _person.IsFemale.ToFlag() * distanceLog);
alternative.AddUtilityTerm(15, _person.IsStudentAge.ToFlag() * distanceFromSchool);
alternative.AddUtilityTerm(16, _person.IsFulltimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(17, _person.IsPartTimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(18, _person.IsNotFullOrPartTimeWorker.ToFlag() * aggregateLogsum);
alternative.AddUtilityTerm(19, parcelParkingDensity);
alternative.AddUtilityTerm(20, c34Ratio);
//Neighborhood
alternative.AddUtilityTerm(21, _person.Household.HasValidIncome.ToFlag() * serviceBuffer);
alternative.AddUtilityTerm(22, _person.Household.HasValidIncome.ToFlag() * educationBuffer);
alternative.AddUtilityTerm(23, _person.Household.HasValidIncome.ToFlag() * foodBuffer);
alternative.AddUtilityTerm(24, _person.Household.HasValidIncome.ToFlag() * governmentBuffer);
alternative.AddUtilityTerm(25, _person.Household.HasValidIncome.ToFlag() * officeBuffer);
alternative.AddUtilityTerm(26, _person.Household.HasValidIncome.ToFlag() * medicalBuffer);
alternative.AddUtilityTerm(27, _person.Household.HasValidIncome.ToFlag() * householdsBuffer);
alternative.AddUtilityTerm(28, _person.Household.HasValidIncome.ToFlag() * studentsUniversityBuffer);
alternative.AddUtilityTerm(29, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * studentsK12Buffer);
alternative.AddUtilityTerm(30, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * studentsUniversityBuffer);
alternative.AddUtilityTerm(31, _person.Household.HasValidIncome.ToFlag() * _person.IsPartTimeWorker.ToFlag() * industrialAgricultureConstructionBuffer);
alternative.AddUtilityTerm(32, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * foodBuffer);
alternative.AddUtilityTerm(33, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * medicalBuffer);
alternative.AddUtilityTerm(34, _person.IsFulltimeWorker.ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * employmentTotalBuffer);
alternative.AddUtilityTerm(35, _person.IsNotFullOrPartTimeWorker.ToFlag() * _person.Household.HasIncomeUnder50K.ToFlag() * governmentBuffer);
alternative.AddUtilityTerm(36, _person.IsNotFullOrPartTimeWorker.ToFlag() * _person.Household.HasIncomeUnder50K.ToFlag() * employmentTotalBuffer);
//PSRC specific district constants for model calibration
if (Global.Configuration.PSRC == true)
{
var homedist = _person.Household.ResidenceParcel.District;
var zonedist = destinationParcel.District;
var homeSKitWorkTRP = homedist == 11 && (zonedist == 8 || zonedist == 10 || zonedist == 7) ? 1 : 0;
var homeKitWorkTRP = homedist == 9 && (zonedist == 8 || zonedist == 10 || zonedist == 7) ? 1 : 0;
var homeEastWorkCBD = homedist == 6 && zonedist == 4 ? 1 : 0;
var homeKitWorkCBD = homedist == 9 && (zonedist == 4) ? 1 : 0;
var homeTacWorkKit = homedist == 8 && (zonedist == 9 || zonedist == 11) ? 1 : 0;
var homeEvWorkEv = homedist == 2 && zonedist == 2 ? 1 : 0;
var homeWSWorkEast = homedist == 5 && zonedist == 6 ? 1 : 0;
alternative.AddUtilityTerm(37, homeTacWorkKit);
alternative.AddUtilityTerm(38, homeEvWorkEv);
alternative.AddUtilityTerm(39, homeWSWorkEast);
alternative.AddUtilityTerm(40, homeSKitWorkTRP);
alternative.AddUtilityTerm(45, homeKitWorkTRP);
alternative.AddUtilityTerm(47, homeKitWorkTRP);
alternative.AddUtilityTerm(48, homeEastWorkCBD);
alternative.AddUtilityTerm(49, homeKitWorkCBD);
}
//Size
alternative.AddUtilityTerm(51, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(52, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(53, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(54, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(55, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(56, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(57, _person.Household.HasValidIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(58, _person.Household.HasValidIncome.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(59, _person.Household.HasValidIncome.ToFlag() * destinationParcel.StudentsUniversity);
alternative.AddUtilityTerm(60, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(61, _person.Household.HasValidIncome.ToFlag() * _person.IsFulltimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(62, _person.Household.HasValidIncome.ToFlag() * _person.IsPartTimeWorker.ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(63, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(64, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentFood);
alternative.AddUtilityTerm(65, _person.Household.HasValidIncome.ToFlag() * _person.IsNotFullOrPartTimeWorker.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(66, _person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(67, _person.Household.HasIncomeUnder50K.ToFlag() * destinationParcel.EmploymentService);
alternative.AddUtilityTerm(68, _person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(69, _person.Household.Has50To75KIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(70, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentEducation);
alternative.AddUtilityTerm(71, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(72, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(73, _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(74, _person.IsFulltimeWorker.ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentGovernment);
alternative.AddUtilityTerm(75, _person.IsFulltimeWorker.ToFlag() * (!_person.Household.Has75KPlusIncome).ToFlag() * employmentIndustrialAgricultureConstruction);
alternative.AddUtilityTerm(76, _person.IsPartTimeWorker.ToFlag() * (!_person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentMedical);
alternative.AddUtilityTerm(77, (!_person.IsFulltimeWorker).ToFlag() * _person.Household.Has75KPlusIncome.ToFlag() * destinationParcel.EmploymentOffice);
alternative.AddUtilityTerm(78, _person.IsNotFullOrPartTimeWorker.ToFlag() * (!_person.Household.HasIncomeUnder50K).ToFlag() * destinationParcel.EmploymentRetail);
alternative.AddUtilityTerm(79, _person.Household.HasMissingIncome.ToFlag() * destinationParcel.EmploymentTotal);
alternative.AddUtilityTerm(80, _person.Household.HasMissingIncome.ToFlag() * destinationParcel.StudentsUniversity);
// set shadow price depending on persontype and add it to utility
// we are using the sampling adjustment factor assuming that it is 1
alternative.AddUtilityTerm(1, destinationParcel.ShadowPriceForEmployment);
//remove nesting for estimation of conditional MNL
alternative.AddNestedAlternative(_sampleSize + 2, 0, 98);
}
