***************************************************************
**Syntax producing the episode and aggregate version of the  **
**USA 1965-66 time use survey, created by Dr. Kimberly Fisher**
**Centre for Time Use Research, University of Oxford, UK     **
**20 August 2010, updated  24 January 2012                   **
***************************************************************.

*extract additional variables.

GET FILE='C:\Users\Kimberly Fisher\Documents\MTUS-archive\USA\1965-66\us65-66turp.sav'.
recode v245 v246 v518 v199 (sysmis=-8).
freq v245 v246 v518.
cro v245 by v246 v518.

Do if (v245>1 and v245<6) or ((v245=1 and (v246>0 and v246<4)) or (v245=0 and (v246=2 or v246=3))).
   compute vehicle=4.
Else if v245=1 or v245=6 or (v245<1 and v246=1).
   compute vehicle=3.
Else if v245=0.
   compute vehicle=0.
Else.
   compute vehicle=-8.
End if.
cro vehicle by v245 v246.
cro vehicle by v199.
if vehicle=-8 and v199>0 vehicle=3.
freq vehicle.

Do if id<1000.
     compute year=1965.
Else if id>1999.
     compute year=1966.
End if. 

rename variables (id=pid).
compute survey=1.
if urbrural>0 survey=0.
recode urbrural (sysmis=-8).
val lab survey
     0 "1965-66 USA Szalai Jackson sample"
     1 "1965-66 USA national sample".
cro survey by urbrural.

sort cases by survey pid.

save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa1965addveh.sav'
   /keep survey pid vehicle.

*create background variables.

Get File='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Questionnaire\USA65_66quest.sav'.

sort cases by survey pid.

match files files=* 
      files='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa1965addveh.sav'
   /BY survey pid.
select if sex>0.
recode vehicle (sysmis=-8).
freq vehicle.

rename variables (pid hhtype urban ageyngst wagelm survey= 
       persid hhtypeo urbano agekid2 empinclm surveyo).

compute hldid=0.
compute computer=-9.
compute citizen=-9.
compute sector=-9.
compute rushed=-9.
compute health=-9.
recode surveyo (0=1) (1=2) into msamp.
freq rushed health sex msamp.

Do if under5>0.
   compute agekidx=1.
Else if under18>0.
   compute agekidx=2.
Else if under18=-8 and (hhtypeo=1 or hhtypeo=3 or hhtypeo=8).
   compute agekidx=-8.
Else.
   compute agekidx=-7.
End if.
freq agekidx.
cro agekidx by under5 under18.

Temp.
select if (hhtypeo=1 or hhtypeo=3) and under18=0.
freq age.

do if hhtypeo=2.
   compute hhtype = 2.
else if hhtypeo=1.
   compute hhtype = 3.
else if under18<1 and (hhtypeo=6 or hhtypeo=7).
   compute hhtype = 1.
else if hhtypeo>0.
   compute hhtype = 4.
else.
   compute hhtype = -8. 
end if.

compute test=0.
if hhtype=4 and civstat=1 test=1.
freq test.
If test=1 hhtype=3.
cro hhtype by hhtypeo civstat.

compute nchild = under18.
Do if under18>0 and nadult>0.
   compute hhldsize=under18+nadult.
Else if nadult>0.
   compute hhldsize=nadult.
Else.
   compute hhldsize=-8.
End if.
*correct a number of inconsistencies and errors - 
*1 person households, unknown number of children 
*and the like.
Do if hhldsize=-8 and hhtype=1.
   compute hhldsize=1.
   compute nchild=0.
   compute agekidx=-7.
Else if hhldsize=-8 and hhtype=2.
   compute hhldsize=2. 
   compute nchild=0.
   compute agekidx=-7.
End if.

compute test=0.
if hhldsize<2 and hhtype=4 test=1.
freq test.

Do if test=1.
   compute hhldsize=1.
   compute hhtype=1.
   compute nchild=0.
End if.

*visually check the 14 cases of hhldsize=-8 and hhtype=3.
temp.
select if hhldsize=-8.
freq age.
*these are all cases of people missing age, so they will be 
*0-weighted. No action taken.

compute test=0.
if hhldsize=2 and hhtype=3 test=1.
freq test.

Do if test=1 and nadult=2 and under5>0.
   compute hhldsize=nadult+under5.
Else if test=1 and nadult=2 and under5=0.
   compute hhtype=2.
Else if test=1 and nadult=1 and under5=1.
   compute hhldsize=3.
Else if test=1.
   compute hhldsize=3.
End if.

if hhldsize>2 and hhtype=2 hhtype=3.
freq hhldsize.
cro hhldsize by hhtype nchild.

recode civstat (1=1) (else=2). 
cro civstat by hhtype hhtypeo.
cro famstat by nchild agekidx.

recode incomeqt (1=1) (2,3=2) (4=3) (else=-8) into income.
If incomeqt=1 incorig=50.
If incomeqt=2 incorig=51.
If incomeqt=3 incorig=52.
If incomeqt=4 incorig=53.
If incomeqt=-8 incorig=99.
cro incorig by income incomeqt.

