open output ts2.out
**************************************************************************
env batch
cal 54 1 4
all 0 93:4
**************************************************************************
** set location of data sets here
open data c:\allq.rat
data(format=rats) / emp2 ur2 upr2 lfp2 pall p1617 p1619 p2024 afp1619 $
                    kaitz kaitzw enr1619 enr1619b empm2554 rwage union $
                    minr cpi le6hm coverw coveru
close data

open data c:\sexrace.rat 
data(format=rats) / empm empf empw em emw p1617m p1617f p1617w $
                    p1619m p1619f p1619w p1617b p1619b empb
close data

**************************************************************************
******create base variables
**************************************************************************
set empnw = (em-emw)/(p1619-p1619w)
set empm = .001*empm
set empf = .001*empf
set empw = .001*empw
set empb = .001*empb

set p1617nw = p1617-p1617w
set p1619nw = p1619-p1619w

set synw = p1617nw/p1619nw
set syw = p1617w/p1619w
set sym = p1617m/p1619m
set syf = p1617f/p1619f
set syb = p1617b/p1619b
set sy = p1617/p1619

set empam = .001*empm2554
set pop = p1619/pall
set popa = p2024/pall
set afp = afp1619/p1619
set emp = .001*emp2
set upr = .001*upr2

set lsym = log(sym)
set lsyf = log(syf)
set lsyw = log(syw)
set lsynw = log(synw)
set lsyb = log(syb)
set lempm = log(empm)
set lempf = log(empf)
set lempw = log(empw)
set lempnw = log(empnw)
set lempb = log(empb)

******combine AW's variables with ours
set enp 54:1 79:4 = enr1619b
set enp 80:1 93:4 = enr1619
set yk 54:1 86:4 = kaitzw
set yk 87:1 93:4 = kaitz
set cover 54:1 86:4 = coverw
set cover 87:1 93:4 = coveru

******create other variables
set trend = t
set trendsq = t**2
seasonal seasons
set lyk = log(yk)
set lemp = log(emp)
set lempam = log(empam)
set lupr = log(upr)
set lpop = log(pop)
set lpopa = log(popa)
set lminr = log(minr)
set tq2 = trend*seasons{-2}
set tq3 = trend*seasons{-1}
set tq4 = trend*seasons
set t2q2 = trendsq*seasons{-2}
set t2q3 = trendsq*seasons{-1}
set t2q4 = trendsq*seasons
set lyksq = lyk**2

set lcover = log(cover)
set wmanr = le6hm/(cpi*.01)
set lwmanr = log(le6hm/cpi)
set clminr = cover*lminr
set csqlminr = (cover**2)*lminr

***fix ups for messed up minimum wage
set yk2 = yk
set yk2 90:01 90:01 = yk*3.35/3.25
set minr2 = minr
set minr2 90:01 90:01 = minr*3.35/3.25
set lyk2 = log(yk2)
set lminr2 = log(minr2)
set minchk = minr*cpi/144.4583
set minchk2 = minr2*cpi/144.4583

print / minchk2 cover yk2 wmanr
print / emp empam empw empnw empm empf empb


table 54:01 93:04 emp empm empf empw empnw empb empam
table 54:01 93:04 sy sym syf syw synw syb
table 54:01 93:04 yk2 cover wmanr minchk2

*******  Correlations

cross emp empam / 0 0
cross lemp lempam / 0 0
cross lemp lupr / 0 0
cross lemp lempf / 0 0
cross lemp lempb / 0 0
cross lemp lempw / 0 0
cross lemp lyk2 / 0 0


******* Wellington's basic spec. log version
***ar(1) version

ar1(method=maxl) lemp
 #constant lyk2 lupr sy afp lpop trend trendsq seasons{-2 to 0} $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

*** add adult male epop
ar1(method=maxl) lemp
 #constant lyk2 lupr lempam sy afp lpop trend trendsq seasons{-2 to 0} $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

***************************************************************************
*******FEMALES
*******AR1 -- with and without adult male epop -- log model

ar1(method=maxl) lempf
 #constant lyk2 lupr syf afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

ar1(method=maxl) lempf
 #constant lyk2 lupr lempam syf afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4


*******MALES
*******AR1 -- with and without adult male epop -- log model

ar1(method=maxl) lempm
 #constant lyk2 lupr sym afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

ar1(method=maxl) lempm
 #constant lyk2 lupr lempam sym afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4


******* WHITES
*******AR1 -- with and without adult male epop -- log model

ar1(method=maxl) lempw
 #constant lyk2 lupr syw afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

ar1(method=maxl) lempw
 #constant lyk2 lupr lempam syw afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4



******* NONWHITES
*******AR1 -- with and without adult male epop -- log model

ar1(method=maxl) lempnw
 #constant lyk2 lupr synw afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

ar1(method=maxl) lempnw
 #constant lyk2 lupr lempam synw afp lpop trend trendsq seasons{-2 to 0}  $
  tq2 tq3 tq4 t2q2 t2q3 t2q4

end