* Reproduces results in Lake, HIR, Tables 5.3 and A.2
* use "Lake_HIR_ICB.dta" 

* Descriptive statistics for ICD data, Table A.2

sum us_SH1995 us_indally us_indexmilperpop1995 us_EH1995 us_indexexchangerate1995 us_indexreltradep1995 duration geo_salience violence usinv

* MULTINOMIAL LOGIT TEST
*drop b*
		*mlogit usinv us_SH1995 us_EH1995 duration violence geo_salience if usactor~=2, base(1) robust
		*estimates store a1
		estsimp mlogit usinv us_SH1995 us_EH1995 duration violence geo_salience if after_week==1, base(1) robust
		*estimates store a2
		
setx mean
setx us_SH1995 0 
simqi, prval(2) genpr(start)
setx us_SH1995 75
simqi, prval(2) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff

		 setx mean
setx us_SH1995 0 
simqi, prval(3) genpr(start)
setx us_SH1995 75
simqi, prval(3) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff


setx mean
setx us_SH1995 min
simqi, prval(2) genpr(pr1)
setx us_SH1995 p75
simqi, prval(2) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

setx mean
setx us_SH1995 min
simqi, prval(3) genpr(pr1)
setx us_SH1995 p75
simqi, prval(3) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

drop b*
		estsimp mlogit usinv us_indally us_indexmilperpop1995 us_indexexchangerate1995 us_indexreltradep1995 duration violence geo_salience if after_week==1, base(1) robust

		*estimates store a1

		 setx mean
setx us_indally 0 
simqi, prval(2) genpr(start)
setx us_indally 75
simqi, prval(2) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff

		 setx mean
setx us_indally 0 
simqi, prval(3) genpr(start)
setx us_indally 75
simqi, prval(3) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff

setx mean
setx us_indally min
simqi, prval(2) genpr(pr1)
setx us_indally p75
simqi, prval(2) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

setx mean
setx us_indally min
simqi, prval(3) genpr(pr1)
setx us_indally p75
simqi, prval(3) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

		 setx mean
setx us_indexreltradep1995 0 
simqi, prval(2) genpr(start)
setx us_indexreltradep1995 75
simqi, prval(2) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff

		 setx mean
setx us_indexreltradep1995 0 
simqi, prval(3) genpr(start)
setx us_indexreltradep1995 75
simqi, prval(3) genpr(end)
gen diff=end-start
sumqi diff, level(95)

drop start end diff


setx mean
setx us_indexreltradep1995 min
simqi, prval(2) genpr(pr1)
setx us_indexreltradep1995 p75
simqi, prval(2) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

setx mean
setx us_indexreltradep1995 min
simqi, prval(3) genpr(pr1)
setx us_indexreltradep1995 p75
simqi, prval(3) genpr(pr2)
gen relrisk = pr2/pr1
sumqi relrisk

drop pr1 pr2 relrisk

drop b*

		

		
*estout a*,  cells(b(star fmt(3)) se(par))   starlevels(* 0.10 ** 0.05 *** 0.01)   stardetach   stats(N, fmt(0) labels("Number of Observations"))   varwidth(18)   modelwidth(8)   mlabels("Model 1" "Model 2" "Model 3" "Model 4" "Model 5" "Model 6")  order(us_SH1995 us_indexreltradep1995 duration geo_salience number_actors coldwar)   varlabels(us_SH1995 "SH Score" duration "Duration" geo_salience "Salience" number_actors "Number of Actors" violence "Violence" us_indexreltradep1995 "Rel. Trade Dep." coldwar "Coldwar"_cons "Constant")  legend   prehead(@title "(Multinomial Logit Estimation, 2 = Diplomatic or Economic Activity, 3 = Indirect of Direct Military Activity)" "" @hline)  posthead(@hline)  prefoot("") postfoot(@hline)  style(fixed)  title("Security and Economic Hierarchy and the Probability of Joining a Crisis, 1950-2000") 
*estout b*,  cells(b(star fmt(3)) se(par))   starlevels(* 0.10 ** 0.05 *** 0.01)   stardetach   stats(N, fmt(0) labels("Number of Observations"))   varwidth(18)   modelwidth(8)   mlabels("Model 1" "Model 2" "Model 3" "Model 4" "Model 5" "Model 6")  order(us_SH1995 us_EH1995 duration geo_salience number_actors coldwar)   varlabels(us_SH1995 "SH Score" duration "Duration" geo_salience "Salience" number_actors "Number of Actors" violence "Violence" us_EH1995 "EH Score" coldwar "Coldwar"_cons "Constant")  legend   prehead(@title "(Multinomial Logit Estimation, 2 = Diplomatic or Economic Activity, 3 = Indirect of Direct Military Activity)" "" @hline)  posthead(@hline)  prefoot("") postfoot(@hline)  style(fixed)  title("Security and Economic Hierarchy and the Probability of Joining a Crisis, 1950-2000")