recode urbano (1=1) (0=2) (else=-8) into urban.
cro urbano by urban.

compute relrefp=1.
compute parntid1=-9.
compute parntid2=-9.
compute partid=-9.
compute cphome=-9.

do if hhtypeo=3 or (hhtypeo=5 & under18>0).
   compute singpar = 1.
else if hhtypeo>0.
   compute singpar = 0.
else.
   compute singpar=-8.
end if.
freq singpar.
cro singpar by sex hhtypeo.

freq wkhrs.
cro wkhrs by fulltime parttime noemploy.
compute workhrs=wkhrs.
If (fulltime=1 or parttime=1) and workhrs=-4 workhrs=5.
If noemploy=1 and workhrs=-4 workhrs=-7.
execute.
recode workhrs (-4=-8).
means workhrs by fulltime parttime
   /cells=min max mean.

do if fulltime=1.
   compute empstat=1.
else if parttime=1.
   compute empstat=2.
else if workhrs>0.
   compute empstat=3.
else if noemploy=1 or wkhrs=0.
   compute empstat=4.
else.
   compute empstat=-8.
end if.
cro empstat by fulltime parttime.

recode empstat (1 thru 3=1) (4=0) (else=-8) into emp.

freq cohab unemp retired disab student.
*cohab, unemp, retired, disab & student 
*already have the values/names for the MTUS.
if age<35 and retired=-8 retired=0.
means age by retired by empstat
   /cells min max mean.

compute educa=educ.
if (educa >0 and educa <3) edcat = 1.
if educa=3 edcat = 2.
if educa>3 edcat = 3.
if educa<0 edcat = -8.

apply dictionary from 
   'C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\MTUS database\schema&etc\MTUShafschema.sav'.

freq hhtype nchild incorig civstat cohab empstat emp unemp student retired edcat singpar disab workhrs educa edcat.

means workhrs by unemp student retired empstat.
cro incorig by income.
cro urbano by urban.
if singpar=1 and hhldsize=1 singpar=0.
cro hhtype by singpar.
cro empstat by fulltime parttime noemploy. 
cro emp by empstat.
cro unemp student by empstat.
cro retired by empstat. 
cro educ by edcat.
compute empsp=-9.
freq occup.
recode occup (2,6=2) (3,4=3) (5=4) (7=5) (8=6) (9=7) 
       (10=8) (11=10) (12,15=11) (13=9) (14=12) 
       (16=13) (17=14) (else=copy).
val lab educa 
   1 "grades 0-8 only"
   2 "grades 9-11, not high school graduate"
   3 high school graduate
   4 some college
   5 college graduate
   6 postgraduate.

sort cases by msamp hldid persid.

save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa65ind.sav'
   /keep msamp hldid persid parntid1 parntid2 partid hhtype hhldsize nchild agekidx agekid2 incorig income vehicle
    ownhome urban computer famstat cphome singpar relrefp civstat cohab citizen empstat emp empinclm
    unemp student retired empsp workhrs occup sector educa edcat disab health rushed.

*pull out sex and age for the episode file.
save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\sexage.sav'
   /keep msamp hldid persid sex age.

*time use variables and episode file.
GET FILE='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Diary\diary1965\USA1965hfep.sav'.
rename variables (pid=persid).
recode survey (0=1) (1=2) into msamp.

match files files=*
      table='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\sexage.sav'
   /BY msamp persid.
desc age.

compute countrya=38.
compute survey=1965.
compute swave=0.
compute id=1.
compute day=diaryday.
*note this dataset already contains the variables 
*month and year in mtus format.
compute diary=wave.
vector avx(41).
   loop i=1 to 41.
      compute avx(i)=0.
   end loop.
execute.
vector mainx(69) secx(69).
   loop i=1 to 69.
      compute mainx(i)=0.
      compute secx(i)=69.
   end loop.
compute mainx=0.
compute secx=69.
compute av=0.
compute carer=0.
if main=40 carer=1.
freq carer.

apply dictionary from 
   'C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\MTUS database\schema&etc\MTUShafschema.sav'.

cro day by diaryday.
cro main by carer.

***********************************
*computation of time use variables*
***********************************.

*fill in AV variables.
compute residual=0.
do if main=8 or main=10 or main=13 or main=14.
   compute avx1=time.
   compute av=1.
else if main=11.
   compute avx2=time.
   compute av=2.
else if main=12 or main=15.
   compute avx3=time.
   compute av=3.
else if main=16 or main=18 or main=19.
   compute avx4=time.
   compute av=4.
else if main=92 or main=93 or main=94.
   compute avx5=time.
   compute av=5.
else if main=20 or main=21.
   compute avx6=time.
   compute av=6.
else if main=22 or main=23.
   compute avx7=time.
   compute av=7.
else if main=24 or main=25 or main=40 or main=68.
   compute avx8=time.
   compute av=8.
