freduse MEHOINUSMAA672N LES1252881600Q CUURA103SA0 if daten>=mdy(1,1,2000), clear
gen m = mofd(daten)
tsset m
tssmooth ma cpi_boston = CUURA103SA0, window(1 1 1)
gen t = qofd(daten)
collapse (mean) MEHOINUSMAA672N LES1252881600Q cpi_boston, by(t)
format %tq t
tsset t
gen delta_muwre = log(LES1252881600Q) - log(L4.LES1252881600Q)
gen delta_boston_income = log(MEHOINUSMAA672N) - log(L4.MEHOINUSMAA672N)
tssmooth ma delta_boston_income2 = delta_boston_income, window(3 1 3)
gen delta_cpi_boston = log(cpi_boston) - log(L4.cpi_boston)
newey delta_boston_income2 delta_muwre if t>=yq(2010,1), lag(3)
predict fitted_boston_income
quietly summarize fitted_boston_income if yofd(dofq(t))==2018
di `r(mean)'