Changeset 3304

Timestamp:

05/05/08 14:18:27 (2 months ago)

Author:

bhylee

Message:

(bhylee) updated files

Location:

trunk/inprocess/bhylee

Files:

• hlcm_parcel_estimation.py (modified) (2 diffs)
• hlcm_parcel_specification.py (modified) (6 diffs)

trunk/inprocess/bhylee/hlcm_parcel_estimation.py

@@ -14 +14 @@
 
 from urbansim.configs.hlcm_estimation_config import HLCMEstimationConfig
+from psrc_parcel.configs.baseline_estimation import BaselineEstimation
 from opus_core.session_configuration import SessionConfiguration
 from opus_core.store.attribute_cache import AttributeCache
 from my_estimation_config import my_configuration
 
-class HlcmParcelEstimation(HLCMEstimationConfig):
+class HlcmParcelEstimation(BaselineEstimation): # comment out for urbansim.configs.hlcm_estimation_config
+#class HlcmParcelEstimation(HLCMEstimationConfig): # comment out for psrc_parcel.configs.baseline_estimation
     def update_config(self):
-        HLCMEstimationConfig.update_config(self)
+#        HLCMEstimationConfig.update_config(self) # comment out for psrc_parcel.configs.baseline_estimation
 #        
         self.replace(my_configuration)
         estimate_config = {}
-        estimate_config["export_estimation_data"]=False
-        estimate_config["estimation_data_file_name"]="/tmp/HLCM_building_estimate_data"
-        estimate_config["use_biogeme_data_format"]=True
-        estimate_config["weights_for_estimation_string"]=  "has_eg_1_units=building.residential_units>=1" #"psrc.parcel.residential_units_when_has_eg_1_surveyed_households_and_is_in_county_033"
+#        estimate_config["export_estimation_data"]=False
+#        estimate_config["estimation_data_file_name"]="/tmp/HLCM_building_estimate_data"
+#        estimate_config["use_biogeme_data_format"]=True
+#        estimate_config["weights_for_estimation_string"]=  "has_eg_1_units=building.residential_units>=1" #"psrc.parcel.residential_units_when_has_eg_1_surveyed_households_and_is_in_county_033"
         #"sampling_filter=(building.disaggregate(building_type.building_type_name)=='single_family_residential') + (building.disaggregate(building_type.building_type_name)=='multi_family_residential') + (building.disaggregate(building_type.building_type_name)=='condo_residential')"
         #"has_eg_1_units=urbansim.building.residential_units>=1"
             
-        estimate_config["stratum"] = "psrc.parcel.is_in_city_seattle" #"psrc.parcel.stratify_by_is_in_city_seattle_and_is_single_family_unit"
-        estimate_config["sample_size_from_each_stratum"] = 5
-        estimate_config["sample_size_from_chosen_stratum"] = 4
-        estimate_config["include_chosen_choice"] = True
+#        estimate_config["stratum"] = "psrc.parcel.is_in_city_seattle" #"psrc.parcel.stratify_by_is_in_city_seattle_and_is_single_family_unit"
+#        estimate_config["sample_size_from_each_stratum"] = 5
+#        estimate_config["sample_size_from_chosen_stratum"] = 4
+#        estimate_config["include_chosen_choice"] = True
+        estimate_config['wesml_sampling_correction_variable'] = 'psrc_parcel.building.wesml_sampling_correction_variable'
         #estimate_config['submodel_string'] = "None"
-        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['sample_size_locations'] = 30
+#        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['sample_size_locations'] = 30
         self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['sampler']="'opus_core.samplers.weighted_sampler'"#"'opus_core.samplers.stratified_sampler'"   #
-        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["estimate_config"]= \
-                            estimate_config
+        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["estimate_config"] = estimate_config
+        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["estimation_weight_string"] = "'has_eg_1_units=building.residential_units>=1'"
+#        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["estimation_weight_string"] = "'urbansim_parcel.building.vacant_residential_units'"
+#        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["estimation_weight_string"] = "'psrc_parcel.building.residential_units'"
                            #{"weights_for_estimation_string":"psrc.parcel.residential_units_when_has_eg_1_surveyed_households_and_is_in_county_033"}
         self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["location_set"] = "building"
         #self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["location_id_string"] = "'household.parcel_id'"
