************************************************************** **********************CLEANING AND STATS********************** /*variables and descriptive statistics*/ import excel "2024-12-23_fleet_pilot.xlsx", sheet("2024-12-23_fleet_pilot") firstrow clear save "fleet_pilot.dta", replace use fleet_pilot.dta sum ourID *drop time rename QextID_time QextID_me rename Qend_time Qend_me drop *time drop *tim drop *ti rename QextID_me QextID_time rename Qend_me Qend_time *drop known TESTs drop if ourID==255 | ourID==256 | ourID==179 | ourID<=169 | ourID==191 *br QextID *drop nedokončené dotazníky drop if Qend==. gen zamer = (QIVQ11==1 | QIVQ12==1 | QIVQ13==1) tab zamer tab zamer QIVQ9, missing *zjistit dokončené experimenty foreach i in 1 2 3 4 5 6 7 8 { generate počet`i'=QDCERES`i'cvp+QDCERES`i'cvd+QDCERES`i'phev+QDCERES`i'bev+QDCERES`i'optout } *drop nedokončené experimenty foreach i in 1 2 3 4 5 6 7 8 { drop if počet`i' ==. } * mediany, speedery gen day= substr(QextID_time,9,2) destring day, replace gen started_hour= substr(QextID_time,12,2) gen started_minute= substr(QextID_time,15,2) gen started_second= substr(QextID_time,18,2) gen finished_hour= substr(Qend_time,12,2) gen finished_minute= substr(Qend_time,15,2) gen finished_second= substr(Qend_time,18,2) destring started_hour, replace destring started_minute, replace destring started_second, replace destring finished_hour, replace destring finished_minute, replace destring finished_second, replace gen time=((finished_hour*60*60+finished_minute*60+finished_second)-(started_hour*60*60+started_minute*60+started_second))/60 replace time=0 if time<=0 sum time, d gen speeder_cz=1 if time<0.48*18.25 replace speeder_cz=0 if speeder_cz==. br speeder_cz *time length of questionnaire // Step 1: Convert string timestamps to Stata datetime format gen double datetime1 = clock(Qend_time, "YMDhms") gen double datetime2 = clock(QextID_time, "YMDhms") // Step 2: Calculate the time difference in seconds gen double time_diff_seconds = abs(datetime1 - datetime2) / 1000 // Step 3: Extract hours, minutes, and seconds gen int hours = floor(time_diff_seconds / 3600) gen int minutes = floor(mod(time_diff_seconds, 3600) / 60) gen int seconds = mod(time_diff_seconds, 60) // Step 4: Combine hours, minutes, and seconds into the desired format gen time_diff_display = string(hours) + ":" + /// string(minutes, "%02.0f") + ":" + /// string(seconds, "%02.0f") // Step 5: Verify the results list Qend_time QextID_time time_diff_display if _n <= 10 br time_diff_display *TADY SE POČÍTÁ KOLIK ZAMÝŠLÍ KOUPIT VOZIDEL* * Replace missing values in QIVQ2used_j_7 and QIVQ3used_j_7 with 0 gen QIVQ2used_j_7_tmp = cond(missing(QIVQ2used_j_7), 0.0, QIVQ2used_j_7) gen QIVQ3used_j_7_tmp = cond(missing(QIVQ3used_j_7), 0.0, QIVQ3used_j_7) * Replace missing values in QIVQ2_j_7 and QIVQ3_j_7 with 0 (if needed) gen QIVQ2_j_7_tmp = cond(missing(QIVQ2_j_7), 0.0, QIVQ2_j_7) gen QIVQ3_j_7_tmp = cond(missing(QIVQ3_j_7), 0.0, QIVQ3_j_7) * Compute differences using temporary variables gen maxIVQ2_diff = QIVQ2_j_7_tmp - QIVQ2used_j_7_tmp gen maxIVQ3_diff = QIVQ3_j_7_tmp - QIVQ3used_j_7_tmp * Initialize result gen zamysleny = . * Calculate result based on conditions replace zamysleny = maxIVQ2_diff if QIVQ2_j_7_tmp > 0 & maxIVQ2_diff > QIVQ2used_j_7_tmp replace zamysleny = QIVQ2used_j_7_tmp if QIVQ2_j_7_tmp > 0 & maxIVQ2_diff <= QIVQ2used_j_7_tmp replace zamysleny = maxIVQ3_diff if QIVQ2_j_7_tmp <= 0 & maxIVQ3_diff > QIVQ3used_j_7_tmp replace zamysleny = QIVQ3used_j_7_tmp if QIVQ2_j_7_tmp <= 0 & maxIVQ3_diff <= QIVQ3used_j_7_tmp * Clean up temporary variables drop QIVQ2used_j_7_tmp QIVQ3used_j_7_tmp QIVQ2_j_7_tmp QIVQ3_j_7_tmp * Check the results list QIVQ2_j_7 QIVQ2used_j_7 maxIVQ2_diff QIVQ3_j_7 QIVQ3used_j_7 maxIVQ3_diff zamysleny list počet1 zamysleny generate zk5=1 if počet1!=zamysleny | počet2!=zamysleny | počet3!=zamysleny | počet4!=zamysleny | počet5!=zamysleny | počet6!=zamysleny | počet7!=zamysleny | počet8!=zamysleny list zk5 QVQ3 tab QVQ3 /* 1 Nabízené ceny byly jiné, než jaké jsem měl/a na mysli, takže jsem upravil/a počet 2 Uvažuji o jiném typu (technologii) vozidel 3 Změnil/a jsem názor, reálně jich potřebuji méně než jsem uvedl/a 4 Nejsem schopný/á se rozhodnout jakou technologii mají mít zbývající auta 5 Už jsem byl/a unaven/a 6 Nevšiml/a jsem si, že součet nesedí a chci počet upravit */ generate zk4=QVQ2 if QDCERES1optout==1 | QDCERES2optout==1 | QDCERES3optout==1 | QDCERES4optout==1 | QDCERES5optout==1 | QDCERES6optout==1 | QDCERES7optout==1 | QDCERES8optout==1 list zk4 QVQ2 tab QVQ2 /*1 Žádná z nabízených alternativ neodpovídala mým preferencím. 2 Cena nabízených alternativ byla příliš vysoká. 3 Mám jiné preference ohledně typu vozidla (např. preferuji jiný typ pohonu). 4 Nemám momentálně potřebu ani zájem pořídit nové vozidlo. 5 Chyběly mi další důležité informace k rozhodnutí (např. detaily o výkonu, spotřebě). */ sum time if zamysleny==1, d sum time if zamysleny>1, d ******COMP****** ******PRICE***** substring *CVp price *QdceCompCvdPrice1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvpPrice_cl_num`i' = real(subinstr(subinstr(QdceCompCvpPrice`i', " Kč", "", .), " ", "", .)) } *CVd price *QdceCompCvdPrice1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvdPrice_cl_num`i' = real(subinstr(subinstr(QdceCompCvdPrice`i', " Kč", "", .), " ", "", .)) } *phev price *QdceCompPhevPrice1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevPrice_cl_num`i' = real(subinstr(subinstr(QdceCompPhevPrice`i', " Kč", "", .), " ", "", .)) } *bev price *QdceCompBevPrice1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevPrice_cl_num`i' = real(subinstr(subinstr(QdceCompBevPrice`i', " Kč", "", .), " ", "", .)) } br QdceCompCvdPrice_cl_num1 ******SUBSIDY***** substring *phev subsidy *QdceCompPhevSubsidy1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevSubsidy_cl_num`i' = real(subinstr(subinstr(QdceCompPhevSubsidy`i', " Kč", "", .), " ", "", .)) } *bev subsidy *QdceCompBevSubsidy1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevSubsidy_cl_num`i' = real(subinstr(subinstr(QdceCompBevSubsidy`i', " Kč", "", .), " ", "", .)) } br QdceCompBevSubsidy_cl_num1 ******PRICE WITH SUBSIDY***** substring *phev price w subsidy *QdceCompPhevPriceWithSubsidy1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevPriceSub_cl_num`i' = real(subinstr(subinstr(QdceCompPhevPriceWithSubsidy`i', " Kč", "", .), " ", "", .)) } *bev price w subsidy *QdceCompBevPriceWithSubsidy1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevPriceSub_cl_num`i' = real(subinstr(subinstr(QdceCompBevPriceWithSubsidy`i', " Kč", "", .), " ", "", .)) } *br QdceCompBevPriceSub_cl_num1 ******COST***** substring *CVp cost *QdceCompCvpCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvpCost_cl_num`i' = real(subinstr(subinstr(QdceCompCvpCost`i', " Kč/km", "", .), ",", ".", .)) } *CVd cost *QdceCompCvdCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvdCost_cl_num`i' = real(subinstr(subinstr(QdceCompCvdCost`i', " Kč/km", "", .), ",", ".", .)) } *phev cost *QdceCompPhevCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevCost_cl_num`i' = real(subinstr(subinstr(QdceCompPhevCost`i', " Kč/km", "", .), ",", ".", .)) } *bev cost *QdceCompBevCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevCost_cl_num`i' = real(subinstr(subinstr(QdceCompBevCost`i', " Kč/km", "", .), ",", ".", .)) } ******RANGE***** substring *CVp range *QdceCompCvpRange1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvpRange_cl_num`i' = real(subinstr(subinstr(QdceCompCvpRange`i', " km", "", .), " ", "", .)) } *CVd range *QdceCompCvdRange1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvdRange_cl_num`i' = real(subinstr(subinstr(QdceCompCvdRange`i', " km", "", .), " ", "", .)) } *phev range *QdceCompPhevRange1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevElRange_cl_num`i' = real(regexs(1)) if regexm(QdceCompPhevRange`i', "^([0-9]+) km") } foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevCvRange_cl_num`i' = real(regexs(1)) if regexm(QdceCompPhevRange`i', "([0-9]+) km$") } *bev range *QdceCompBevRange1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevRange_cl_num`i' = real(subinstr(subinstr(QdceCompBevRange`i', " km", "", .), " ", "", .)) } ******YEARLY_COST***** substring *CVp price *QdceCompCvdPriceYearlyCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvpPYearCost_cl_num`i' = real(subinstr(subinstr(QdceCompCvpPriceYearlyCost`i', " Kč/rok", "", .), " ", "", .)) } *CVd price *QdceCompCvdPriceYearlyCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompCvdPYearCost_cl_num`i' = real(subinstr(subinstr(QdceCompCvdPriceYearlyCost`i', " Kč/rok", "", .), " ", "", .)) } *phev price *QdceCompPhevPriceYearlyCost1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompPhevPYearCost_cl_num`i' = real(subinstr(subinstr(QdceCompPhevPriceYearlyCost`i', " Kč/rok", "", .), " ", "", .)) } *bev price *QdceCompBevPrice1 foreach i in 1 2 3 4 5 6 7 8 { generate QdceCompBevPYearCost_cl_num`i' = real(subinstr(subinstr(QdceCompBevPriceYearlyCost`i', " Kč/rok", "", .), " ", "", .)) } save "fleet_pilot.dta", replace use "fleet_pilot.dta", clear export excel "fleet_pilot.xlsx", firstrow(variables) replace ******************************************************************* *************************DATA PREPARATION************************** *Preparing the data /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear ** first experiment, 8 choice sets in total, QDCE_RUN==8 ** * DR denotes total driving range in KM, while DRe keeps range driven from a baterry. It means we need to generate DRoil=DR-DRe *****first experiment, 4 alternatives + opt-out***** **** first choice set ***** ******** CVp ********** generate pair=1 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num1 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num1 gen narok = 0 gen cost = QdceCompCvpCost_cl_num1 gen DR = QdceCompCvpRange_cl_num1 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES1max = max(QDCERES1bev, QDCERES1cvd, QDCERES1cvp, QDCERES1optout, QDCERES1phev) *br QDCERES1max generate choice = 0 replace choice = 1 if QDCERES1max==QDCERES1cvp *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES1cvp *br QDCERES1cvp * Count how many alternatives have the maximum utility generate num_max1 = (QDCERES1bev == QDCERES1max) + /// (QDCERES1cvd == QDCERES1max) + /// (QDCERES1cvp == QDCERES1max) + /// (QDCERES1optout == QDCERES1max) + /// (QDCERES1phev == QDCERES1max) // Drop observations where num_max is greater than 1 for the current task drop if num_max1 > 1 save pair11.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=1 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num1 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num1 gen narok = 0 gen cost = QdceCompCvdCost_cl_num1 gen DR = QdceCompCvdRange_cl_num1 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES1max = max(QDCERES1bev, QDCERES1cvd, QDCERES1cvp, QDCERES1optout, QDCERES1phev) *br QDCERES1max generate choice = 0 replace choice = 1 if QDCERES1max==QDCERES1cvd *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES1cvd *br QDCERES1cvp * Count how many alternatives have the maximum utility generate num_max1 = (QDCERES1bev == QDCERES1max) + /// (QDCERES1cvd == QDCERES1max) + /// (QDCERES1cvp == QDCERES1max) + /// (QDCERES1optout == QDCERES1max) + /// (QDCERES1phev == QDCERES1max) // Drop observations where num_max is greater than 1 for the current task drop if num_max1 > 1 save pair12.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=1 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num1 gen price_subs = QdceCompBevSubsidy_cl_num1 gen price = QdceCompBevPriceSub_cl_num1 gen narok = QdceBevNarok1 gen cost = QdceCompBevCost_cl_num1 gen DR = 0 gen DRe = QdceCompBevRange_cl_num1 gen pubcharge = QdceBevCharge1 gen polv2g = QdceBevV2g1 gen polwallbox = QdceBevSubskw1 generate QDCERES1max = max(QDCERES1bev, QDCERES1cvd, QDCERES1cvp, QDCERES1optout, QDCERES1phev) *br QDCERES1max generate choice = 0 replace choice = 1 if QDCERES1max==QDCERES1bev *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES1bev *br QDCERES1cvp * Count how many alternatives have the maximum utility generate num_max1 = (QDCERES1bev == QDCERES1max) + /// (QDCERES1cvd == QDCERES1max) + /// (QDCERES1cvp == QDCERES1max) + /// (QDCERES1optout == QDCERES1max) + /// (QDCERES1phev == QDCERES1max) // Drop observations where num_max is greater than 1 for the current task drop if num_max1 > 1 save pair13.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=1 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num1 gen price_subs = QdceCompPhevSubsidy_cl_num1 gen price = QdceCompPhevPriceSub_cl_num1 gen narok = QdcePhevNarok1 gen cost = QdceCompPhevCost_cl_num1 gen DR = QdceCompPhevCvRange_cl_num1 gen DRe = QdceCompPhevElRange_cl_num1 gen pubcharge = QdcePhevCharge1 gen polv2g = 0 gen polwallbox = 0 generate QDCERES1max = max(QDCERES1bev, QDCERES1cvd, QDCERES1cvp, QDCERES1optout, QDCERES1phev) *br QDCERES1max generate choice = 0 replace choice = 1 if QDCERES1max==QDCERES1phev *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES1phev *br QDCERES1cvp * Count how many alternatives have the maximum utility generate num_max1 = (QDCERES1bev == QDCERES1max) + /// (QDCERES1cvd == QDCERES1max) + /// (QDCERES1cvp == QDCERES1max) + /// (QDCERES1optout == QDCERES1max) + /// (QDCERES1phev == QDCERES1max) // Drop observations where num_max is greater than 1 for the current task drop if num_max1 > 1 save pair14.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=1 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES1max = max(QDCERES1bev, QDCERES1cvd, QDCERES1cvp, QDCERES1optout, QDCERES1phev) *br QDCERES1max generate choice = 0 replace choice = 1 if QDCERES1max==QDCERES1optout *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES1optout *br QDCERES1cvp * Count how many alternatives have the maximum utility generate num_max1 = (QDCERES1bev == QDCERES1max) + /// (QDCERES1cvd == QDCERES1max) + /// (QDCERES1cvp == QDCERES1max) + /// (QDCERES1optout == QDCERES1max) + /// (QDCERES1phev == QDCERES1max) // Drop observations where num_max is greater than 1 for the current task drop if num_max1 > 1 save pair15.dta br num_max1 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie1 = (num_max1 > 1) br has_tie1 * Optional: List the choice sets with ties for verification list QDCERES1bev QDCERES1cvd QDCERES1cvp QDCERES1optout QDCERES1phev has_tie1 if has_tie1==1 *list QDCERES1bev QDCERES1cvd QDCERES1cvp QDCERES1optout QDCERES1phev has_tie1 if has_tie1==0 **** second choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=2 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num2 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num2 gen narok = 0 gen cost = QdceCompCvpCost_cl_num2 gen DR = QdceCompCvpRange_cl_num2 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES2max = max(QDCERES2bev, QDCERES2cvd, QDCERES2cvp, QDCERES2optout, QDCERES2phev) *br QDCERES2max generate choice = 0 replace choice = 1 if QDCERES2max==QDCERES2cvp *br choice generate choice_bhat = 0 replace choice_bhat = QDCERES2cvp * Count how many alternatives have the maximum utility generate num_max2 = (QDCERES2bev == QDCERES2max) + /// (QDCERES2cvd == QDCERES2max) + /// (QDCERES2cvp == QDCERES2max) + /// (QDCERES2optout == QDCERES2max) + /// (QDCERES2phev == QDCERES2max) // Drop observations where num_max is greater than 1 for the current task drop if num_max2 > 1 save pair21.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=2 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num2 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num2 gen narok = 0 gen cost = QdceCompCvdCost_cl_num2 gen DR = QdceCompCvdRange_cl_num2 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES2max = max(QDCERES2bev, QDCERES2cvd, QDCERES2cvp, QDCERES2optout, QDCERES2phev) *br QDCERES2max generate choice = 0 replace choice = 1 if QDCERES2max==QDCERES2cvd br choice generate choice_bhat = 0 replace choice_bhat = QDCERES2cvd * Count how many alternatives have the maximum utility generate num_max2 = (QDCERES2bev == QDCERES2max) + /// (QDCERES2cvd == QDCERES2max) + /// (QDCERES2cvp == QDCERES2max) + /// (QDCERES2optout == QDCERES2max) + /// (QDCERES2phev == QDCERES2max) // Drop observations where num_max is greater than 1 for the current task drop if num_max2 > 1 save pair22.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=2 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num2 gen price_subs = QdceCompBevSubsidy_cl_num2 gen price = QdceCompBevPriceSub_cl_num2 gen narok = QdceBevNarok2 gen cost = QdceCompBevCost_cl_num2 gen DR = 0 gen DRe = QdceCompBevRange_cl_num2 gen pubcharge = QdceBevCharge2 gen polv2g = QdceBevV2g2 gen polwallbox = QdceBevSubskw2 generate QDCERES2max = max(QDCERES2bev, QDCERES2cvd, QDCERES2cvp, QDCERES2optout, QDCERES2phev) *br QDCERES2max generate choice = 0 replace choice = 1 if QDCERES2max==QDCERES2bev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES2bev * Count how many alternatives have the maximum utility generate num_max2 = (QDCERES2bev == QDCERES2max) + /// (QDCERES2cvd == QDCERES2max) + /// (QDCERES2cvp == QDCERES2max) + /// (QDCERES2optout == QDCERES2max) + /// (QDCERES2phev == QDCERES2max) // Drop observations where num_max is greater than 1 for the current task drop if num_max2 > 1 save pair23.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=2 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num2 gen price_subs = QdceCompPhevSubsidy_cl_num2 gen price = QdceCompPhevPriceSub_cl_num2 gen narok = QdcePhevNarok2 gen cost = QdceCompPhevCost_cl_num2 gen DR = QdceCompPhevCvRange_cl_num2 gen DRe = QdceCompPhevElRange_cl_num2 gen pubcharge = QdcePhevCharge2 gen polv2g = 0 gen polwallbox = 0 generate QDCERES2max = max(QDCERES2bev, QDCERES2cvd, QDCERES2cvp, QDCERES2optout, QDCERES2phev) *br QDCERES2max generate choice = 0 replace choice = 1 if QDCERES2max==QDCERES2phev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES2phev * Count how many alternatives have the maximum utility generate num_max2 = (QDCERES2bev == QDCERES2max) + /// (QDCERES2cvd == QDCERES2max) + /// (QDCERES2cvp == QDCERES2max) + /// (QDCERES2optout == QDCERES2max) + /// (QDCERES2phev == QDCERES2max) // Drop observations where num_max is greater than 1 for the current task drop if num_max2 > 1 save pair24.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=2 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES2max = max(QDCERES2bev, QDCERES2cvd, QDCERES2cvp, QDCERES2optout, QDCERES2phev) *br QDCERES2max generate choice = 0 replace choice = 1 if QDCERES2max==QDCERES2optout *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES2optout * Count how many alternatives have the maximum utility generate num_max2 = (QDCERES2bev == QDCERES2max) + /// (QDCERES2cvd == QDCERES2max) + /// (QDCERES2cvp == QDCERES2max) + /// (QDCERES2optout == QDCERES2max) + /// (QDCERES2phev == QDCERES2max) // Drop observations where num_max is greater than 1 for the current task drop if num_max2 > 1 save pair25.dta br num_max2 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie2 = (num_max2 > 1) br has_tie2 * Optional: List the choice sets with ties for verification list ourID QDCERES2bev QDCERES2cvd QDCERES2cvp QDCERES2optout QDCERES2phev has_tie2 if has_tie2==1 *list QDCERES2bev QDCERES2cvd QDCERES2cvp QDCERES2optout QDCERES2phev has_tie2 if has_tie2==0 **** third choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=3 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num3 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num3 gen narok = 0 gen cost = QdceCompCvpCost_cl_num3 gen DR = QdceCompCvpRange_cl_num3 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES3max = max(QDCERES3bev, QDCERES3cvd, QDCERES3cvp, QDCERES3optout, QDCERES3phev) *br QDCERES3max generate choice = 0 replace choice = 1 if QDCERES3max==QDCERES3cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES3cvp * Count how many alternatives have the maximum utility generate num_max3 = (QDCERES3bev == QDCERES3max) + /// (QDCERES3cvd == QDCERES3max) + /// (QDCERES3cvp == QDCERES3max) + /// (QDCERES3optout == QDCERES3max) + /// (QDCERES3phev == QDCERES3max) // Drop observations where num_max is greater than 1 for the current task drop if num_max3 > 1 save pair31.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=3 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num3 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num3 gen narok = 0 gen cost = QdceCompCvdCost_cl_num3 gen DR = QdceCompCvdRange_cl_num3 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES3max = max(QDCERES3bev, QDCERES3cvd, QDCERES3cvp, QDCERES3optout, QDCERES3phev) *br QDCERES3max generate choice = 0 replace choice = 1 if QDCERES3max==QDCERES3cvd br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES3cvd * Count how many alternatives have the maximum utility generate num_max3 = (QDCERES3bev == QDCERES3max) + /// (QDCERES3cvd == QDCERES3max) + /// (QDCERES3cvp == QDCERES3max) + /// (QDCERES3optout == QDCERES3max) + /// (QDCERES3phev == QDCERES3max) // Drop observations where num_max is greater than 1 for the current task drop if num_max3 > 1 save pair32.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=3 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num3 gen price_subs = QdceCompBevSubsidy_cl_num3 gen price = QdceCompBevPriceSub_cl_num3 gen narok = QdceBevNarok3 gen cost = QdceCompBevCost_cl_num3 gen DR = 0 gen DRe = QdceCompBevRange_cl_num3 gen pubcharge = QdceBevCharge3 gen polv2g = QdceBevV2g3 gen polwallbox = QdceBevSubskw3 generate QDCERES3max = max(QDCERES3bev, QDCERES3cvd, QDCERES3cvp, QDCERES3optout, QDCERES3phev) *br QDCERES3max generate choice = 0 replace choice = 1 if QDCERES3max==QDCERES3bev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES3bev * Count how many alternatives have the maximum utility generate num_max3 = (QDCERES3bev == QDCERES3max) + /// (QDCERES3cvd == QDCERES3max) + /// (QDCERES3cvp == QDCERES3max) + /// (QDCERES3optout == QDCERES3max) + /// (QDCERES3phev == QDCERES3max) // Drop observations where num_max is greater than 1 for the current task drop if num_max3 > 1 save pair33.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=3 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num3 gen price_subs = QdceCompPhevSubsidy_cl_num3 gen price = QdceCompPhevPriceSub_cl_num3 gen narok = QdcePhevNarok3 gen cost = QdceCompPhevCost_cl_num3 gen DR = QdceCompPhevCvRange_cl_num3 gen DRe = QdceCompPhevElRange_cl_num3 gen pubcharge = QdcePhevCharge3 gen polv2g = 0 gen polwallbox = 0 generate QDCERES3max = max(QDCERES3bev, QDCERES3cvd, QDCERES3cvp, QDCERES3optout, QDCERES3phev) *br QDCERES3max generate choice = 0 replace choice = 1 if QDCERES3max==QDCERES3phev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES3phev * Count how many alternatives have the maximum utility generate num_max3 = (QDCERES3bev == QDCERES3max) + /// (QDCERES3cvd == QDCERES3max) + /// (QDCERES3cvp == QDCERES3max) + /// (QDCERES3optout == QDCERES3max) + /// (QDCERES3phev == QDCERES3max) // Drop observations where num_max is greater than 1 for the current task drop if num_max3 > 1 save pair34.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=3 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES3max = max(QDCERES3bev, QDCERES3cvd, QDCERES3cvp, QDCERES3optout, QDCERES3phev) *br QDCERES3max generate choice = 0 replace choice = 1 if QDCERES3max==QDCERES3optout br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES3optout * Count how many alternatives have the maximum utility generate num_max3 = (QDCERES3bev == QDCERES3max) + /// (QDCERES3cvd == QDCERES3max) + /// (QDCERES3cvp == QDCERES3max) + /// (QDCERES3optout == QDCERES3max) + /// (QDCERES3phev == QDCERES3max) // Drop observations where num_max is greater than 1 for the current task drop if num_max3 > 1 save pair35.dta br num_max3 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie3 = (num_max3 > 1) br has_tie3 * Optional: List the choice sets with ties for verification list ourID QDCERES3bev QDCERES3cvd QDCERES3cvp QDCERES3optout QDCERES3phev has_tie3 if has_tie3==1 *list QDCERES3bev QDCERES3cvd QDCERES3cvp QDCERES3optout QDCERES3phev has_tie3 if has_tie3==0 **** fourth choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=4 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num4 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num4 gen narok = 0 gen cost = QdceCompCvpCost_cl_num4 gen DR = QdceCompCvpRange_cl_num4 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES4max = max(QDCERES4bev, QDCERES4cvd, QDCERES4cvp, QDCERES4optout, QDCERES4phev) *br QDCERES4max generate choice = 0 replace choice = 1 if QDCERES4max==QDCERES4cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES4cvp * Count how many alternatives have the maximum utility generate num_max4 = (QDCERES4bev == QDCERES4max) + /// (QDCERES4cvd == QDCERES4max) + /// (QDCERES4cvp == QDCERES4max) + /// (QDCERES4optout == QDCERES4max) + /// (QDCERES4phev == QDCERES4max) // Drop observations where num_max is greater than 1 for the current task drop if num_max4 > 1 save pair41.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=4 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num4 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num4 gen narok = 0 gen cost = QdceCompCvdCost_cl_num4 gen DR = QdceCompCvdRange_cl_num4 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES4max = max(QDCERES4bev, QDCERES4cvd, QDCERES4cvp, QDCERES4optout, QDCERES4phev) *br QDCERES4max generate choice = 0 replace choice = 1 if QDCERES4max==QDCERES4cvd br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES4cvd * Count how many alternatives have the maximum utility generate num_max4 = (QDCERES4bev == QDCERES4max) + /// (QDCERES4cvd == QDCERES4max) + /// (QDCERES4cvp == QDCERES4max) + /// (QDCERES4optout == QDCERES4max) + /// (QDCERES4phev == QDCERES4max) // Drop observations where num_max is greater than 1 for the current task drop if num_max4 > 1 save pair42.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=4 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num4 gen price_subs = QdceCompBevSubsidy_cl_num4 gen price = QdceCompBevPriceSub_cl_num4 gen narok = QdceBevNarok4 gen cost = QdceCompBevCost_cl_num4 gen DR = 0 gen DRe = QdceCompBevRange_cl_num4 gen pubcharge = QdceBevCharge4 gen polv2g = QdceBevV2g4 gen polwallbox = QdceBevSubskw4 generate QDCERES4max = max(QDCERES4bev, QDCERES4cvd, QDCERES4cvp, QDCERES4optout, QDCERES4phev) *br QDCERES4max generate choice = 0 replace choice = 1 if QDCERES4max==QDCERES4bev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES4bev * Count how many alternatives have the maximum utility generate num_max4 = (QDCERES4bev == QDCERES4max) + /// (QDCERES4cvd == QDCERES4max) + /// (QDCERES4cvp == QDCERES4max) + /// (QDCERES4optout == QDCERES4max) + /// (QDCERES4phev == QDCERES4max) // Drop observations where num_max is greater than 1 for the current task drop if num_max4 > 1 save pair43.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=4 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num4 gen price_subs = QdceCompPhevSubsidy_cl_num4 gen price = QdceCompPhevPriceSub_cl_num4 gen narok = QdcePhevNarok4 gen cost = QdceCompPhevCost_cl_num4 gen DR = QdceCompPhevCvRange_cl_num4 gen DRe = QdceCompPhevElRange_cl_num4 gen pubcharge = QdcePhevCharge4 gen polv2g = 0 gen polwallbox = 0 generate QDCERES4max = max(QDCERES4bev, QDCERES4cvd, QDCERES4cvp, QDCERES4optout, QDCERES4phev) *br QDCERES4max generate choice = 0 replace choice = 1 if QDCERES4max==QDCERES4phev br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES4phev * Count how many alternatives have the maximum utility generate num_max4 = (QDCERES4bev == QDCERES4max) + /// (QDCERES4cvd == QDCERES4max) + /// (QDCERES4cvp == QDCERES4max) + /// (QDCERES4optout == QDCERES4max) + /// (QDCERES4phev == QDCERES4max) // Drop observations where num_max is greater than 1 for the current task drop if num_max4 > 1 save pair44.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=4 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES4max = max(QDCERES4bev, QDCERES4cvd, QDCERES4cvp, QDCERES4optout, QDCERES4phev) *br QDCERES4max generate choice = 0 replace choice = 1 if QDCERES4max==QDCERES4optout br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES4optout * Count how many alternatives have the maximum utility generate num_max4 = (QDCERES4bev == QDCERES4max) + /// (QDCERES4cvd == QDCERES4max) + /// (QDCERES4cvp == QDCERES4max) + /// (QDCERES4optout == QDCERES4max) + /// (QDCERES4phev == QDCERES4max) // Drop observations where num_max is greater than 1 for the current task drop if num_max4 > 1 save pair45.dta br num_max4 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie4 = (num_max4 > 1) br has_tie4 * Optional: List the choice sets with ties for verification list QDCERES4bev QDCERES4cvd QDCERES4cvp QDCERES4optout QDCERES4phev has_tie4 if has_tie4==1 *list QDCERES4bev QDCERES4cvd QDCERES4cvp QDCERES4optout QDCERES4phev has_tie4 if has_tie4==0 **** fifth choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=5 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num5 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num5 gen narok = 0 gen cost = QdceCompCvpCost_cl_num5 gen DR = QdceCompCvpRange_cl_num5 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES5max = max(QDCERES5bev, QDCERES5cvd, QDCERES5cvp, QDCERES5optout, QDCERES5phev) *br QDCERES5max generate choice = 0 replace choice = 1 if QDCERES5max==QDCERES5cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES5cvp * Count how many alternatives have the maximum utility generate num_max5 = (QDCERES5bev == QDCERES5max) + /// (QDCERES5cvd == QDCERES5max) + /// (QDCERES5cvp == QDCERES5max) + /// (QDCERES5optout == QDCERES5max) + /// (QDCERES5phev == QDCERES5max) // Drop observations where num_max is greater than 1 for the current task drop if num_max5 > 1 save pair51.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=5 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num5 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num5 gen narok = 0 gen cost = QdceCompCvdCost_cl_num5 gen DR = QdceCompCvdRange_cl_num5 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES5max = max(QDCERES5bev, QDCERES5cvd, QDCERES5cvp, QDCERES5optout, QDCERES5phev) *br QDCERES5max generate choice = 0 replace choice = 1 if QDCERES5max==QDCERES5cvd *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES5cvd * Count how many alternatives have the maximum utility generate num_max5 = (QDCERES5bev == QDCERES5max) + /// (QDCERES5cvd == QDCERES5max) + /// (QDCERES5cvp == QDCERES5max) + /// (QDCERES5optout == QDCERES5max) + /// (QDCERES5phev == QDCERES5max) // Drop observations where num_max is greater than 1 for the current task drop if num_max5 > 1 save pair52.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=5 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num5 gen price_subs = QdceCompBevSubsidy_cl_num5 gen price = QdceCompBevPriceSub_cl_num5 gen narok = QdceBevNarok5 gen cost = QdceCompBevCost_cl_num5 gen DR = 0 gen DRe = QdceCompBevRange_cl_num5 gen pubcharge = QdceBevCharge5 gen polv2g = QdceBevV2g5 gen polwallbox = QdceBevSubskw5 generate QDCERES5max = max(QDCERES5bev, QDCERES5cvd, QDCERES5cvp, QDCERES5optout, QDCERES5phev) *br QDCERES5max generate choice = 0 replace choice = 1 if QDCERES5max==QDCERES5bev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES5bev * Count how many alternatives have the maximum utility generate num_max5 = (QDCERES5bev == QDCERES5max) + /// (QDCERES5cvd == QDCERES5max) + /// (QDCERES5cvp == QDCERES5max) + /// (QDCERES5optout == QDCERES5max) + /// (QDCERES5phev == QDCERES5max) // Drop observations where num_max is greater than 1 for the current task drop if num_max5 > 1 save pair53.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=5 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num5 gen price_subs = QdceCompPhevSubsidy_cl_num5 gen price = QdceCompPhevPriceSub_cl_num5 gen narok = QdcePhevNarok5 gen cost = QdceCompPhevCost_cl_num5 gen DR = QdceCompPhevCvRange_cl_num5 gen DRe = QdceCompPhevElRange_cl_num5 gen pubcharge = QdcePhevCharge5 gen polv2g = 0 gen polwallbox = 0 generate QDCERES5max = max(QDCERES5bev, QDCERES5cvd, QDCERES5cvp, QDCERES5optout, QDCERES5phev) *br QDCERES5max generate choice = 0 replace choice = 1 if QDCERES5max==QDCERES5phev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES5phev * Count how many alternatives have the maximum utility generate num_max5 = (QDCERES5bev == QDCERES5max) + /// (QDCERES5cvd == QDCERES5max) + /// (QDCERES5cvp == QDCERES5max) + /// (QDCERES5optout == QDCERES5max) + /// (QDCERES5phev == QDCERES5max) // Drop observations where num_max is greater than 1 for the current task drop if num_max5 > 1 save pair54.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=5 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES5max = max(QDCERES5bev, QDCERES5cvd, QDCERES5cvp, QDCERES5optout, QDCERES5phev) *br QDCERES5max generate choice = 0 replace choice = 1 if QDCERES5max==QDCERES5optout *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES5optout * Count how many alternatives have the maximum utility generate num_max5 = (QDCERES5bev == QDCERES5max) + /// (QDCERES5cvd == QDCERES5max) + /// (QDCERES5cvp == QDCERES5max) + /// (QDCERES5optout == QDCERES5max) + /// (QDCERES5phev == QDCERES5max) // Drop observations where num_max is greater than 1 for the current task drop if num_max5 > 1 save pair55.dta br num_max5 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie5 = (num_max5 > 1) br has_tie5 * Optional: List the choice sets with ties for verification list QDCERES5bev QDCERES5cvd QDCERES5cvp QDCERES5optout QDCERES5phev has_tie5 if has_tie5==1 *list QDCERES5bev QDCERES5cvd QDCERES5cvp QDCERES5optout QDCERES5phev has_tie5 if has_tie5==0 **** sixth choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=6 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num6 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num6 gen narok = 0 gen cost = QdceCompCvpCost_cl_num6 gen DR = QdceCompCvpRange_cl_num6 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES6max = max(QDCERES6bev, QDCERES6cvd, QDCERES6cvp, QDCERES6optout, QDCERES6phev) *br QDCERES6max generate choice = 0 replace choice = 1 if QDCERES6max==QDCERES6cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES6cvp * Count how many alternatives have the maximum utility generate num_max6 = (QDCERES6bev == QDCERES6max) + /// (QDCERES6cvd == QDCERES6max) + /// (QDCERES6cvp == QDCERES6max) + /// (QDCERES6optout == QDCERES6max) + /// (QDCERES6phev == QDCERES6max) // Drop observations where num_max is greater than 1 for the current task drop if num_max6 > 1 save pair61.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=6 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num6 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num6 gen narok = 0 gen cost = QdceCompCvdCost_cl_num6 gen DR = QdceCompCvdRange_cl_num6 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES6max = max(QDCERES6bev, QDCERES6cvd, QDCERES6cvp, QDCERES6optout, QDCERES6phev) *br QDCERES6max generate choice = 0 replace choice = 1 if QDCERES6max==QDCERES6cvd *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES6cvd * Count how many alternatives have the maximum utility generate num_max6 = (QDCERES6bev == QDCERES6max) + /// (QDCERES6cvd == QDCERES6max) + /// (QDCERES6cvp == QDCERES6max) + /// (QDCERES6optout == QDCERES6max) + /// (QDCERES6phev == QDCERES6max) // Drop observations where num_max is greater than 1 for the current task drop if num_max6 > 1 save pair62.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=6 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num6 gen price_subs = QdceCompBevSubsidy_cl_num6 gen price = QdceCompBevPriceSub_cl_num6 gen narok = QdceBevNarok6 gen cost = QdceCompBevCost_cl_num6 gen DR = 0 gen DRe = QdceCompBevRange_cl_num6 gen pubcharge = QdceBevCharge6 gen polv2g = QdceBevV2g6 gen polwallbox = QdceBevSubskw6 generate QDCERES6max = max(QDCERES6bev, QDCERES6cvd, QDCERES6cvp, QDCERES6optout, QDCERES6phev) *br QDCERES6max generate choice = 0 replace choice = 1 if QDCERES6max==QDCERES6bev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES6bev * Count how many alternatives have the maximum utility generate num_max6 = (QDCERES6bev == QDCERES6max) + /// (QDCERES6cvd == QDCERES6max) + /// (QDCERES6cvp == QDCERES6max) + /// (QDCERES6optout == QDCERES6max) + /// (QDCERES6phev == QDCERES6max) // Drop observations where num_max is greater than 1 for the current task drop if num_max6 > 1 save pair63.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=6 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num6 gen price_subs = QdceCompPhevSubsidy_cl_num6 gen price = QdceCompPhevPriceSub_cl_num6 gen narok = QdcePhevNarok6 gen cost = QdceCompPhevCost_cl_num6 gen DR = QdceCompPhevCvRange_cl_num6 gen DRe = QdceCompPhevElRange_cl_num6 gen pubcharge = QdcePhevCharge6 gen polv2g = 0 gen polwallbox = 0 generate QDCERES6max = max(QDCERES6bev, QDCERES6cvd, QDCERES6cvp, QDCERES6optout, QDCERES6phev) *br QDCERES6max generate choice = 0 replace choice = 1 if QDCERES6max==QDCERES6phev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES6phev * Count how many alternatives have the maximum utility generate num_max6 = (QDCERES6bev == QDCERES6max) + /// (QDCERES6cvd == QDCERES6max) + /// (QDCERES6cvp == QDCERES6max) + /// (QDCERES6optout == QDCERES6max) + /// (QDCERES6phev == QDCERES6max) // Drop observations where num_max is greater than 1 for the current task drop if num_max6 > 1 save pair64.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=6 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES6max = max(QDCERES6bev, QDCERES6cvd, QDCERES6cvp, QDCERES6optout, QDCERES6phev) *br QDCERES6max generate choice = 0 replace choice = 1 if QDCERES6max==QDCERES6optout *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES6optout * Count how many alternatives have the maximum utility generate num_max6 = (QDCERES6bev == QDCERES6max) + /// (QDCERES6cvd == QDCERES6max) + /// (QDCERES6cvp == QDCERES6max) + /// (QDCERES6optout == QDCERES6max) + /// (QDCERES6phev == QDCERES6max) // Drop observations where num_max is greater than 1 for the current task drop if num_max6 > 1 save pair65.dta br num_max6 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie6 = (num_max6 > 1) br has_tie6 * Optional: List the choice sets with ties for verification list QDCERES6bev QDCERES6cvd QDCERES6cvp QDCERES6optout QDCERES6phev has_tie6 if has_tie6==1 *list QDCERES6bev QDCERES6cvd QDCERES6cvp QDCERES6optout QDCERES6phev has_tie6 if has_tie6==0 **** seventh choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=7 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num7 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num7 gen narok = 0 gen cost = QdceCompCvpCost_cl_num7 gen DR = QdceCompCvpRange_cl_num7 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES7max = max(QDCERES7bev, QDCERES7cvd, QDCERES7cvp, QDCERES7optout, QDCERES7phev) *br QDCERES7max generate choice = 0 replace choice = 1 if QDCERES7max==QDCERES7cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES7cvp * Count how many alternatives have the maximum utility generate num_max7 = (QDCERES7bev == QDCERES7max) + /// (QDCERES7cvd == QDCERES7max) + /// (QDCERES7cvp == QDCERES7max) + /// (QDCERES7optout == QDCERES7max) + /// (QDCERES7phev == QDCERES7max) // Drop observations where num_max is greater than 1 for the current task drop if num_max7 > 1 save pair71.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=7 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num7 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num7 gen narok = 0 gen cost = QdceCompCvdCost_cl_num7 gen DR = QdceCompCvdRange_cl_num7 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES7max = max(QDCERES7bev, QDCERES7cvd, QDCERES7cvp, QDCERES7optout, QDCERES7phev) *br QDCERES7max generate choice = 0 replace choice = 1 if QDCERES7max==QDCERES7cvd *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES7cvd * Count how many alternatives have the maximum utility generate num_max7 = (QDCERES7bev == QDCERES7max) + /// (QDCERES7cvd == QDCERES7max) + /// (QDCERES7cvp == QDCERES7max) + /// (QDCERES7optout == QDCERES7max) + /// (QDCERES7phev == QDCERES7max) // Drop observations where num_max is greater than 1 for the current task drop if num_max7 > 1 save pair72.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=7 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num7 gen price_subs = QdceCompBevSubsidy_cl_num7 gen price = QdceCompBevPriceSub_cl_num7 gen narok = QdceBevNarok7 gen cost = QdceCompBevCost_cl_num7 gen DR = 0 gen DRe = QdceCompBevRange_cl_num7 gen pubcharge = QdceBevCharge7 gen polv2g = QdceBevV2g7 gen polwallbox = QdceBevSubskw7 generate QDCERES7max = max(QDCERES7bev, QDCERES7cvd, QDCERES7cvp, QDCERES7optout, QDCERES7phev) *br QDCERES7max generate choice = 0 replace choice = 1 if QDCERES7max==QDCERES7bev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES7bev * Count how many alternatives have the maximum utility generate num_max7 = (QDCERES7bev == QDCERES7max) + /// (QDCERES7cvd == QDCERES7max) + /// (QDCERES7cvp == QDCERES7max) + /// (QDCERES7optout == QDCERES7max) + /// (QDCERES7phev == QDCERES7max) // Drop observations where num_max is greater than 1 for the current task drop if num_max7 > 1 save pair73.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=7 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num7 gen price_subs = QdceCompPhevSubsidy_cl_num7 gen price = QdceCompPhevPriceSub_cl_num7 gen narok = QdcePhevNarok7 gen cost = QdceCompPhevCost_cl_num7 gen DR = QdceCompPhevCvRange_cl_num7 gen DRe = QdceCompPhevElRange_cl_num7 gen pubcharge = QdcePhevCharge7 gen polv2g = 0 gen polwallbox = 0 generate QDCERES7max = max(QDCERES7bev, QDCERES7cvd, QDCERES7cvp, QDCERES7optout, QDCERES7phev) *br QDCERES7max generate choice = 0 replace choice = 1 if QDCERES7max==QDCERES7phev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES7phev * Count how many alternatives have the maximum utility generate num_max7 = (QDCERES7bev == QDCERES7max) + /// (QDCERES7cvd == QDCERES7max) + /// (QDCERES7cvp == QDCERES7max) + /// (QDCERES7optout == QDCERES7max) + /// (QDCERES7phev == QDCERES7max) // Drop observations where num_max is greater than 1 for the current task drop if num_max7 > 1 save pair74.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=7 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES7max = max(QDCERES7bev, QDCERES7cvd, QDCERES7cvp, QDCERES7optout, QDCERES7phev) *br QDCERES7max generate choice = 0 replace choice = 1 if QDCERES7max==QDCERES7optout *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES7optout * Count how many alternatives have the maximum utility generate num_max7 = (QDCERES7bev == QDCERES7max) + /// (QDCERES7cvd == QDCERES7max) + /// (QDCERES7cvp == QDCERES7max) + /// (QDCERES7optout == QDCERES7max) + /// (QDCERES7phev == QDCERES7max) // Drop observations where num_max is greater than 1 for the current task drop if num_max7 > 1 save pair75.dta br num_max7 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie7 = (num_max7 > 1) br has_tie7 * Optional: List the choice sets with ties for verification list QDCERES7bev QDCERES7cvd QDCERES7cvp QDCERES7optout QDCERES7phev has_tie7 if has_tie7==1 *list QDCERES7bev QDCERES7cvd QDCERES7cvp QDCERES7optout QDCERES7phev has_tie7 if has_tie7==0 **** eighth choice set ***** ******** CVp ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=8 generate altern=1 gen price_pre = QdceCompCvpPrice_cl_num8 gen price_subs = 0 gen price = QdceCompCvpPrice_cl_num8 gen narok = 0 gen cost = QdceCompCvpCost_cl_num8 gen DR = QdceCompCvpRange_cl_num8 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES8max = max(QDCERES8bev, QDCERES8cvd, QDCERES8cvp, QDCERES8optout, QDCERES8phev) *br QDCERES8max generate choice = 0 replace choice = 1 if QDCERES8max==QDCERES8cvp *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES8cvp * Count how many alternatives have the maximum utility generate num_max8 = (QDCERES8bev == QDCERES8max) + /// (QDCERES8cvd == QDCERES8max) + /// (QDCERES8cvp == QDCERES8max) + /// (QDCERES8optout == QDCERES8max) + /// (QDCERES8phev == QDCERES8max) // Drop observations where num_max is greater than 1 for the current task drop if num_max8 > 1 save pair81.dta ******** CVd ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=8 generate altern=2 gen price_pre = QdceCompCvdPrice_cl_num8 gen price_subs = 0 gen price = QdceCompCvdPrice_cl_num8 gen narok = 0 gen cost = QdceCompCvdCost_cl_num8 gen DR = QdceCompCvdRange_cl_num8 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES8max = max(QDCERES8bev, QDCERES8cvd, QDCERES8cvp, QDCERES8optout, QDCERES8phev) *br QDCERES8max generate choice = 0 replace choice = 1 if QDCERES8max==QDCERES8cvd *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES8cvd * Count how many alternatives have the maximum utility generate num_max8 = (QDCERES8bev == QDCERES8max) + /// (QDCERES8cvd == QDCERES8max) + /// (QDCERES8cvp == QDCERES8max) + /// (QDCERES8optout == QDCERES8max) + /// (QDCERES8phev == QDCERES8max) // Drop observations where num_max is greater than 1 for the current task drop if num_max8 > 1 save pair82.dta ******** BEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=8 generate altern=3 gen price_pre = QdceCompBevPrice_cl_num8 gen price_subs = QdceCompBevSubsidy_cl_num8 gen price = QdceCompBevPriceSub_cl_num8 gen narok = QdceBevNarok8 gen cost = QdceCompBevCost_cl_num8 gen DR = 0 gen DRe = QdceCompBevRange_cl_num8 gen pubcharge = QdceBevCharge8 gen polv2g = QdceBevV2g8 gen polwallbox = QdceBevSubskw8 generate QDCERES8max = max(QDCERES8bev, QDCERES8cvd, QDCERES8cvp, QDCERES8optout, QDCERES8phev) *br QDCERES8max generate choice = 0 replace choice = 1 if QDCERES8max==QDCERES8bev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES8bev * Count how many alternatives have the maximum utility generate num_max8 = (QDCERES8bev == QDCERES8max) + /// (QDCERES8cvd == QDCERES8max) + /// (QDCERES8cvp == QDCERES8max) + /// (QDCERES8optout == QDCERES8max) + /// (QDCERES8phev == QDCERES8max) // Drop observations where num_max is greater than 1 for the current task drop if num_max8 > 1 save pair83.dta ******** PHEV ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=8 generate altern=4 gen price_pre = QdceCompPhevPrice_cl_num8 gen price_subs = QdceCompPhevSubsidy_cl_num8 gen price = QdceCompPhevPriceSub_cl_num8 gen narok = QdcePhevNarok8 gen cost = QdceCompPhevCost_cl_num8 gen DR = QdceCompPhevCvRange_cl_num8 gen DRe = QdceCompPhevElRange_cl_num8 gen pubcharge = QdcePhevCharge8 gen polv2g = 0 gen polwallbox = 0 generate QDCERES8max = max(QDCERES8bev, QDCERES8cvd, QDCERES8cvp, QDCERES8optout, QDCERES8phev) *br QDCERES8max generate choice = 0 replace choice = 1 if QDCERES8max==QDCERES8phev *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES8phev * Count how many alternatives have the maximum utility generate num_max8 = (QDCERES8bev == QDCERES8max) + /// (QDCERES8cvd == QDCERES8max) + /// (QDCERES8cvp == QDCERES8max) + /// (QDCERES8optout == QDCERES8max) + /// (QDCERES8phev == QDCERES8max) // Drop observations where num_max is greater than 1 for the current task drop if num_max8 > 1 save pair84.dta ******** OPT OUT ********** /*use "pilot_clean.dta", clear*/ use "fleet_pilot.dta", clear generate pair=8 generate altern=5 gen price_pre = 0 gen price_subs = 0 gen price = 0 gen narok = 0 gen cost = 0 gen DR = 0 gen DRe = 0 gen pubcharge = 0 gen polv2g = 0 gen polwallbox = 0 generate QDCERES8max = max(QDCERES8bev, QDCERES8cvd, QDCERES8cvp, QDCERES8optout, QDCERES8phev) br QDCERES8max generate choice = 0 replace choice = 1 if QDCERES8max==QDCERES8optout *br choice *bhat generate choice_bhat = 0 replace choice_bhat = QDCERES8optout * Count how many alternatives have the maximum utility generate num_max8 = (QDCERES8bev == QDCERES8max) + /// (QDCERES8cvd == QDCERES8max) + /// (QDCERES8cvp == QDCERES8max) + /// (QDCERES8optout == QDCERES8max) + /// (QDCERES8phev == QDCERES8max) // Drop observations where num_max is greater than 1 for the current task drop if num_max8 > 1 save pair85.dta br num_max8 * Mark choice sets with ties (where more than one alternative has max utility) generate has_tie8 = (num_max8 > 1) br has_tie8 * Optional: List the choice sets with ties for verification list QDCERES8bev QDCERES8cvd QDCERES8cvp QDCERES8optout QDCERES8phev has_tie8 if has_tie8==1 list QDCERES8bev QDCERES8cvd QDCERES8cvp QDCERES8optout QDCERES8phev has_tie8 if has_tie8==0 ******************* *****merge files***** save dce_cz_mainall.dta, replace use "pair11" /*1st experiment, 4+1 */ append using pair12 append using pair13 append using pair14 append using pair15 append using pair21 append using pair22 append using pair23 append using pair24 append using pair25 append using pair31 append using pair32 append using pair33 append using pair34 append using pair35 append using pair41 append using pair42 append using pair43 append using pair44 append using pair45 append using pair51 append using pair52 append using pair53 append using pair54 append using pair55 append using pair61 append using pair62 append using pair63 append using pair64 append using pair65 append using pair71 append using pair72 append using pair73 append using pair74 append using pair75 append using pair81 append using pair82 append using pair83 append using pair84 append using pair85 /*create pid*/ generate pid=pair*10000+ourID sort ourID pair altern gen DRtot = DR + DRe ****constants (asc)***** gen asc_CVp = 0 replace asc_CVp = 1 if altern==1 gen asc_CVd = 0 replace asc_CVd = 1 if altern==2 gen asc_PHEV = 0 replace asc_PHEV = 1 if altern==4 gen asc_BEV = 0 replace asc_BEV = 1 if altern==3 gen asc_nocar= 0 replace asc_nocar=1 if altern==5 /*drop extra entries I do not need*/ *drop if QDCE_RUN==4 & (pair==15 | pair==16 | pair==17 | pair==18) *drop if QDCE_RUN==8 & (pair==21 | pair==22 | pair==23 | pair==24 | pair==25 | pair==26) save dce_cz_mainall.dta, replace * to check the data br ourID pair altern pid choice choice_bhat QDCE_RUN DR price price_pre Qend asc_CVp asc_CVd asc_BEV asc_PHEV asc_nocar save dce_cz_mainall.dta, replace ****************************************************** *********************ESTIMATION*********************** ***Estimation Poisson*** use "fleet_pilot.dta", clear * data * generate owner = (QIIQ1 == 1) // vlastník a ostatní tab owner/*vlastník*/ tab QIIQ2/*pravomoc*/ tab QIVQ4/*nájezd*/ summarize QIVQ4 tab QVIIQ1/*provozfve*/ tab QVIIQ8/*ESG*/ tab QVIIIQ1/*obrat*/ tab QVIIQ1 QVIIQ8 /*provozfve*/ /*ESG*/ **PLÁNUJÍ** * bev ************* summarize QIVQ2_e_7, detail summarize QIVQ3_e_7, detail tab QIVQ2_e_7 tab QIVQ3_e_7 * Fit a Poisson regression model poisson QIVQ2_e_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estimates store poisson_model2 estat ic estat gof poisson QIVQ3_e_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estimates store poisson_model3 estat ic estat gof * Fit begative binomial reg model nbreg QIVQ2_e_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estimates store nb_model2 estat ic nbreg QIVQ3_e_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estimates store nb_model3 estat ic * phev************* tab QIVQ2_d_7 tab QIVQ3_d_7 * Fit a Poisson regression model poisson QIVQ2_d_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof poisson QIVQ3_d_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof * Fit begative binomial reg model nbreg QIVQ2_d_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic nbreg QIVQ3_d_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic * ev=bev+phev********** gen bev_phev2 = QIVQ2_e_7 + QIVQ2_d_7 gen bev_phev3 = QIVQ3_e_7 + QIVQ3_d_7 tab bev_phev2 tab bev_phev3 * Fit a Poisson regression model poisson bev_phev2/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof poisson bev_phev3/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof * Fit begative binomial reg model nbreg bev_phev2/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic nbreg bev_phev3/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic **MAJÍ** tab QIIIQ2_e_7 //bev tab QIIIQ3_e_7 //bev tab QIIIQ2_d_7 //phev tab QIIIQ3_d_7 //phev * Fit a Poisson regression model poisson QIIIQ2_e_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof poisson QIIIQ3_e_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof * Fit begative binomial reg model nbreg QIIIQ2_e_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof nbreg QIIIQ3_e_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic * phev************* tab QIIIQ2_d_7 tab QIIIQ3_d_7 * Fit a Poisson regression model poisson QIIIQ2_d_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic poisson QIIIQ3_d_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic * Fit begative binomial reg model nbreg QIIIQ2_d_7/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic nbreg QIIIQ3_d_7/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic * ev=bev+phev********** gen bev_phev2_maji = QIIIQ2_e_7 + QIIIQ2_d_7 gen bev_phev3_maji = QIIIQ3_e_7 + QIIIQ3_d_7 tab bev_phev2_maji tab bev_phev3_maji * Fit a Poisson regression model poisson bev_phev2_maji/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof poisson bev_phev3_maji/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic estat gof * Fit begative binomial reg model nbreg bev_phev2_maji/*own*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic nbreg bev_phev3_maji/*lease*/ owner/*vlastník*/ QIIQ2/*pravomoc*/ QIVQ4/*nájezd*/ QVIIQ1/*provozfve*/ QVIIQ8/*ESG*/ QVIIIQ1/*obrat*/, vce(robust) estat ic ****Estimation Logits**** use dce_cz_mainall.dta, clear br clogit choice asc_CVp asc_CVd asc_PHEV asc_BEV price_pre price_subs narok cost DR DRe pubcharge polv2g polwallbox if speeder==0, group(pid) cluster(ourID) /*all price_pre*/ estat ic clogit choice asc_CVp asc_CVd asc_PHEV asc_BEV price narok cost DR DRe pubcharge polv2g polwallbox if speeder==0, group(pid) cluster(ourID) /*price*/ estat ic clogit choice asc_CVp asc_CVd asc_PHEV asc_BEV price_pre narok cost DR DRe pubcharge polv2g polwallbox if speeder==0, group(pid) cluster(ourID) /*1price_pre*/ *install mixlogit* ssc install mixlogit mixlogit choice price_pre price_subs price /*if speeder==0*/, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok cost DR DRe pubcharge polv2g polwallbox) /*all*/ estat ic mixlogit choice price_pre if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV price_subs narok cost DR DRe pubcharge polv2g polwallbox) /*1price*/ estat ic mixlogit choice price_pre/*cena bez dotace*/ price_subs/*dotace*/ if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok cost DR DRe pubcharge polv2g polwallbox) /*-price main*/ estat ic mixlogit choice cost price_pre/*cena bez dotace*/ price_subs/*dotace*/ if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok DR DRe pubcharge polv2g polwallbox) /*COST*/ estat ic *interactions gen intPV=0 replace intPV=1 if QVIIQ1==1 gen cost_intPV = cost * intPV tab cost_intPV tab QVIIQ1 /*FVE current*/ mixlogit choice price_pre/*cena bez dotace*/ price_subs/*dotace*/ if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok cost cost_intPV DR DRe pubcharge polv2g polwallbox ) /**/ estat ic mixlogit choice cost cost_intPV price_pre/*cena bez dotace*/ price_subs/*dotace*/ if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok DR DRe pubcharge polv2g polwallbox ) /*cost_intPV*/ estat ic testnl _b[price_pre]- _b[price_subs] = 0 *interactions - insignificant gen ascbev_intPV = asc_BEV * intPV tab ascbev_intPV tab QVIIQ1 /*FVE current*/ mixlogit choice cost price_pre/*cena bez dotace*/ price_subs/*dotace*/ if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV ascbev_intPV narok DR DRe pubcharge polv2g polwallbox) /*ascbev_intPV*/ estat ic ***1Price**** mixlogit choice cost price if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok DR DRe pubcharge polv2g polwallbox) /*COST*/ estat ic mixlogit choice cost cost_intPV price if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok DR DRe pubcharge polv2g polwallbox) /*COST*/ estat ic mixlogit choice price if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV narok cost DR DRe pubcharge polv2g polwallbox) /*main price*/ estat ic mixlogit choice price if speeder==0, group(pid) id(ourID) nrep(50) /// rand(asc_CVp asc_CVd asc_PHEV asc_BEV ascbev_intPV narok cost DR DRe pubcharge polv2g polwallbox) /*ascbev_intPV*/ * WTP nlcom -_b[asc_CVp]/_b[price_pre] nlcom -_b[asc_CVd]/_b[price_pre] nlcom -_b[asc_PHEV]/_b[price_pre] nlcom -_b[asc_BEV]/_b[price_pre] nlcom -_b[narok]/_b[price_pre] nlcom -_b[cost]/_b[price_pre] nlcom -_b[DRtot]/_b[price_pre] nlcom -_b[DRe]/_b[price_pre] nlcom -_b[pubcharge]/_b[price_pre] nlcom -_b[polv2g]/_b[price_pre] nlcom -_b[polwallbox]/_b[price_pre] /* **************************************************** *********************QUESTIONS********************** import excel "2024-12-23_fleet_pilot.xlsx", sheet("2024-12-23_fleet_pilot") firstrow clear save "fleet_pilot.dta", replace use fleet_pilot.dta sum ourID *drop time rename QextID_time QextID_me rename Qend_time Qend_me drop *time drop *tim drop *ti rename QextID_me QextID_time rename Qend_me Qend_time *drop known TESTs drop if ourID==255 | ourID==256 | ourID==179 | ourID<=169 | ourID==191 *br QextID *drop nedokončené dotazníky drop if Qend==. *mediany, speedery gen day= substr(QextID_time,9,2) destring day, replace gen started_hour= substr(QextID_time,12,2) gen started_minute= substr(QextID_time,15,2) gen started_second= substr(QextID_time,18,2) gen finished_hour= substr(Qend_time,12,2) gen finished_minute= substr(Qend_time,15,2) gen finished_second= substr(Qend_time,18,2) destring started_hour, replace destring started_minute, replace destring started_second, replace destring finished_hour, replace destring finished_minute, replace destring finished_second, replace gen time=((finished_hour*60*60+finished_minute*60+finished_second)-(started_hour*60*60+started_minute*60+started_second))/60 replace time=0 if time<=0 sum time, d gen speeder_cz=1 if time<0.48*18.25 replace speeder_cz=0 if speeder_cz==. */ *VI tab QVIQ1 tab QVIQ2 tab QVIQ3 tab QVIQ4 tab QVIQ5 *VII tab QVIIQ1 tab QVIIQ2 tab QVIIQ3 tab QVIIQ4 tab QVIIQ5 tab QVIIQ6 *VIII tab QVIIIQ1 tab QVIIIQ2 tab QIIQ1 * data * generate owner = (QIIQ1 == 1) // vlastník a ostatní tab owner/*vlastník*/ tab QIIQ2/*pravomoc*/ tab QIVQ4/*nájezd*/ summarize QIVQ4 tab QVIIQ1/*provozfve*/ tab QVIIQ8/*ESG*/ tab QVIIIQ1/*obrat*/ tab QVIIQ1 QVIIQ8 /*provozfve*/ /*ESG*/ *******TABULKY******** /*NEW*/ /*vlastnictví*/ /*benzín*/ tab QIVQ2_a_1 /*mini*/ tab QIVQ2_a_2 /*malá*/ tab QIVQ2_a_3 /*střední*/ tab QIVQ2_a_4 /*velká*/ tab QIVQ2_a_5 /*SUV*/ tab QIVQ2_a_6 /*MPV*/ tab QIVQ2_a_7 /*CELKEM*/ /*nafta*/ tab QIVQ2_b_1 /*mini*/ tab QIVQ2_b_2 /*malá*/ tab QIVQ2_b_3 /*střední*/ tab QIVQ2_b_4 /*velká*/ tab QIVQ2_b_5 /*SUV*/ tab QIVQ2_b_6 /*MPV*/ tab QIVQ2_b_7 /*CELKEM*/ /*hybridní*/ tab QIVQ2_c_1 /*mini*/ tab QIVQ2_c_2 /*malá*/ tab QIVQ2_c_3 /*střední*/ tab QIVQ2_c_4 /*velká*/ tab QIVQ2_c_5 /*SUV*/ tab QIVQ2_c_6 /*MPV*/ tab QIVQ2_c_7 /*CELKEM*/ /*PHEV*/ tab QIVQ2_d_1 /*mini*/ tab QIVQ2_d_2 /*malá*/ tab QIVQ2_d_3 /*střední*/ tab QIVQ2_d_4 /*velká*/ tab QIVQ2_d_5 /*SUV*/ tab QIVQ2_d_6 /*MPV*/ tab QIVQ2_d_7 /*CELKEM*/ /*BEV*/ tab QIVQ2_e_1 /*mini*/ tab QIVQ2_e_2 /*malá*/ tab QIVQ2_e_3 /*střední*/ tab QIVQ2_e_4 /*velká*/ tab QIVQ2_e_5 /*SUV*/ tab QIVQ2_e_6 /*MPV*/ tab QIVQ2_e_7 /*CELKEM*/ /*LPG*/ tab QIVQ2_f_1 /*mini*/ tab QIVQ2_f_2 /*malá*/ tab QIVQ2_f_3 /*střední*/ tab QIVQ2_f_4 /*velká*/ tab QIVQ2_f_5 /*SUV*/ tab QIVQ2_f_6 /*MPV*/ tab QIVQ2_f_7 /*CELKEM*/ /*CNG*/ tab QIVQ2_g_1 /*mini*/ tab QIVQ2_g_2 /*malá*/ tab QIVQ2_g_3 /*střední*/ tab QIVQ2_g_4 /*velká*/ tab QIVQ2_g_5 /*SUV*/ tab QIVQ2_g_6 /*MPV*/ tab QIVQ2_g_7 /*CELKEM*/ /*vodík*/ tab QIVQ2_h_1 /*mini*/ tab QIVQ2_h_2 /*malá*/ tab QIVQ2_h_3 /*střední*/ tab QIVQ2_h_4 /*velká*/ tab QIVQ2_h_5 /*SUV*/ tab QIVQ2_h_6 /*MPV*/ tab QIVQ2_h_7 /*CELKEM*/ /*jiné*/ tab QIVQ2_i_1 /*mini*/ tab QIVQ2_i_2 /*malá*/ tab QIVQ2_i_3 /*střední*/ tab QIVQ2_i_4 /*velká*/ tab QIVQ2_i_5 /*SUV*/ tab QIVQ2_i_6 /*MPV*/ tab QIVQ2_i_7 /*CELKEM*/ /*leasing*/ /*benzín*/ tab QIVQ3_a_1 /*mini*/ tab QIVQ3_a_2 /*malá*/ tab QIVQ3_a_3 /*střední*/ tab QIVQ3_a_4 /*velká*/ tab QIVQ3_a_5 /*SUV*/ tab QIVQ3_a_6 /*MPV*/ tab QIVQ3_a_7 /*CELKEM*/ /*nafta*/ tab QIVQ3_b_1 /*mini*/ tab QIVQ3_b_2 /*malá*/ tab QIVQ3_b_3 /*střední*/ tab QIVQ3_b_4 /*velká*/ tab QIVQ3_b_5 /*SUV*/ tab QIVQ3_b_6 /*MPV*/ tab QIVQ3_b_7 /*CELKEM*/ /*hybridní*/ tab QIVQ3_c_1 /*mini*/ tab QIVQ3_c_2 /*malá*/ tab QIVQ3_c_3 /*střední*/ tab QIVQ3_c_4 /*velká*/ tab QIVQ3_c_5 /*SUV*/ tab QIVQ3_c_6 /*MPV*/ tab QIVQ3_c_7 /*CELKEM*/ /*PHEV*/ tab QIVQ3_d_1 /*mini*/ tab QIVQ3_d_2 /*malá*/ tab QIVQ3_d_3 /*střední*/ tab QIVQ3_d_4 /*velká*/ tab QIVQ3_d_5 /*SUV*/ tab QIVQ3_d_6 /*MPV*/ tab QIVQ3_d_7 /*CELKEM*/ /*BEV*/ tab QIVQ3_e_1 /*mini*/ tab QIVQ3_e_2 /*malá*/ tab QIVQ3_e_3 /*střední*/ tab QIVQ3_e_4 /*velká*/ tab QIVQ3_e_5 /*SUV*/ tab QIVQ3_e_6 /*MPV*/ tab QIVQ3_e_7 /*CELKEM*/ /*LPG*/ tab QIVQ3_f_1 /*mini*/ tab QIVQ3_f_2 /*malá*/ tab QIVQ3_f_3 /*střední*/ tab QIVQ3_f_4 /*velká*/ tab QIVQ3_f_5 /*SUV*/ tab QIVQ3_f_6 /*MPV*/ tab QIVQ3_f_7 /*CELKEM*/ /*CNG*/ tab QIVQ3_g_1 /*mini*/ tab QIVQ3_g_2 /*malá*/ tab QIVQ3_g_3 /*střední*/ tab QIVQ3_g_4 /*velká*/ tab QIVQ3_g_5 /*SUV*/ tab QIVQ3_g_6 /*MPV*/ tab QIVQ3_g_7 /*CELKEM*/ /*vodík*/ tab QIVQ3_h_1 /*mini*/ tab QIVQ3_h_2 /*malá*/ tab QIVQ3_h_3 /*střední*/ tab QIVQ3_h_4 /*velká*/ tab QIVQ3_h_5 /*SUV*/ tab QIVQ3_h_6 /*MPV*/ tab QIVQ3_h_7 /*CELKEM*/ /*jiné*/ tab QIVQ3_i_1 /*mini*/ tab QIVQ3_i_2 /*malá*/ tab QIVQ3_i_3 /*střední*/ tab QIVQ3_i_4 /*velká*/ tab QIVQ3_i_5 /*SUV*/ tab QIVQ3_i_6 /*MPV*/ tab QIVQ3_i_7 /*CELKEM*/ /*USED*/ /*vlastnictví*/ /*benzín*/ tab QIVQ2used_a_1 /*mini*/ tab QIVQ2used_a_2 /*malá*/ tab QIVQ2used_a_3 /*střední*/ tab QIVQ2used_a_4 /*velká*/ tab QIVQ2used_a_5 /*SUV*/ tab QIVQ2used_a_6 /*MPV*/ tab QIVQ2used_a_7 /*CELKEM*/ /*nafta*/ tab QIVQ2used_b_1 /*mini*/ tab QIVQ2used_b_2 /*malá*/ tab QIVQ2used_b_3 /*střední*/ tab QIVQ2used_b_4 /*velká*/ tab QIVQ2used_b_5 /*SUV*/ tab QIVQ2used_b_6 /*MPV*/ tab QIVQ2used_b_7 /*CELKEM*/ /*hybridní*/ tab QIVQ2used_c_1 /*mini*/ tab QIVQ2used_c_2 /*malá*/ tab QIVQ2used_c_3 /*střední*/ tab QIVQ2used_c_4 /*velká*/ tab QIVQ2used_c_5 /*SUV*/ tab QIVQ2used_c_6 /*MPV*/ tab QIVQ2used_c_7 /*CELKEM*/ /*PHEV*/ tab QIVQ2used_d_1 /*mini*/ tab QIVQ2used_d_2 /*malá*/ tab QIVQ2used_d_3 /*střední*/ tab QIVQ2used_d_4 /*velká*/ tab QIVQ2used_d_5 /*SUV*/ tab QIVQ2used_d_6 /*MPV*/ tab QIVQ2used_d_7 /*CELKEM*/ /*BEV*/ tab QIVQ2used_e_1 /*mini*/ tab QIVQ2used_e_2 /*malá*/ tab QIVQ2used_e_3 /*střední*/ tab QIVQ2used_e_4 /*velká*/ tab QIVQ2used_e_5 /*SUV*/ tab QIVQ2used_e_6 /*MPV*/ tab QIVQ2used_e_7 /*CELKEM*/ /*LPG*/ tab QIVQ2used_f_1 /*mini*/ tab QIVQ2used_f_2 /*malá*/ tab QIVQ2used_f_3 /*střední*/ tab QIVQ2used_f_4 /*velká*/ tab QIVQ2used_f_5 /*SUV*/ tab QIVQ2used_f_6 /*MPV*/ tab QIVQ2used_f_7 /*CELKEM*/ /*CNG*/ tab QIVQ2used_g_1 /*mini*/ tab QIVQ2used_g_2 /*malá*/ tab QIVQ2used_g_3 /*střední*/ tab QIVQ2used_g_4 /*velká*/ tab QIVQ2used_g_5 /*SUV*/ tab QIVQ2used_g_6 /*MPV*/ tab QIVQ2used_g_7 /*CELKEM*/ /*vodík*/ tab QIVQ2used_h_1 /*mini*/ tab QIVQ2used_h_2 /*malá*/ tab QIVQ2used_h_3 /*střední*/ tab QIVQ2used_h_4 /*velká*/ tab QIVQ2used_h_5 /*SUV*/ tab QIVQ2used_h_6 /*MPV*/ tab QIVQ2used_h_7 /*CELKEM*/ /*jiné*/ tab QIVQ2used_i_1 /*mini*/ tab QIVQ2used_i_2 /*malá*/ tab QIVQ2used_i_3 /*střední*/ tab QIVQ2used_i_4 /*velká*/ tab QIVQ2used_i_5 /*SUV*/ tab QIVQ2used_i_6 /*MPV*/ tab QIVQ2used_i_7 /*CELKEM*/ /*leasing*/ /*benzín*/ tab QIVQ3used_a_1 /*mini*/ tab QIVQ3used_a_2 /*malá*/ tab QIVQ3used_a_3 /*střední*/ tab QIVQ3used_a_4 /*velká*/ tab QIVQ3used_a_5 /*SUV*/ tab QIVQ3used_a_6 /*MPV*/ tab QIVQ3used_a_7 /*CELKEM*/ /*nafta*/ tab QIVQ3used_b_1 /*mini*/ tab QIVQ3used_b_2 /*malá*/ tab QIVQ3used_b_3 /*střední*/ tab QIVQ3used_b_4 /*velká*/ tab QIVQ3used_b_5 /*SUV*/ tab QIVQ3used_b_6 /*MPV*/ tab QIVQ3used_b_7 /*CELKEM*/ /*hybridní*/ tab QIVQ3used_c_1 /*mini*/ tab QIVQ3used_c_2 /*malá*/ tab QIVQ3used_c_3 /*střední*/ tab QIVQ3used_c_4 /*velká*/ tab QIVQ3used_c_5 /*SUV*/ tab QIVQ3used_c_6 /*MPV*/ tab QIVQ3used_c_7 /*CELKEM*/ /*PHEV*/ tab QIVQ3used_d_1 /*mini*/ tab QIVQ3used_d_2 /*malá*/ tab QIVQ3used_d_3 /*střední*/ tab QIVQ3used_d_4 /*velká*/ tab QIVQ3used_d_5 /*SUV*/ tab QIVQ3used_d_6 /*MPV*/ tab QIVQ3used_d_7 /*CELKEM*/ /*BEV*/ tab QIVQ3used_e_1 /*mini*/ tab QIVQ3used_e_2 /*malá*/ tab QIVQ3used_e_3 /*střední*/ tab QIVQ3used_e_4 /*velká*/ tab QIVQ3used_e_5 /*SUV*/ tab QIVQ3used_e_6 /*MPV*/ tab QIVQ3used_e_7 /*CELKEM*/ /*LPG*/ tab QIVQ3used_f_1 /*mini*/ tab QIVQ3used_f_2 /*malá*/ tab QIVQ3used_f_3 /*střední*/ tab QIVQ3used_f_4 /*velká*/ tab QIVQ3used_f_5 /*SUV*/ tab QIVQ3used_f_6 /*MPV*/ tab QIVQ3used_f_7 /*CELKEM*/ /*CNG*/ tab QIVQ3used_g_1 /*mini*/ tab QIVQ3used_g_2 /*malá*/ tab QIVQ3used_g_3 /*střední*/ tab QIVQ3used_g_4 /*velká*/ tab QIVQ3used_g_5 /*SUV*/ tab QIVQ3used_g_6 /*MPV*/ tab QIVQ3used_g_7 /*CELKEM*/ /*vodík*/ tab QIVQ3used_h_1 /*mini*/ tab QIVQ3used_h_2 /*malá*/ tab QIVQ3used_h_3 /*střední*/ tab QIVQ3used_h_4 /*velká*/ tab QIVQ3used_h_5 /*SUV*/ tab QIVQ3used_h_6 /*MPV*/ tab QIVQ3used_h_7 /*CELKEM*/ /*jiné*/ tab QIVQ3used_i_1 /*mini*/ tab QIVQ3used_i_2 /*malá*/ tab QIVQ3used_i_3 /*střední*/ tab QIVQ3used_i_4 /*velká*/ tab QIVQ3used_i_5 /*SUV*/ tab QIVQ3used_i_6 /*MPV*/ tab QIVQ3used_i_7 /*CELKEM*/ **************************************************** *********************PICTURES*********************** use dce_cz_mainall.dta, clear // Inicializace proměnných gen maxIVQ2 = QIVQ2_j_7 gen maxIVQ3 = QIVQ3_j_7 // Pokud jsou hodnoty missing, nastavíme je na 0 replace maxIVQ2 = 0 if missing(maxIVQ2) replace maxIVQ3 = 0 if missing(maxIVQ3) gen ownership = "" // Create a new variable to store the result replace ownership = "own" if maxIVQ2 > 0 // Assign 'own' if maxIVQ2 > 0 replace ownership = "rent" if maxIVQ3 > 0 & maxIVQ2 <= 0 // Assign 'rent' if maxIVQ3 > 0 and maxIVQ2 is not greater than 0 *drop if ownership == "rent" // Step 1: Create a separate variable for each alternative gen price_1_cvp = price if altern == 1 & ownership == "rent" //CVp gen price_2_cvd = price if altern == 2 & ownership == "rent" //CVd gen price_3_bev = price if altern == 3 & ownership == "rent" //BEV gen price_4_phev = price if altern == 4 & ownership == "rent" //PHEV gen price_5_optout = price if altern == 5 & ownership == "rent" // Optional: Opt-out alternative if needed twoway (histogram price_1_cvp, frequency color(orange%50) lwidth(vthin) title("CVp")) /// (histogram price_2_cvd, frequency color(yellow%50) lwidth(vthin) title("CVd")) /// (histogram price_3_bev, frequency color(green%50) lwidth(vthin) title("BEV")) /// (histogram price_4_phev, frequency color(blue%50) lwidth(vthin) title("PHEV")) /// , title("Price Distributions for Vehicle Alternatives") /// subtitle("Comparing CVp, CVd, BEV, and PHEV") /// ytitle("Frequency") /// xtitle("Price (CZK)") /// ylabel(, format(%9.1f)) /// legend(order(1 "CVp" 2 "CVd" 3 "BEV" 4 "PHEV") /// rows(1) size(small) position(6)) /// graphregion(color(white)) bgcolor(white) *means summarize price_1_cvp summarize price_2_cvd summarize price_3_bev summarize price_4_phev use "fleet_pilot.dta", clear // Tabulate the frequencies of QDCE_SIZE tabulate QDCE_SIZE, generate(size_) br QDCE_SIZE // Create a bar chart using actual frequencies graph bar (count), over(QDCE_SIZE, label(angle(45))) /// bar(1, color(blue%60)) bar(2, color(orange%60)) bar(3, color(green%60)) bar(4, color(red%60)) /// ytitle("Frequency") /// title("Vehicle Size Distribution") /// subtitle("Counts of 'small', 'medium', 'large', and 'utility' vehicles") /// graphregion(color(white)) bgcolor(white)