///Deklarasi data panel
xtset ID Year, yearly

///Transformasi log
gen lnPOV=log(POV)
gen lnCapitalExp=log(CapitalExp)
gen lnMandatoryExp=log(MandatoryExp)

///Statistik Deskriptif
sum lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19

///Estimasi FEM
xtreg lnPOV l1.lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, fe robust

///Estimasi OLS
reg lnPOV l1.lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19

///Estimasi One Step First Difference GMM
xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2) robust
///Abond test  
estat abond
///Sargan Test
xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2) 
estat sargan

///Estimasi Two Step First Difference GMM
xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2) robust
///Abond test
estat abond
///Sargan Test
xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2)
estat sargan

///Estimasi One Step System GMM
xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2) vce(robust)
///Abond test
estat abond
///Sargan Test
xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2)
estat sargan

///Estimasi Two Step System GMM
xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2) vce(robust)
///Abond test
estat abond
///Sargan Test
xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2)
estat sargan


eststo FEM: xtreg lnPOV l1.lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, fe robust
eststo OLS: reg lnPOV l1.lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19
eststo OneStep_FD: xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2) robust
eststo TwoStep_FD: xtabond lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2) robust
eststo OneStep_Sys: xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, maxldep(1) maxlags(1) artest(2) vce(robust)
eststo TwoStep_Sys: xtdpdsys lnPOV LE PrimarySchooling SecondarySchooling lnCapitalExp lnMandatoryExp Unemployment Water DCov19, twostep maxldep(1) maxlags(1) artest(2) vce(robust)

///
estimates table FEM OLS OneStep_FD TwoStep_FD OneStep_Sys TwoStep_Sys, star stats (N r2 r2_a)

esttab FEM OLS OneStep_FD TwoStep_FD OneStep_Sys TwoStep_Sys, star stats (N r2 r2_a)

esttab FEM OLS OneStep_FD TwoStep_FD OneStep_Sys TwoStep_Sys using Estimasi_Model.rtf, b(3) se scalars("N Obs" "r2 R-sq" "r2_a Adj. R-sq") star(* 0.10 ** 0.05 *** 0.01) nonumbers nonotes mtitles replace title("Hasil Estimasi Regresi") addnote ("***p < 0.01, **p < 0.05, *p < 0.1") 

//////////////////////////////////////