Wrote statistical codes on STATA for research purposes.
*Setting up STATA*
clear all
cd "
capture log close
log using abc.log, replace text
set more off (or, set more off, perm)
set memory 100m
pause on
log using pwd, replace
dis "`c(username)'"
dis "`c(machine_type)'"
dis "`c(os)'"
*Understanding what we're working with
codebook
describe
list in 1/10
summarize (var)
misstable summarize
misstable patterns
ds, has(type string) // describing string variables
*Understanding the merged variable.
su merge_1_2_3 if [merge_1_2_3 ==1]
su merge_1_2_3 if [merge_1_2_3 ==2]
su merge_1_2_3 if [merge_1_2_3 ==3]
su merge_1_2 if [merge_1_3 ==1]
su merge_1_2 if [merge_1_3 ==2]
su merge_1_2 if [merge_1_3 ==3]
su merge_2_3 if [merge_2_3 ==1]
su merge_2_3 if [merge_2_3 ==2]
su merge_2_3 if [merge_2_3 ==3]
*Drop merged variable (Not Matched)*
drop if [merge_1_2_3== 1]
drop if [merge_1_2_3== 2]
drop merge_2_3 merge_1_2
*DESCRIPTIVE STATS*
*Ratio on Settlement between Host and Refugee
*Demographic Variables*
sum ls02_age if merge_1_2_3 ==3, detail
tab (ls05_sex) (settlement_type) if merge_1_2_3 == 3
tab (upazila) (settlement_type) if merge_1_2_3 == 3
tab (ls09_religion) (settlement_type) if merge_1_2_3 == 3
tab (ls10_marital_status) (settlement_type) if merge_1_2_3 == 3
tab (ls06_relationship) (settlement_type) if merge_1_2_3 == 3
tab (hh_size_baseline_roster) (settlement_type) if merge_1_2_3 == 3
tab (ls18_school) (settlement_type) if merge_1_2_3 == 3
tab (ls16_read) (settlement_type) if merge_1_2_3 == 3
tab (ls17_write) (settlement_type) if merge_1_2_3 == 3
tab (num_working_mem) (settlement_type) if merge_1_2_3 ==3
tab (upazila) (settlement_type) if merge_1_2_3 == 3
tab (ls12_country_born) (settlement_type) if merge_1_2_3 == 3
*Labor force participation & employment status*
tab (l05_a_isic) (settlement_type) if merge_1_2_3 == 3
tab (l05_a_isco) (settlement_type) if merge_1_2_3 == 3
tab (ls25_working) (settlement_type) if merge_1_2_3 == 3
tab (ls26_available_work) (settlement_type) if merge_1_2_3 == 3
tab (ls27_lookfor_work) (settlement_type) if merge_1_2_3 == 3
tab (lm_work_7_days_r3) (settlement_type) if merge_1_2_3 ==3
tab (l08_agri_or_not) (settlement_type) if merge_1_2_3 ==3
tab (l06a_whynotsearch) (settlement_type) if merge_1_2_3 ==3
tab (l03_lastwork) (settlement_type) if merge_1_2_3 ==3
tab (l10_org_type) (settlement_type) if merge_1_2_3 ==3
tab (l10a_num_employees) (settlement_type) if merge_1_2_3 ==3
tab (l22_ent_chnge) (settlement_type) if merge_1_2_3 ==3
tab (l23_ent_problems) (settlement_type) if merge_1_2_3 ==3
tab (l18_ent_operations) (settlement_type) if merge_1_2_3 ==3
tab (l19_ent_ownership) (settlement_type) if merge_1_2_3 ==3
tab (l21_ent_place_operations) (settlement_type) if merge_1_2_3 ==3
*Earning & Income*
sum (l12_month_rate) if merge_1_2_3 ==3, detail
sum l14_7dayrate, if merge_1_2_3 ==3, detail
sum l15_7daylastrate if merge_1_2_3 ==3, detail
tabstat (l12_month_rate) (l14_7dayrate)(l15_7daylastrate) if merge_1_2_3==3, stat(mean median sd min max) col(stat)
tab (l09_wage_or_not) (settlement_type) if merge_1_2_3 ==3
tab (l04_lastpay) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_1) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_2) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_3) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_4) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_5) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_6) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_7) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_8) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src_9) (settlement_type) if merge_1_2_3 ==3
tab (l00a_current_inc_src__96) (settlement_type) if merge_1_2_3 ==3
*Working Hours & Job Condition*
sum l15a_7daylasthour if merge_1_2_3 ==3, detail
sum l14a_7dayhour if merge_1_2_3 ==3, detail
tabstat (l15a_7daylasthour) (l14a_7dayhour) if merge_1_2_3==3, stat(mean median sd min max) col(stat)
tab (l11_salary_freq) (settlement_type) if merge_1_2_3 ==3
tab (l01a_temp_absent) (settlement_type) if merge_1_2_3 ==3
*Job Search & Unemployment*
tab (l06_jobsearch) (settlement_type) if merge_1_2_3 ==3
tab (l06b_jobexpect) (settlement_type) if merge_1_2_3 ==3
*Ratio of Distances for spillover*
tab (ls05_sex) (stratum) if merge_1_2_3 == 3
tab (ls25_working) (stratum) if merge_1_2_3 == 3
tab (ls18_school) (stratum) if merge_1_2_3 == 3
tab (ls12_country_born) (stratum) if merge_1_2_3 == 3
*Applying weights*
svyset cluster [pw=aa_weight_r1], strata(stratum) singleunit(scaled)
svyset cluster [pw=aa_weight_r2], strata(stratum) singleunit(scaled)
svyset cluster [pw=aa_weight_r3], strata(stratum) singleunit(scaled)
*weighted descriptive statistics*
svy: tabulate ls25_working, missing
svy: tabulate ls25_working settlement_type
svy: tabulate ls05_sex settlement_type
svy: tab (ls18_school) (settlement_type) if merge_1_2_3 == 3
svy: tab (ls16_read) (settlement_type) if merge_1_2_3 == 3
svy: tab (ls17_write) (settlement_type) if merge_1_2_3 == 3
svy: tab (ls25_working) (settlement_type) if merge_1_2_3 == 3
svy: tab (ls26_available_work) (settlement_type) if merge_1_2_3 == 3
svy: tab (ls27_lookfor_work) (settlement_type) if merge_1_2_3 == 3
svy: tab (job_code_r2) (settlement_type) if merge_1_2_3 == 3
svy: tab (job_code_r3) (settlement_type) if merge_1_2_3 == 3
svy: tab (adult_surveydone) (settlement_type) if merge_1_2_3 == 3
svy: tab upazila settlement_type if merge_1_2_3 == 3
*Weighted measure on spillover*
svy: tab (ls05_sex) (stratum) if merge_1_2_3 == 3
svy: tab (ls25_working) (stratum) if merge_1_2_3 == 3
svy: tab (ls18_school) (stratum) if merge_1_2_3 == 3
svy: tab (ls16_read) (stratum) if merge_1_2_3 == 3
svy: tab (ls17_write) (stratum) if merge_1_2_3 == 3
svy: tab (ls26_available_work) (stratum) if merge_1_2_3 == 3
svy: tab (ls27_lookfor_work) (stratum) if merge_1_2_3 == 3
svy: tab (ls12_country_born) (stratum) if merge_1_2_3 == 3