else if main=67.
   compute avx9=time.
   compute av=9.
else if main=26 or main=27 or main=30 or main=31.
   compute avx10=time.
   compute av=10.
else if main=65 or (main>32 and main<40).
   compute avx11=time.
   compute av=11.
else if main=91 or (main>94 and main<98).
   compute avx12=time.
   compute av=12.
else if main=1 or main=2 or main=6.
   compute avx13=time.
   compute av=13.
else if main=7 or main=28 or main=29 or main=32.
   compute avx14=time.
   compute av=14.
else if main=9 and eloc~=6.
   compute avx15=time.
   compute av=15.
else if main=3 or main=4 or main=5.
   compute avx16=time.
   compute av=16.
else if main=90 or main=98.
   compute avx17=time.
   compute av=17.
else if main=50 or main=54 or main=55 or main=64.
   compute avx18=time.
   compute av=18.
else if main=60 or main=63 or main=66.
   compute avx19=time.
   compute av=19.
else if main=51.
   compute avx20=time.
   compute av=20.
else if main=62.
   compute avx21=time.
   compute av=21.
else if main=49.
   compute avx22=time.
   compute av=22.
else if (main>40 and main<47) or main=48.
   compute avx23=time.
   compute av=23.
else if main=52 or main=53.
   compute avx24=time.
   compute av=24.
else if main=57.
   compute avx25=time.
   compute av=25.
else if main=56.
   compute avx27=time.
   compute av=27.
else if main=9 and eloc=6.
   compute avx28=time.
   compute av=28.
else if (main=71 or main=72) and eloc>1.
   compute avx29=time.
   compute av=29.
else if main=85.
   compute avx30=time.
   compute av=30.
else if main=86.
   compute avx31=time.
   compute av=31.
else if main=84.
   compute avx32=time.
   compute av=32.
else if main=17.
   compute avx33=time.
   compute av=33.
else if main=81.
   compute avx34=time.
   compute av=34.
else if main=82 or main=83.
   compute avx35=time.
   compute av=35.
else if main=78.
   compute avx36=time.
   compute av=36.
else if main=88.
   compute avx37=time.
   compute av=37.
else if (main=71 or main=72) and eloc<2.
   compute avx38=time.
   compute av=38.
else if main=76.
   compute avx39=time.
   compute av=39.
else if main=70 or main=73 or main=74 or main=75 or main=77 or main=87 or main=89.
   compute avx40=time.
   compute av=40.
else if main=-8 or main=58.
   compute avx41=time.
   compute av=41.
else.
   compute residual=main.
end if.
freq residual av.

temp.
select if av<8.
cro main by av.

temp.
select if av>7 and av<16.
cro main by av.

temp.
select if av>15 and av<24.
cro main by av.

temp.
select if av>23 and av<32.
cro main by av.

temp.
select if av>31.
cro main by av.

*check if can separately identify dog walking.
temp.
select if main=68.
freq mtrav sec.
*9 episodes potentially could count as dog walking, 
*not enough to code separately.

*fill in main variables.
compute residual=0.
Do if main=2.
   compute mainx1=time.
   compute mainx=1.
Else if main=4.
   compute mainx3=time.
   compute mainx=3.
Else if main=3 or main=5.
   compute mainx2=time.
   compute mainx=2.
Else if main=1 or main=6 or main=7.
   compute mainx4=time.
   compute mainx=4.
Else if main=8.
   compute mainx5=time.
   compute mainx=5.
Else if main=9.
   compute mainx6=time.
   compute mainx=6.
Else if main=10.
   compute mainx7=time.
   compute mainx=7.
Else if main=11.
   compute mainx8=time.
   compute mainx=8.
Else if main=12.
   compute mainx9=time.
   compute mainx=9.
Else if main=92.
   compute mainx11=time.
   compute mainx=11.
Else if main=13.
   compute mainx12=time.
   compute mainx=12.
Else if main=14.
   compute mainx13=time.
   compute mainx=13.
Else if main=15.
   compute mainx14=time.
   compute mainx=14.
Else if main=16.
   compute mainx15=time.
   compute mainx=15.
Else if main=17.
   compute mainx16=time.
   compute mainx=16.
Else if main=18 or main=19.
   compute mainx17=time.
   compute mainx=17.
Else if main=20.
   compute mainx18=time.
   compute mainx=18.
Else if main=21.
   compute mainx19=time.
   compute mainx=19.
Else if main=22.
   compute mainx20=time.
   compute mainx=20.
Else if main=23.
   compute mainx21=time.
   compute mainx=21.
Else if main=24.
   compute mainx22=time.
   compute mainx=22.
Else if main=25.
   compute mainx23=time.
   compute mainx=23.
Else if main=26 or main=27.
   compute mainx24=time.
   compute mainx=24.
Else if main=28 or main=29.
   compute mainx25=time.
   compute mainx=25.
