Loading...
墨滴

Shiyuan Li

2021/09/29  阅读:55  主题:橙心

Simulation Method for Cross Model Comparison

1. For Categorical Groups

sysuse nlsw88,clear
capture program drop myboot_category
program define myboot_category, eclass
 preserve
 global X c.age##c.age c.ttl_exp##c.ttl_exp i.race i.married i.c_city i.south i.smsa i.occupation i.union

 gen lnwage=ln(wage)

 logit collgrad $X,asis
 predict pscore

 gen strata=1 if pscore>=0 & pscore<0.1
 replace strata=2 if pscore>=0.1 & pscore<0.2
 replace strata=3 if pscore>=0.2 & pscore<0.3
 replace strata=4 if pscore>=0.3 & pscore<0.6
 replace strata=5 if pscore>=0.6 & pscore<=1

 reg lnwage collgrad $X if strata==1
 scalar b1=_b[collgrad]
 reg lnwage collgrad $X if strata==2
 scalar b2=_b[collgrad]
 reg lnwage collgrad $X if strata==3
 scalar b3=_b[collgrad]
 reg lnwage collgrad $X if strata==4
 scalar b4=_b[collgrad]
 reg lnwage collgrad $X if strata==5
 scalar b5=_b[collgrad]
end

cd "D:"  //Note: This code specifies the place in which you want to save the simulation results 

bootstrap b1 b2 b3 b4 b5, saving(bsample.dta) nowarn reps(999) seed(101): myboot_category
test _bs_1=_bs_2=_bs_3 =_bs_4=_bs_5

2. For Ordinal Groups

*---Preparation--*
ssc install runby 

capture program drop sub_reg
 program define sub_reg
 reg lnwage collgrad $X
 gen b1= _b[collgrad]
end

global X ///
c.age##c.age c.ttl_exp##c.ttl_exp i.race i.married i.c_city ///
i.south i.smsa i.occupation i.union

*---Conducting the Test--*
sysuse nlsw88,clear

capture program drop myboot_ordinal
 program define myboot_ordinal, eclass
 preserve

 gen lnwage=ln(wage)

 logit collgrad $X,asis
 predict pscore

 gen strata=1 if pscore>=0 & pscore<0.1
 replace strata=2 if pscore>=0.1 & pscore<0.2
 replace strata=3 if pscore>=0.2 & pscore<0.3
 replace strata=4 if pscore>=0.3 & pscore<0.6
 replace strata=5 if pscore>=0.6 & pscore<=1

 runby sub_reg,by(strata) 
 collapse (mean) b1_j=b1 (percent) percent_j=b1, by(strata)

 reg b1_j strata [aweight=percent_j],r
 scalar slope=_b[strata]
end

bootstrap slope, nowarn reps(999) seed(10101): myboot_ordinal

Shiyuan Li

2021/09/29  阅读:55  主题:橙心

作者介绍

Shiyuan Li

PhD student GSM,PKU