----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  C:\Users\chouse\Dropbox\basu-house\analysis\code\../temp/entry2quart_ipd_cumtight.log
  log type:  text
 opened on:  18 Jan 2017, 15:01:30

. /**************************************************************************************/
> 
> /*** NSLY DATA ON USER COST ***/
> /* Do the MATLAB interpolattion */
> 
> capture confirm file `myinput'entry_quarterly_interpolated_ipd.csv;

. if _rc == 0 {;
.     di "File existed and is now gone.";
File existed and is now gone.
.     erase `myinput'entry_quarterly_interpolated_ipd.csv;
. };

. if _rc != 0 {;
.     di "File did not exist.";
. };

. !matlab -r run('./entry2quart_ipd_cumtight.m');

. sleep 10000;

.                              /* sleep 1,000 = 1 second  */
> capture confirm file `myinput'entry_quarterly_interpolated_ipd_cumtight.csv;

. while _rc {;
  2.                             /* file does not exist! */
>            sleep 1000;
  3.            capture confirm file `myinput'entry_quarterly_interpolated_ipd_cumtight.csv;
  4.            };

.            /* Get output from UC2QUART.M */
> import delimited using `myinput'entry_quarterly_interpolated_ipd_cumtight, clear;
(3 vars, 104 obs)

. /* merge in CHOW data to extrapolate */
> gen dateq = yq(year,quarter);

. format dateq %tq;

. merge 1:1 dateq using `myinput'fred4uc2quart_data;
(note: variable year was int, now float to accommodate using data's values)
(note: variable quarter was byte, now float to accommodate using data's values)

    Result                           # of obs.
    -----------------------------------------
    not matched                           173
        from master                         0  (_merge==1)
        from using                        173  (_merge==2)

    matched                               104  (_merge==3)
    -----------------------------------------

. gen lrgdp = log(rgdp);

. gen lempl = log(empl);

. gen lcomp = log(comp);

. reg lentry_wageq ur lrgdp lempl lcomp;

      Source |       SS           df       MS      Number of obs   =       104
-------------+----------------------------------   F(4, 99)        =   1776.00
       Model |    41.30417         4  10.3260425   Prob > F        =    0.0000
    Residual |  .575605829        99    .0058142   R-squared       =    0.9863
-------------+----------------------------------   Adj R-squared   =    0.9857
       Total |  41.8797758       103  .406599765   Root MSE        =    .07625

------------------------------------------------------------------------------
lentry_wageq |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          ur |   .0181153   .0079127     2.29   0.024     .0024148    .0338158
       lrgdp |   3.115966   .4006961     7.78   0.000     2.320898    3.911034
       lempl |  -2.832952     .57532    -4.92   0.000    -3.974512   -1.691392
       lcomp |   1.293074   .4666182     2.77   0.007     .3672026    2.218946
       _cons |   1.230506   4.388836     0.28   0.780    -7.477898     9.93891
------------------------------------------------------------------------------

. predict lentry_wage;
(option xb assumed; fitted values)
(4 missing values generated)

. replace lentry_wage = lentry_wageq if !missing(lentry_wageq);
(104 real changes made)

. tsset dateq;
        time variable:  dateq, 1947q1 to 2016q1
                delta:  1 quarter

. tsline lentry_wage, lc(black) lw(1.1) xtitle("") ytitle("Log(Real Dollars)")
> ylabel(,angle(h));

.  rename lentry_wage lentry_wage_ipd_cumtight;

. la var lentry_wage_ipd_cumtight "Log(entry wage) (IPD deflated, cumtight)";

. rename lentry_wageq lentry_wage_ipd_interp_cumtight;

. /*** merge all the data together and output ***/
> /* original data used by the MATLAB program */
> merge m:1 year using `myinput'fixed_effects_ipd_cumtight_entry_wage, nogen;

    Result                           # of obs.
    -----------------------------------------
    not matched                           173
        from master                       173  
        from using                          0  

    matched                               104  
    -----------------------------------------

. rename lentry_wage lentry_wage_ipd_fes_cumtight;

. tsline lentry_wage_ipd_cumtight lentry_wage_ipd_interp_cumtight lentry_wage_ipd_fes_cumtight if tin(1978q1,2012q4), lw(0.8 0.8 0.8) lc(black red mint)
> legend(c(3) label(1 "Final Series") label(2 "Interpolated only") label(3 "Time FEs") region(lstyle(none)));

. /* export excel year quarter lentry_wage_ipd_cumtight lentry_wage_ipd_interp_cumtight lentry_wage_ipd_fes_cumtight using `myoutput'lentry_quarterly_ipd_all, firstrow(variables)
>  replace; */
> 
> /*** output more "final" dataset ***/
> keep year quarter dateq lentry_wage_ipd_cumtight;

. save `myinput'lentry_quarterly_ipd_cumtight, replace;
file ../input/lentry_quarterly_ipd_cumtight.dta saved

. 
end of do-file

. do user_cost_var.do;

. /**************************************************************
> Program: user_cost_var.do
> Purpose: CEE VAR with special treatment given to user cost.
> 
> Date Started:  7 September 2015
> Date Revised: 12 July      2016
> **************************************************************/
. #delimit ;
delimiter now ;
. clear all;

. set mo off;

. set scheme s1mono;

. capture log close;