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\bs_se"
gl out 			"$root\DATA\CODE\v4_analysis_tests\idb_portal\pearson\bs_se"
gl core 		"$root\DATA\CODE\v4_analysis_tests\idb_portal\core"

log using "$code/bayley_bayley_pearson_12mths.log", replace 

gl col_vars r_b_st_cog_kd r_b_st_lr_kd r_b_st_le_kd r_b_st_mf_kd r_b_st_mg_kd  		// variables Bayley
gl row_vars r_b_st_cog_kd r_b_st_lr_kd r_b_st_le_kd r_b_st_mf_kd r_b_st_mg_kd  		// variables Bayley


********Definir dataset de trabajo*******
	use "$directory", replace 																//open dataset
	
	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

	
//****************Bayley vs Bayley: 6-18 months
	u `dataset', clear 																	

	foreach y of glo col_vars {
	foreach x of glo row_vars { 
		bootstrap corr =r(rho), reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corr `y' `x' if b_age_mth>=6 & b_age_mth<=18
		mat `y'_`x' =r(table)'
		}
	}
	
	foreach y of glo col_vars {
		mat b_`y' =(`y'_r_b_st_cog_kd \ `y'_r_b_st_lr_kd \ `y'_r_b_st_le_kd \ `y'_r_b_st_mf_kd \ `y'_r_b_st_mg_kd )
	}
	
	mat b_618 =(b_r_b_st_cog_kd \ b_r_b_st_lr_kd \ b_r_b_st_le_kd \ b_r_b_st_mf_kd \ b_r_b_st_mg_kd)
	
	di "****************6-18 months number of observations**************
	di e(N)

//****************Bayley vs Bayley: 19-30 months
	u `dataset', clear 		
	
	foreach y of glo col_vars {
	foreach x of glo row_vars { 
		bootstrap corr =r(rho), reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corr `y' `x' if b_age_mth>=19 & b_age_mth<=30
		mat `y'_`x' =r(table)'
		}
	}
	
	foreach y of glo col_vars {
		mat b_`y' =(`y'_r_b_st_cog_kd \ `y'_r_b_st_lr_kd \ `y'_r_b_st_le_kd \ `y'_r_b_st_mf_kd \ `y'_r_b_st_mg_kd)
	}
	
	mat b_1930 =(b_r_b_st_cog_kd \ b_r_b_st_lr_kd \ b_r_b_st_le_kd \ b_r_b_st_mf_kd \ b_r_b_st_mg_kd)
	
	di "****************19-30 months number of observations**************
	di e(N)

	
//****************Bayley vs Bayley: 31-42 months
	u `dataset', clear 																	

	foreach y of glo col_vars {
	foreach x of glo row_vars { 
		bootstrap corr =r(rho), reps(1000) seed(1978) level(95) nodots nowarn strata(strata): corr `y' `x' if b_age_mth>=31 & b_age_mth<=42
		mat `y'_`x' =r(table)'
		}
	}
	
	foreach y of glo col_vars {
		mat b_`y' =(`y'_r_b_st_cog_kd \ `y'_r_b_st_lr_kd \ `y'_r_b_st_le_kd \ `y'_r_b_st_mf_kd \ `y'_r_b_st_mg_kd)
	}
	
	mat b_3142 =(b_r_b_st_cog_kd \ b_r_b_st_lr_kd \ b_r_b_st_le_kd \ b_r_b_st_mf_kd \ b_r_b_st_mg_kd)

	di "****************31-42 months number of observations**************
	di e(N)


//*****Convert Matrices to Datasets
foreach x in 618 1930 3142 {
	svmat b_`x'
	keep b_`x'1-b_`x'6
	drop if _n >=40
	
	ren b_`x'1 b
	ren b_`x'2 se
	ren b_`x'3 z
	ren b_`x'4 pvalue
	ren b_`x'5 ll
	ren b_`x'6 ul
	
	save "$out\bayley_`x'.dta", replace
}