clear
cap log close
cap clear matrix
set more off
set mem 200m
set matsize 1100
program drop _all
version 11.0

gl root 		"C:\DROPBOX\Dropbox\TESTS_COLOMBIA"

gl testdata 	"$root\DATA\CLEANDATA\v4_tests\idb_portal\tests_bogota.dta"
gl out 			"$root\DATA\CODE\v4_analysis_tests\idb_portal\descriptives"

log using "$out/sample charact battery a.log", replace

use "$testdata"
	keep if battery == 2
	drop if b_flag_outlier == 1
    sum b_tot_cog b_tot_lr b_tot_le b_tot_mf b_tot_mg 

	*final analysis sample 1305 observations

	* Age Categories
	local c = 12
	g age_cat_`c'mth =.
	lab var age_cat_`c'mth "Months Category, `c' months intervals"
	local j = 1
		forvalues i=7(`c')37 {
			local aux = `i'+`c'-1
			replace age_cat_`c'mth = `j' if b_age_mth>=`i' & b_age_mth<=`aux'
			local j = `j'+1
			}
	replace age_cat_`c'mth = 1 if b_age_mth==6 | b_age_mth==5

	tempfile batta
	save `batta', replace

***********Bayley************
u `batta', clear

foreach m in cog lr le mf mg {
		label var b_tot_`m' "All Children"		
}
	
forvalues i=1(1)3 {
foreach m in cog lr le mf mg {
		gen a_b_tot_`m'_`i' = b_tot_`m' if age_cat_`c'mth == `i'
		}
}

local j = 1
foreach cat in " 6 - 18" " 19 - 30" " 31 - 42" {
foreach m in cog lr le mf mg {
		label var a_b_tot_`m'_`j' "`cat'"
		}
	local j = `j'+1
}

foreach m in cog lr le mf mg {
	local r b_tot_`m'
	local a_`m' 	a_`r'_1 a_`r'_2 a_`r'_3 b_tot_`m'
}

