*----------------------------------------------------------------------------* *Continuacion del Ingreso laboral *1. Obtenemos los promedios por edad ponderados con los Ns cd "/Users/lrosero/Desktop/ElBono1" use lcd-icv06,clear gen unos=1 collapse (rawsum) unos (sum) factor_f (mean) yle yls yls_raw [w= factor_f],by(edad) ren unos nobs ren factor pob * 2. Suavizamiento *Replicamos para en el suavizamiento tomar en cuenta el N. de observaciones en cada edad sum nobs gen obs10 = round(nobs/10) sum obs10 *Se crean 5473 observaciones replicadas = obs10 expand obs10,gen(duplicada) *tres suavizamientos lowess yle edad if edad>=10,bw(.2) gen(yle_s2) lowess yle edad if edad>=10,bw(.1) gen(yle_s1) lowess yle edad if edad>=10,bw(.07) gen(yle_s07) line yle* edad if dup==0 *El primero parece mejor ajuste ren yle_s2 yle_s drop yle_s07 yle_s1 *valores negativos y missing a cero recode yle_s -200/0 .=0 lowess yls edad if edad>=10,bw(.2) gen(yls_s2) lowess yls edad if edad>=10,bw(.1) gen(yls_s1) lowess yls edad if edad>=10,bw(.06) gen(yls_s06) line yls* edad if dup==0 *El primero parece adecuado ren yls_s2 yls_s drop yls_s1 yls_s06 recode yls_s -200/0 .=0 *Suavizamiento similar en los valores raw lowess yls_raw edad if edad>=10,bw(.2) gen(yls_raws) recode yls_raws -200/0 .=0 *Nos deshacemos de los duplicados drop if dup==1 *etiquetado lab var yle "Ingreso laboral asalariados observado" lab var yls "Ingreso laboral independiente observado" lab var yle_s "Ingreso laboral asalariados suavizado" lab var yls_s "Ingreso laboral independiente suavizado" lab var yls_raws "Ingreso laboral independiente CRUDO suavizado" save tempo,replace *----------------------------------------------------------------------------* * 3 AJUSTE POR CUENTAS NACIONALES *controles macro *YLE 15542 millones de dolares *YLS 9169 gen yle_macro = 15542 gen yls_macro = 9169 *agregados crudos egen double yle_tot = total(pob/1000000 * yle_s) egen double yls_tot = total(pob/1000000 * yls_s) *Factores de correccion gen fyle = yle_macro / yle_tot gen fyls = yls_macro / yls_tot *examinamos los factores de correccion summ f* *Luego usare estos factores para ajustar tambien los micro-datos sort edad save factores,replace *Valores finales gen YLE = fyle * yle_s gen YLS = fyls * yls_s * recien ahora tengo la curva de ingreso laboral gen YL = YLE + YLS *Etiquetas lab var YLE "Labor Earnings" lab var YLS "Self-Employment income" lab var YL "Labor Income" *gr‡ficas *Graficas ilustrativas de suavizamiento y de correccion line yle yle_s edad,legend(off) lw(medium thick) subtit("Suavizamiento del ingreso asalariado Ecuador 2006") ytit("D—lares per c‡pita anual") ysize(5) xsize(6) xla(0(10)90) xtit("A–os de edad") gr export fig-suaviza-YLE2006.pdf,replace line yls_raws yls_s edad, legend(label(1 "Observado") label(2 "Corregido") col(1) pos(1) ring(0) ) lw(medium thick) subtit("Correcci—n del ingreso aut—nomo Ecuador 2006") ytit("D—lares per c‡pita anual") ysize(5) xsize(6) xla(0(10)90) xtit("A–os de edad") gr export fig-corrige-YLS2006.pdf,replace *Grafica de los dos YL como areas *Ingreso laboral twoway area YL YLE edad,legend(label(1 "Aut—nomo") label(2 "Asalariado") col(1) pos(1) ring(0)) subtit("Ingreso laboral per c‡pita Ecuador 2006") ytit("D—lares per c‡pita anual") ysize(5) xsize(6) xla(0(10)90) xtit("A–os de edad") gr export fig-inglab2006.pdf,replace keep edad pob YLE YLS YL sort edad save LCD_Ecua2006,replace ****************************************** * 4 CHEQUEOS DE VALIDEZ *Comparaci—n con otros pa’ses *Para ello hay que normalizar los resultados *Normalizar: todos los pa’ses en una misma unidad monetaria: 1=una unidad de ingreso * = promedio simple del YL de 30 a 49 annos * Calculamos esta unidad de ingreso use LCD_Ecua2006,clear sort edad egen YL3049 = mean(YL) if edad>=30 & edad<=49 replace YL3049 = YL3049[35] if YL3049==. *examinar que c3049 sea una constante summ YL3049 *Cada unidad de ingreso aqui equivale a $3.009 rename YL3049 incunit *calculamos el YL normalizado gen YL_nor = YL/incunit *Queremos comparar con otros paises ren edad age ren pob POP gen year=2006 gen str15 country = "Ecuador" append using LCD-americas-age1 *borramos lo que habia para ecuador drop if codcoun == 218 *asignamos este codigo a lo recien calculado (que esta en missing) recode codcoun . = 218 *Podemos comparar ahora nuestra estimacion de YL *Es mas facil si hacemos un archivo "wide"con solo YL_nor keep age YL_nor codcount rename YL_nor YL reshape wide YL,i(age) j(codcount) rename YL218 YL_Ecuador order age YL_Ecuador *Grafico con la linea gruesa de Ecuador twoway line YL* age if age<90,lwidth(thick) tit("Ingreso laboral en pa’ses de Amˇrica") ytit("Unidades de ingreso*") ylabel(#6) legend(off ) ytit("") ysize(5) xsize(6.5) xla(0(10)90) xtit("Edad") note("* Promedio del ingreso laboral, edades 30-49. Ecuador 2006 es la l’nea gruesa") gr export Amˇrica-YLnor.pdf,replace ****************************************** * 5. es conveniente ajustar los valores micro a los controles: use lcd-icv06,clear OJO METO A MANO LOS DOS FACTORES gen yle_aj = yle * fyle gen yls_aj = yls * fyls drop nobs - cgx_total lab var yle_aj "Ingreso laboral asalariados ajustado a macros" lab var yls_aj "Ingreso laboral independiente ajustado a macros" sort idhog persona save lcd-icv06, replace