Else if main=30 or main=31 or main=32.
   compute mainx26=time.
   compute mainx=26.
Else if main=68.
   compute mainx27=time.
   compute mainx=27.
Else if main=33 or main=34 or main=35.
   compute mainx28=time.
   compute mainx=28.
Else if main=37.
   compute mainx29=time.
   compute mainx=29.
Else if main=36 or main=38 or main=65.
   compute mainx30=time.
   compute mainx=30.
Else if main=39.
   compute mainx31=time.
   compute mainx=31.
Else if main=40.
   compute mainx32=time.
   compute mainx=32.
Else if (main>40 and main<47) or main=48.
   compute mainx33=time.
   compute mainx=33.
Else if main=49.
   compute mainx34=time.
   compute mainx=34.
Else if main=50.
   compute mainx35=time.
   compute mainx=35.
Else if main=51.
   compute mainx36=time.
   compute mainx=36.
Else if main=52 or main=53.
   compute mainx37=time.
   compute mainx=37.
Else if main=54 or main=55.
   compute mainx38=time.
   compute mainx=38.
Else if main=56 or (main=9 and eloc=6).
   compute mainx39=time.
   compute mainx=39.
Else if main=57.
   compute mainx40=time.
   compute mainx=40.
Else if main=58.
   compute mainx41=time.
   compute mainx=41.
Else if main=60.
   compute mainx42=time.
   compute mainx=42.
Else if main=62.
   compute mainx43=time.
   compute mainx=43.
Else if main=63.
   compute mainx44=time.
   compute mainx=44.
Else if main=64.
   compute mainx45=time.
   compute mainx=45.
Else if main=66 or main=67.
   compute mainx46=time.
   compute mainx=46.
Else if main=70.
   compute mainx51=time.
   compute mainx=51.
Else if main=72.
   compute mainx48=time.
   compute mainx=48.
Else if main=88.
   compute mainx49=time.
   compute mainx=49.
Else if main=71 or main=73.
   compute mainx50=time.
   compute mainx=50.
Else if main=74 or main=75.
   compute mainx52=time.
   compute mainx=52.
Else if main=87.
   compute mainx53=time.
   compute mainx=53.
Else if main=76 or main=77.
   compute mainx54=time.
   compute mainx=54.
Else if main=78.
   compute mainx55=time.
   compute mainx=55.
Else if main=81 or main=82 or main=83.
   compute mainx56=time.
   compute mainx=56.
Else if main=84.
   compute mainx57=time.
   compute mainx=57.
Else if main=85.
   compute mainx58=time.
   compute mainx=58.
Else if main=86.
   compute mainx59=time.
   compute mainx=59.
Else if main=89.
   compute mainx61=time.
   compute mainx=61.
Else if main=90.
   compute mainx62=time.
   compute mainx=62.
Else if main=93.
   compute mainx63=time.
   compute mainx=63.
Else if main=94.
   compute mainx64=time.
   compute mainx=64.
Else if main=97.
   compute mainx65=time.
   compute mainx=65.
Else if main=96.
   compute mainx66=time.
   compute mainx=66.
Else if main=91 or main=95.
   compute mainx67=time.
   compute mainx=67.
Else if main=98.
   compute mainx68=time.
   compute mainx=68.
Else if main=-8.
   compute mainx69=time.
   compute mainx=69.
else.
   compute residual=main.
end if.
freq residual mainx.

temp.
select if mainx<11.
cro main by mainx.

temp.
select if mainx>10 and mainx<21.
cro main by mainx.

temp.
select if mainx>20 and mainx<31.
cro main by mainx.

temp.
select if mainx>30 and mainx<41.
cro main by mainx.

temp.
select if mainx>40 and mainx<51.
cro main by mainx.

temp.
select if mainx>50 and mainx<61.
cro main by mainx.

temp.
select if mainx>60.
cro main by mainx.

*fill in secondary activity variables.
compute residual=0.
Do if sec=2.
   compute secx1=time.
   compute secx=1.
Else if sec=4.
   compute secx3=time.
   compute secx=3.
Else if sec=3 or sec=5.
   compute secx2=time.
   compute secx=2.
Else if sec=1 or sec=6 or sec=7.
   compute secx4=time.
   compute secx=4.
Else if sec=8.
   compute secx5=time.
   compute secx=5.
Else if sec=9.
   compute secx6=time.
   compute secx=6.
Else if sec=10.
   compute secx7=time.
   compute secx=7.
Else if sec=11.
   compute secx8=time.
   compute secx=8.
Else if sec=12.
   compute secx9=time.
   compute secx=9.
Else if sec=92.
   compute secx11=time.
   compute secx=11.
Else if sec=13.
   compute secx12=time.
   compute secx=12.
Else if sec=14.
   compute secx13=time.
   compute secx=13.
Else if sec=15.
   compute secx14=time.
   compute secx=14.
Else if sec=16.
   compute secx15=time.
   compute secx=15.
Else if sec=17.
   compute secx16=time.
   compute secx=16.