estpost sum `a_cog' `a_lr' `a_le' `a_mf' `a_mg', d
est store A

	esttab A using "$out\desc_battery A_`c'm.csv", replace wide plain ///
	title(Table: Raw Score Battery B, All Children and by 6 Months Age Groups) ///
	refcat(a_b_tot_cog_1 "Cognitive" a_b_tot_lr_1 "Language Receptive" a_b_tot_le_1 "Language Expressive" a_b_tot_mf_1 "Fine Motor" ///
	a_b_tot_mg_1 "Gross Motor" , nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

		label var b_tot_cog "Cognitive"		
		label var b_tot_lr "Receptive Language"		
		label var b_tot_le "Expressive Language"		
		label var b_tot_mf "Fine Motor"		
		label var b_tot_mg "Gross Motor"		

estpost sum b_tot_cog b_tot_lr b_tot_le b_tot_mf b_tot_mg, d
est store B

	esttab B using "$out\desc_battery A_all.csv", replace wide plain ///
	title(Table: Raw Score Battery B, All Children) ///
	refcat(b_tot_cog "Bayley-III", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

***********ASQ************
u `batta', clear

drop if asq_flag_outlier == 1
keep if flag_all_b_asq ==1	

foreach m in cog lang mf mg soc tot {
		label var asq_`m' "All Children"		
}
	
forvalues i=1(1)3 {
foreach m in cog lang mf mg soc tot {
		gen a_asq_`m'_`i' = asq_`m' if age_cat_`c'mth == `i'
		}
}

local j = 1
foreach cat in " 6 - 18" " 19 - 30" " 31 - 42" {
foreach m in cog lang mf mg soc tot {
		label var a_asq_`m'_`j' "`cat'"
		}
	local j = `j'+1
}

foreach m in cog lang mf mg soc tot {
	local r asq_`m'
	local a_`m' 	a_`r'_1 a_`r'_2 a_`r'_3 asq_`m'
}

estpost sum `a_cog' `a_lang' `a_mf' `a_mg' `a_soc' `a_tot', d
est store A

	esttab A using "$out\desc_battery A_`c'm.csv", append wide plain ///
	refcat(a_asq_cog_1 "Cognitive" a_asq_lang_1 "Communication" a_asq_mf_1 "Fine Motor" a_asq_mg_1 "Gross Motor" a_asq_soc_1 "Personal-Social" ///
	a_asq_total_1 "Total Score", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs
	
		label var asq_cog "Cognitive"		
		label var asq_lang "Communication"		
		label var asq_mf "Fine Motor"		
		label var asq_mg "Gross Motor"
		label var asq_soc "Personal-Social"		
		label var asq_tot "Total Score"

estpost sum asq_cog asq_lang asq_mf asq_mg asq_soc asq_tot, d
est store B

	esttab B using "$out\desc_battery A_all.csv", append wide plain ///
	refcat(asq_cog "ASQ-3", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

***********Denver************
u `batta', clear

drop if den_flag_outlier == 1
keep if flag_all_b_den ==1	

foreach m in lang mf mg soc tot {
		label var den_`m' "All Children"		
}

forvalues i=1(1)3 {
foreach m in lang mf mg soc tot {
		gen a_den_`m'_`i' = den_`m' if age_cat_`c'mth == `i'
		}
}

local j = 1
foreach cat in " 6 - 18" " 19 - 30" " 31 - 42" {
foreach m in lang mf mg soc tot {
		label var a_den_`m'_`j' "`cat'"
		}
	local j = `j'+1
}

foreach m in lang mf mg soc tot {
	local r den_`m'
	local a_`m' 	a_`r'_1 a_`r'_2 a_`r'_3 den_`m'
}

estpost sum `a_lang' `a_mf' `a_mg' `a_soc' `a_tot', detail
est store A

	esttab A using "$out\desc_battery A_`c'm.csv", append wide plain ///
	refcat(a_den_lang_1 "Language" a_den_mf_1 "Fine Motor" a_den_mg_1 "Gross Motor" a_den_soc_1 "Personal-Social" ///
	a_den_tot_1 "Total Score", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

		label var den_lang "Communication"		
		label var den_mf "Fine Motor"		
		label var den_mg "Gross Motor"
		label var den_soc "Personal-Social"		
		label var den_tot "Total Score"

estpost sum den_lang den_mf den_mg den_soc den_tot, d
est store B

	esttab B using "$out\desc_battery A_all.csv", append wide plain ///
	refcat(den_lang "Denver-II", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

***********Mac-Arthur************
u `batta', clear

drop if mac_flag_outlier == 1
keep if flag_all_b_mac ==1	

foreach m in lr le {
		label var mac_`m' "All Children"		
}


forvalues i=1(1)3 {
	foreach m in lr le {
		gen a_mac_`m'_`i' = mac_`m' if age_cat_`c'mth == `i' 
		}
}

local j = 1
foreach cat in " 8 - 18" " 19 - 30" " 31 - 42" {
	foreach m in lr le {
		label var a_mac_`m'_`j' "`cat'"
		}
	local j = `j'+1
}


foreach m in lr le {
	local r mac_`m'
	local a_`m' 	a_`r'_1 a_`r'_2 a_`r'_3 mac_`m'
}

estpost sum a_mac_lr_1 mac_lr a_mac_le_1 a_mac_le_2 mac_le if b_age_mth>=8 & b_age_mth<=30, detail
est store A

	esttab A using "$out\desc_battery A_`c'm.csv", append wide plain ///
	refcat(a_mac_lr_1 "Receptive Language" a_mac_le_1 "Expressive Language", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs

		label var mac_lr  "Receptive Language"		
		label var mac_le  "Expressive Language"		
		*label var mac_le_12 "Expressive Language"		
		*label var mac_le_34 "Expressive Language"		

estpost sum mac_lr mac_le if b_age_mth>=8 & b_age_mth<=30, d
est store B

	esttab B using "$out\desc_battery A_all.csv", append wide plain ///
	refcat(mac_lr "SFI (MacArthur)" mac_le34 "SFII (MacArthur)", nolabel) ///
	nomtitle collabel("N" "Mean" "Median" "SD" "Min" "Max") ///
	cells ("count(fmt(0)) mean(fmt(3)) p50(fmt(1)) sd(fmt(3)) min(fmt(2)) max(fmt(2))") label lines nonum noobs


log close