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_cross 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
*gl flag_row asq_flag_outlier																		//define row_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 			

cap program drop corrs_cross
program define corrs_cross, rclass
	
		*Denver 6-18 mths
		corr r_b_st_cog_kd r_den_st_mf_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_cog_cog_618 =r(rho)
		
		corr r_b_st_cog_kd r_den_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_cog_soc_618 =r(rho)

		corr r_b_st_lr_kd r_den_st_lang_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_lr_lang_618 =r(rho)
		
		corr r_b_st_lr_kd r_den_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_lr_soc_618 =r(rho)

		corr r_b_st_le_kd r_den_st_lang_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_le_lang_618 =r(rho)
		
		corr r_b_st_le_kd r_den_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_den_le_soc_618 =r(rho)

		*ASQ 6-18 mths
		corr r_b_st_cog_kd r_asq_st_cog_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_cog_cog_618 =r(rho)
		
		corr r_b_st_cog_kd r_asq_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_cog_soc_618 =r(rho)

		corr r_b_st_lr_kd r_asq_st_lang_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_lr_lang_618 =r(rho)
		
		corr r_b_st_lr_kd r_asq_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_lr_soc_618 =r(rho)

		corr r_b_st_le_kd r_asq_st_lang_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_le_lang_618 =r(rho)
		
		corr r_b_st_le_kd r_asq_st_soc_kd if b_age_mth>=6 & b_age_mth<=18
		return scalar corr_b_asq_le_soc_618 =r(rho)
		
		*WHO-Motor vs. Denver 6-18 mths
		corr r_b_st_cog_kd r_den_st_mf_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_den_cog_cog_615 =r(rho)
		
		corr r_b_st_cog_kd r_mil_st_tot_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_mil_cog_mg_615 =r(rho)

		corr r_b_st_lr_kd r_den_st_lang_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_den_lr_lang_615 =r(rho)
		
		corr r_b_st_lr_kd r_mil_st_tot_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_mil_lr_mg_615 =r(rho)
		
		corr r_b_st_le_kd r_den_st_lang_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_den_le_lang_615 =r(rho)
		
		corr r_b_st_le_kd r_mil_st_tot_kd if b_age_mth>=6 & b_age_mth<=15
		return scalar corr_b_mil_le_mg_615 =r(rho)

		*Language vs. Cognition 19 and older
		corr r_b_st_cog_kd r_bat_st_cog_kd if b_age_mth>=19 & b_age_mth<=30
		return scalar corr_b_bat_cog_cog_1930 =r(rho)
		
		corr r_b_st_cog_kd r_bat_st_lang_kd if b_age_mth>=19 & b_age_mth<=30
		return scalar corr_b_bat_cog_lang_1930 =r(rho)
		
		corr r_b_st_cog_kd r_den_st_mf_kd if b_age_mth>=31 & b_age_mth<=42
		return scalar corr_b_den_cog_cog_3142 =r(rho)
		
		corr r_b_st_cog_kd r_den_st_lang_kd if b_age_mth>=31 & b_age_mth<=42
		return scalar corr_b_den_cog_lang_3142 =r(rho)
end

bootstrap ///
(r(corr_b_den_cog_cog_618) - r(corr_b_den_cog_soc_618)) (r(corr_b_den_lr_lang_618) - r(corr_b_den_lr_soc_618)) (r(corr_b_den_le_lang_618) - r(corr_b_den_le_soc_618)) ///
(r(corr_b_asq_cog_cog_618) - r(corr_b_asq_cog_soc_618)) (r(corr_b_asq_lr_lang_618) - r(corr_b_asq_lr_soc_618)) (r(corr_b_asq_le_lang_618) - r(corr_b_asq_le_soc_618)) ///
(r(corr_b_den_cog_cog_615) - r(corr_b_mil_cog_mg_615))  (r(corr_b_den_lr_lang_615) - r(corr_b_mil_lr_mg_615)) (r(corr_b_den_le_lang_615) - r(corr_b_mil_le_mg_615)) ///
(r(corr_b_bat_cog_cog_1930) - r(corr_b_bat_cog_lang_1930)) (r(corr_b_den_cog_cog_3142) - r(corr_b_den_cog_lang_3142)) ///
, reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corrs_cross


log close