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_denver_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_den_st_lang_kd r_den_st_mf_kd r_den_st_mg_kd r_den_st_soc_kd r_den_st_tot_kd		

gl flag_all flag_all_b_den																			//define flag with all data between tests
gl flag_col b_flag_outlier																			//define col_vars outlier
gl flag_row den_flag_outlier																		//define row_vars outlier

gl col_test Bayley-III																					//Title test col_vars
gl row_test Denver-II

********Definir dataset de trabajo*******
	use "$directory", replace 																//open dataset
	keep if $flag_all == 1																	//keep useful data, deleting outliers and missed
	*keep if $flag_row == 0
	keep if $flag_col == 0
	keep if battery == 2		//replace with corresponding battery (1=batB, 2=batA)
	
	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 Denver: 6-18 months
	u `dataset', clear 																	//open dataset

	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 den_`y' =(`y'_r_den_st_lang_kd \ `y'_r_den_st_mf_kd \ `y'_r_den_st_mg_kd \ `y'_r_den_st_soc_kd \ `y'_r_den_st_tot_kd)
	}

	mat denver_618 =(den_r_b_st_cog_kd \ den_r_b_st_lr_kd \ den_r_b_st_le_kd \ den_r_b_st_mf_kd \ den_r_b_st_mg_kd)
	
	di "****************6-18 months number of observations**************
	di e(N)

//****************Bayley vs Denver: 19-30 months
	u `dataset', clear 																	//open dataset

	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 den_`y' =(`y'_r_den_st_lang_kd \ `y'_r_den_st_mf_kd \ `y'_r_den_st_mg_kd \ `y'_r_den_st_soc_kd \ `y'_r_den_st_tot_kd)
	}

	mat denver_1930 =(den_r_b_st_cog_kd \ den_r_b_st_lr_kd \ den_r_b_st_le_kd \ den_r_b_st_mf_kd \ den_r_b_st_mg_kd)
	
	di "****************19-30 months number of observations**************
	di e(N)

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

	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 den_`y' =(`y'_r_den_st_lang_kd \ `y'_r_den_st_mf_kd \ `y'_r_den_st_mg_kd \ `y'_r_den_st_soc_kd \ `y'_r_den_st_tot_kd)
	}

	mat denver_3142 =(den_r_b_st_cog_kd \ den_r_b_st_lr_kd \ den_r_b_st_le_kd \ den_r_b_st_mf_kd \ den_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 denver_`x'
	keep denver_`x'1-denver_`x'6
	drop if _n >=40
	
	ren denver_`x'1 b
	ren denver_`x'2 se
	ren denver_`x'3 z
	ren denver_`x'4 pvalue
	ren denver_`x'5 ll
	ren denver_`x'6 ul
	
	save "$out\denver_`x'.dta", replace
}