Else if sec=18 or sec=19.
   compute secx17=time.
   compute secx=17.
Else if sec=20.
   compute secx18=time.
   compute secx=18.
Else if sec=21.
   compute secx19=time.
   compute secx=19.
Else if sec=22.
   compute secx20=time.
   compute secx=20.
Else if sec=23.
   compute secx21=time.
   compute secx=21.
Else if sec=24.
   compute secx22=time.
   compute secx=22.
Else if sec=25.
   compute secx23=time.
   compute secx=23.
Else if sec=26 or sec=27.
   compute secx24=time.
   compute secx=24.
Else if sec=28 or sec=29.
   compute secx25=time.
   compute secx=25.
Else if sec=30 or sec=31 or sec=32.
   compute secx26=time.
   compute secx=26.
Else if sec=68.
   compute secx27=time.
   compute secx=27.
Else if sec=33 or sec=34 or sec=35.
   compute secx28=time.
   compute secx=28.
Else if sec=37.
   compute secx29=time.
   compute secx=29.
Else if sec=36 or sec=38 or sec=65.
   compute secx30=time.
   compute secx=30.
Else if sec=39.
   compute secx31=time.
   compute secx=31.
Else if sec=40.
   compute secx32=time.
   compute secx=32.
Else if (sec>40 and sec<47) or sec=48.
   compute secx33=time.
   compute secx=33.
Else if sec=49.
   compute secx34=time.
   compute secx=34.
Else if sec=50.
   compute secx35=time.
   compute secx=35.
Else if sec=51.
   compute secx36=time.
   compute secx=36.
Else if sec=52 or sec=53.
   compute secx37=time.
   compute secx=37.
Else if sec=54 or sec=55.
   compute secx38=time.
   compute secx=38.
Else if sec=56 or (sec=9 and eloc=6).
   compute secx39=time.
   compute secx=39.
Else if sec=57.
   compute secx40=time.
   compute secx=40.
Else if sec=58.
   compute secx41=time.
   compute secx=41.
Else if sec=60.
   compute secx42=time.
   compute secx=42.
Else if sec=62.
   compute secx43=time.
   compute secx=43.
Else if sec=63.
   compute secx44=time.
   compute secx=44.
Else if sec=64.
   compute secx45=time.
   compute secx=45.
Else if sec=66 or sec=67.
   compute secx46=time.
   compute secx=46.
Else if sec=70.
   compute secx51=time.
   compute secx=51.
Else if sec=72.
   compute secx48=time.
   compute secx=48.
Else if sec=88.
   compute secx49=time.
   compute secx=49.
Else if sec=71 or sec=73.
   compute secx50=time.
   compute secx=50.
Else if sec=74 or sec=75.
   compute secx52=time.
   compute secx=52.
Else if sec=87.
   compute secx53=time.
   compute secx=53.
Else if sec=76 or sec=77.
   compute secx54=time.
   compute secx=54.
Else if sec=78.
   compute secx55=time.
   compute secx=55.
Else if sec=81 or sec=82 or sec=83.
   compute secx56=time.
   compute secx=56.
Else if sec=84.
   compute secx57=time.
   compute secx=57.
Else if sec=85.
   compute secx58=time.
   compute secx=58.
Else if sec=86.
   compute secx59=time.
   compute secx=59.
Else if sec=89.
   compute secx61=time.
   compute secx=61.
Else if sec=90.
   compute secx62=time.
   compute secx=62.
Else if sec=93.
   compute secx63=time.
   compute secx=63.
Else if sec=94.
   compute secx64=time.
   compute secx=64.
Else if sec=97.
   compute secx65=time.
   compute secx=65.
Else if sec=96.
   compute secx66=time.
   compute secx=66.
Else if sec=91 or sec=95.
   compute secx67=time.
   compute secx=67.
Else if sec=98.
   compute secx68=time.
   compute secx=68.
else.
   compute residual=sec.
end if.
freq residual secx.

*cday, time, start, end, epnum, eloc are in the MTUS W6 format.
compute clockst=clockst/100.
recode inout (1=2) (2=1) (else=copy).
compute ict=-9.
freq ict mtrav inout alone sppart child.

*alone, sppart, child already in MTUS format.
Do if otherp=1 or shoprof=1 or sppart=1.
   compute oad=1.
Else if hhadult=-7.
   compute oad=-7.
Else if cowork=1 or hhadult=1.
   compute oad=1.
Else.
   compute oad=0.
End if.
cro oad by otherp cowork shoprof hhadult sppart.

*create sppart time for W58.
compute sptime=0.
if sppart=1 sptime=time.
desc sptime.

cro baddem by lowqual by day.
*27 episodes missing day of week.
compute eat=0.
compute trv=0.
compute pcr=0.
compute slp=0.
If main=8 or main=9 or main=20 or main=21 or main=56 or 
   sec=2 or sec=8 or sec=9 or sec=20 or sec=21 eat=1. 