-        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['submodel_string'] = "'psrc.household.number_of_nonhome_based_workers'"
-        
+#        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['submodel_string'] = "'psrc.household.number_of_nonhome_based_workers'"
+
+        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]["variable_package"]="'urbansim_parcel'"
+
         self["models_configuration"]["household_location_choice_model"]["controller"]["prepare_for_estimate"]["arguments"]["join_datasets"] = 'True'
         self["models_configuration"]["household_location_choice_model"]["controller"]["prepare_for_estimate"]["arguments"]["index_to_unplace"] = 'None'
-        self["models_configuration"]["household_location_choice_model"]["controller"]["prepare_for_estimate"]["arguments"]["filter"] = None #"'psrc.household.customized_filter'"
+        self["models_configuration"]["household_location_choice_model"]["controller"]["prepare_for_estimate"]["arguments"]["filter"] = "'household.move == 1'"#None #"'psrc.household.customized_filter'"
+        self["models_configuration"]["household_location_choice_model"]["controller"]["init"]["arguments"]['filter'] = None
         
         self["models"] = [
-        #    {"household_relocation_model": ["run"]},
+#            {"household_relocation_model": ["run"]},
             {"household_location_choice_model": ["estimate"]}
         ]
@@ -66 +74 @@
                             run_configuration, "household_location_choice_model",
                             "inprocess.bhylee.hlcm_parcel_specification")
-    run_configuration.replace(my_configuration)
     er = Estimator(run_configuration, save_estimation_results=False)
     
     er.estimate()
+    er.create_prediction_success_table()
+    er.create_prediction_success_table(choice_geography_id="area_type_id=building.disaggregate(zone.area_type_id, intermediates=[parcel])" )
+    er.create_prediction_success_table(choice_geography_id="building_type_id=building.building_type_id" )
+#    er.create_prediction_success_table(choice_geography_id="large_area_id=building.disaggregate(faz.large_area_id, intermediates=[zone, parcel])" )
 #    er.reestimate("hlcm_parcel_specification")

trunk/inprocess/bhylee/hlcm_parcel_specification.py

