/*------------------------------------------------------------------------------ Stata do file for Kim, M-K. and Jakus, P.M. (2019) Wildfire, national park visitation, and changes in regional economic activity, Journal of Outdoor Recreation and Tourism 26:34-42 FIGURE 1 Monthly visitation to NPS in Utah TABLE 1 Visitation and wildfire data by NP and radius of wildfire zone TABLE 2 Corrected AICs TABLE 3 Best relative NPS visitation models (estimation results) ------------------------------------------------------------------------------*/ clear all cd "E:\Box Sync\Box Sync\01research\wildfire\NPS_visit\stata" use "NPS_visitors", clear /*----- NPS_visitors.dta includes monthly data for gas_p : gasoline price cpi : consumer price index v_arch : number of visitors in Arches NP v_bryce : number of visitors in Bryce Canyon NP v_canyon : number of visitors in Canyonlands NP v_capitol: number of visitors in Capitol Reef NP v_zion : number of visitors in Zion NP t : trend 1 to 273 recession: recession dummy eventmonth: year month r_income: real income mighty 5: mighty 5 dummy */ * Declare data to be time series tsset eventmonth * Call wildfire data and merge merge 1:1 eventmonth using "m_wildfire_AZ CO NM NV UT_data.dta" drop _merge /* m_wildfire_AZ CO NM NV UT_data.dta includes monthly wildfire data with various distances to NPs, for example wf080_arch is the "ACRES" burned by wildfire with radious of 80 km */ * missing data or no fire replace wf080_arch = 0 if missing(wf080_arch) replace wf080_bryc = 0 if missing(wf080_bryc) replace wf080_cany = 0 if missing(wf080_cany) replace wf080_capi = 0 if missing(wf080_capi) replace wf080_zion = 0 if missing(wf080_zion) replace wf160_arch = 0 if missing(wf160_arch) replace wf160_bryc = 0 if missing(wf160_bryc) replace wf160_cany = 0 if missing(wf160_cany) replace wf160_capi = 0 if missing(wf160_capi) replace wf160_zion = 0 if missing(wf160_zion) replace wf320_arch = 0 if missing(wf320_arch) replace wf320_bryc = 0 if missing(wf320_bryc) replace wf320_cany = 0 if missing(wf320_cany) replace wf320_capi = 0 if missing(wf320_capi) replace wf320_zion = 0 if missing(wf320_zion) * Acres to hectare: 1 acre = 0.404686 hectare *-----> Area burned in 1000 hectare * 080 --> wildfire <= 80 km gen wf080_arch000h = wf080_arch*0.404686/1000 gen wf080_bryc000h = wf080_bryc*0.404686/1000 gen wf080_cany000h = wf080_cany*0.404686/1000 gen wf080_capi000h = wf080_capi*0.404686/1000 gen wf080_zion000h = wf080_zion*0.404686/1000 * 160 --> 80 km < wildifre <= 160 km gen wf160_arch000h = wf160_arch*0.404686/1000 gen wf160_bryc000h = wf160_bryc*0.404686/1000 gen wf160_cany000h = wf160_cany*0.404686/1000 gen wf160_capi000h = wf160_capi*0.404686/1000 gen wf160_zion000h = wf160_zion*0.404686/1000 * 320 --> 160 km < wildifre <= 320 km gen wf320_arch000h = wf320_arch*0.404686/1000 gen wf320_bryc000h = wf320_bryc*0.404686/1000 gen wf320_cany000h = wf320_cany*0.404686/1000 gen wf320_capi000h = wf320_capi*0.404686/1000 gen wf320_zion000h = wf320_zion*0.404686/1000 * Different zoning * wildfire <= 160 km gen wfin160_arch000h = wf080_arch000h + wf160_arch000h gen wfin160_bryc000h = wf080_bryc000h + wf160_bryc000h gen wfin160_cany000h = wf080_cany000h + wf160_cany000h gen wfin160_capi000h = wf080_capi000h + wf160_capi000h gen wfin160_zion000h = wf080_zion000h + wf160_zion000h * wildfire <= 320 km gen wfin320_arch000h = wf080_arch000h + wf160_arch000h + wf320_arch000h gen wfin320_bryc000h = wf080_bryc000h + wf160_bryc000h + wf320_bryc000h gen wfin320_cany000h = wf080_cany000h + wf160_cany000h + wf320_cany000h gen wfin320_capi000h = wf080_capi000h + wf160_capi000h + wf320_capi000h gen wfin320_zion000h = wf080_zion000h + wf160_zion000h + wf320_zion000h * Check first five obs list eventmonth v_zion wf080_zion000h wf160_zion000h wf320_zion000 /// wfin160_zion000h wfin320_zion000h in 1/5 * Data manupulation ------------------------------------- * Logged visitors gen ln_v_arch = ln(v_arch) gen ln_v_bryc = ln(v_bryce) gen ln_v_cany = ln(v_canyon) gen ln_v_capi = ln(v_capitol) gen ln_v_zion = ln(v_zion) * income in 1000$ gen r_inc000 = r_income/1000 * real gas price gen r_gas_p = (gas_p*100/cpi) * income adjusted gas price gen adj_r_gas_p = r_gas_p/r_inc000 * month and year gen m = month(dofm(eventmonth)) gen yr = year(dofm(eventmonth)) * visitors in 1000 persons for Fig 1 gen v_arch000 = v_arch/1000 gen v_bryc000 = v_bryc/1000 gen v_cany000 = v_cany/1000 gen v_capi000 = v_capi/1000 gen v_zion000 = v_zion/1000 *-------------------- * FIGURE 1 *-------------------- twoway (tsline v_arch000, lwidth(medthick) lcolor(gs5)) /// (qfit v_arch000 eventmonth, lcolor(gs0)) /// , graphregion(fcolor(white)) ytitle("000 Persons") xtitle("") legend(off) /// ylabel(0(100)300) title("Arches NP", size(med)) xlabel(400(54)671, format(%tmCY)) /// saving(v_arch,replace) twoway (tsline v_bryc000, lwidth(medthick) lcolor(gs5)) /// (qfit v_bryc000 eventmonth, lcolor(gs0)) /// , graphregion(fcolor(white)) ytitle("000 Persons") xtitle("") legend(off) /// title("Bryce Canyon NP", size(med)) ylabel(0(100)400) xlabel(400(54)671, format(%tmCY)) /// saving(v_bryc,replace) twoway (tsline v_cany000, lwidth(medthick) lcolor(gs5)) /// (qfit v_cany000 eventmonth, lcolor(gs0)) /// , graphregion(fcolor(white)) ytitle("000 Persons") xtitle("") legend(off) /// title("Canyonlands NP", size(med)) ylabel(0(50)200) xlabel(400(54)671, format(%tmCY)) /// saving(v_cany,replace) twoway (tsline v_capi000, lwidth(medthick) lcolor(gs5)) /// (qfit v_capi000 eventmonth, lcolor(gs0)) /// , graphregion(fcolor(white)) ytitle("000 Persons") xtitle("") legend(off) /// title("Capitol Reef NP", size(med)) ylabel(0(50)200) xlabel(400(54)671, format(%tmCY)) /// saving(v_capi,replace) twoway (tsline v_zion000, lwidth(medthick) lcolor(gs5)) /// (qfit v_zion000 eventmonth, lcolor(gs0)) /// , graphregion(fcolor(white)) ytitle("000 Persons") xtitle("") legend(off) /// title("Zion NP", size(med)) ylabel(0(100)500) xlabel(400(54)671, format(%tmCY)) /// saving(v_zion,replace) gr combine v_zion.gph v_bryc.gph v_arch.gph v_capi.gph v_cany.gph /// , graphregion(fcolor(white)) saving(figure_1,replace) *-------------------------- * TABLE 1: summary statistics *--------------------------- * Arches NP sum v_arch, detail * 80 km radius sum wf080_arch000h if m == 5, d sum wf080_arch000h if m == 6, d sum wf080_arch000h if m == 7, d sum wf080_arch000h if m == 8, d sum wf080_arch000h if m == 9, d * 160 km radius sum wfin160_arch000h if m == 5, d sum wfin160_arch000h if m == 6, d sum wfin160_arch000h if m == 7, d sum wfin160_arch000h if m == 8, d sum wfin160_arch000h if m == 9, d * 320 km radius sum wfin320_arch000h if m == 5, d sum wfin320_arch000h if m == 6, d sum wfin320_arch000h if m == 7, d sum wfin320_arch000h if m == 8, d sum wfin320_arch000h if m == 9, d * Bryce Canyon NP sum v_bryc, detail * 80 km radius sum wf080_bryc000h if m == 5, d sum wf080_bryc000h if m == 6, d sum wf080_bryc000h if m == 7, d sum wf080_bryc000h if m == 8, d sum wf080_bryc000h if m == 9, d * 160 km radius sum wfin160_bryc000h if m == 5, d sum wfin160_bryc000h if m == 6, d sum wfin160_bryc000h if m == 7, d sum wfin160_bryc000h if m == 8, d sum wfin160_bryc000h if m == 9, d * 320 km radius sum wfin320_bryc000h if m == 5, d sum wfin320_bryc000h if m == 6, d sum wfin320_bryc000h if m == 7, d sum wfin320_bryc000h if m == 8, d sum wfin320_bryc000h if m == 9, d * Canyonlands NP sum v_cany, detail * 80 km radius sum wf080_cany000h if m == 5, d sum wf080_cany000h if m == 6, d sum wf080_cany000h if m == 7, d sum wf080_cany000h if m == 8, d sum wf080_cany000h if m == 9, d * 160 km radius sum wfin160_cany000h if m == 5, d sum wfin160_cany000h if m == 6, d sum wfin160_cany000h if m == 7, d sum wfin160_cany000h if m == 8, d sum wfin160_cany000h if m == 9, d * 320 km radius sum wfin320_cany000h if m == 5, d sum wfin320_cany000h if m == 6, d sum wfin320_cany000h if m == 7, d sum wfin320_cany000h if m == 8, d sum wfin320_cany000h if m == 9, d * Capitol Reef NP sum v_capi, detail * 80 km radius sum wf080_capi000h if m == 5, d sum wf080_capi000h if m == 6, d sum wf080_capi000h if m == 7, d sum wf080_capi000h if m == 8, d sum wf080_capi000h if m == 9, d * 160 km radius sum wfin160_capi000h if m == 5, d sum wfin160_capi000h if m == 6, d sum wfin160_capi000h if m == 7, d sum wfin160_capi000h if m == 8, d sum wfin160_capi000h if m == 9, d * 320 km radius sum wfin320_capi000h if m == 5, d sum wfin320_capi000h if m == 6, d sum wfin320_capi000h if m == 7, d sum wfin320_capi000h if m == 8, d sum wfin320_capi000h if m == 9, d * Zion NP sum v_zion, detail * 80 km radius sum wf080_zion000h if m == 5, d sum wf080_zion000h if m == 6, d sum wf080_zion000h if m == 7, d sum wf080_zion000h if m == 8, d sum wf080_zion000h if m == 9, d * 160 km radius sum wfin160_zion000h if m == 5, d sum wfin160_zion000h if m == 6, d sum wfin160_zion000h if m == 7, d sum wfin160_zion000h if m == 8, d sum wfin160_zion000h if m == 9, d * 320 km radius sum wfin320_zion000h if m == 5, d sum wfin320_zion000h if m == 6, d sum wfin320_zion000h if m == 7, d sum wfin320_zion000h if m == 8, d sum wfin320_zion000h if m == 9, d /*============================================================================ ESTIMATION First obs: 1993 m4 t = 1 Last obs: 2015 m12 t = 273 wildfire activity areas burned in 80 km (50 miles) to the visitor centers wf080_****000h wildfire **** park name in 000 hetare Three functional forms Semi-log form that is linear in wildfire size Semi-log form that is quadratic in wildfire size Log-log from that uses the natural log of wildifre And three distances, 80km, 160km and 320km ============================================================================== */ *========================= * ARCHES *========================= * additional data manupluation for nonlinear terms gen wf080_arch000h_sq = wf080_arch000h*wf080_arch000h gen ln_wf080_arch000h = ln(wf080_arch000h + 1) gen wfin160_arch000h_sq = wfin160_arch000h*wfin160_arch000h gen ln_wfin160_arch000h = ln(wfin160_arch000h + 1) gen wfin320_arch000h_sq = wfin320_arch000h*wfin320_arch000h gen ln_wfin320_arch000h = ln(wfin320_arch000h + 1) * ONE LAG * 80 km * semi-log qui reg ln_v_arch L(0/1).wf080_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/1).wf080_arch000h L(0/1).wf080_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log reg ln_v_arch L(0/1).ln_wf080_arch000h adj_r_gas_p mighty5 t recession b12.m estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_arch L(0/1).wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/1).wfin160_arch000h L(0/1).wfin160_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/1).ln_wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_arch L(0/1).wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/1).wfin320_arch000h L(0/1).wfin320_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/1).ln_wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * TWO LAGS * 80 km * semi-log qui reg ln_v_arch L(0/2).wf080_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/2).wf080_arch000h L(0/2).wf080_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/2).ln_wf080_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_arch L(0/2).wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/2).wfin160_arch000h L(0/2).wfin160_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/2).ln_wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_arch L(0/2).wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/2).wfin320_arch000h L(0/2).wfin320_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/2).ln_wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * THREE LAGS * 80 km * semi-log qui reg ln_v_arch L(0/3).wf080_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/3).wf080_arch000h L(0/3).wf080_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/3).ln_wf080_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_arch L(0/3).wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/3).wfin160_arch000h L(0/3).wfin160_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/3).ln_wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_arch L(0/3).wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_arch L(0/3).wfin320_arch000h L(0/3).wfin320_arch000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_arch L(0/3).ln_wfin320_arch000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) *========================== * BRYCE *========================= gen wf080_bryc000h_sq = wf080_bryc000h*wf080_bryc000h gen ln_wf080_bryc000h = ln(wf080_bryc000h + 1) gen wfin160_bryc000h_sq = wfin160_bryc000h*wfin160_bryc000h gen ln_wfin160_bryc000h = ln(wfin160_bryc000h + 1) gen wfin320_bryc000h_sq = wfin320_bryc000h*wfin320_bryc000h gen ln_wfin320_bryc000h = ln(wfin320_bryc000h + 1) * ONE LAG *--------- * 80 km * semi-log qui reg ln_v_bryc L(0/1).wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/1).wf080_bryc000h L(0/1).wf080_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log reg ln_v_bryc L(0/1).ln_wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km - ONE lag * semi-log qui reg ln_v_bryc L(0/1).wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/1).wfin160_bryc000h L(0/1).wfin160_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/1).ln_wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km - ONE lag * semi-log qui reg ln_v_bryc L(0/1).wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/1).wfin320_bryc000h L(0/1).wfin320_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/1).ln_wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * TWO LAGS * 80 km * semi-log qui reg ln_v_bryc L(0/2).wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/2).wf080_bryc000h L(0/2).wf080_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/2).ln_wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_bryc L(0/2).wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/2).wfin160_bryc000h L(0/2).wfin160_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/2).ln_wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_bryc L(0/2).wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/2).wfin320_bryc000h L(0/2).wfin320_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/2).ln_wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * THREE LAGS * 80 km * semi-log qui reg ln_v_bryc L(0/3).wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/3).wf080_bryc000h L(0/3).wf080_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/3).ln_wf080_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_bryc L(0/3).wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/3).wfin160_bryc000h L(0/3).wfin160_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/3).ln_wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_bryc L(0/3).wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_bryc L(0/3).wfin320_bryc000h L(0/3).wfin320_bryc000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_bryc L(0/3).ln_wfin320_bryc000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) *========================== * CANYONLANDS *========================= gen wf080_cany000h_sq = wf080_cany000h*wf080_cany000h gen ln_wf080_cany000h = ln(wf080_cany000h + 1) gen wfin160_cany000h_sq = wfin160_cany000h*wfin160_cany000h gen ln_wfin160_cany000h = ln(wfin160_cany000h + 1) gen wfin320_cany000h_sq = wfin320_cany000h*wfin320_cany000h gen ln_wfin320_cany000h = ln(wfin320_cany000h + 1) * ONE LAG * 80 km * semi-log qui reg ln_v_cany L(0/1).wf080_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/1).wf080_cany000h L(0/1).wf080_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log reg ln_v_cany L(0/1).ln_wf080_cany000h adj_r_gas_p mighty5 t recession b12.m estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km - ONE lag * semi-log qui reg ln_v_cany L(0/1).wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/1).wfin160_cany000h L(0/1).wfin160_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/1).ln_wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km - ONE lag * semi-log qui reg ln_v_cany L(0/1).wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/1).wfin320_cany000h L(0/1).wfin320_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/1).ln_wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * TWO LAGS * 80 km * semi-log qui reg ln_v_cany L(0/2).wf080_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/2).wf080_cany000h L(0/2).wf080_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/2).ln_wf080_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_cany L(0/2).wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/2).wfin160_cany000h L(0/2).wfin160_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/2).ln_wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_cany L(0/2).wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/2).wfin320_cany000h L(0/2).wfin320_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/2).ln_wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * THREE LAGS * 80 km * semi-log qui reg ln_v_cany L(0/3).wf080_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/3).wf080_cany000h L(0/3).wf080_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/3).ln_wf080_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_cany L(0/3).wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/3).wfin160_cany000h L(0/3).wfin160_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/3).ln_wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_cany L(0/3).wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_cany L(0/3).wfin320_cany000h L(0/3).wfin320_cany000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_cany L(0/3).ln_wfin320_cany000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) *========================== * CAPITOL REEF *========================= gen wf080_capi000h_sq = wf080_capi000h*wf080_capi000h gen ln_wf080_capi000h = ln(wf080_capi000h + 1) gen wfin160_capi000h_sq = wfin160_capi000h*wfin160_capi000h gen ln_wfin160_capi000h = ln(wfin160_capi000h + 1) gen wfin320_capi000h_sq = wfin320_capi000h*wfin320_capi000h gen ln_wfin320_capi000h = ln(wfin320_capi000h + 1) * ONE LAG * 80 km * semi-log qui reg ln_v_capi L(0/1).wf080_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/1).wf080_capi000h L(0/1).wf080_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log reg ln_v_capi L(0/1).ln_wf080_capi000h adj_r_gas_p mighty5 t recession b12.m estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km - ONE lag * semi-log qui reg ln_v_capi L(0/1).wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/1).wfin160_capi000h L(0/1).wfin160_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/1).ln_wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km - ONE lag * semi-log qui reg ln_v_capi L(0/1).wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/1).wfin320_capi000h L(0/1).wfin320_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/1).ln_wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * TWO LAGS * 80 km * semi-log qui reg ln_v_capi L(0/2).wf080_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/2).wf080_capi000h L(0/2).wf080_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/2).ln_wf080_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_capi L(0/2).wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/2).wfin160_capi000h L(0/2).wfin160_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/2).ln_wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_capi L(0/2).wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/2).wfin320_capi000h L(0/2).wfin320_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/2).ln_wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * THREE LAGS * 80 km * semi-log qui reg ln_v_capi L(0/3).wf080_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/3).wf080_capi000h L(0/3).wf080_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/3).ln_wf080_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_capi L(0/3).wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/3).wfin160_capi000h L(0/3).wfin160_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/3).ln_wfin160_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_capi L(0/3).wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_capi L(0/3).wfin320_capi000h L(0/3).wfin320_capi000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_capi L(0/3).ln_wfin320_capi000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) *========================== * ZION *========================= gen wf080_zion000h_sq = wf080_zion000h*wf080_zion000h gen ln_wf080_zion000h = ln(wf080_zion000h + 1) gen wfin160_zion000h_sq = wfin160_zion000h*wfin160_zion000h gen ln_wfin160_zion000h = ln(wfin160_zion000h + 1) gen wfin320_zion000h_sq = wfin320_zion000h*wfin320_zion000h gen ln_wfin320_zion000h = ln(wfin320_zion000h + 1) * ONE LAG * 80 km * semi-log qui reg ln_v_zion L(0/1).wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/1).wf080_zion000h L(0/1).wf080_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/1).ln_wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km - ONE lag * semi-log qui reg ln_v_zion L(0/1).wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/1).wfin160_zion000h L(0/1).wfin160_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/1).ln_wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km - ONE lag * semi-log qui reg ln_v_zion L(0/1).wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/1).wfin320_zion000h L(0/1).wfin320_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/1).ln_wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * TWO LAGS * 80 km * semi-log qui reg ln_v_zion L(0/2).wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/2).wf080_zion000h L(0/2).wf080_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/2).ln_wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_zion L(0/2).wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/2).wfin160_zion000h L(0/2).wfin160_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/2).ln_wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_zion L(0/2).wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/2).wfin320_zion000h L(0/2).wfin320_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/2).ln_wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * THREE LAGS * 80 km * semi-log qui reg ln_v_zion L(0/3).wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/3).wf080_zion000h L(0/3).wf080_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/3).ln_wf080_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 160 km * semi-log qui reg ln_v_zion L(0/3).wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/3).wfin160_zion000h L(0/3).wfin160_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/3).ln_wfin160_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * 320 km * semi-log qui reg ln_v_zion L(0/3).wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * semi-log Q qui reg ln_v_zion L(0/3).wfin320_zion000h L(0/3).wfin320_zion000h_sq adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) * log-log qui reg ln_v_zion L(0/3).ln_wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m qui estat ic display -2*e(ll) + 2*e(rank) + (2*e(rank)*(e(rank)+1))/(e(N)-e(rank)-1) *------------------------------------- * TABLE 3. Best NPS Visitation Models *------------------------------------- * Arches semi-log linear, 160 km, one lag newey ln_v_arch L(0/1).wfin160_arch000h adj_r_gas_p mighty5 t recession b12.m, lag(12) * Bryce semi-log linear, 160 km, one lag newey ln_v_bryc L(0/1).wfin160_bryc000h adj_r_gas_p mighty5 t recession b12.m, lag(12) * Canyonlands semi-log linear, 160 km, one lag newey ln_v_cany L(0/1).wfin160_cany000h adj_r_gas_p mighty5 t recession b12.m, lag(12) * Capitol Reef semi-log linear, 80 km, one lag newey ln_v_capi L(0/1).wf080_capi000h adj_r_gas_p mighty5 t recession b12.m, lag(12) * Zion log-log, 320 km, one lag newey ln_v_zion L(0/2).ln_wfin320_zion000h adj_r_gas_p mighty5 t recession b12.m , lag(12) *------------------------------------------------ * TABLE 4. Fixed Effects Panel (Pooled) Model *------------------------------------------------ * See 2_Kim_Jakus_2019_JORT26_Pooled.do * End of do file ===============================================================