If main=3 or main=4 or main=5 or main=78 or sec=3 or 
   sec=5 or sec=78 slp=1.
If main=1 or main=2 or main=6 or main=7 or main=29 or 
   main=30 or sec=1 or sec=6 or sec=7 or sec=29 or sec=30 pcr=1.
If (main>59 and main<70) or (main>89 and main<99) or (sec>59
   and sec<70) or (sec>89 and sec<99) trv=1.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid persid
  /miss=SUM(mainx69) /maxep=MAX(epnum) /eatdrk=MAX(eat)
  /sleep=MAX(slp) /pcare=MAX(pcr) /trav=MAX(trv).
freq miss maxep.
compute nep=0.
if maxep<7 nep=1.
compute dobasic=eatdrk+sleep+pcare+trav.
freq dobasic nep.
compute misbasic=0.
if lowqual=1 and dobasic<3 misbasic=1.
cro misbasic by nep maxep.

Do If baddem=0 and lowqual=0.
   compute badcase=0.
Else if baddem=1 and lowqual=0 and day>0.
   compute badcase=1.
Else if baddem=0 and lowqual=0 and day=-8.
   compute badcase=2.
Else if baddem=0 and day>-8 and nep=0 and miss=1 and misbasic=0.
   compute badcase=3.
Else if baddem=0 and day>-8 and nep=1 and miss=0 and misbasic=0.
   compute badcase=4.
Else if baddem=0 and day>-8 and nep=0 and miss=0 and misbasic=1.
   compute badcase=5.
Else if baddem=1 and day=-8 and nep=0 and miss=0 and misbasic=0.
   compute badcase=6.
Else if baddem=1 and day>-8 and nep=0 and miss=1 and misbasic=0.
   compute badcase=7.
Else if baddem=1 and day>-8 and nep=1 and miss=0 and misbasic=0.
   compute badcase=8.
Else if baddem=1 and day>-8 and nep=0 and miss=0 and misbasic=1.
   compute badcase=9.
Else if baddem=0 and day=-8 and nep=0 and miss=1 and misbasic=0.
   compute badcase=10.
Else if baddem=0 and day=-8 and nep=1 and miss=0 and misbasic=0.
   compute badcase=11.
Else if baddem=0 and day=-8 and nep=0 and miss=0 and misbasic=1.
   compute badcase=12.
Else if baddem=0 and day>-8 and nep=1 and miss=1 and misbasic=0.
   compute badcase=13.
Else if baddem=0 and day>-8 and nep=0 and miss=1 and misbasic=1.
   compute badcase=14.
Else if baddem=0 and day>-8 and nep=1 and miss=0 and misbasic=1.
   compute badcase=15.
Else if baddem=1 and day=-8 and nep=0 and miss=1 and misbasic=0.
   compute badcase=16.
Else if baddem=1 and day=-8 and nep=1 and miss=0 and misbasic=0.
   compute badcase=17.
Else if baddem=1 and day=-8 and nep=0 and miss=0 and misbasic=1.
   compute badcase=18.
Else if baddem=1 and day>-8 and nep=1 and miss=1 and misbasic=0.
   compute badcase=19.
Else if baddem=1 and day>-8 and nep=0 and miss=1 and misbasic=1.
   compute badcase=20.
Else if baddem=1 and day>-8 and nep=1 and miss=0 and misbasic=1.
   compute badcase=21.
Else if baddem=0 and day=-8 and nep=1 and miss=1 and misbasic=0.
   compute badcase=22.
Else if baddem=0 and day=-8 and nep=0 and miss=1 and misbasic=1.
   compute badcase=23.
Else if baddem=0 and day=-8 and nep=1 and miss=0 and misbasic=1.
   compute badcase=24.
Else if baddem=0 and day>-8 and nep=1 and miss=1 and misbasic=1.
   compute badcase=25.
Else if baddem=1 and day=-8 and nep=1 and miss=1 and misbasic=0.
   compute badcase=26.
Else if baddem=1 and day=-8 and nep=0 and miss=1 and misbasic=1.
   compute badcase=27.
Else if baddem=1 and day=-8 and nep=1 and miss=0 and misbasic=1.
   compute badcase=28.
Else if baddem=1 and day>-8 and nep=1 and miss=1 and misbasic=1.
   compute badcase=29.
Else if baddem=0 and day=-8 and nep=1 and miss=1 and misbasic=1.
   compute badcase=30.
Else.
   compute badcase=31.
