version 13 clear all macro drop _all set linesize 80 cap log close set more off *============= PEARSON CORRELATION =========== gl root "C:\DROPBOX\Dropbox\TESTS_COLOMBIA" gl directory "$root\DATA\CLEANDATA\v4_tests\idb_portal\tests_bogota.dta" gl code "$root\DATA\CODE\v4_analysis_tests\idb_portal\pearson\tests_corrs" log using "$code/test significancy level correlations_same function.log", replace gl flag_all flag_all_b_asq //define flag with all data between tests gl flag_col b_flag_outlier //define col_vars outlier *Same function scales (matching domains) *Set up Dataset use "$directory", replace keep if $flag_col == 0 g age_cat = 1 if b_age_mth >=6 & b_age_mth <=14 replace age_cat = 2 if b_age_mth >=15 & b_age_mth <=23 replace age_cat = 3 if b_age_mth >=24 & b_age_mth <=32 replace age_cat = 4 if b_age_mth >=33 & b_age_mth <=42 g estrato_flag =estrato replace estrato_flag =3 if estrato ==4 g strata =. foreach x in 1 2 3 4 { foreach y in 1 2 3 { qui replace strata =`y'`x' if age_cat ==`x' & estrato_flag ==`y' } } tempfile dataset save `dataset', replace ***********************Cognition********************** u `dataset', replace ren r_den_st_mf_kd r_den_st_cog_kd cap program drop corrs_cog program define corrs_cog, rclass foreach x in bat asq den { corr r_b_st_cog_kd r_`x'_st_cog_kd if b_age_mth>=6 & b_age_mth<=18 return scalar corr_b_`x'_cog_618 =r(rho) corr r_b_st_cog_kd r_`x'_st_cog_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_cog_1930 =r(rho) corr r_b_st_cog_kd r_`x'_st_cog_kd if b_age_mth>=31 & b_age_mth<=42 return scalar corr_b_`x'_cog_3142 =r(rho) } end bootstrap /// (r(corr_b_den_cog_618) - r(corr_b_asq_cog_618)) (r(corr_b_den_cog_618) - r(corr_b_bat_cog_618)) (r(corr_b_asq_cog_618) - r(corr_b_bat_cog_618)) /// (r(corr_b_den_cog_1930) - r(corr_b_asq_cog_1930)) (r(corr_b_den_cog_1930) - r(corr_b_bat_cog_1930)) (r(corr_b_asq_cog_1930) - r(corr_b_bat_cog_1930)) /// (r(corr_b_den_cog_3142) - r(corr_b_asq_cog_3142)) (r(corr_b_den_cog_3142) - r(corr_b_bat_cog_3142)) (r(corr_b_asq_cog_3142) - r(corr_b_bat_cog_3142)) /// , reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_cog ***********************Receptive Language********************** u `dataset', replace ren r_asq_st_lang_kd r_asq_st_lr_kd ren r_den_st_lang_kd r_den_st_lr_kd ren r_bat_st_lang_kd r_bat_st_lr_kd cap program drop corrs_lr program define corrs_lr, rclass foreach x in bat asq den { corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=6 & b_age_mth<=18 return scalar corr_b_`x'_lr_618 =r(rho) corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lr_818 =r(rho) corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_lr_1930 =r(rho) corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=31 & b_age_mth<=42 return scalar corr_b_`x'_lr_3142 =r(rho) } foreach x in mac { corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lrr_818 =r(rho) corr r_b_st_le_kd r_`x'_st_lr_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_ler_818 =r(rho) corr r_b_st_lr_kd r_`x'_st_le_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lre_818 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lee_818 =r(rho) corr r_b_st_lr_kd r_`x'_st_le_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_lre_1930 =r(rho) } end bootstrap /// (r(corr_b_mac_lre_818) - r(corr_b_mac_lrr_818)) (r(corr_b_mac_lee_818) - r(corr_b_mac_ler_818)) /// (r(corr_b_den_lr_618) - r(corr_b_mac_lrr_818)) (r(corr_b_den_lr_618) - r(corr_b_mac_lre_818)) /// (r(corr_b_bat_lr_618) - r(corr_b_mac_lrr_818)) (r(corr_b_bat_lr_618) - r(corr_b_mac_lre_818)) /// (r(corr_b_asq_lr_618) - r(corr_b_mac_lrr_818)) (r(corr_b_asq_lr_618) - r(corr_b_mac_lre_818)) /// (r(corr_b_den_lr_618) - r(corr_b_bat_lr_618)) (r(corr_b_den_lr_618) - r(corr_b_asq_lr_618)) (r(corr_b_asq_lr_618) - r(corr_b_bat_lr_618)) /// (r(corr_b_den_lr_1930) - r(corr_b_bat_lr_1930)) (r(corr_b_den_lr_1930) - r(corr_b_asq_lr_1930)) (r(corr_b_asq_lr_1930) - r(corr_b_bat_lr_1930)) /// (r(corr_b_den_lr_3142) - r(corr_b_bat_lr_3142)) (r(corr_b_den_lr_3142) - r(corr_b_asq_lr_3142)) (r(corr_b_asq_lr_3142) - r(corr_b_bat_lr_3142)) /// (r(corr_b_den_lr_1930) - r(corr_b_mac_lre_1930)) (r(corr_b_bat_lr_1930) - r(corr_b_mac_lre_1930)) (r(corr_b_asq_lr_1930) - r(corr_b_mac_lre_1930)) /// (r(corr_b_den_lr_818) - r(corr_b_mac_lrr_818)) (r(corr_b_den_lr_818) - r(corr_b_mac_lre_818)) /// (r(corr_b_bat_lr_818) - r(corr_b_mac_lrr_818)) (r(corr_b_bat_lr_818) - r(corr_b_mac_lre_818)) /// (r(corr_b_asq_lr_818) - r(corr_b_mac_lrr_818)) (r(corr_b_asq_lr_818) - r(corr_b_mac_lre_818)) /// , reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_lr ***********************Expressive Language********************** u `dataset', replace ren r_asq_st_lang_kd r_asq_st_le_kd ren r_den_st_lang_kd r_den_st_le_kd ren r_bat_st_lang_kd r_bat_st_le_kd cap program drop corrs_le program define corrs_le, rclass foreach x in bat asq den { corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=6 & b_age_mth<=18 return scalar corr_b_`x'_le_618 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_le_818 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_le_1930 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=31 & b_age_mth<=42 return scalar corr_b_`x'_le_3142 =r(rho) } foreach x in mac { corr r_b_st_lr_kd r_`x'_st_lr_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lrr_818 =r(rho) corr r_b_st_le_kd r_`x'_st_lr_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_ler_818 =r(rho) corr r_b_st_lr_kd r_`x'_st_le_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lre_818 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=8 & b_age_mth<=18 return scalar corr_b_`x'_lee_818 =r(rho) corr r_b_st_le_kd r_`x'_st_le_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_lee_1930 =r(rho) } end bootstrap /// (r(corr_b_den_le_618) - r(corr_b_mac_ler_818)) (r(corr_b_den_le_618) - r(corr_b_mac_lee_818)) /// (r(corr_b_bat_le_618) - r(corr_b_mac_ler_818)) (r(corr_b_bat_le_618) - r(corr_b_mac_lee_818)) /// (r(corr_b_asq_le_618) - r(corr_b_mac_ler_818)) (r(corr_b_asq_le_618) - r(corr_b_mac_lee_818)) /// (r(corr_b_den_le_618) - r(corr_b_bat_le_618)) (r(corr_b_den_le_618) - r(corr_b_asq_le_618)) (r(corr_b_asq_le_618) - r(corr_b_bat_le_618)) /// (r(corr_b_den_le_1930) - r(corr_b_bat_le_1930)) (r(corr_b_den_le_1930) - r(corr_b_asq_le_1930)) (r(corr_b_asq_le_1930) - r(corr_b_bat_le_1930)) /// (r(corr_b_den_le_3142) - r(corr_b_bat_le_3142)) (r(corr_b_den_le_3142) - r(corr_b_asq_le_3142)) (r(corr_b_asq_le_3142) - r(corr_b_bat_le_3142)) /// (r(corr_b_den_le_1930) - r(corr_b_mac_lee_1930)) (r(corr_b_bat_le_1930) - r(corr_b_mac_lee_1930)) (r(corr_b_asq_le_1930) - r(corr_b_mac_lee_1930)) /// (r(corr_b_den_le_818) - r(corr_b_mac_ler_818)) (r(corr_b_den_le_818) - r(corr_b_mac_lee_818)) /// (r(corr_b_bat_le_818) - r(corr_b_mac_ler_818)) (r(corr_b_bat_le_818) - r(corr_b_mac_lee_818)) /// (r(corr_b_asq_le_818) - r(corr_b_mac_ler_818)) (r(corr_b_asq_le_818) - r(corr_b_mac_lee_818)) /// , reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_le ***********************Fine Motor********************** u `dataset', replace ren r_bat_st_mot_kd r_bat_st_mf_kd cap program drop corrs_mf program define corrs_mf, rclass foreach x in bat asq den { corr r_b_st_mf_kd r_`x'_st_mf_kd if b_age_mth>=6 & b_age_mth<=18 return scalar corr_b_`x'_mf_618 =r(rho) corr r_b_st_mf_kd r_`x'_st_mf_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_mf_1930 =r(rho) corr r_b_st_mf_kd r_`x'_st_mf_kd if b_age_mth>=31 & b_age_mth<=42 return scalar corr_b_`x'_mf_3142 =r(rho) } end bootstrap /// (r(corr_b_den_mf_618) - r(corr_b_asq_mf_618)) (r(corr_b_den_mf_618) - r(corr_b_bat_mf_618)) (r(corr_b_asq_mf_618) - r(corr_b_bat_mf_618)) /// (r(corr_b_den_mf_1930) - r(corr_b_asq_mf_1930)) (r(corr_b_den_mf_1930) - r(corr_b_bat_mf_1930)) (r(corr_b_asq_mf_1930) - r(corr_b_bat_mf_1930)) /// (r(corr_b_den_mf_3142) - r(corr_b_asq_mf_3142)) (r(corr_b_den_mf_3142) - r(corr_b_bat_mf_3142)) (r(corr_b_asq_mf_3142) - r(corr_b_bat_mf_3142)) /// , reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_mf ***********************Gross Motor********************** u `dataset', replace ren r_bat_st_mot_kd r_bat_st_mg_kd ren r_mil_st_tot_kd r_mil_st_mg_kd cap program drop corrs_mg program define corrs_mg, rclass foreach x in bat asq den { corr r_b_st_mg_kd r_`x'_st_mg_kd if b_age_mth>=6 & b_age_mth<=18 return scalar corr_b_`x'_mg_618 =r(rho) corr r_b_st_mg_kd r_`x'_st_mg_kd if b_age_mth>=6 & b_age_mth<=15 return scalar corr_b_`x'_mg_615 =r(rho) corr r_b_st_mg_kd r_`x'_st_mg_kd if b_age_mth>=19 & b_age_mth<=30 return scalar corr_b_`x'_mg_1930 =r(rho) corr r_b_st_mg_kd r_`x'_st_mg_kd if b_age_mth>=31 & b_age_mth<=42 return scalar corr_b_`x'_mg_3142 =r(rho) } foreach x in mil { corr r_b_st_mg_kd r_`x'_st_mg_kd if b_age_mth>=6 & b_age_mth<=15 return scalar corr_b_`x'_mg_615 =r(rho) } end bootstrap /// (r(corr_b_mil_mg_615) - r(corr_b_den_mg_618)) (r(corr_b_mil_mg_615) - r(corr_b_bat_mg_618)) (r(corr_b_mil_mg_615) - r(corr_b_asq_mg_618)) /// (r(corr_b_mil_mg_615) - r(corr_b_den_mg_615)) (r(corr_b_mil_mg_615) - r(corr_b_bat_mg_615)) (r(corr_b_mil_mg_615) - r(corr_b_asq_mg_615)) /// (r(corr_b_den_mg_618) - r(corr_b_asq_mg_618)) (r(corr_b_den_mg_618) - r(corr_b_bat_mg_618)) (r(corr_b_asq_mg_618) - r(corr_b_bat_mg_618)) /// (r(corr_b_den_mg_1930) - r(corr_b_asq_mg_1930)) (r(corr_b_den_mg_1930) - r(corr_b_bat_mg_1930)) (r(corr_b_asq_mg_1930) - r(corr_b_bat_mg_1930)) /// (r(corr_b_den_mg_3142) - r(corr_b_asq_mg_3142)) (r(corr_b_den_mg_3142) - r(corr_b_bat_mg_3142)) (r(corr_b_asq_mg_3142) - r(corr_b_bat_mg_3142)) /// , reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_mg log close