* II. EDUCATION************************************************** use 02_enighur11_personas_ingresos,clear gen edad = p03 gen edad5 =int(edad/5)*5 ****p15a: nivel de enseñanza al que asiste Kinder, primaria, etc)**** ****p17: tipo de centro al que asiste (privado, publico....)**** * Obviaremos que unos pocos matriculados (<1%) no asisten tab p15a p17,miss gen student=p17 recode student 1/4=1 .=0 lab var student "Es estudiante" gen public=p17 recode public 1 4=1 2 3=0 *tres niveles de educacion: 1=basica, 2=secundaria y 3=terciaria gen levedu=p15a recode levedu 2/5=1 6 7=2 8/10=3 1=. tab p15a levedu,miss tab levedu publ,miss gen edu=levedu + 3*(public) *Hacemos dummies para cada nivel-tipo (seis dummies tab edu,gen(ed) lab var ed1 "Basic-private" lab var ed2 "Secondary-private" lab var ed3 "Tertiary-private" lab var ed4 "Basic-public" lab var ed5 "Secondary-public" lab var ed6 "Tertiary-public" ****Find the number of enrolled students of each edad by household**** gen sedad=edad if student==1 recode sedad 29/90=29 tab sedad,gen(n) *Notar que creamos 29 dummies n1.... * n1 1 si es estudiante de edad 5 etc... ****Colapsar información de personas a hogar**** sort identif_hog collapse (sum) student public ed1 - ed6 n1 - n25,by(identif_hog) sort identif_hog *notar que las dummies se convirtieron numero de personas en el hogar que asisten a escuela en esa edad save educhh,replace *este es un archivo temporal con info de educ en cada hogar *Cargo el archivo con info de hogares de la encuesta use 10_enighur11_hogares_agregados, clear sort identif_hog merge identif_hog using educhh tab _mer drop _merge ****d10: gasto mensual en educación**** *probaremos dos metodos. * el metodo 2 es el sugerido por el manual * en este el gastodepende de la edad solamente ****Método 1 (propuesto por nosotros): El gasto depende del nivel educativo**** * Sigue el metodo 1 * d10 es el gasto en educaci—n del hogar recode d10 .=0 *La regresion sin constante reg d10 ed1 ed2 ed3 ed4 ed5 ed6 [w=fexp],robust noconstant *notar que R2=0.41 *los coeficientes de regresion matrix bes = get(_b) svmat bes, name(be) predict cfepred1 ***********************2**************************** *Método 2 (sugerido por el manual): Educación depende de la edad reg d10 n1 -n25 [w=fexp],robust noconstant *notar que R2=0.17: predice mejor la regresion previa matrix bns = get(_b) svmat bns, name(bn) predict cfepred2 ************************2*********************************** summ cfepred1 cfepred2 d10 replace cfepred1=0 if cfepred1<0 ****Assign coefficients to all enrolled by edad**** sort identif_hog keep identif_hog d10 be1 - cfepred2 fexp order identif_hog d10 cfepred1 cfepred2 fexp sort identif_hog count save educhh,replace *dejamos solo los coeficientes de las regresiones use educhh,clear keep in 1 drop identif_hog - fexp gen id=1 reshape long be,i(id) j(edu) drop id reshape long bn,i(edu) j(sedad) replace sedad=sedad+4 *para no tener negativos sumamos 2 a be replace be = be+2 sort edu sedad save educoef,replace ***Volvemos a file de personas*** use 02_enighur11_personas_ingresos,clear gen edad = p03 gen student=p17 recode student 1/4=1 .=0 gen public=p17 recode public 1 4=0 2 3=1 gen levedu=p15a recode levedu 2/5=1 6 7=2 8/10=3 1=. gen edu=levedu + 3*(public) gen sedad=edad if student==1 recode sedad 29/90=29 sort edu sedad merge edu sedad using educoef,nokeep tab _mer drop _merge sort identif_hog merge identif_hog using educhh keep identif_hog identif_per12 fexp_cen2010 edad d10 edu cfepred1 cfepred2 sedad edu be bn ****Total Estimated average household expenditures on education**** sort identif_hog by identif_hog: egen pred1=sum(be) by identif_hog: egen pred2=sum(bn) *Proportion of estimated household expenditures on each enrolled individual *times the household educational expenditure* *regression coeffienct divided by TOTAL household expenditure times the household expenditure on education* gen cfei1= be*d10/pred1 gen cfei2= bn*d10/pred2 recode cfei1 .=0 recode cfei2 .=0 graph bar cfei1 cfei2 if edad<40,over(edad) *dejamos solamente la estimacion por el primer metodo keep identif_hog identif_per12 edad cfei1 ren cfei1 cfei *anualizado replace cfei = cfei *12 sort identif_per12 save cfei,replace