End if.
val lab badcase
   0 good case
   1 miss age or sex only
   2 miss day of week only
   3 "miss 91+ minutes only"
   4 "<7 episodes only"
   5 "miss 2+ basic acts only"
   6 "miss age/sex + day of week"
   7 "miss age/sex + 91+ min"
   8 "miss age/sex + <7 episodes"
   9 "miss age/sex + 2+ basic acts"
  10 "miss day of week + 91+ min"
  11 "miss day of week + <7 episodes"
  12 "miss day of week + 2+ basic acts"
  13 "miss 91+ min + <7 episodes"
  14 "miss 91+ min + 2+ basic acts"
  15 "<7 ep & miss 2+ basic acts"
  16 "miss age/sex + day or week + 91+ min"
  17 "miss age/sex + day or week + <7 episodes"
  18 "miss age/sex + day or week + 2+ basic acts"
  19 "miss age/sex + 91+ min + <7 eps"
  20 "miss age/sex + 91+ min + 2+ basic acts"
  21 "miss age/sex + 2+ basic acts + <7 eps"
  22 "miss day or week + 91+ min + <7 eps"
  23 "miss day or week + 91+ min + 2+ basic acts"
  24 "miss day or week + 2+ basic acts + <7 eps"
  25 "miss 91+ min + 2+ basic acts + <7 eps"
  26 "miss age/sex + day or week + 91+ min + <7 episodes"
  27 "miss age/sex + day or week + 91+ min + 2+ basic acts"
  28 "miss age/sex + day or week + <7 episodes + 2+ basic acts"
  29 "miss age/sex + 91+ min + <7 episodes + 2+ basic acts"
  30 "day or week + 91+ min + <7 episodes + 2+ basic acts"
  31 bad on all five points.
freq badcase.

sort cases by countrya survey swave msamp hldid persid id epnum.  

SAVE OUTFILE='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa1965hef.sav'
  /keep countrya survey swave msamp hldid persid id day cday month year diary badcase sex age 
   time clockst start end epnum mainx secx av inout eloc ict mtrav alone child sppart oad.

aggregate
  /outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\sum.sav'
  /break=countrya survey swave msamp hldid persid id
  /day=max(day) /month=max(month) /year=max(year) /diary=max(diary) /badcase=max(badcase) 
  /sex=MAX(sex) /age=MAX(age) /carer=max(carer) 
  /av1 = sum(avx1) /av2 = sum(avx2)  /av3 = sum(avx3) /av4 = sum(avx4)  
  /av5 = sum(avx5) /av6 = sum(avx6)  /av7 = sum(avx7) /av8 = sum(avx8)  
  /av9 = sum(avx9) /av10= sum(avx10) /av11=sum(avx11) /av12=sum(avx12)  
  /av13=sum(avx13) /av14= sum(avx14) /av15=sum(avx15) /av16=sum(avx16)  
  /av17=sum(avx17) /av18= sum(avx18) /av19=sum(avx19) /av20=sum(avx20)  
  /av21=sum(avx21) /av22= sum(avx22) /av23=sum(avx23) /av24=sum(avx24)  
  /av25=sum(avx25) /av26= max(avx26) /av27=sum(avx27) /av28=sum(avx28)  
  /av29=sum(avx29) /av30= sum(avx30) /av31=sum(avx31) /av32=sum(avx32)  
  /av33=sum(avx33) /av34= sum(avx34) /av35=sum(avx35) /av36=sum(avx36)  
  /av37=sum(avx37) /av38= sum(avx38) /av39=sum(avx39) /av40=sum(avx40)  
  /av41=sum(avx41) /main1=sum(mainx1) /main2=sum(mainx2) /main3=SUM(mainx3)
  /main4 =sum(mainx4 ) /main5 =SUM(mainx5 ) /main6 =SUM(mainx6 ) /main7 =SUM(mainx7 )
  /main8 =SUM(mainx8 ) /main9 =SUM(mainx9 ) /main10=SUM(mainx10) /main11=SUM(mainx11)
  /main12=SUM(mainx12) /main13=SUM(mainx13) /main14=SUM(mainx14) /main15=SUM(mainx15)
  /main16=SUM(mainx16) /main17=SUM(mainx17) /main18=SUM(mainx18) /main19=SUM(mainx19)
  /main20=SUM(mainx20) /main21=SUM(mainx21) /main22=SUM(mainx22) /main23=SUM(mainx23)
  /main24=SUM(mainx24) /main25=SUM(mainx25) /main26=SUM(mainx26) /main27=SUM(mainx27)
  /main28=SUM(mainx28) /main29=SUM(mainx29) /main30=SUM(mainx30) /main31=SUM(mainx31)
  /main32=SUM(mainx32) /main33=SUM(mainx33) /main34=SUM(mainx34) /main35=SUM(mainx35)
  /main36=SUM(mainx36) /main37=SUM(mainx37) /main38=SUM(mainx38) /main39=SUM(mainx39)
  /main40=SUM(mainx40) /main41=SUM(mainx41) /main42=SUM(mainx42) /main43=SUM(mainx43)
  /main44=SUM(mainx44) /main45=SUM(mainx45) /main46=SUM(mainx46) /main47=SUM(mainx47)
  /main48=SUM(mainx48) /main49=SUM(mainx49) /main50=SUM(mainx50) /main51=SUM(mainx51)
  /main52=SUM(mainx52) /main53=SUM(mainx53) /main54=SUM(mainx54) /main55=SUM(mainx55)
  /main56=SUM(mainx56) /main57=SUM(mainx57) /main58=SUM(mainx58) /main59=SUM(mainx59)
  /main60=SUM(mainx60) /main61=SUM(mainx61) /main62=SUM(mainx62) /main63=SUM(mainx63)
  /main64=SUM(mainx64) /main65=SUM(mainx65) /main66=SUM(mainx66) /main67=SUM(mainx67)
  /main68=SUM(mainx68) /main69=SUM(mainx69) /sppart=SUM(sptime)  /ocombwt=max(origwght) 
  /propwt= max(recwght).

