# Hoja de Ejercicios 1 - ejercicio 1
open "...\USmacro.gdt"

#(a)
setobs 1 1947 --time-series

setinfo wkearns -d "Ingresos semanales"
setinfo wkhours -d "Horas de Trabajo"
setinfo outphr -d "Productividad" -n "Productividad"
setinfo time -n "tiempo"

#(b) Variable year
genr year = time + 1946
setinfo year -d "Year"

#(c) Productividad a lo largo del tiempo
gnuplot outphr time --output=display

#(d) Salario promedio para cada año
genr sal = wkearns/wkhours
setinfo sal -d "Salario Promedio"
setinfo year --discrete
summary sal --simple --by=year

#(e) Tasas de crecimiento: diferencias logaritmos
genr lsal = log(sal)
genr g_sal = lsal - lsal(-1)
setinfo g_sal -d "Tasa de crecimiento del salario"

#Alternativamente:
diff lsal # first difference of log(sal)
#o:
ldiff sal # first log difference of sal

# Borramos las repetidas
delete d_lsal ld_sal

#Para el producto usamos ldiff
ldiff outphr
rename ld_outphr g_outphr 
setinfo g_outphr -d "Tasa de crecimiento del producto por hora"
gnuplot g_outphr g_sal time --with-lines --time-series --output=display

#(f) Correlaciones

corr sal outphr
corr g_sal g_outphr
gnuplot g_outphr g_sal time --with-lines --time-series --output=display
gnuplot outphr time --output=display

# 1947 - 1980
smpl year<1981 --restrict
corr sal outphr
corr g_sal g_outphr
gnuplot g_outphr g_sal time --with-lines --time-series --output=display
gnuplot outphr time --output=display
smpl --full

# Guardar base
store "USmacro_new.gdt"