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

. /*** NSLY DATA ON USER COST ***/
> /* Do the MATLAB interpolattion */
> /* !del \\lsa-research03.m.storage.umich.edu\lsa-research03\richryan\nlsy\analysis\input\entry_quarterly_interpolated_ipd.csv; */
> 
> capture confirm file `myinput'entry_quarterly_interpolated_ipd.csv;

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

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

. /* !matlab -r run('\\lsa-research03.m.storage.umich.edu\lsa-research03\richryan\nlsy\analysis\code\entry2quart_ipd.m'); */
> !matlab -r run('./entry2quart_ipd.m');

. sleep 10000;

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

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

.            /* Get output from UC2QUART.M */
> import delimited using `myinput'entry_quarterly_interpolated_ipd, 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)        =   1834.65
       Model |  40.3858913         4  10.0964728   Prob > F        =    0.0000
    Residual |  .544818279        99  .005503215   R-squared       =    0.9867
-------------+----------------------------------   Adj R-squared   =    0.9862
       Total |  40.9307095       103  .397385529   Root MSE        =    .07418

------------------------------------------------------------------------------
lentry_wageq |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          ur |   .0140256   .0076982     1.82   0.071    -.0012493    .0293004
       lrgdp |   3.337793   .3898328     8.56   0.000      2.56428    4.111306
       lempl |  -3.203123   .5597224    -5.72   0.000    -4.313734   -2.092513
       lcomp |   1.116343   .4539676     2.46   0.016     .2155728    2.017113
       _cons |   4.316892    4.26985     1.01   0.314    -4.155417     12.7892
------------------------------------------------------------------------------

. 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;

. la var lentry_wage_ipd "Log(entry wage) (IPD deflated)";

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

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

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

. tsline lentry_wage_ipd lentry_wageq lentry_wage 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)));

. rename lentry_wageq lentry_wage_ipd_interp;

. rename lentry_wage lentry_wage_ipd_fes;

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

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

. 
end of do-file

. do entry2quart_ipd_cumtight.do;

. /********************************************************************
> Program: entry2quart_ipd_cumtight.do
> Purpose: Take annual measure of entry wages and interpolate and
> extrapolate to get quarterly series of entry wages.
> 
> Date Started: 29 December  2015
> Date Revised: 17 July      2016
> ********************************************************************/
. #delimit ;
delimiter now ;
. #delimit ;
delimiter now ;
. clear all;

. set mo off;

. set scheme s1mono;

. capture log close;