*put aggregate file together.

match files
      files='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\sum.sav'
      table='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa65ind.sav'
     /by msamp hldid persid.
execute.

compute tottime = sum(av1 to av41).
compute tot2 = sum(main1 to main69).
freq tottime tot2.

desc av41 main69.
desc av1 to av40.
desc main1 to main68.
*recode time use variables which could not be create 
*for this study to -9.
recode av26 main1 main3 main10 main11 main14 main41 main44 main45 
       main47 main60 main61 main62 (0=-9).

apply dictionary from 'C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\MTUS database\schema&etc\MTUShafschema.sav'.

save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\USA1965haf.sav'
   /keep countrya survey swave msamp hldid persid id parntid1 parntid2 partid day 
    month year diary badcase hhtype hhldsize nchild agekidx agekid2 incorig income 
    ownhome urban computer vehicle sex age famstat cphome singpar relrefp civstat 
    cohab citizen empstat emp unemp student retired empsp workhrs empinclm occup 
    sector educa edcat rushed health carer disab av1 av2 av3 av4 av5 av6 av7 av8 
    av9 av10 av11 av12 av13 av14 av15 av16 av17 av18 av19 av20 av21 av22 av23 av24 
    av25 av26 av27 av28 av29 av30 av31 av32 av33 av34 av35 av36 av37 av38 av39 av40 
    av41 main1 main2 main3 main4 main5 main6 main7 main8 main9 main10 main11 main12 
    main13 main14 main15 main16 main17 main18 main19 main20 main21 main22 main23 
    main24 main25 main26 main27 main28 main29 main30 main31 main32 main33 main34 
    main35 main36 main37 main38 main39 main40 main41 main42 main43 main44 main45 
    main46 main47 main48 main49 main50 main51 main52 main53 main54 main55 main56 
    main57 main58 main59 main60 main61 main62 main63 main64 main65 main66 main67 
    main68 main69 sppart ocombwt propwt.

*final tidy of aggregate file.

get file='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\USA1965haf.sav'.

desc hldid persid id incorig age workhrs empinclm av1 av2 av3 av4 av5 av6 av7 av8 av9 av10 av11 
    av12 av13 av14 av15 av16 av17 av18 av19 av20 av21 av22 av23 av24 av25 av26 av27 av28 av29 
    av30 av31 av32 av33 av34 av35 av36 av37 av38 av39 av40 av41 ocombwt propwt.

desc age main1 to main69 sppart.
freq badcase.

AGGREGATE
  /OUTFILE=*
  MODE=ADDVARIABLES
  /BREAK=countrya
  /wtsum = SUM(propwt).
freq wtsum.

freq countrya survey swave msamp parntid1 parntid2 partid day month year diary badcase 
    hhtype hhldsize nchild agekidx agekid2 income ownhome urban computer vehicle sex age famstat 
    cphome singpar relrefp civstat cohab citizen empstat emp unemp student retired empsp workhrs 
    occup sector educa edcat rushed health carer disab incorig.

means age hhldsize nchild by famstat
  /cells min max.

means age hhldsize by nchild retired
  /cells min max median.

cro famstat by nchild civstat hhtype agekidx cphome.
cro agekidx by nchild.
cro empstat by emp unemp retired student disab.
cro empsp by civstat relrefp partid.
cro civstat by hhtype relrefp partid cohab cphome.
*note no relrefp=2 should appear if civstat=2.
cro hhtype by hhldsize.
cro parntid1 by parntid2 cphome.

weight by propwt.
freq day badcase.
cro day by sex
   /cells= count column.
weight off.

save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\USA1965haf.sav'
   /drop wtsum.

*final tidy of episode file.

get file='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa1965hef.sav'.

rename variables (mainx secx=main sec).

apply dictionary from 
   'C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\MTUS database\schema&etc\MTUShefschema.sav'.

*check labels and distrubutions.
freq countrya survey swave msamp day cday month year diary badcase sex age 
     clockst epnum main sec av inout eloc ict mtrav alone child sppart oad.
if mtrav=3 and inout=3 inout=2.
if mtrav=5 and inout=1 inout=3.
cro mtrav by inout.
desc hldid persid id age epnum time start end.

save outfile='C:\Documents and Settings\Kimberly Fisher\My Documents\Kimberly\time\mtus\mtus data\W553\USA\1965-66\usa1965hef.sav'.