@@ -17 +17 @@
             [
 
+#             ("urbansim_parcel.household_x_building.ln_sampling_probability_for_bias_correction_mnl_residential_units", "bias", 1),
+#             ("urbansim_parcel.household_x_building.ln_sampling_probability_for_bias_correction_mnl_vacant_residential_units", "bias", 1),
+
+
 #        ALTERNATIVE
+             "A_is_in_same_area_type = (household.previous_area_type_id == building.disaggregate(zone.area_type_id))",
+#             "A_is_in_same_district = (household.previous_district_id == building.disaggregate(zone.district_id))",
+#             "A_is_in_same_faz = (household.previous_faz_id == building.disaggregate(zone.faz_id))",
+             "A_is_in_same_large_area = (household.previous_large_area_id == building.disaggregate(faz.large_area_id))",
+#             "A_is_in_same_zone = (household.previous_zone_id == building.disaggregate(zone.zone_id))",
 #             "A_is_new_construction = (urbansim_parcel.building.age_masked < 2)",
 #             "A_is_pre1945 = (urbansim_parcel.building.age_masked > 60)",
@@ -26 +35 @@
 #             "A_ln_price_per_sqft = ln(urbansim_parcel.building.unit_price)",
 #             "A_ln_price_per_unit = ln((urbansim_parcel.building.unit_price) * urbansim_parcel.building.building_sqft_per_unit)",
-#             "A_ln_residential_units = ln(psrc_parcel.building.residential_units)",
+             "A_ln_residential_units = ln(psrc_parcel.building.residential_units)",
+#             "A_ln_vacant_residential_units = ln(urbansim_parcel.building.vacant_residential_units)",
 #             "A_unit_price = urbansim_parcel.building.unit_price * urbansim_parcel.building.building_sqft_per_unit",
+
                           
 #        NEIGHBOURHOOD
@@ -38 +49 @@
 #             "B_number_of_high_income_households_within_walking_distance = building.disaggregate(psrc.parcel.number_of_high_income_households_within_walking_distance)",
 #             "B_percent_high_income_households_within_walking_distance = building.disaggregate(psrc.parcel.percent_high_income_households_within_walking_distance)",
+#             "B_percent_home_owner_households_within_walking_distance = building.disaggregate(psrc.parcel.percent_home_owners_within_walking_distance)",
 #             "B_percent_low_income_households_within_walking_distance = building.disaggregate(psrc.parcel.percent_low_income_households_within_walking_distance)",
 #             "B_percent_med_income_households_within_walking_distance = building.disaggregate(psrc.parcel.percent_mid_income_households_within_walking_distance)",
 #             "B_percent_young_households_within_walking_distance = building.disaggregate(psrc.parcel.percent_young_households_within_walking_distance)",
-#             "urbansim_parcel.household_x_building.ln_sampling_probability_for_bias_correction_mnl_vacant_residential_units",
 
 
@@ -58 +69 @@
 #             "C_is_low_income_x_ln_zonal_pop_den = urbansim.household.is_low_income * (ln(building.disaggregate(urbansim_parcel.zone.population_per_acre))).astype(float32)",
 #             "C_is_mid_income_x_ln_average_income_within_walking_distance = urbansim.household.is_mid_income * ln(building.disaggregate(psrc.parcel.average_income_within_walking_distance))",
+             "C_is_renter_x_is_multi_family_residential = (household.tenure == 3) * urbansim.building.is_multi_family_residential",
              "C_is_single_person_x_is_not_single_family_residential = (household.persons < 2) * numpy.logical_not(urbansim.building.is_single_family_residential)",
 #             "C_is_young_x_percent_young_households_within_walking_distance = urbansim.household.is_young * building.disaggregate(psrc.parcel.percent_young_households_within_walking_distance)",
@@ -68 +80 @@
 #             "C_ln_income_x_ln_average_zonal_income = ln(household.income) * ln(building.disaggregate(urbansim_parcel.zone.average_income))",
              "C_ln_income_x_ln_price_per_sqft = ln(household.income) * ln(urbansim_parcel.building.unit_price)",
-#             "C_ln_income_x_ln_zonal_pop_den = ln(household.income) * (ln(building.disaggregate(urbansim_parcel.zone.population_per_acre))).astype(float32)",
+             "C_ln_income_x_ln_zonal_pop_den = ln(household.income) * (ln(building.disaggregate(urbansim_parcel.zone.population_per_acre))).astype(float32)",
 
 #             "C_hh_size_x_ln_zonal_pop_den = household.persons * (ln(building.disaggregate(urbansim_parcel.zone.population_per_acre))).astype(float32)",
 #             "C_has_children_x_ln_zonal_pop_den = (household.children > 0) * (ln(building.disaggregate(urbansim_parcel.zone.population_per_acre))).astype(float32)",
-             
+
                           
 #        ACCESSIBILITY
@@ -84 +96 @@
 #             "D_logsum_income_break3_x_trip_weighted_average_logsum_hbw_am_income_3 = (psrc.household.logsum_income_break==3) * building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_3)",
 #             "D_logsum_income_break4_x_trip_weighted_average_logsum_hbw_am_income_4 = (psrc.household.logsum_income_break==4) * building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_4)",
-#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone = ln(urbansim.building.ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone)",
-#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_transit_walk = ln(urbansim.building.ln_employment_within_30_minutes_travel_time_hbw_am_transit_walk)",
+#             "D_logsum_income_break1_x_exp_trip_weighted_average_logsum_hbw_am_income_1 = (psrc.household.logsum_income_break==1) * exp(building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_1))",
+#             "D_logsum_income_break2_x_exp_trip_weighted_average_logsum_hbw_am_income_2 = (psrc.household.logsum_income_break==2) * exp(building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_2))",
+#             "D_logsum_income_break3_x_exp_trip_weighted_average_logsum_hbw_am_income_3 = (psrc.household.logsum_income_break==3) * exp(building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_3))",
+#             "D_logsum_income_break4_x_exp_trip_weighted_average_logsum_hbw_am_income_4 = (psrc.household.logsum_income_break==4) * exp(building.disaggregate(psrc.zone.trip_weighted_average_logsum_hbw_am_income_4))",
+#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone = building.disaggregate(psrc.zone.ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone)",
+#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone2 = building.disaggregate(ln_bounded(urbansim_parcel.zone.employment_within_30_minutes_travel_time_hbw_am_drive_alone))",
+#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone3 = urbansim.building.ln_employment_within_30_minutes_travel_time_hbw_am_drive_alone",
+#             "D_ln_employment_within_30_minutes_travel_time_hbw_am_transit_walk = urbansim.building.ln_employment_within_30_minutes_travel_time_hbw_am_transit_walk",
 #             "D_ln_generalized_cost_weighted_access_to_employment_hbw_am_drive_alone = ln(building.disaggregate(psrc.zone.generalized_cost_weighted_access_to_employment_hbw_am_drive_alone))",
 #             "D_ln_generalized_cost_weighted_access_to_employment_hbw_am_transit_walk = ln(building.disaggregate(psrc.zone.generalized_cost_weighted_access_to_employment_hbw_am_transit_walk))",