version 12

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

gl namefile 	"$out\znonpar_bayley_battelle_12mths.xlsx"		

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

gl excelsheet 	"Bayley-Battelle-12Mths"																	// sheet excel donde se guardan tablas

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_bat_st_cog_kd r_bat_st_lang_kd r_bat_st_mot_kd r_bat_st_soc_kd r_bat_st_adp_kd r_bat_st_tot_kd			

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

gl col_test Bayley-III																					//Title test col_vars
gl row_test BSID-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 == 1																//replace with corresponding battery (1=batB, 2=batA)
	tempfile dataset
	save `dataset', replace

//******Bayley vs Battelle: 6-42 months
	u `dataset', clear 																	//open dataset
	keep if b_age_mth>=6
	
	do "$core\pvaluesmat.do"
	
	correlate $row_vars $col_vars

	gl typetest "Pearson Correlation $col_test VS $row_test 6 - 42 months"				// define title of table

	matrix A = r(C)																			// matrix with coef. matrix
	matrix B = P																			// matrix with p-values

	do "$core\matstring.do"

	export excel using "$namefile", sheet("$excelsheet") firstrow(varlabels) cell(A1) sheetreplace	// export matrix stored as dataset in Excel

	local numtest = wordcount("$row_vars")+5								// count variables to put the next matrix bellow in the same sheet
																								

*Repeat for each age category by adding age conditions (theses matrices are appended bellow the last one)		
		
//****************Bayley vs Battelle: 6-18 months
	u `dataset', clear 																	//open dataset
	keep if b_age_mth>=6 & b_age_mth<=18 
	
	do "$core\pvaluesmat.do"															// generate matrix with p-values
	
	correlate $row_vars $col_vars

	gl typetest "Pearson Correlation $col_test VS $row_test 6 - 18 months"				// define title of table

	matrix A = r(C)																		// matrix with coef. matrix
	matrix B = P																		// matrix with p-values (got from pvaluesmat.do)

	do "$core\matstring.do"																

	export excel using "$namefile", sheet("$excelsheet") firstrow(varlabels) cell(A`numtest') sheetmodify	

	local numtest = `numtest'+wordcount("$row_vars")+4															// actualizo numtest para pegar el siguiente	

	
//***************Bayley vs Battelle: 19-30 months
	u `dataset', clear 																	//open dataset
	keep if b_age_mth>=19 & b_age_mth<=30 
	
	do "$core\pvaluesmat.do"
	
	correlate $row_vars $col_vars

	gl typetest "Pearson Correlation $col_test VS $row_test 19 - 30 months"				// define title of table

	matrix A = r(C)																			// matrix with coef. matrix
	matrix B = P																			// matrix with p-values

	do "$core\matstring.do"

	export excel using "$namefile", sheet("$excelsheet") firstrow(varlabels) cell(A`numtest') sheetmodify	

	local numtest = `numtest'+wordcount("$row_vars")+4													// actualizo numtest para pegar el siguiente

//***************Bayley vs Battelle: 31-42 months
	u `dataset', clear 																	//open dataset
	keep if b_age_mth>=31 & b_age_mth<=42 
	
	do "$core\pvaluesmat.do"
	
	correlate $row_vars $col_vars

	gl typetest "Pearson Correlation $col_test VS $row_test 31 - 42 months"				// define title of table

	matrix A = r(C)																			// matrix with coef. matrix
	matrix B = P																			// matrix with p-values

	do "$core\matstring.do"

	export excel using "$namefile", sheet("$excelsheet") firstrow(varlabels) cell(A`numtest') sheetmodify	

	local numtest = `numtest'+wordcount("$row_vars")+4													// actualizo numtest para pegar el siguiente