********************************************************
**Syntax producing the episode and aggregate versions **
**of the USA 1985 time use survey, created by Dr.     **
**Kimberly Fisher, Centre for Time Use Research,      **
**University of Oxford, UK, revised 24 January 2013   **
********************************************************.

*create background variables and construct household matrix.

Get File='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Questionnaire\USA85quest.sav'.

ADD FILES /FILE=*
  /RENAME (empflag income03 income65 income75 kidflag marrflag nkidflag wkhrflag=d0 d1 d2 d3 d4 d5 d6 d7)
  /FILE='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Supplementary files\1985\usa1985childsum.sav'
  /RENAME (athome atwrksc baddem caremflg cday diaryday elsewhr infltwt inside inveh locunk lowqual 
    lunk month numep origwght outside owghtflg recwght t0pcare t1paid t2ed t3unpaid t4acvol t5outhm 
    t6exerc t7inhm t8media t9trav tmain1 tmain10 tmain11 tmain12 tmain13 tmain14 tmain15 tmain16 
    tmain17 tmain18 tmain19 tmain2 tmain20 tmain21 tmain22 tmain23 tmain24 tmain25 tmain26 tmain27 
    tmain28 tmain29 tmain3 tmain30 tmain31 tmain32 tmain33 tmain34 tmain35 tmain36 tmain37 tmain38 
    tmain39 tmain4 tmain40 tmain41 tmain42 tmain43 tmain44 tmain45 tmain46 tmain48 tmain49 tmain5 
    tmain50 tmain51 tmain52 tmain53 tmain54 tmain55 tmain56 tmain57 tmain58 tmain6 tmain60 tmain62 
    tmain63 tmain64 tmain65 tmain66 tmain67 tmain68 tmain7 tmain70 tmain71 tmain72 tmain73 tmain74 
    tmain75 tmain76 tmain77 tmain78 tmain8 tmain81 tmain82 tmain83 tmain84 tmain85 tmain86 tmain87 
    tmain88 tmain89 tmain9 tmain90 tmain91 tmain92 tmain93 tmain94 tmain95 tmain96 tmain97 tmain98 
    tmiss tsc1 tsc10 tsc11 tsc12 tsc13 tsc14 tsc15 tsc16 tsc17 tsc18 tsc19 tsc2 tsc20 tsc21 tsc22 tsc23 
    tsc24 tsc25 tsc26 tsc27 tsc28 tsc29 tsc3 tsc30 tsc31 tsc32 tsc33 tsc34 tsc35 tsc36 tsc37 tsc38 
    tsc39 tsc4 tsc40 tsc41 tsc42 tsc43 tsc44 tsc45 tsc46 tsc48 tsc49 tsc5 tsc50 tsc51 tsc52 tsc53 tsc54 
    tsc55 tsc56 tsc57 tsc58 tsc6 tsc60 tsc62 tsc63 tsc64 tsc65 tsc66 tsc67 tsc68 tsc7 tsc70 tsc71 tsc72 
    tsc73 tsc74 tsc75 tsc76 tsc77 tsc78 tsc8 tsc81 tsc82 tsc83 tsc84 tsc85 tsc86 tsc87 tsc88 tsc89 tsc9 
    tsc90 tsc91 tsc92 tsc93 tsc94 tsc95 tsc96 tsc97 tsc98 walone wchild wclsfam withunk wother wsppart 
    xtimewt year=d8 d9 d10 d11 d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 d26 d27 d28 d29 
    d30 d31 d32 d33 d34 d35 d36 d37 d38 d39 d40 d41 d42 d43 d44 d45 d46 d47 d48 d49 d50 d51 d52 d53 d54 
    d55 d56 d57 d58 d59 d60 d61 d62 d63 d64 d65 d66 d67 d68 d69 d70 d71 d72 d73 d74 d75 d76 d77 d78 d79 
    d80 d81 d82 d83 d84 d85 d86 d87 d88 d89 d90 d91 d92 d93 d94 d95 d96 d97 d98 d99 d100 d101 d102 d103 
    d104 d105 d106 d107 d108 d109 d110 d111 d112 d113 d114 d115 d116 d117 d118 d119 d120 d121 d122 d123 
    d124 d125 d126 d127 d128 d129 d130 d131 d132 d133 d134 d135 d136 d137 d138 d139 d140 d141 d142 d143 
    d144 d145 d146 d147 d148 d149 d150 d151 d152 d153 d154 d155 d156 d157 d158 d159 d160 d161 d162 d163 
    d164 d165 d166 d167 d168 d169 d170 d171 d172 d173 d174 d175 d176 d177 d178 d179 d180 d181 d182 d183 
    d184 d185 d186 d187 d188 d189 d190 d191 d192 d193 d194 d195 d196 d197 d198 d199 d200 d201 d202 d203 
    d204 d205 d206 d207 d208 d209 d210 d211 d212 d213 d214 d215 d216 d217 d218 d219 d220 d221 d222 d223 
    d224 d225 d226 d227 d228 d229)
  /DROP=d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11 d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 
    d25 d26 d27 d28 d29 d30 d31 d32 d33 d34 d35 d36 d37 d38 d39 d40 d41 d42 d43 d44 d45 d46 d47 d48 d49 
    d50 d51 d52 d53 d54 d55 d56 d57 d58 d59 d60 d61 d62 d63 d64 d65 d66 d67 d68 d69 d70 d71 d72 d73 d74 
    d75 d76 d77 d78 d79 d80 d81 d82 d83 d84 d85 d86 d87 d88 d89 d90 d91 d92 d93 d94 d95 d96 d97 d98 d99 
    d100 d101 d102 d103 d104 d105 d106 d107 d108 d109 d110 d111 d112 d113 d114 d115 d116 d117 d118 d119 
    d120 d121 d122 d123 d124 d125 d126 d127 d128 d129 d130 d131 d132 d133 d134 d135 d136 d137 d138 d139 
    d140 d141 d142 d143 d144 d145 d146 d147 d148 d149 d150 d151 d152 d153 d154 d155 d156 d157 d158 d159 
    d160 d161 d162 d163 d164 d165 d166 d167 d168 d169 d170 d171 d172 d173 d174 d175 d176 d177 d178 d179 
    d180 d181 d182 d183 d184 d185 d186 d187 d188 d189 d190 d191 d192 d193 d194 d195 d196 d197 d198 d199 
    d200 d201 d202 d203 d204 d205 d206 d207 d208 d209 d210 d211 d212 d213 d214 d215 d216 d217 d218 d219 
    d220 d221 d222 d223 d224 d225 d226 d227 d228 d229.
EXECUTE.

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

sort cases by hldid persid.
*
*match files files=* /in=base
*      files='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\addrshhlth.sav' /in=add**
*   /BY hldid persid.
*cro base by add.

*remove non-diary cases.
*select if base=1.

compute computer=-9.
compute citizen=-9.
compute sector=-9.
compute vehicle=-9.
compute rushed=-9.
compute health=-9.
*recode wagerate (sysmis=-5).
*freq rushed health wagerate sex.
freq rushed health empinclm sex.

recode agekid2 (-8=-8) (-7=-7) (3=1) (9 thru 12=2)
   (13 thru 17=3) into agekidx.
means agekid2 by agekidx
   /cells = min max.

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.
*the two cases not yet coded are the two identical twins from the child file.
recode hhtype (sysmis=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.
freq hhldsize.
cro hhldsize by hhtype.
compute test=0.
if hhtype=4 and hhldsize=1 test=1.
freq test.
*correct inconsistencies between the household size and type variables.
if hhtype=4 and hhldsize=1 hhtype=1.
if hhtype=3 and hhldsize=2 hhldsize=3.
if hhtype=2 and hhldsize<2 hhldsize=2.
if hhtype=2 and hhldsize>2 hhtype=3.
cro hhldsize civstat by hhtype.
compute test=0.
if civstat>1 and hhtype=2 test=1.
freq test.

Do if hldid=152 or hldid=305 or hldid=1481 or hldid=7047.
   compute hhtype=3.
   compute hhldsize=3.
Else if hldid=883 or hldid=989 or hldid=6149.
   compute hhtype=4.
End if.
cro hhldsize civstat by hhtype.
recode civstat (1=1) (2 thru high=2) (else=-8). 

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.

*Do if hldid=1466 and persid=5.
*   compute wagerate=0.
*   compute incorig=51.
*Else if hldid=1547 and persid=2.
*   compute wagerate=0.
*   compute incorig=99.
*End if.
*cro incorig by incomeqt.
*
*compute empinclm=-8.
*if empstat=3 empinclm=-7.
*if (empstat=1 or empstat=2) and wkhrs>0 and wagerate>0 empinclm=wagerate* wkhrs*4.
*desc empinclm.
*freq empinclm.

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

temp.
select if hhtype=2.
freq persid civstat.
*all people in couple only households have person ids 1 or 2.

Do if hhtype=2 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hhtype=2 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome=0.
   compute singpar=0.
Else if hhtype=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome=0.
   compute singpar=0.
Else if civstat=2.
   compute relrefp=-5.
   compute parntid1=-5.
   compute parntid2=-5.
   compute partid=-7.
   compute cphome=-5.
   compute singpar=-5.
Else.
   compute relrefp=-5.
   compute parntid1=-5.
   compute parntid2=-5.
   compute partid=-5.
   compute cphome=-5.
   compute singpar=-5.
End if.
freq relrefp parntid1 parntid2 partid cphome singpar.

compute coupsex1=0.
compute coupsex2=99.
if civstat=1 coupsex1=sex.
if civstat=1 coupsex2=sex.
compute coup=0.
if civstat=1 coup=1.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid
  /mincsex=MIN(coupsex2) /maxcsex=MAX(coupsex1) /ncoup=SUM(coup).
freq ncoup.
cro ncoup by hhtype.

temp.
select if hhtype=3 and ncoup=2 and civstat=1.
cro mincsex by maxcsex.

*deal with 4 person households. Assume that when a large
*age disparity, that the people closer in age are in the
*same couple, that person 1 is usually the hh reference 
*persons, that family is listed after parents followed
*by in-laws.
Do if hldid=1319 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1319 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1319 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=4.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1319 and persid=4.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1319 and persid=5.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7418 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7418 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7418 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=4.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7418 and persid=4.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7594 and persid=1.
   compute relrefp=1.
   compute parntid1=5.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7594 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7594 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7594 and persid=4.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7594 and persid=5.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7594 and persid=6.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7594 and persid=7.
   compute relrefp=6.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=6.
   compute cphome =0.
   compute singpar=0.
End if.
freq partid relrefp.
*if people are in a single person hh and missing 
*civstat, code as not in couple.
if hhtype=1 and civstat=-8 civstat=2.
cro ncoup by partid by civstat.

Do if hldid=517 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=4.
   compute cphome =0.
   compute singpar=0.
Else if hldid=517 and persid=2.
   compute relrefp=10.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=517 and persid=3.
   compute relrefp=10.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=517 and persid=4.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=549 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-8.
   compute cphome =0.
   compute singpar=0.
Else if hldid=549 and persid=2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-8.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=549 and persid=3.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=549 and persid=4.
   compute relrefp=9.
   compute parntid1=2.
   compute parntid2=3.
   compute cphome =1.
   compute singpar=0.
Else if hldid=549 and persid=5.
   compute relrefp=9.
   compute parntid1=2.
   compute parntid2=3.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1448 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-8.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1448 and persid=2.
   compute relrefp=5.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1448 and persid=3.
   compute relrefp=8.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6081 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6081 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6081 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute civstat=2.
   compute singpar=0.
Else if hldid=6180 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6180 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6180 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=4.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6180 and persid=4.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=7085 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7085 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7085 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute civstat=2.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7112 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7112 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7112 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute civstat=2.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7200 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7200 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7200 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute civstat=2.
   compute singpar=0. 
Else if hldid=7496 and persid=2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7496 and persid=3.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7496 and persid=4.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=5.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=7496 and persid=5.
   compute relrefp=6.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=4.
   compute cphome =0.
   compute singpar=0. 
   compute civstat=1.
End if.
cro ncoup by partid by civstat.

Do if hldid=138 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=138 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=152 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=152 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=168 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0. 
Else if hldid=168 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=194 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
   compute civstat=2.
Else if hldid=194 and persid>1.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=203 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=203 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=232 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-8.
   compute cphome =0.
   compute singpar=0.
Else if hldid=232 and persid>1.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
End if. 
cro ncoup by partid by civstat.

do if hldid=255 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=255 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=255 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=265 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=265 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=305 and persid=1.
   compute relrefp=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=305 and persid=2.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
Else if hldid=312 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=312 and persid>2.
   compute relrefp=1.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=399 and persid=1.
   compute relrefp=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute hhtype=4.
Else if hldid=399 and persid=2.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
   compute civstat=2.
   compute hhtype=4.
Else if hldid=399 and persid=3.
   compute relrefp=5.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute hhtype=4.
Else if hldid=513 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=513 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=524 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=524 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=794 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=794 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=807 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=807 and persid=2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=3.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=961 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
   compute civstat=2.
   compute hhtype=4.
Else if hldid=961 and (persid=2 or persid=3).
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=1.
   compute hhtype=4.
Else if hldid=961 and persid=4.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute hhtype=4.
Else if hldid=1091 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1091 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1142 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1142 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1161 and persid=1.
   compute relrefp=1.
   compute parntid1=2.
   compute parntid2=3.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1161 and persid=2.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1161 and persid=3.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=1161 and persid=4.
   compute relrefp=5.
   compute parntid1=2.
   compute parntid2=3.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1336 and persid=1.
   compute relrefp=1.
   compute parntid1=2.
   compute parntid2=3.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1336 and persid=2.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1343 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1343 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=1343 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1481 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1481 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1626 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1626 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1638 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1638 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1665 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=1665 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1705 and persid=1.
   compute relrefp=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=1705 and persid=2.
   compute relrefp=4.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
   compute civstat=2.
Else if hldid=1705 and persid>2.
   compute relrefp=5.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=6003 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=4.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6003 and (persid=2 or persid=3).
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=4.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=6003 and persid=4.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=6164 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6164 and persid=4.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7007 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7007 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7030 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
   compute civstat=2.
Else if hldid=7030 and persid>1.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7047 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7047 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7064 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
   compute hhtype=1.
   compute hhldsize=1.
Else if hldid=7067 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7067 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7103 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7103 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7114 and (persid>1 and persid<7).
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7114 and persid=7.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7283 and persid=1.
   compute relrefp=1.
   compute parntid1=3.
   compute parntid2=4.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7283 and persid=2.
   compute relrefp=5.
   compute parntid1=3.
   compute parntid2=4.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7297 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7297 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=7297 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7302 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7302 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7305 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7305 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7385 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute hhldsize=4.
Else if hldid=7385 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute hhldsize=4.
Else if hldid=7503 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7503 and persid=3.
   compute relrefp=5.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7736 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7736 and persid=4.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if ncoup=1.
   compute partid=-8.
End if.
if persid=1 relrefp=1.
cro ncoup by partid by civstat.
freq relrefp partid.

compute test=0.
if partid=-5 and ncoup=0 test=1.
freq test.

Do if hldid=825 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
Else if hldid=825 and persid>1.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-8.
   compute cphome =0.
   compute singpar=0.
Else if hldid=825 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6056 and persid=1.
   compute relrefp=1.
   compute parntid1=5.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6056 and (persid>1 and persid<5).
   compute relrefp=5.
   compute parntid1=5.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6056 and persid=5.
   compute relrefp=4.
   compute parntid1=7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
Else if hldid=6056 and persid=6.
   compute relrefp=5.
   compute parntid1=5.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
Else if hldid=6056 and persid=7.
   compute relrefp=9.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7042 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=7042 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7074 and persid=1.
   compute relrefp=1.
   compute parntid1=3.
   compute parntid2=4.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7074 and persid=2.
   compute relrefp=5.
   compute parntid1=3.
   compute parntid2=4.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7123 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
   compute civstat=1.
Else if hldid=7123 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
Else if hldid=7672 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=1.
Else if hldid=7672 and persid=2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=1.
   compute civstat=2.
Else if hldid=7672 and persid>2.
   compute relrefp=9.
   compute parntid1=2.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
End if.
if civstat=-8 and ncoup=0 partid=-8.
if civstat=2 partid=-7.
cro ncoup by partid by civstat.
freq relrefp partid.

Do if hldid=810 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=810 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=810 and persid>2.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
Else if hldid=999 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=999 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=999 and persid=3.
   compute relrefp=10.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=-7.
   compute cphome =0.
   compute singpar=0.
   compute civstat=2.
Else if hldid=7262 and persid=1.
   compute relrefp=1.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7262 and persid=2.
   compute relrefp=2.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute cphome =0.
   compute singpar=0.
Else if hldid=7262 and persid=3.
   compute relrefp=3.
   compute parntid1=1.
   compute parntid2=2.
   compute partid=-7.
   compute cphome =1.
   compute singpar=0.
   compute civstat=2.
End if.
cro partid by civstat.

temp.
select if partid=-5.
cro mincsex by maxcsex.

compute cpid1=0.
compute cpid2=99.
if partid=-5 cpid1=persid.
if partid=-5 cpid2=persid.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid
  /lowcpid=MIN(cpid2) /hicpid=MAX(cpid1) /ndiary=n.
cro lowcpid by hicpid.

if partid=-5 and persid=lowcpid partid=hicpid.
if partid=-5 and persid=hicpid partid=lowcpid.
freq partid hhtype.
cro partid by civstat.

if partid=1 and relrefp=-5 relrefp=2.
cro relrefp by hhtype persid civstat ndiary.

compute test=0.
if ndiary=1 and relrefp=-5 test=1.
freq test.
cro civstat by persid by test.
Do if ndiary=1 and civstat=1 and persid=2 and relrefp=-5.
   compute relrefp=2.
   compute partid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
End if.
Do if hldid=823.
   compute relrefp=3.
   compute partid=-7.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
Else if hldid=1422.
   compute relrefp=3.
   compute partid=-7.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
Else if hldid=1688.
   compute relrefp=2.
   compute partid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
End if.

compute test=0.
if (age>-8 and age<18) and (civstat=1 or ndiary=2) and relrefp=-5 test=1.
freq test.

Do if (hldid=2 or hldid=116 or hldid=144 or hldid=156 or hldid=183 or hldid=259 or hldid=490
   or hldid=529 or hldid=604 or hldid=694 or hldid=1207 or hldid=1576 or hldid=1711 or hldid=6053
   or hldid=6134 or hldid=7317 or hldid=7332 or hldid=7405 or hldid=7493 or hldid=7546 or 
   hldid=7629) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=2 or hldid=116 or hldid=144 or hldid=156 or hldid=183 or hldid=259 or hldid=490
   or hldid=529 or hldid=604 or hldid=694 or hldid=1207 or hldid=1576 or hldid=1711 or hldid=6053
   or hldid=6134 or hldid=7317 or hldid=7332 or hldid=7405 or hldid=7493 or hldid=7546 or 
   hldid=7629) and persid=2.
   compute relrefp=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
Else if (hldid=502 or hldid=1384) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=502 or hldid=1384) and persid=2.
   compute relrefp=5.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1149 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
Else if hldid=1149 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=6112 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute relrefp=2.
Else if (hldid=6112 or hldid=7516) and (persid=3 or persid=4).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
   compute partid=-7.
   compute civstat=2.
End if.
freq relrefp.

temp.
select if hhtype=3 and relrefp=-5 and (age>-8 and age<18).
freq persid ndiary civstat.

compute test=0.
if hhtype=3 and relrefp=-5 and (age>-8 and age<18) and persid=2 test=1.
freq test.

Do if hldid=454 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=454 and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=688 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
Else if hldid=688 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
   compute relrefp=5.
Else if hldid=688 and (persid=4 or persid=5).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1261 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
Else if hldid=1261 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=1261 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1712 and (persid=1 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1712 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=3.
   compute relrefp=3.
End if.

Do if hldid=6014 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=5.
Else if hldid=6014 and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=5.
   compute relrefp=5.
Else if hldid=6014 and (persid=3 or persid=5).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=6068 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=6068 and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=3.
   compute relrefp=3.
Else if hldid=6068 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7510 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
Else if hldid=7510 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=7510 and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7553 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
Else if hldid=7553 and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
   compute relrefp=5.
Else if hldid=7624 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
   compute relrefp=1.
Else if hldid=7624 and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7651 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=5.
Else if hldid=7651 and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=5.
   compute relrefp=5.
End if.
compute test=0.
if civstat=1 and relrefp=-5 and (ndiary=3) test=1.
freq test.
cro test by persid.

Do if (hldid=0 or hldid=148 or hldid=268 or hldid=418 or hldid=474 or hldid=632 
   or hldid=1158 or hldid=1171 or hldid=1208 or hldid=1234 or hldid=1337 or hldid=1510 
   or hldid=1518 or hldid=6155 or hldid=7062 or hldid=7714) and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
Else if (hldid=0 or hldid=148 or hldid=268 or hldid=418 or hldid=474 or hldid=632 
   or hldid=1158 or hldid=1171 or hldid=1208 or hldid=1234 or hldid=1337 or hldid=1510 
   or hldid=1518 or hldid=6155 or hldid=7062 or hldid=7714) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=161 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=161 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=7708 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7708 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7708 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=6.
End if.
compute test=0.
if civstat=1 and relrefp=-5 and (ndiary=2 or ndiary=6) test=1.
freq test.
cro test by persid.

do if hldid=1509 or hldid=7119 or (hldid=7237 and (persid=2 or persid=6)).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7237 and persid=1.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=6.
Else if hldid=7237 and persid=5.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=6.
   compute relrefp=5.
Else if hldid=7237 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7237 and persid=4.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=5.
   compute parntid2=-7.
   compute relrefp=9.
End if.
compute test=0.
if civstat=1 and relrefp=-5 and (ndiary=4 or ndiary=5) test=1.
freq test.
cro test by persid.

Do if hldid=149 and persid=1.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=149 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=149 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=6.
Else if hldid=149 and persid>3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=425 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
Else if hldid=425 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=4.
   compute parntid2=5.
   compute relrefp=5.
Else if hldid=425 and (persid=4 or persid=5).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=749 or hldid=1126) and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
Else if (hldid=749 or hldid=1126) and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
   compute relrefp=5.
Else if (hldid=749 or hldid=1126) and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1256 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
Else if hldid=1256 and persid=4.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
   compute relrefp=5.
Else if hldid=1256 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1261 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
Else if hldid=1261 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=1261 and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1326 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
Else if hldid=1326 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=1326 and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1407 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
Else if hldid=1407 and (persid=4 or persid=5).
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
   compute relrefp=5.
Else if hldid=1407 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=1465 or hldid=1547 or hldid=7015 or hldid=7290) and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
Else if (hldid=1465 or hldid=1547 or hldid=7015 or hldid=7290) and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=4.
   compute relrefp=5.
Else if (hldid=1465 or hldid=1547 or hldid=7015 or hldid=7290) and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
End if.

Do if hldid=1582 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
Else if hldid=1582 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=4.
   compute relrefp=5.
Else if hldid=1582 and (persid=2 or persid=4).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=6018 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
Else if hldid=6018 and persid=4.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=3.
   compute relrefp=5.
Else if hldid=6018 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7314 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7314 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=7314 and persid=4.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=7607 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=2.
   compute parntid2=3.
Else if hldid=7607 and (persid=2 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7607 and persid=4.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=2.
   compute parntid2=3.
   compute relrefp=5.
End if.

temp.
select if relrefp=-5.
freq age civstat hhtype.

compute test=0.
if relrefp=-5 and age>59 test=1.
freq test.

Do if hldid=109.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   if persid=3 relrefp=5.
Else If (hldid=204 or hldid=387 or hldid=413 or hldid=527) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=2.
   compute parntid2=-7.
Else If ((hldid=204 or hldid=387 or hldid=413 or hldid=527) and persid=2) or 
   ((hldid=279 or hldid=608) and persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=279 or hldid=608) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=3.
   compute parntid2=-7.
Else if hldid=279 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=558 or hldid=1373 or hldid=1446 or hldid=1630 
   or hldid=7398 or hldid=7757) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=558 or hldid=608 or hldid=1373 or hldid=1446 or 
   hldid=1630 or hldid=7398 or hldid=7757) and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=5.
Else if (hldid=1106 or hldid=1276 or hldid=1642 or hldid=6057 or hldid=6144 or 
   hldid=7275 or hldid=7424) and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if (hldid=1106 or hldid=1276 or hldid=1642 or hldid=6057 or hldid=6144 or 
   hldid=7275 or hldid=7424) and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=6162 or hldid=7570) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=6162 or hldid=7570) and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=7638 and persid=1. 
   compute singpar=0.
   compute cphome=0.
   compute parntid1=5.
   compute parntid2=-7.
Else if hldid=7638 and persid=2. 
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7638 and (persid=3 or persid=4).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
Else if hldid=7638 and persid=5. 
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
End if.

Do if hldid=757 and persid<3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=757 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=759 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
Else if hldid=759 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=759 and persid=3.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=4.
End if.

temp.
select if relrefp=-5.
freq age ndiary.

compute test=0.
if relrefp=-5 and ndiary=6 test=1.
freq test.

Do if (hldid=154 or hldid=227 or hldid=1218 or hldid=1405 or hldid=1655 or hldid=7643) and persid<3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=154 or hldid=227 or hldid=1218 or hldid=1405 or hldid=1655 or hldid=7643) and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
End if.
Do if hldid=522 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute hhtype=3.
   compute civstat=1.
   compute partid=2.
Else if hldid=522 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute hhtype=3.
   compute civstat=1.
   compute partid=1.
   compute relrefp=2.
Else if hldid=522 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
   compute hhtype=3.
   compute civstat=2.
   compute partid=-7.
   compute relrefp=3.
Else if hldid=7299 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7299 and persid=2.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7299 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=9.
End if.

Do if hldid=164 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=3.
Else if hldid=164 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=3.
   compute partid=-7.
   compute relrefp=3.
End if.

temp.
select if relrefp=-5.
cro ndiary by hhtype by persid.

compute test=0.
if relrefp=-5 and hhtype=4 test=1.
freq test.

Do if (hldid=5 or hldid=22 or hldid=36 or hldid=146 or hldid=205 or hldid=218 or
   hldid=467 or hldid=634 or hldid=668 or hldid=685 or hldid=793) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=5 or hldid=22 or hldid=36 or hldid=146 or hldid=205 or hldid=218 or
   hldid=467 or hldid=634 or hldid=668 or hldid=685 or hldid=793) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=54 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=54 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=54 and persid>2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
else if hldid=193 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
else if hldid=193 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
   compute relrefp=5.
else if hldid=193 and persid=3.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
else if hldid=195 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=6.
   compute parntid2=-7.
else if hldid=195 and persid<6.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=6.
   compute parntid2=-7.
   compute relrefp=5.
else if hldid=195 and persid<6.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
else if hldid=607 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
else if hldid=607 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=614 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=614 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=614 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid<900 and test=1 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid<900 and test=1 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
End if.

Do if (hldid=1021 or hldid=1123 or hldid=1160 or hldid=1220) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=1021 or hldid=1123 or hldid=1160) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=1220 and persid=2.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=1220 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=1080 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1080 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=1125 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=5.
   compute parntid2=-7.
Else if hldid=1125 and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=5.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=1137 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=1137 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1137 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid<1233 and test=1 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid<1233 and test=1 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
End if.

Do if hldid=1233 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=1233 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1246 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1246 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=1247 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1247 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if (hldid=1259 or hldid=1304)and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=1259 or hldid=1304) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if (hldid=1291 or hldid=1298) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1291 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=1298 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
End if.
Execute.

Do if hldid=1309 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1309 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=1317 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1317 and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if (hldid=1342 or hldid=1347 or hldid=1372 or hldid=1475 or hldid=1514 or hldid=1519 
   or hldid=1563 or hldid=1608) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=1342 or hldid=1347 or hldid=1372 or hldid=1475 or hldid=1514 or hldid=1519 
   or hldid=1563 or hldid=1608) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=1479 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=1479 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1479 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=1506 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1506 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if (hldid=1625 or hldid=1641) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=1625 or hldid=1641) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
End if.
Execute.

Do if (hldid=1647 or hldid=1657 or hldid=1678 or hldid=1693 or hldid=1696 or hldid=1722
   or hldid=5004 or hldid=6027 or hldid=6046 or hldid=6047 or hldid=6048) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=1647 or hldid=1657 or hldid=1678 or hldid=1693 or hldid=1696 or hldid=1722
   or hldid=5004 or hldid=6027 or hldid=6046 or hldid=6047 or hldid=6048) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=1682 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=1682 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=1735 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=2.
   compute civstat=1.
Else if hldid=1735 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute partid=1.
   compute civstat=1.
   compute relrefp=2.
Else if hldid=1735 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
Else if hldid=6050 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=6050 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=6050 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=6072 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=6072 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
End if.
Execute.

Do if hldid=6087 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
Else if hldid=6087 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=6087 and persid=3.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=6138 or hldid=6157 or hldid=6173 or hldid=6182 or hldid=7000
   or hldid=7048 or hldid=7056 or hldid=7078 or hldid=7097 or hldid=7131 
   or hldid=7151 or hldid=7170 or hldid=7190 or hldid=7207 or hldid=7222) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=6138 or hldid=6157 or hldid=6173 or hldid=6182 or hldid=7000
   or hldid=7048 or hldid=7056 or hldid=7078 or hldid=7097 or hldid=7131
   or hldid=7151 or hldid=7170 or hldid=7190 or hldid=7207 or hldid=7222) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if (hldid=6143 or hldid=7069 or hldid=7084 or hldid=7171) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=6143 or hldid=7069 or hldid=7084 or hldid=7171) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
End if.
Execute.

Do if hldid=6087 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
Else if hldid=6087 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
   compute relrefp=5.
Else if hldid=6087 and persid=3.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if (hldid=7235 or hldid=7340 or hldid=7354 or hldid=7357 or hldid=7364 
   or hldid=7376 or hldid=7401 or hldid=7417 or hldid=7495 or hldid=7507) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=7235 or hldid=7340 or hldid=7354 or hldid=7357 or hldid=7364
   or hldid=7376 or hldid=7401 or hldid=7417 or hldid=7495 or hldid=7507) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=7250 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7250 and persid=2.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=7250 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7254 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7254 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7322 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=7322 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7352 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7352 and persid=2.
   compute singpar=1.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7352 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=7430 and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7430 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=7430 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=10.
Else if hldid=7436 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7436 and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
End if.
Execute.

Do if hldid=7548 and persid=1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
Else if hldid=7548 and persid=2.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
Else if hldid=7548 and persid=3.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=2.
   compute parntid2=-7.
   compute relrefp=5.
Else if (hldid=7578 or hldid=7630 or hldid=7640 or hldid=7647 or hldid=7723
   or hldid=7729) and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=7578 or hldid=7630 or hldid=7640 or hldid=7647 or hldid=7723
   or hldid=7729) and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if (hldid=7622 or hldid=7742) and persid=1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if (hldid=7622 or hldid=7742) and persid>1.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=10.
End if.
Do if hldid=7356 and persid=1.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7356 and persid>1.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=-7.
   compute relrefp=3.
Else if hldid=7733 and persid=2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=3.
   compute parntid2=-7.
   compute relrefp=1.
Else if hldid=7733 and persid=3.
   compute singpar=1.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=4.
End if.
Execute.

compute test=0.
if relrefp=-5 and hhtype=4 test=1.
freq test.

compute test=0.
if relrefp=-5 and age<26 and persid>2 test=1.
freq test.

Do if hldid=289 and (persid=1 or persid=4).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=289.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=4.
   compute relrefp=3.
Else if hldid=1549 and (persid=1 or persid=5).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1549.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=5.
   compute relrefp=3.
Else if hldid=1656 and (persid=1 or persid=5).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=1656.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=5.
   compute relrefp=3.
Else if hldid=7748 and (persid=1 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=7748.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=3.
   compute relrefp=3.
Else if relrefp=-5 and age<26 and persid>2.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
End if.

compute test=0.
if relrefp=-5 test=1.
freq test.

Do if (hldid=451 or hldid=816 or hldid=7741 or hldid=7748) and (persid=1 or persid=3).
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
Else if hldid=451 or hldid=816 or hldid=7741 or hldid=7748.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=3.
   compute relrefp=3.
Else if (hldid=686 or hldid=1254 or hldid=1441) and persid=4.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if hldid=1606 and persid=3.
   compute singpar=0.
   compute cphome=0.
   compute parntid1=-7.
   compute parntid2=-7.
   compute relrefp=9.
Else if relrefp=-5.
   compute singpar=0.
   compute cphome=1.
   compute parntid1=1.
   compute parntid2=2.
   compute relrefp=3.
End if.
freq relrefp.

recode cphome singpar (-5=0).
recode parntid1 parntid2 (-5=-7).
cro relrefp by parntid1.
Do if (age>-8 and age<18) and parntid1>0.
   compute famstat=4.
Else if age>-8 and age<18.
   compute famstat=5.
End if.
freq famstat.

cro famstat by nchild agekidx.

recode wkhrs (-9=-7) (else=copy) into workhrs.
means workhrs by fulltime parttime
   /cells=min max.

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.

*cohab, unemp, retired, disab & student 
*already have the values/names for the MTUS.
recode disab (-3=-9) (else=copy).
recode unemp student retired (-3=0) (else=copy).
means age by retired by empstat
   /cells min max mean.
if age>69 retired=1.
freq empinclm.

cro empstat by empinclm.

temp.
select if empinclm<0.
cro empstat by empinclm.

temp.
select if age>-8 and age<18.
freq empinclm.

if empinclm=-7 and empstat<4 and age>17 empinclm=-8.

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:\Users\Kimberly Fisher\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.
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.

cro civstat by cphome.
if civstat=1 and cphome=1 cphome=0.
compute test=0.
if civstat<1 and cphome=1 test=1.
freq test.
if test=1 civstat=2.
cro civstat by cphome.

sort cases by hldid persid.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85ind.sav'.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\spemp.sav'
   /keep hldid partid empstat civstat.

get file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\spemp.sav'.
freq partid civstat.
select if civstat=1 and partid>0.
freq partid civstat.
rename variables (partid empstat=persid empsp2).
sort cases by hldid persid.
save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\spemp.sav'
   /drop civstat.

match files
      file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85ind.sav'
      file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\spemp.sav'
   /BY hldid persid.
recode empsp2 (sysmis=-7).
recode sex (sysmis=-5).
freq sex.
select if sex>0.
freq sex.
if civstat=1 and empsp2=-7 empsp2=-8.
cro empsp by empsp2.
if empsp<0 and empsp2>-8 empsp=empsp2.
cro empsp by empsp2 civstat.
*3309 (99%) cases match up exactly. In 1 case, the spouse did not complete a diary, but the respondent reported the
*employment status of the spouse. In the other 29 cases, the respondent did not report the employment status 
*of her or his spouse, but the spouses reported their own employment status.

compute coup=0.
if civstat=1 coup=1.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid
  /ncoup2=SUM(coup).
cro ncoup2 by hhtype.
if ncoup2>0 and hhtype=4 hhtype=3.

temp.
select if ncoup2>2.
cro hldid by ncoup2.

*amend agekidx and agekid2.
compute kidage=999.
compute kidp=0.
if relrefp=3 kidage=age.
if relrefp=3 kidp=1.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid
  /minkidage=MIN(kidage) /kidpres=MAX(kidp).
cro agekidx by minkidage kidpres.

Do if agekidx=-7 and (minkidage>0 and minkidage<999).
   compute agekidx=4.
   compute agekid2=minkidage.
End if.
freq agekidx agekid2.
cro agekidx by kidpres.

cro famstat by agekidx nchild.
cro empsp by partid by civstat.
Do if civstat=2 and partid=-8.
   compute partid=-7.
   compute empsp=-7.
End if.
cro empsp by partid by civstat.

cro parntid1 by parntid2.

sort cases by hldid persid.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85ind.sav'
   /keep 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.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85match.sav'
   /keep hldid persid parntid1 parntid2 partid.

*pull out sex and age for the W6 file.
Get File='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Questionnaire\USA85quest.sav'
   /keep hhid pid age sex.
rename variables (hhid pid=hldid persid).
recode age (80 thru hi=80).
freq age.
save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\sexage.sav'.

*time use variables and W6 file - need to combine AHTUS adult and child files.
GET FILE='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Diary\diary1985\USA1985hfep.sav'.
rename variables (hhid pid=hldid persid).

match files files=*
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\sexage.sav'
   /BY hldid persid.
desc age.

ADD FILES /FILE=*
  /FILE='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\AHTUS final files\Harmonised files\Supplementary files\1985\usa1985childep.sav'
  /RENAME (ageyngst civstat cohab disab educ empsp empstat ethnic ethnic2 famstat fulltime 
    hhtype hisp homemakr incomeqt nadult noemploy nwork occup ownhome parttime regionc regione retired 
    state student under18 under5 unemp urban wagelm wkhrs=d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11 d12 
    d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 d25 d26 d27 d28 d29 d30 d31)
  /DROP=d0 d1 d2 d3 d4 d5 d6 d7 d8 d9 d10 d11 d12 d13 d14 d15 d16 d17 d18 d19 d20 d21 d22 d23 d24 
    d25 d26 d27 d28 d29 d30 d31.
desc age.
recode hldid persid (sysmis=-8).

compute countrya=38.
compute survey=1985.
compute swave=0.
compute msamp=0.
Do if hldid=-8.
   compute hldid=hhid.
End if.
Do if persid=-8.
   compute persid=pid.
End if.
execute.
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:\Users\Kimberly Fisher\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.
*4 episodes 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 main 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=0.
if main=89 ict=1.
cro main by ict.

*alone, sppart, child already in MTUS format.
compute oad=-9.

cro baddem by lowqual by day.
*no missing day of week.
compute eat=0.
compute trv=0.
compute pcr=0.
compute slp=0.
If main=8 or main=9 or main=56 or sec=2 or sec=8 or sec=9 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 sec=1 or sec=6 or sec=7 pcr=1.
If main=62 or main=63 or (main>89 and main<99) or sec=62 or sec=63 or (sec>89 and sec<99) trv=1.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hhid pid
  /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:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa1985hef.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:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\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=MAX(sppart)  /ocombwt=max(origwght) 
  /propwt= max(recwght).

Match files
      file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa1985hef.sav'
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85match.sav'
   /BY hldid persid.
freq sppart partid.
select if partid>0.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp1.sav'
   /keep hldid persid start epnum.
save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp2.sav'
   /keep hldid partid start.

get file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp1.sav'.
compute partmark=1.
save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp1.sav'.

get file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp2.sav'.
rename variables (partid=persid).
sort cases by hldid persid.
compute partmark=1.
save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp2.sav'.

GET file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\1985\USA85withrate.SAV'.

match files files=*
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp2.sav'
  /BY hldid persid start.
execute.

recode partmark (sysmis=0).
execute.

AGGREGATE
  /OUTFILE=* MODE=ADDVARIABLES
  /BREAK= hldid persid
  /partner=MAX(partmark).
freq partner.

select if partner=1.
execute.

vector actasp(288) locsp(288).
Do repeat a=actasp1 to locsp288.
   compute a=0.
End repeat.

Do if start=0.
   compute actasp1=actcode.
   compute locsp1= where.
Else if start<5.
   compute actasp2=actcode.
   compute locsp2= where.
Else if start<10.
   compute actasp3=actcode.
   compute locsp3= where.
Else if start<15.
   compute actasp4=actcode.
   compute locsp4= where.
Else if start<20.
   compute actasp5=actcode.
   compute locsp5= where.
Else if start<25.
   compute actasp6=actcode.
   compute locsp6= where.
Else if start<30.
   compute actasp7=actcode.
   compute locsp7= where.
Else if start<35.
   compute actasp8=actcode.
   compute locsp8= where.
Else if start<40.
   compute actasp9=actcode.
   compute locsp9= where.
Else if start<45.
   compute actasp10=actcode.
   compute locsp10= where.
Else if start<50.
   compute actasp11=actcode.
   compute locsp11= where.
Else if start<55.
   compute actasp12=actcode.
   compute locsp12= where.
Else if start<60.
   compute actasp13=actcode.
   compute locsp13= where.
Else if start<65.
   compute actasp14=actcode.
   compute locsp14= where.
Else if start<70.
   compute actasp15=actcode.
   compute locsp15= where.
Else if start<75.
   compute actasp16=actcode.
   compute locsp16= where.
Else if start<80.
   compute actasp17=actcode.
   compute locsp17= where.
Else if start<85.
   compute actasp18=actcode.
   compute locsp18= where.
Else if start<90.
   compute actasp19=actcode.
   compute locsp19= where.
Else if start<95.
   compute actasp20=actcode.
   compute locsp20= where.
Else if start<100.
   compute actasp21=actcode.
   compute locsp21= where.
Else if start<105.
   compute actasp22=actcode.
   compute locsp22= where.
Else if start<110.
   compute actasp23=actcode.
   compute locsp23= where.
Else if start<115.
   compute actasp24=actcode.
   compute locsp24= where.
Else if start<120.
   compute actasp25=actcode.
   compute locsp25= where.
Else if start<125.
   compute actasp26=actcode.
   compute locsp26= where.
Else if start<130.
   compute actasp27=actcode.
   compute locsp27= where.
Else if start<135.
   compute actasp28=actcode.
   compute locsp28= where.
Else if start<140.
   compute actasp29=actcode.
   compute locsp29= where.
Else if start<145.
   compute actasp30=actcode.
   compute locsp30= where.
Else if start<150.
   compute actasp30=actcode.
   compute locsp30= where.
Else if start<155.
   compute actasp31=actcode.
   compute locsp31= where.
Else if start<160.
   compute actasp32=actcode.
   compute locsp32= where.
Else if start<165.
   compute actasp33=actcode.
   compute locsp33= where.
Else if start<170.
   compute actasp34=actcode.
   compute locsp34= where.
Else if start<175.
   compute actasp35=actcode.
   compute locsp35= where.
Else if start<180.
   compute actasp36=actcode.
   compute locsp36= where.
Else if start<185.
   compute actasp37=actcode.
   compute locsp37= where.
Else if start<190.
   compute actasp38=actcode.
   compute locsp38= where.
Else if start<195.
   compute actasp39=actcode.
   compute locsp39= where.
Else if start<200.
   compute actasp40=actcode.
   compute locsp40= where.
Else if start<205.
   compute actasp41=actcode.
   compute locsp41= where.
Else if start<210.
   compute actasp42=actcode.
   compute locsp42= where.
Else if start<215.
   compute actasp43=actcode.
   compute locsp43= where.
Else if start<220.
   compute actasp44=actcode.
   compute locsp44= where.
Else if start<225.
   compute actasp45=actcode.
   compute locsp45= where.
Else if start<230.
   compute actasp46=actcode.
   compute locsp46= where.
Else if start<235.
   compute actasp47=actcode.
   compute locsp47= where.
Else if start<240.
   compute actasp48=actcode.
   compute locsp48= where.
Else if start<245.
   compute actasp49=actcode.
   compute locsp49= where.
Else if start<250.
   compute actasp50=actcode.
   compute locsp50= where.
Else if start<255.
   compute actasp51=actcode.
   compute locsp51= where.
Else if start<260.
   compute actasp52=actcode.
   compute locsp52= where.
Else if start<265.
   compute actasp53=actcode.
   compute locsp53= where.
Else if start<270.
   compute actasp54=actcode.
   compute locsp54= where.
Else if start<275.
   compute actasp55=actcode.
   compute locsp55= where.
Else if start<280.
   compute actasp56=actcode.
   compute locsp56= where.
Else if start<285.
   compute actasp57=actcode.
   compute locsp57= where.
Else if start<290.
   compute actasp58=actcode.
   compute locsp58= where.
Else if start<295.
   compute actasp59=actcode.
   compute locsp59= where.
Else if start<300.
   compute actasp60=actcode.
   compute locsp60= where.
Else if start<305.
   compute actasp61=actcode.
   compute locsp61= where.
Else if start<310.
   compute actasp62=actcode.
   compute locsp62= where.
Else if start<315.
   compute actasp63=actcode.
   compute locsp63= where.
Else if start<320.
   compute actasp64=actcode.
   compute  locsp64=where.
Else if start<325.
   compute actasp65=actcode.
   compute  locsp65=where.
Else if start<330.
   compute actasp66=actcode.
   compute  locsp66=where.
Else if start<335.
   compute actasp67=actcode.
   compute  locsp67=where.
Else if start<340.
   compute actasp68=actcode.
   compute  locsp68=where.
Else if start<345.
   compute actasp69=actcode.
   compute  locsp69=where.
Else if start<350.
   compute actasp70=actcode.
   compute  locsp70=where.
Else if start<355.
   compute actasp71=actcode.
   compute  locsp71=where.
Else if start<360.
   compute actasp72=actcode.
   compute  locsp72=where.
Else if start<365.
   compute actasp73=actcode.
   compute  locsp73=where.
Else if start<370.
   compute actasp74=actcode.
   compute  locsp74=where.
Else if start<375.
   compute actasp75=actcode.
   compute  locsp75=where.
Else if start<380.
   compute actasp76=actcode.
   compute  locsp76=where.
Else if start<385.
   compute actasp77=actcode.
   compute  locsp77=where.
Else if start<390.
   compute actasp78=actcode.
   compute  locsp78=where.
Else if start<395.
   compute actasp79=actcode.
   compute  locsp79=where.
Else if start<400.
   compute actasp80=actcode.
   compute  locsp80=where.
Else if start<405.
   compute actasp81=actcode.
   compute  locsp81=where.
Else if start<410.
   compute actasp82=actcode.
   compute  locsp82=where.
Else if start<415.
   compute actasp83=actcode.
   compute  locsp83=where.
Else if start<420.
   compute actasp84=actcode.
   compute  locsp84=where.
Else if start<425.
   compute actasp85=actcode.
   compute  locsp85=where.
Else if start<430.
   compute actasp86=actcode.
   compute  locsp86=where.
Else if start<435.
   compute actasp87=actcode.
   compute  locsp87=where.
Else if start<440.
   compute actasp88=actcode.
   compute  locsp88=where.
Else if start<445.
   compute actasp89=actcode.
   compute  locsp89=where.
Else if start<450.
   compute actasp90=actcode.
   compute  locsp90=where.
Else if start<455.
   compute actasp91=actcode.
   compute  locsp91=where.
Else if start<460.
   compute actasp92=actcode.
   compute  locsp92=where.
Else if start<465.
   compute actasp93=actcode.
   compute  locsp93=where.
Else if start<470.
   compute actasp94=actcode.
   compute  locsp94=where.
Else if start<475.
   compute actasp95=actcode.
   compute  locsp95=where.
Else if start<480.
   compute actasp96=actcode.
   compute  locsp96=where.
Else if start<485.
   compute actasp97=actcode.
   compute  locsp97=where.
Else if start<490.
   compute actasp98=actcode.
   compute  locsp98=where.
Else if start<495.
   compute actasp99=actcode.
   compute  locsp99=where.
Else if start<500.
   compute actasp100=actcode.
   compute  locsp100=where.
Else if start<505.
   compute actasp101=actcode.
   compute  locsp101=where.
Else if start<510.
   compute actasp102=actcode.
   compute  locsp102=where.
Else if start<515.
   compute actasp103=actcode.
   compute  locsp103=where.
Else if start<520.
   compute actasp104=actcode.
   compute  locsp104=where.
Else if start<525.
   compute actasp105=actcode.
   compute  locsp105=where.
Else if start<530.
   compute actasp106=actcode.
   compute  locsp106=where.
Else if start<535.
   compute actasp107=actcode.
   compute  locsp107=where.
Else if start<540.
   compute actasp108=actcode.
   compute  locsp108=where.
Else if start<545.
   compute actasp109=actcode.
   compute  locsp109=where.
Else if start<550.
   compute actasp110=actcode.
   compute  locsp110=where.
Else if start<555.
   compute actasp111=actcode.
   compute  locsp111=where.
Else if start<560.
   compute actasp112=actcode.
   compute  locsp112=where.
Else if start<565.
   compute actasp113=actcode.
   compute  locsp113=where.
Else if start<570.
   compute actasp114=actcode.
   compute  locsp114=where.
Else if start<575.
   compute actasp115=actcode.
   compute  locsp115=where.
Else if start<580.
   compute actasp116=actcode.
   compute  locsp116=where.
Else if start<585.
   compute actasp117=actcode.
   compute  locsp117=where.
Else if start<590.
   compute actasp118=actcode.
   compute  locsp118=where.
Else if start<595.
   compute actasp119=actcode.
   compute  locsp119=where.
Else if start<600.
   compute actasp120=actcode.
   compute  locsp120=where.
Else if start<605.
   compute actasp121=actcode.
   compute  locsp121=where.
Else if start<610.
   compute actasp122=actcode.
   compute  locsp122=where.
Else if start<615.
   compute actasp123=actcode.
   compute  locsp123=where.
Else if start<620.
   compute actasp124=actcode.
   compute  locsp124=where.
Else if start<625.
   compute actasp125=actcode.
   compute  locsp125=where.
Else if start<630.
   compute actasp126=actcode.
   compute  locsp126=where.
Else if start<635.
   compute actasp127=actcode.
   compute  locsp127=where.
Else if start<640.
   compute actasp128=actcode.
   compute  locsp128=where.
Else if start<645.
   compute actasp129=actcode.
   compute  locsp129=where.
Else if start<650.
   compute actasp130=actcode.
   compute  locsp130=where.
Else if start<655.
   compute actasp131=actcode.
   compute  locsp131=where.
Else if start<660.
   compute actasp132=actcode.
   compute  locsp132=where.
Else if start<665.
   compute actasp133=actcode.
   compute  locsp133=where.
Else if start<670.
   compute actasp134=actcode.
   compute  locsp134=where.
Else if start<675.
   compute actasp135=actcode.
   compute  locsp135=where.
Else if start<680.
   compute actasp136=actcode.
   compute  locsp136=where.
Else if start<685.
   compute actasp137=actcode.
   compute  locsp137=where.
Else if start<690.
   compute actasp138=actcode.
   compute  locsp138=where.
Else if start<695.
   compute actasp139=actcode.
   compute  locsp139=where.
Else if start<700.
   compute actasp140=actcode.
   compute  locsp140=where.
Else if start<705.
   compute actasp141=actcode.
   compute  locsp141=where.
Else if start<710.
   compute actasp142=actcode.
   compute  locsp142=where.
Else if start<715.
   compute actasp143=actcode.
   compute  locsp143=where.
Else if start<720.
   compute actasp144=actcode.
   compute  locsp144=where.
Else if start<725.
   compute actasp145=actcode.
   compute  locsp145=where.
Else if start<730.
   compute actasp146=actcode.
   compute  locsp146=where.
Else if start<735.
   compute actasp147=actcode.
   compute  locsp147=where.
Else if start<740.
   compute actasp148=actcode.
   compute  locsp148=where.
Else if start<745.
   compute actasp149=actcode.
   compute  locsp149=where.
Else if start<750.
   compute actasp150=actcode.
   compute  locsp150=where.
Else if start<755.
   compute actasp151=actcode.
   compute  locsp151=where.
Else if start<760.
   compute actasp152=actcode.
   compute  locsp152=where.
Else if start<765.
   compute actasp153=actcode.
   compute  locsp153=where.
Else if start<770.
   compute actasp154=actcode.
   compute  locsp154=where.
Else if start<775.
   compute actasp155=actcode.
   compute  locsp155=where.
Else if start<780.
   compute actasp156=actcode.
   compute  locsp156=where.
Else if start<785.
   compute actasp157=actcode.
   compute  locsp157=where.
Else if start<790.
   compute actasp158=actcode.
   compute  locsp158=where.
Else if start<795.
   compute actasp159=actcode.
   compute  locsp159=where.
Else if start<800.
   compute actasp160=actcode.
   compute  locsp160=where.
Else if start<805.
   compute actasp161=actcode.
   compute  locsp161=where.
Else if start<810.
   compute actasp162=actcode.
   compute  locsp162=where.
Else if start<815.
   compute actasp163=actcode.
   compute  locsp163=where.
Else if start<820.
   compute actasp164=actcode.
   compute  locsp164=where.
Else if start<825.
   compute actasp165=actcode.
   compute  locsp165=where.
Else if start<830.
   compute actasp166=actcode.
   compute  locsp166=where.
Else if start<835.
   compute actasp167=actcode.
   compute  locsp167=where.
Else if start<840.
   compute actasp168=actcode.
   compute  locsp168=where.
Else if start<845.
   compute actasp169=actcode.
   compute  locsp169=where.
Else if start<850.
   compute actasp170=actcode.
   compute  locsp170=where.
Else if start<855.
   compute actasp171=actcode.
   compute  locsp171=where.
Else if start<860.
   compute actasp172=actcode.
   compute  locsp172=where.
Else if start<865.
   compute actasp173=actcode.
   compute  locsp173=where.
Else if start<870.
   compute actasp174=actcode.
   compute  locsp174=where.
Else if start<875.
   compute actasp175=actcode.
   compute  locsp175=where.
Else if start<880.
   compute actasp176=actcode.
   compute  locsp176=where.
Else if start<885.
   compute actasp177=actcode.
   compute  locsp177=where.
Else if start<890.
   compute actasp178=actcode.
   compute  locsp178=where.
Else if start<895.
   compute actasp179=actcode.
   compute  locsp179=where.
Else if start<900.
   compute actasp180=actcode.
   compute  locsp180=where.
Else if start<905.
   compute actasp181=actcode.
   compute  locsp181=where.
Else if start<910.
   compute actasp182=actcode.
   compute  locsp182=where.
Else if start<915.
   compute actasp183=actcode.
   compute  locsp183=where.
Else if start<920.
   compute actasp184=actcode.
   compute  locsp184=where.
Else if start<925.
   compute actasp185=actcode.
   compute  locsp185=where.
Else if start<930.
   compute actasp186=actcode.
   compute  locsp186=where.
Else if start<935.
   compute actasp187=actcode.
   compute  locsp187=where.
Else if start<940.
   compute actasp188=actcode.
   compute  locsp188=where.
Else if start<945.
   compute actasp189=actcode.
   compute  locsp189=where.
Else if start<950.
   compute actasp190=actcode.
   compute  locsp190=where.
Else if start<955.
   compute actasp191=actcode.
   compute  locsp191=where.
Else if start<960.
   compute actasp192=actcode.
   compute  locsp192=where.
Else if start<965.
   compute actasp193=actcode.
   compute  locsp193=where.
Else if start<970.
   compute actasp194=actcode.
   compute  locsp194=where.
Else if start<975.
   compute actasp195=actcode.
   compute  locsp195=where.
Else if start<980.
   compute actasp196=actcode.
   compute  locsp196=where.
Else if start<985.
   compute actasp197=actcode.
   compute  locsp197=where.
Else if start<990.
   compute actasp198=actcode.
   compute  locsp198=where.
Else if start<995.
   compute actasp199=actcode.
   compute  locsp199=where.
Else if start<1000.
   compute actasp200=actcode.
   compute  locsp200=where.
Else if start<1005.
   compute actasp201=actcode.
   compute  locsp201=where.
Else if start<1010.
   compute actasp202=actcode.
   compute  locsp202=where.
Else if start<1015.
   compute actasp203=actcode.
   compute  locsp203=where.
Else if start<1020.
   compute actasp204=actcode.
   compute  locsp204=where.
Else if start<1025.
   compute actasp205=actcode.
   compute  locsp205=where.
Else if start<1030.
   compute actasp206=actcode.
   compute  locsp206=where.
Else if start<1035.
   compute actasp207=actcode.
   compute  locsp207=where.
Else if start<1040.
   compute actasp208=actcode.
   compute  locsp208=where.
Else if start<1045.
   compute actasp209=actcode.
   compute  locsp209=where.
Else if start<1050.
   compute actasp210=actcode.
   compute  locsp210=where.
Else if start<1055.
   compute actasp211=actcode.
   compute  locsp211=where.
Else if start<1060.
   compute actasp212=actcode.
   compute  locsp212=where.
Else if start<1065.
   compute actasp213=actcode.
   compute  locsp213=where.
Else if start<1070.
   compute actasp214=actcode.
   compute  locsp214=where.
Else if start<1075.
   compute actasp215=actcode.
   compute  locsp215=where.
Else if start<1080.
   compute actasp216=actcode.
   compute  locsp216=where.
Else if start<1085.
   compute actasp217=actcode.
   compute  locsp217=where.
Else if start<1090.
   compute actasp218=actcode.
   compute  locsp218=where.
Else if start<1095.
   compute actasp219=actcode.
   compute  locsp219=where.
Else if start<1100.
   compute actasp220=actcode.
   compute  locsp220=where.
Else if start<1105.
   compute actasp221=actcode.
   compute  locsp221=where.
Else if start<1110.
   compute actasp222=actcode.
   compute  locsp222=where.
Else if start<1115.
   compute actasp223=actcode.
   compute  locsp223=where.
Else if start<1120.
   compute actasp224=actcode.
   compute  locsp224=where.
Else if start<1125.
   compute actasp225=actcode.
   compute  locsp225=where.
Else if start<1130.
   compute actasp226=actcode.
   compute  locsp226=where.
Else if start<1135.
   compute actasp227=actcode.
   compute  locsp227=where.
Else if start<1140.
   compute actasp228=actcode.
   compute  locsp228=where.
Else if start<1145.
   compute actasp229=actcode.
   compute  locsp229=where.
Else if start<1150.
   compute actasp230=actcode.
   compute  locsp230=where.
Else if start<1155.
   compute actasp231=actcode.
   compute  locsp231=where.
Else if start<1160.
   compute actasp232=actcode.
   compute  locsp232=where.
Else if start<1165.
   compute actasp233=actcode.
   compute  locsp233=where.
Else if start<1170.
   compute actasp234=actcode.
   compute  locsp234=where.
Else if start<1175.
   compute actasp235=actcode.
   compute  locsp235=where.
Else if start<1180.
   compute actasp236=actcode.
   compute  locsp236=where.
Else if start<1185.
   compute actasp237=actcode.
   compute  locsp237=where.
Else if start<1190.
   compute actasp238=actcode.
   compute  locsp238=where.
Else if start<1195.
   compute actasp239=actcode.
   compute  locsp239=where.
Else if start<1200.
   compute actasp240=actcode.
   compute  locsp240=where.
Else if start<1205.
   compute actasp241=actcode.
   compute  locsp241=where.
Else if start<1210.
   compute actasp242=actcode.
   compute  locsp242=where.
Else if start<1215.
   compute actasp243=actcode.
   compute  locsp243=where.
Else if start<1220.
   compute actasp244=actcode.
   compute  locsp244=where.
Else if start<1225.
   compute actasp245=actcode.
   compute  locsp245=where.
Else if start<1230.
   compute actasp246=actcode.
   compute  locsp246=where.
Else if start<1235.
   compute actasp247=actcode.
   compute  locsp247=where.
Else if start<1240.
   compute actasp248=actcode.
   compute  locsp248=where.
Else if start<1245.
   compute actasp249=actcode.
   compute  locsp249=where.
Else if start<1250.
   compute actasp250=actcode.
   compute  locsp250=where.
Else if start<1255.
   compute actasp251=actcode.
   compute  locsp251=where.
Else if start<1260.
   compute actasp252=actcode.
   compute  locsp252=where.
Else if start<1265.
   compute actasp253=actcode.
   compute  locsp253=where.
Else if start<1270.
   compute actasp254=actcode.
   compute  locsp254=where.
Else if start<1275.
   compute actasp255=actcode.
   compute  locsp255=where.
Else if start<1280.
   compute actasp256=actcode.
   compute  locsp256=where.
Else if start<1285.
   compute actasp257=actcode.
   compute  locsp257=where.
Else if start<1290.
   compute actasp258=actcode.
   compute  locsp258=where.
Else if start<1295.
   compute actasp259=actcode.
   compute  locsp259=where.
Else if start<1300.
   compute actasp260=actcode.
   compute  locsp260=where.
Else if start<1305.
   compute actasp261=actcode.
   compute  locsp261=where.
Else if start<1310.
   compute actasp262=actcode.
   compute  locsp262=where.
Else if start<1315.
   compute actasp263=actcode.
   compute  locsp263=where.
Else if start<1320.
   compute actasp264=actcode.
   compute  locsp264=where.
Else if start<1325.
   compute actasp265=actcode.
   compute  locsp265=where.
Else if start<1330.
   compute actasp266=actcode.
   compute  locsp266=where.
Else if start<1335.
   compute actasp267=actcode.
   compute  locsp267=where.
Else if start<1340.
   compute actasp268=actcode.
   compute  locsp268=where.
Else if start<1345.
   compute actasp269=actcode.
   compute  locsp269=where.
Else if start<1350.
   compute actasp270=actcode.
   compute  locsp270=where.
Else if start<1355.
   compute actasp271=actcode.
   compute  locsp271=where.
Else if start<1360.
   compute actasp272=actcode.
   compute  locsp272=where.
Else if start<1365.
   compute actasp273=actcode.
   compute  locsp273=where.
Else if start<1370.
   compute actasp274=actcode.
   compute  locsp274=where.
Else if start<1375.
   compute actasp275=actcode.
   compute  locsp275=where.
Else if start<1380.
   compute actasp276=actcode.
   compute  locsp276=where.
Else if start<1385.
   compute actasp277=actcode.
   compute  locsp277=where.
Else if start<1390.
   compute actasp278=actcode.
   compute  locsp278=where.
Else if start<1395.
   compute actasp279=actcode.
   compute  locsp279=where.
Else if start<1400.
   compute actasp280=actcode.
   compute  locsp280=where.
Else if start<1405.
   compute actasp281=actcode.
   compute  locsp281=where.
Else if start<1410.
   compute actasp282=actcode.
   compute  locsp282=where.
Else if start<1415.
   compute actasp283=actcode.
   compute  locsp283=where.
Else if start<1420.
   compute actasp284=actcode.
   compute  locsp284=where.
Else if start<1425.
   compute actasp285=actcode.
   compute  locsp285=where.
Else if start<1430.
   compute actasp286=actcode.
   compute  locsp286=where.
Else if start<1435.
   compute actasp287=actcode.
   compute  locsp287=where.
Else.
   compute actasp288=actcode.
   compute  locsp288=where.
End if.
execute.

sort cases hldid(a) persid(a) epnum(d).
Do if hldid=lag(hldid) and persid=lag(persid).
   Do repeat a=actasp1 to locsp288.
      if lag(a)>0 a=-5.
   End repeat.
End if.
freq actasp1 locsp1.

sort cases hldid(a) persid(a) epnum(a).

vector actasp=actasp1 to actasp288 /locsp=locsp1 to locsp288.
loop i=1 to 287.
compute j=i+1.
   Do if actasp(i)>0 and actasp(j)=0.
      compute actasp(j)=actasp(i).
   End if.
   Do if locsp(i)>0 and locsp(j)=0.
      compute locsp(j)=locsp(i).
   End if.
End loop.
execute.

Aggregate
   outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\1985\USA85spousevals.SAV'
   /break=hldid persid
   /actasp1=MAX(actasp1) /actasp2=MAX(actasp2) /actasp3=MAX(actasp3) /actasp4=MAX(actasp4) /actasp5=MAX(actasp5) 
   /actasp6=MAX(actasp6) /actasp7=MAX(actasp7) /actasp8=MAX(actasp8) /actasp9=MAX(actasp9) /actasp10=MAX(actasp10) 
   /actasp11=MAX(actasp11) /actasp12=MAX(actasp12) /actasp13=MAX(actasp13) /actasp14=MAX(actasp14) /actasp15=MAX(actasp15) 
   /actasp16=MAX(actasp16) /actasp17=MAX(actasp17) /actasp18=MAX(actasp18) /actasp19=MAX(actasp19) /actasp20=MAX(actasp20) 
   /actasp21=MAX(actasp21) /actasp22=MAX(actasp22) /actasp23=MAX(actasp23) /actasp24=MAX(actasp24) /actasp25=MAX(actasp25) 
   /actasp26=MAX(actasp26) /actasp27=MAX(actasp27) /actasp28=MAX(actasp28) /actasp29=MAX(actasp29) /actasp30=MAX(actasp30) 
   /actasp31=MAX(actasp31) /actasp32=MAX(actasp32) /actasp33=MAX(actasp33) /actasp34=MAX(actasp34) /actasp35=MAX(actasp35) 
   /actasp36=MAX(actasp36) /actasp37=MAX(actasp37) /actasp38=MAX(actasp38) /actasp39=MAX(actasp39) /actasp40=MAX(actasp40) 
   /actasp41=MAX(actasp41) /actasp42=MAX(actasp42) /actasp43=MAX(actasp43) /actasp44=MAX(actasp44) /actasp45=MAX(actasp45) 
   /actasp46=MAX(actasp46) /actasp47=MAX(actasp47) /actasp48=MAX(actasp48) /actasp49=MAX(actasp49) /actasp50=MAX(actasp50) 
   /actasp51=MAX(actasp51) /actasp52=MAX(actasp52) /actasp53=MAX(actasp53) /actasp54=MAX(actasp54) /actasp55=MAX(actasp55) 
   /actasp56=MAX(actasp56) /actasp57=MAX(actasp57) /actasp58=MAX(actasp58) /actasp59=MAX(actasp59) /actasp60=MAX(actasp60) 
   /actasp61=MAX(actasp61) /actasp62=MAX(actasp62) /actasp63=MAX(actasp63) /actasp64=MAX(actasp64) /actasp65=MAX(actasp65) 
   /actasp66=MAX(actasp66) /actasp67=MAX(actasp67) /actasp68=MAX(actasp68) /actasp69=MAX(actasp69) /actasp70=MAX(actasp70) 
   /actasp71=MAX(actasp71) /actasp72=MAX(actasp72) /actasp73=MAX(actasp73) /actasp74=MAX(actasp74) /actasp75=MAX(actasp75) 
   /actasp76=MAX(actasp76) /actasp77=MAX(actasp77) /actasp78=MAX(actasp78) /actasp79=MAX(actasp79) /actasp80=MAX(actasp80) 
   /actasp81=MAX(actasp81) /actasp82=MAX(actasp82) /actasp83=MAX(actasp83) /actasp84=MAX(actasp84) /actasp85=MAX(actasp85) 
   /actasp86=MAX(actasp86) /actasp87=MAX(actasp87) /actasp88=MAX(actasp88) /actasp89=MAX(actasp89) /actasp90=MAX(actasp90) 
   /actasp91=MAX(actasp91) /actasp92=MAX(actasp92) /actasp93=MAX(actasp93) /actasp94=MAX(actasp94) /actasp95=MAX(actasp95) 
   /actasp96=MAX(actasp96) /actasp97=MAX(actasp97) /actasp98=MAX(actasp98) /actasp99=MAX(actasp99) /actasp100=MAX(actasp100) 
   /actasp101=MAX(actasp101) /actasp102=MAX(actasp102) /actasp103=MAX(actasp103) /actasp104=MAX(actasp104) 
   /actasp105=MAX(actasp105) /actasp106=MAX(actasp106) /actasp107=MAX(actasp107) /actasp108=MAX(actasp108) 
   /actasp109=MAX(actasp109) /actasp110=MAX(actasp110) /actasp111=MAX(actasp111) /actasp112=MAX(actasp112) 
   /actasp113=MAX(actasp113) /actasp114=MAX(actasp114) /actasp115=MAX(actasp115) /actasp116=MAX(actasp116) 
   /actasp117=MAX(actasp117) /actasp118=MAX(actasp118) /actasp119=MAX(actasp119) /actasp120=MAX(actasp120) 
   /actasp121=MAX(actasp121) /actasp122=MAX(actasp122) /actasp123=MAX(actasp123) /actasp124=MAX(actasp124) 
   /actasp125=MAX(actasp125) /actasp126=MAX(actasp126) /actasp127=MAX(actasp127) /actasp128=MAX(actasp128) 
   /actasp129=MAX(actasp129) /actasp130=MAX(actasp130) /actasp131=MAX(actasp131) /actasp132=MAX(actasp132) 
   /actasp133=MAX(actasp133) /actasp134=MAX(actasp134) /actasp135=MAX(actasp135) /actasp136=MAX(actasp136) 
   /actasp137=MAX(actasp137) /actasp138=MAX(actasp138) /actasp139=MAX(actasp139) /actasp140=MAX(actasp140) 
   /actasp141=MAX(actasp141) /actasp142=MAX(actasp142) /actasp143=MAX(actasp143) /actasp144=MAX(actasp144) 
   /actasp145=MAX(actasp145) /actasp146=MAX(actasp146) /actasp147=MAX(actasp147) /actasp148=MAX(actasp148) 
   /actasp149=MAX(actasp149) /actasp150=MAX(actasp150) /actasp151=MAX(actasp151) /actasp152=MAX(actasp152) 
   /actasp153=MAX(actasp153) /actasp154=MAX(actasp154) /actasp155=MAX(actasp155) /actasp156=MAX(actasp156) 
   /actasp157=MAX(actasp157) /actasp158=MAX(actasp158) /actasp159=MAX(actasp159) /actasp160=MAX(actasp160)
   /actasp161=MAX(actasp161) /actasp162=MAX(actasp162) /actasp163=MAX(actasp163) /actasp164=MAX(actasp164) 
   /actasp165=MAX(actasp165) /actasp166=MAX(actasp166) /actasp167=MAX(actasp167) /actasp168=MAX(actasp168) 
   /actasp169=MAX(actasp169) /actasp170=MAX(actasp170) /actasp171=MAX(actasp171) /actasp172=MAX(actasp172) 
   /actasp173=MAX(actasp173) /actasp174=MAX(actasp174) /actasp175=MAX(actasp175) /actasp176=MAX(actasp176) 
   /actasp177=MAX(actasp177) /actasp178=MAX(actasp178) /actasp179=MAX(actasp179) /actasp180=MAX(actasp180) 
   /actasp181=MAX(actasp181) /actasp182=MAX(actasp182) /actasp183=MAX(actasp183) /actasp184=MAX(actasp184) 
   /actasp185=MAX(actasp185) /actasp186=MAX(actasp186) /actasp187=MAX(actasp187) /actasp188=MAX(actasp188) 
   /actasp189=MAX(actasp189) /actasp190=MAX(actasp190) /actasp191=MAX(actasp191) /actasp192=MAX(actasp192) 
   /actasp193=MAX(actasp193) /actasp194=MAX(actasp194) /actasp195=MAX(actasp195) /actasp196=MAX(actasp196) 
   /actasp197=MAX(actasp197) /actasp198=MAX(actasp198) /actasp199=MAX(actasp199) /actasp200=MAX(actasp200) 
   /actasp201=MAX(actasp201) /actasp202=MAX(actasp202) /actasp203=MAX(actasp203) /actasp204=MAX(actasp204) 
   /actasp205=MAX(actasp205) /actasp206=MAX(actasp206) /actasp207=MAX(actasp207) /actasp208=MAX(actasp208) 
   /actasp209=MAX(actasp209) /actasp210=MAX(actasp210) /actasp211=MAX(actasp211) /actasp212=MAX(actasp212) 
   /actasp213=MAX(actasp213) /actasp214=MAX(actasp214) /actasp215=MAX(actasp215) /actasp216=MAX(actasp216) 
   /actasp217=MAX(actasp217) /actasp218=MAX(actasp218) /actasp219=MAX(actasp219) /actasp220=MAX(actasp220) 
   /actasp221=MAX(actasp221) /actasp222=MAX(actasp222) /actasp223=MAX(actasp223) /actasp224=MAX(actasp224) 
   /actasp225=MAX(actasp225) /actasp226=MAX(actasp226) /actasp227=MAX(actasp227) /actasp228=MAX(actasp228) 
   /actasp229=MAX(actasp229) /actasp230=MAX(actasp230) /actasp231=MAX(actasp231) /actasp232=MAX(actasp232) 
   /actasp233=MAX(actasp233) /actasp234=MAX(actasp234) /actasp235=MAX(actasp235) /actasp236=MAX(actasp236) 
   /actasp237=MAX(actasp237) /actasp238=MAX(actasp238) /actasp239=MAX(actasp239) /actasp240=MAX(actasp240) 
   /actasp241=MAX(actasp241) /actasp242=MAX(actasp242) /actasp243=MAX(actasp243) /actasp244=MAX(actasp244) 
   /actasp245=MAX(actasp245) /actasp246=MAX(actasp246) /actasp247=MAX(actasp247) /actasp248=MAX(actasp248) 
   /actasp249=MAX(actasp249) /actasp250=MAX(actasp250) /actasp251=MAX(actasp251) /actasp252=MAX(actasp252) 
   /actasp253=MAX(actasp253) /actasp254=MAX(actasp254) /actasp255=MAX(actasp255) /actasp256=MAX(actasp256) 
   /actasp257=MAX(actasp257) /actasp258=MAX(actasp258) /actasp259=MAX(actasp259) /actasp260=MAX(actasp260)
   /actasp261=MAX(actasp261) /actasp262=MAX(actasp262) /actasp263=MAX(actasp263) /actasp264=MAX(actasp264) 
   /actasp265=MAX(actasp265) /actasp266=MAX(actasp266) /actasp267=MAX(actasp267) /actasp268=MAX(actasp268) 
   /actasp269=MAX(actasp269) /actasp270=MAX(actasp270) /actasp271=MAX(actasp271) /actasp272=MAX(actasp272) 
   /actasp273=MAX(actasp273) /actasp274=MAX(actasp274) /actasp275=MAX(actasp275) /actasp276=MAX(actasp276) 
   /actasp277=MAX(actasp277) /actasp278=MAX(actasp278) /actasp279=MAX(actasp279) /actasp280=MAX(actasp280) 
   /actasp281=MAX(actasp281) /actasp282=MAX(actasp282) /actasp283=MAX(actasp283) /actasp284=MAX(actasp284) 
   /actasp285=MAX(actasp285) /actasp286=MAX(actasp286) /actasp287=MAX(actasp287) /actasp288=MAX(actasp288)
   /locsp1=MAX(locsp1) /locsp2=MAX(locsp2) /locsp3=MAX(locsp3) /locsp4=MAX(locsp4) /locsp5=MAX(locsp5) 
   /locsp6=MAX(locsp6) /locsp7=MAX(locsp7) /locsp8=MAX(locsp8) /locsp9=MAX(locsp9) /locsp10=MAX(locsp10) 
   /locsp11=MAX(locsp11) /locsp12=MAX(locsp12) /locsp13=MAX(locsp13) /locsp14=MAX(locsp14) /locsp15=MAX(locsp15) 
   /locsp16=MAX(locsp16) /locsp17=MAX(locsp17) /locsp18=MAX(locsp18) /locsp19=MAX(locsp19) /locsp20=MAX(locsp20) 
   /locsp21=MAX(locsp21) /locsp22=MAX(locsp22) /locsp23=MAX(locsp23) /locsp24=MAX(locsp24) /locsp25=MAX(locsp25) 
   /locsp26=MAX(locsp26) /locsp27=MAX(locsp27) /locsp28=MAX(locsp28) /locsp29=MAX(locsp29) /locsp30=MAX(locsp30) 
   /locsp31=MAX(locsp31) /locsp32=MAX(locsp32) /locsp33=MAX(locsp33) /locsp34=MAX(locsp34) /locsp35=MAX(locsp35) 
   /locsp36=MAX(locsp36) /locsp37=MAX(locsp37) /locsp38=MAX(locsp38) /locsp39=MAX(locsp39) /locsp40=MAX(locsp40) 
   /locsp41=MAX(locsp41) /locsp42=MAX(locsp42) /locsp43=MAX(locsp43) /locsp44=MAX(locsp44) /locsp45=MAX(locsp45) 
   /locsp46=MAX(locsp46) /locsp47=MAX(locsp47) /locsp48=MAX(locsp48) /locsp49=MAX(locsp49) /locsp50=MAX(locsp50) 
   /locsp51=MAX(locsp51) /locsp52=MAX(locsp52) /locsp53=MAX(locsp53) /locsp54=MAX(locsp54) /locsp55=MAX(locsp55) 
   /locsp56=MAX(locsp56) /locsp57=MAX(locsp57) /locsp58=MAX(locsp58) /locsp59=MAX(locsp59) /locsp60=MAX(locsp60) 
   /locsp61=MAX(locsp61) /locsp62=MAX(locsp62) /locsp63=MAX(locsp63) /locsp64=MAX(locsp64) /locsp65=MAX(locsp65) 
   /locsp66=MAX(locsp66) /locsp67=MAX(locsp67) /locsp68=MAX(locsp68) /locsp69=MAX(locsp69) /locsp70=MAX(locsp70) 
   /locsp71=MAX(locsp71) /locsp72=MAX(locsp72) /locsp73=MAX(locsp73) /locsp74=MAX(locsp74) /locsp75=MAX(locsp75) 
   /locsp76=MAX(locsp76) /locsp77=MAX(locsp77) /locsp78=MAX(locsp78) /locsp79=MAX(locsp79) /locsp80=MAX(locsp80) 
   /locsp81=MAX(locsp81) /locsp82=MAX(locsp82) /locsp83=MAX(locsp83) /locsp84=MAX(locsp84) /locsp85=MAX(locsp85) 
   /locsp86=MAX(locsp86) /locsp87=MAX(locsp87) /locsp88=MAX(locsp88) /locsp89=MAX(locsp89) /locsp90=MAX(locsp90) 
   /locsp91=MAX(locsp91) /locsp92=MAX(locsp92) /locsp93=MAX(locsp93) /locsp94=MAX(locsp94) /locsp95=MAX(locsp95) 
   /locsp96=MAX(locsp96) /locsp97=MAX(locsp97) /locsp98=MAX(locsp98) /locsp99=MAX(locsp99) /locsp100=MAX(locsp100) 
   /locsp101=MAX(locsp101) /locsp102=MAX(locsp102) /locsp103=MAX(locsp103) /locsp104=MAX(locsp104) 
   /locsp105=MAX(locsp105) /locsp106=MAX(locsp106) /locsp107=MAX(locsp107) /locsp108=MAX(locsp108) 
   /locsp109=MAX(locsp109) /locsp110=MAX(locsp110) /locsp111=MAX(locsp111) /locsp112=MAX(locsp112) 
   /locsp113=MAX(locsp113) /locsp114=MAX(locsp114) /locsp115=MAX(locsp115) /locsp116=MAX(locsp116) 
   /locsp117=MAX(locsp117) /locsp118=MAX(locsp118) /locsp119=MAX(locsp119) /locsp120=MAX(locsp120) 
   /locsp121=MAX(locsp121) /locsp122=MAX(locsp122) /locsp123=MAX(locsp123) /locsp124=MAX(locsp124) 
   /locsp125=MAX(locsp125) /locsp126=MAX(locsp126) /locsp127=MAX(locsp127) /locsp128=MAX(locsp128) 
   /locsp129=MAX(locsp129) /locsp130=MAX(locsp130) /locsp131=MAX(locsp131) /locsp132=MAX(locsp132) 
   /locsp133=MAX(locsp133) /locsp134=MAX(locsp134) /locsp135=MAX(locsp135) /locsp136=MAX(locsp136) 
   /locsp137=MAX(locsp137) /locsp138=MAX(locsp138) /locsp139=MAX(locsp139) /locsp140=MAX(locsp140) 
   /locsp141=MAX(locsp141) /locsp142=MAX(locsp142) /locsp143=MAX(locsp143) /locsp144=MAX(locsp144) 
   /locsp145=MAX(locsp145) /locsp146=MAX(locsp146) /locsp147=MAX(locsp147) /locsp148=MAX(locsp148) 
   /locsp149=MAX(locsp149) /locsp150=MAX(locsp150) /locsp151=MAX(locsp151) /locsp152=MAX(locsp152) 
   /locsp153=MAX(locsp153) /locsp154=MAX(locsp154) /locsp155=MAX(locsp155) /locsp156=MAX(locsp156) 
   /locsp157=MAX(locsp157) /locsp158=MAX(locsp158) /locsp159=MAX(locsp159) /locsp160=MAX(locsp160)
   /locsp161=MAX(locsp161) /locsp162=MAX(locsp162) /locsp163=MAX(locsp163) /locsp164=MAX(locsp164) 
   /locsp165=MAX(locsp165) /locsp166=MAX(locsp166) /locsp167=MAX(locsp167) /locsp168=MAX(locsp168) 
   /locsp169=MAX(locsp169) /locsp170=MAX(locsp170) /locsp171=MAX(locsp171) /locsp172=MAX(locsp172) 
   /locsp173=MAX(locsp173) /locsp174=MAX(locsp174) /locsp175=MAX(locsp175) /locsp176=MAX(locsp176) 
   /locsp177=MAX(locsp177) /locsp178=MAX(locsp178) /locsp179=MAX(locsp179) /locsp180=MAX(locsp180) 
   /locsp181=MAX(locsp181) /locsp182=MAX(locsp182) /locsp183=MAX(locsp183) /locsp184=MAX(locsp184) 
   /locsp185=MAX(locsp185) /locsp186=MAX(locsp186) /locsp187=MAX(locsp187) /locsp188=MAX(locsp188) 
   /locsp189=MAX(locsp189) /locsp190=MAX(locsp190) /locsp191=MAX(locsp191) /locsp192=MAX(locsp192) 
   /locsp193=MAX(locsp193) /locsp194=MAX(locsp194) /locsp195=MAX(locsp195) /locsp196=MAX(locsp196) 
   /locsp197=MAX(locsp197) /locsp198=MAX(locsp198) /locsp199=MAX(locsp199) /locsp200=MAX(locsp200) 
   /locsp201=MAX(locsp201) /locsp202=MAX(locsp202) /locsp203=MAX(locsp203) /locsp204=MAX(locsp204) 
   /locsp205=MAX(locsp205) /locsp206=MAX(locsp206) /locsp207=MAX(locsp207) /locsp208=MAX(locsp208) 
   /locsp209=MAX(locsp209) /locsp210=MAX(locsp210) /locsp211=MAX(locsp211) /locsp212=MAX(locsp212) 
   /locsp213=MAX(locsp213) /locsp214=MAX(locsp214) /locsp215=MAX(locsp215) /locsp216=MAX(locsp216) 
   /locsp217=MAX(locsp217) /locsp218=MAX(locsp218) /locsp219=MAX(locsp219) /locsp220=MAX(locsp220) 
   /locsp221=MAX(locsp221) /locsp222=MAX(locsp222) /locsp223=MAX(locsp223) /locsp224=MAX(locsp224) 
   /locsp225=MAX(locsp225) /locsp226=MAX(locsp226) /locsp227=MAX(locsp227) /locsp228=MAX(locsp228) 
   /locsp229=MAX(locsp229) /locsp230=MAX(locsp230) /locsp231=MAX(locsp231) /locsp232=MAX(locsp232) 
   /locsp233=MAX(locsp233) /locsp234=MAX(locsp234) /locsp235=MAX(locsp235) /locsp236=MAX(locsp236) 
   /locsp237=MAX(locsp237) /locsp238=MAX(locsp238) /locsp239=MAX(locsp239) /locsp240=MAX(locsp240) 
   /locsp241=MAX(locsp241) /locsp242=MAX(locsp242) /locsp243=MAX(locsp243) /locsp244=MAX(locsp244) 
   /locsp245=MAX(locsp245) /locsp246=MAX(locsp246) /locsp247=MAX(locsp247) /locsp248=MAX(locsp248) 
   /locsp249=MAX(locsp249) /locsp250=MAX(locsp250) /locsp251=MAX(locsp251) /locsp252=MAX(locsp252) 
   /locsp253=MAX(locsp253) /locsp254=MAX(locsp254) /locsp255=MAX(locsp255) /locsp256=MAX(locsp256) 
   /locsp257=MAX(locsp257) /locsp258=MAX(locsp258) /locsp259=MAX(locsp259) /locsp260=MAX(locsp260)
   /locsp261=MAX(locsp261) /locsp262=MAX(locsp262) /locsp263=MAX(locsp263) /locsp264=MAX(locsp264) 
   /locsp265=MAX(locsp265) /locsp266=MAX(locsp266) /locsp267=MAX(locsp267) /locsp268=MAX(locsp268) 
   /locsp269=MAX(locsp269) /locsp270=MAX(locsp270) /locsp271=MAX(locsp271) /locsp272=MAX(locsp272) 
   /locsp273=MAX(locsp273) /locsp274=MAX(locsp274) /locsp275=MAX(locsp275) /locsp276=MAX(locsp276) 
   /locsp277=MAX(locsp277) /locsp278=MAX(locsp278) /locsp279=MAX(locsp279) /locsp280=MAX(locsp280) 
   /locsp281=MAX(locsp281) /locsp282=MAX(locsp282) /locsp283=MAX(locsp283) /locsp284=MAX(locsp284) 
   /locsp285=MAX(locsp285) /locsp286=MAX(locsp286) /locsp287=MAX(locsp287) /locsp288=MAX(locsp288).

GET file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\1985\USA85withrate.SAV'.

match files files=*
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp1.sav'
  /BY hldid persid start.
execute.

select if partmark=1.
execute.

vector acta(288) loc(288).
Do repeat a=acta1 to loc288.
   compute a=0.
End repeat.

Do if start=0.
   compute acta1=actcode.
   compute loc1= where.
Else if start<5.
   compute acta2=actcode.
   compute loc2= where.
Else if start<10.
   compute acta3=actcode.
   compute loc3= where.
Else if start<15.
   compute acta4=actcode.
   compute loc4= where.
Else if start<20.
   compute acta5=actcode.
   compute loc5= where.
Else if start<25.
   compute acta6=actcode.
   compute loc6= where.
Else if start<30.
   compute acta7=actcode.
   compute loc7= where.
Else if start<35.
   compute acta8=actcode.
   compute loc8= where.
Else if start<40.
   compute acta9=actcode.
   compute loc9= where.
Else if start<45.
   compute acta10=actcode.
   compute loc10= where.
Else if start<50.
   compute acta11=actcode.
   compute loc11= where.
Else if start<55.
   compute acta12=actcode.
   compute loc12= where.
Else if start<60.
   compute acta13=actcode.
   compute loc13= where.
Else if start<65.
   compute acta14=actcode.
   compute loc14= where.
Else if start<70.
   compute acta15=actcode.
   compute loc15= where.
Else if start<75.
   compute acta16=actcode.
   compute loc16= where.
Else if start<80.
   compute acta17=actcode.
   compute loc17= where.
Else if start<85.
   compute acta18=actcode.
   compute loc18= where.
Else if start<90.
   compute acta19=actcode.
   compute loc19= where.
Else if start<95.
   compute acta20=actcode.
   compute loc20= where.
Else if start<100.
   compute acta21=actcode.
   compute loc21= where.
Else if start<105.
   compute acta22=actcode.
   compute loc22= where.
Else if start<110.
   compute acta23=actcode.
   compute loc23= where.
Else if start<115.
   compute acta24=actcode.
   compute loc24= where.
Else if start<120.
   compute acta25=actcode.
   compute loc25= where.
Else if start<125.
   compute acta26=actcode.
   compute loc26= where.
Else if start<130.
   compute acta27=actcode.
   compute loc27= where.
Else if start<135.
   compute acta28=actcode.
   compute loc28= where.
Else if start<140.
   compute acta29=actcode.
   compute loc29= where.
Else if start<145.
   compute acta30=actcode.
   compute loc30= where.
Else if start<150.
   compute acta30=actcode.
   compute loc30= where.
Else if start<155.
   compute acta31=actcode.
   compute loc31= where.
Else if start<160.
   compute acta32=actcode.
   compute loc32= where.
Else if start<165.
   compute acta33=actcode.
   compute loc33= where.
Else if start<170.
   compute acta34=actcode.
   compute loc34= where.
Else if start<175.
   compute acta35=actcode.
   compute loc35= where.
Else if start<180.
   compute acta36=actcode.
   compute loc36= where.
Else if start<185.
   compute acta37=actcode.
   compute loc37= where.
Else if start<190.
   compute acta38=actcode.
   compute loc38= where.
Else if start<195.
   compute acta39=actcode.
   compute loc39= where.
Else if start<200.
   compute acta40=actcode.
   compute loc40= where.
Else if start<205.
   compute acta41=actcode.
   compute loc41= where.
Else if start<210.
   compute acta42=actcode.
   compute loc42= where.
Else if start<215.
   compute acta43=actcode.
   compute loc43= where.
Else if start<220.
   compute acta44=actcode.
   compute loc44= where.
Else if start<225.
   compute acta45=actcode.
   compute loc45= where.
Else if start<230.
   compute acta46=actcode.
   compute loc46= where.
Else if start<235.
   compute acta47=actcode.
   compute loc47= where.
Else if start<240.
   compute acta48=actcode.
   compute loc48= where.
Else if start<245.
   compute acta49=actcode.
   compute loc49= where.
Else if start<250.
   compute acta50=actcode.
   compute loc50= where.
Else if start<255.
   compute acta51=actcode.
   compute loc51= where.
Else if start<260.
   compute acta52=actcode.
   compute loc52= where.
Else if start<265.
   compute acta53=actcode.
   compute loc53= where.
Else if start<270.
   compute acta54=actcode.
   compute loc54= where.
Else if start<275.
   compute acta55=actcode.
   compute loc55= where.
Else if start<280.
   compute acta56=actcode.
   compute loc56= where.
Else if start<285.
   compute acta57=actcode.
   compute loc57= where.
Else if start<290.
   compute acta58=actcode.
   compute loc58= where.
Else if start<295.
   compute acta59=actcode.
   compute loc59= where.
Else if start<300.
   compute acta60=actcode.
   compute loc60= where.
Else if start<305.
   compute acta61=actcode.
   compute loc61= where.
Else if start<310.
   compute acta62=actcode.
   compute loc62= where.
Else if start<315.
   compute acta63=actcode.
   compute loc63= where.
Else if start<320.
   compute acta64=actcode.
   compute  loc64=where.
Else if start<325.
   compute acta65=actcode.
   compute  loc65=where.
Else if start<330.
   compute acta66=actcode.
   compute  loc66=where.
Else if start<335.
   compute acta67=actcode.
   compute  loc67=where.
Else if start<340.
   compute acta68=actcode.
   compute  loc68=where.
Else if start<345.
   compute acta69=actcode.
   compute  loc69=where.
Else if start<350.
   compute acta70=actcode.
   compute  loc70=where.
Else if start<355.
   compute acta71=actcode.
   compute  loc71=where.
Else if start<360.
   compute acta72=actcode.
   compute  loc72=where.
Else if start<365.
   compute acta73=actcode.
   compute  loc73=where.
Else if start<370.
   compute acta74=actcode.
   compute  loc74=where.
Else if start<375.
   compute acta75=actcode.
   compute  loc75=where.
Else if start<380.
   compute acta76=actcode.
   compute  loc76=where.
Else if start<385.
   compute acta77=actcode.
   compute  loc77=where.
Else if start<390.
   compute acta78=actcode.
   compute  loc78=where.
Else if start<395.
   compute acta79=actcode.
   compute  loc79=where.
Else if start<400.
   compute acta80=actcode.
   compute  loc80=where.
Else if start<405.
   compute acta81=actcode.
   compute  loc81=where.
Else if start<410.
   compute acta82=actcode.
   compute  loc82=where.
Else if start<415.
   compute acta83=actcode.
   compute  loc83=where.
Else if start<420.
   compute acta84=actcode.
   compute  loc84=where.
Else if start<425.
   compute acta85=actcode.
   compute  loc85=where.
Else if start<430.
   compute acta86=actcode.
   compute  loc86=where.
Else if start<435.
   compute acta87=actcode.
   compute  loc87=where.
Else if start<440.
   compute acta88=actcode.
   compute  loc88=where.
Else if start<445.
   compute acta89=actcode.
   compute  loc89=where.
Else if start<450.
   compute acta90=actcode.
   compute  loc90=where.
Else if start<455.
   compute acta91=actcode.
   compute  loc91=where.
Else if start<460.
   compute acta92=actcode.
   compute  loc92=where.
Else if start<465.
   compute acta93=actcode.
   compute  loc93=where.
Else if start<470.
   compute acta94=actcode.
   compute  loc94=where.
Else if start<475.
   compute acta95=actcode.
   compute  loc95=where.
Else if start<480.
   compute acta96=actcode.
   compute  loc96=where.
Else if start<485.
   compute acta97=actcode.
   compute  loc97=where.
Else if start<490.
   compute acta98=actcode.
   compute  loc98=where.
Else if start<495.
   compute acta99=actcode.
   compute  loc99=where.
Else if start<500.
   compute acta100=actcode.
   compute  loc100=where.
Else if start<505.
   compute acta101=actcode.
   compute  loc101=where.
Else if start<510.
   compute acta102=actcode.
   compute  loc102=where.
Else if start<515.
   compute acta103=actcode.
   compute  loc103=where.
Else if start<520.
   compute acta104=actcode.
   compute  loc104=where.
Else if start<525.
   compute acta105=actcode.
   compute  loc105=where.
Else if start<530.
   compute acta106=actcode.
   compute  loc106=where.
Else if start<535.
   compute acta107=actcode.
   compute  loc107=where.
Else if start<540.
   compute acta108=actcode.
   compute  loc108=where.
Else if start<545.
   compute acta109=actcode.
   compute  loc109=where.
Else if start<550.
   compute acta110=actcode.
   compute  loc110=where.
Else if start<555.
   compute acta111=actcode.
   compute  loc111=where.
Else if start<560.
   compute acta112=actcode.
   compute  loc112=where.
Else if start<565.
   compute acta113=actcode.
   compute  loc113=where.
Else if start<570.
   compute acta114=actcode.
   compute  loc114=where.
Else if start<575.
   compute acta115=actcode.
   compute  loc115=where.
Else if start<580.
   compute acta116=actcode.
   compute  loc116=where.
Else if start<585.
   compute acta117=actcode.
   compute  loc117=where.
Else if start<590.
   compute acta118=actcode.
   compute  loc118=where.
Else if start<595.
   compute acta119=actcode.
   compute  loc119=where.
Else if start<600.
   compute acta120=actcode.
   compute  loc120=where.
Else if start<605.
   compute acta121=actcode.
   compute  loc121=where.
Else if start<610.
   compute acta122=actcode.
   compute  loc122=where.
Else if start<615.
   compute acta123=actcode.
   compute  loc123=where.
Else if start<620.
   compute acta124=actcode.
   compute  loc124=where.
Else if start<625.
   compute acta125=actcode.
   compute  loc125=where.
Else if start<630.
   compute acta126=actcode.
   compute  loc126=where.
Else if start<635.
   compute acta127=actcode.
   compute  loc127=where.
Else if start<640.
   compute acta128=actcode.
   compute  loc128=where.
Else if start<645.
   compute acta129=actcode.
   compute  loc129=where.
Else if start<650.
   compute acta130=actcode.
   compute  loc130=where.
Else if start<655.
   compute acta131=actcode.
   compute  loc131=where.
Else if start<660.
   compute acta132=actcode.
   compute  loc132=where.
Else if start<665.
   compute acta133=actcode.
   compute  loc133=where.
Else if start<670.
   compute acta134=actcode.
   compute  loc134=where.
Else if start<675.
   compute acta135=actcode.
   compute  loc135=where.
Else if start<680.
   compute acta136=actcode.
   compute  loc136=where.
Else if start<685.
   compute acta137=actcode.
   compute  loc137=where.
Else if start<690.
   compute acta138=actcode.
   compute  loc138=where.
Else if start<695.
   compute acta139=actcode.
   compute  loc139=where.
Else if start<700.
   compute acta140=actcode.
   compute  loc140=where.
Else if start<705.
   compute acta141=actcode.
   compute  loc141=where.
Else if start<710.
   compute acta142=actcode.
   compute  loc142=where.
Else if start<715.
   compute acta143=actcode.
   compute  loc143=where.
Else if start<720.
   compute acta144=actcode.
   compute  loc144=where.
Else if start<725.
   compute acta145=actcode.
   compute  loc145=where.
Else if start<730.
   compute acta146=actcode.
   compute  loc146=where.
Else if start<735.
   compute acta147=actcode.
   compute  loc147=where.
Else if start<740.
   compute acta148=actcode.
   compute  loc148=where.
Else if start<745.
   compute acta149=actcode.
   compute  loc149=where.
Else if start<750.
   compute acta150=actcode.
   compute  loc150=where.
Else if start<755.
   compute acta151=actcode.
   compute  loc151=where.
Else if start<760.
   compute acta152=actcode.
   compute  loc152=where.
Else if start<765.
   compute acta153=actcode.
   compute  loc153=where.
Else if start<770.
   compute acta154=actcode.
   compute  loc154=where.
Else if start<775.
   compute acta155=actcode.
   compute  loc155=where.
Else if start<780.
   compute acta156=actcode.
   compute  loc156=where.
Else if start<785.
   compute acta157=actcode.
   compute  loc157=where.
Else if start<790.
   compute acta158=actcode.
   compute  loc158=where.
Else if start<795.
   compute acta159=actcode.
   compute  loc159=where.
Else if start<800.
   compute acta160=actcode.
   compute  loc160=where.
Else if start<805.
   compute acta161=actcode.
   compute  loc161=where.
Else if start<810.
   compute acta162=actcode.
   compute  loc162=where.
Else if start<815.
   compute acta163=actcode.
   compute  loc163=where.
Else if start<820.
   compute acta164=actcode.
   compute  loc164=where.
Else if start<825.
   compute acta165=actcode.
   compute  loc165=where.
Else if start<830.
   compute acta166=actcode.
   compute  loc166=where.
Else if start<835.
   compute acta167=actcode.
   compute  loc167=where.
Else if start<840.
   compute acta168=actcode.
   compute  loc168=where.
Else if start<845.
   compute acta169=actcode.
   compute  loc169=where.
Else if start<850.
   compute acta170=actcode.
   compute  loc170=where.
Else if start<855.
   compute acta171=actcode.
   compute  loc171=where.
Else if start<860.
   compute acta172=actcode.
   compute  loc172=where.
Else if start<865.
   compute acta173=actcode.
   compute  loc173=where.
Else if start<870.
   compute acta174=actcode.
   compute  loc174=where.
Else if start<875.
   compute acta175=actcode.
   compute  loc175=where.
Else if start<880.
   compute acta176=actcode.
   compute  loc176=where.
Else if start<885.
   compute acta177=actcode.
   compute  loc177=where.
Else if start<890.
   compute acta178=actcode.
   compute  loc178=where.
Else if start<895.
   compute acta179=actcode.
   compute  loc179=where.
Else if start<900.
   compute acta180=actcode.
   compute  loc180=where.
Else if start<905.
   compute acta181=actcode.
   compute  loc181=where.
Else if start<910.
   compute acta182=actcode.
   compute  loc182=where.
Else if start<915.
   compute acta183=actcode.
   compute  loc183=where.
Else if start<920.
   compute acta184=actcode.
   compute  loc184=where.
Else if start<925.
   compute acta185=actcode.
   compute  loc185=where.
Else if start<930.
   compute acta186=actcode.
   compute  loc186=where.
Else if start<935.
   compute acta187=actcode.
   compute  loc187=where.
Else if start<940.
   compute acta188=actcode.
   compute  loc188=where.
Else if start<945.
   compute acta189=actcode.
   compute  loc189=where.
Else if start<950.
   compute acta190=actcode.
   compute  loc190=where.
Else if start<955.
   compute acta191=actcode.
   compute  loc191=where.
Else if start<960.
   compute acta192=actcode.
   compute  loc192=where.
Else if start<965.
   compute acta193=actcode.
   compute  loc193=where.
Else if start<970.
   compute acta194=actcode.
   compute  loc194=where.
Else if start<975.
   compute acta195=actcode.
   compute  loc195=where.
Else if start<980.
   compute acta196=actcode.
   compute  loc196=where.
Else if start<985.
   compute acta197=actcode.
   compute  loc197=where.
Else if start<990.
   compute acta198=actcode.
   compute  loc198=where.
Else if start<995.
   compute acta199=actcode.
   compute  loc199=where.
Else if start<1000.
   compute acta200=actcode.
   compute  loc200=where.
Else if start<1005.
   compute acta201=actcode.
   compute  loc201=where.
Else if start<1010.
   compute acta202=actcode.
   compute  loc202=where.
Else if start<1015.
   compute acta203=actcode.
   compute  loc203=where.
Else if start<1020.
   compute acta204=actcode.
   compute  loc204=where.
Else if start<1025.
   compute acta205=actcode.
   compute  loc205=where.
Else if start<1030.
   compute acta206=actcode.
   compute  loc206=where.
Else if start<1035.
   compute acta207=actcode.
   compute  loc207=where.
Else if start<1040.
   compute acta208=actcode.
   compute  loc208=where.
Else if start<1045.
   compute acta209=actcode.
   compute  loc209=where.
Else if start<1050.
   compute acta210=actcode.
   compute  loc210=where.
Else if start<1055.
   compute acta211=actcode.
   compute  loc211=where.
Else if start<1060.
   compute acta212=actcode.
   compute  loc212=where.
Else if start<1065.
   compute acta213=actcode.
   compute  loc213=where.
Else if start<1070.
   compute acta214=actcode.
   compute  loc214=where.
Else if start<1075.
   compute acta215=actcode.
   compute  loc215=where.
Else if start<1080.
   compute acta216=actcode.
   compute  loc216=where.
Else if start<1085.
   compute acta217=actcode.
   compute  loc217=where.
Else if start<1090.
   compute acta218=actcode.
   compute  loc218=where.
Else if start<1095.
   compute acta219=actcode.
   compute  loc219=where.
Else if start<1100.
   compute acta220=actcode.
   compute  loc220=where.
Else if start<1105.
   compute acta221=actcode.
   compute  loc221=where.
Else if start<1110.
   compute acta222=actcode.
   compute  loc222=where.
Else if start<1115.
   compute acta223=actcode.
   compute  loc223=where.
Else if start<1120.
   compute acta224=actcode.
   compute  loc224=where.
Else if start<1125.
   compute acta225=actcode.
   compute  loc225=where.
Else if start<1130.
   compute acta226=actcode.
   compute  loc226=where.
Else if start<1135.
   compute acta227=actcode.
   compute  loc227=where.
Else if start<1140.
   compute acta228=actcode.
   compute  loc228=where.
Else if start<1145.
   compute acta229=actcode.
   compute  loc229=where.
Else if start<1150.
   compute acta230=actcode.
   compute  loc230=where.
Else if start<1155.
   compute acta231=actcode.
   compute  loc231=where.
Else if start<1160.
   compute acta232=actcode.
   compute  loc232=where.
Else if start<1165.
   compute acta233=actcode.
   compute  loc233=where.
Else if start<1170.
   compute acta234=actcode.
   compute  loc234=where.
Else if start<1175.
   compute acta235=actcode.
   compute  loc235=where.
Else if start<1180.
   compute acta236=actcode.
   compute  loc236=where.
Else if start<1185.
   compute acta237=actcode.
   compute  loc237=where.
Else if start<1190.
   compute acta238=actcode.
   compute  loc238=where.
Else if start<1195.
   compute acta239=actcode.
   compute  loc239=where.
Else if start<1200.
   compute acta240=actcode.
   compute  loc240=where.
Else if start<1205.
   compute acta241=actcode.
   compute  loc241=where.
Else if start<1210.
   compute acta242=actcode.
   compute  loc242=where.
Else if start<1215.
   compute acta243=actcode.
   compute  loc243=where.
Else if start<1220.
   compute acta244=actcode.
   compute  loc244=where.
Else if start<1225.
   compute acta245=actcode.
   compute  loc245=where.
Else if start<1230.
   compute acta246=actcode.
   compute  loc246=where.
Else if start<1235.
   compute acta247=actcode.
   compute  loc247=where.
Else if start<1240.
   compute acta248=actcode.
   compute  loc248=where.
Else if start<1245.
   compute acta249=actcode.
   compute  loc249=where.
Else if start<1250.
   compute acta250=actcode.
   compute  loc250=where.
Else if start<1255.
   compute acta251=actcode.
   compute  loc251=where.
Else if start<1260.
   compute acta252=actcode.
   compute  loc252=where.
Else if start<1265.
   compute acta253=actcode.
   compute  loc253=where.
Else if start<1270.
   compute acta254=actcode.
   compute  loc254=where.
Else if start<1275.
   compute acta255=actcode.
   compute  loc255=where.
Else if start<1280.
   compute acta256=actcode.
   compute  loc256=where.
Else if start<1285.
   compute acta257=actcode.
   compute  loc257=where.
Else if start<1290.
   compute acta258=actcode.
   compute  loc258=where.
Else if start<1295.
   compute acta259=actcode.
   compute  loc259=where.
Else if start<1300.
   compute acta260=actcode.
   compute  loc260=where.
Else if start<1305.
   compute acta261=actcode.
   compute  loc261=where.
Else if start<1310.
   compute acta262=actcode.
   compute  loc262=where.
Else if start<1315.
   compute acta263=actcode.
   compute  loc263=where.
Else if start<1320.
   compute acta264=actcode.
   compute  loc264=where.
Else if start<1325.
   compute acta265=actcode.
   compute  loc265=where.
Else if start<1330.
   compute acta266=actcode.
   compute  loc266=where.
Else if start<1335.
   compute acta267=actcode.
   compute  loc267=where.
Else if start<1340.
   compute acta268=actcode.
   compute  loc268=where.
Else if start<1345.
   compute acta269=actcode.
   compute  loc269=where.
Else if start<1350.
   compute acta270=actcode.
   compute  loc270=where.
Else if start<1355.
   compute acta271=actcode.
   compute  loc271=where.
Else if start<1360.
   compute acta272=actcode.
   compute  loc272=where.
Else if start<1365.
   compute acta273=actcode.
   compute  loc273=where.
Else if start<1370.
   compute acta274=actcode.
   compute  loc274=where.
Else if start<1375.
   compute acta275=actcode.
   compute  loc275=where.
Else if start<1380.
   compute acta276=actcode.
   compute  loc276=where.
Else if start<1385.
   compute acta277=actcode.
   compute  loc277=where.
Else if start<1390.
   compute acta278=actcode.
   compute  loc278=where.
Else if start<1395.
   compute acta279=actcode.
   compute  loc279=where.
Else if start<1400.
   compute acta280=actcode.
   compute  loc280=where.
Else if start<1405.
   compute acta281=actcode.
   compute  loc281=where.
Else if start<1410.
   compute acta282=actcode.
   compute  loc282=where.
Else if start<1415.
   compute acta283=actcode.
   compute  loc283=where.
Else if start<1420.
   compute acta284=actcode.
   compute  loc284=where.
Else if start<1425.
   compute acta285=actcode.
   compute  loc285=where.
Else if start<1430.
   compute acta286=actcode.
   compute  loc286=where.
Else if start<1435.
   compute acta287=actcode.
   compute  loc287=where.
Else.
   compute acta288=actcode.
   compute  loc288=where.
End if.
execute.

sort cases hldid(a) persid(a) epnum(d).
Do if hldid=lag(hldid) and persid=lag(persid).
   Do repeat a=acta1 to loc288.
      if lag(a)>0 a=-5.
   End repeat.
End if.
freq acta1 loc1.

sort cases hldid(a) persid(a) epnum(a).

vector acta=acta1 to acta288 /loc=loc1 to loc288.
loop i=1 to 287.
compute j=i+1.
   Do if acta(i)>0 and acta(j)=0.
      compute acta(j)=acta(i).
   End if.
   Do if loc(i)>0 and loc(j)=0.
      compute loc(j)=loc(i).
   End if.
End loop.
execute.

match files files=*
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\1985\USA85spousevals.SAV'
   /BY hldid persid.
execute.
freq actasp1.
recode actasp1 to locsp288 (sysmis=-8).

compute sppart=0.

Do repeat a=loc1 to loc288 /b=locsp1 to locsp288 /c=acta1 to acta288 /d=actasp1 to actasp288.
*if both partners in the same place at home another person's home or church and the
*diarist is not sleeping (as who else present not collected during sleep in most surveys)
*code as spouses together.
   if ((a>0 and a<20) or a=40 or a=81) and (c>450 and c<470) and b=a sppart=1.
*If diarist is in other locations other than travel of paid work or education and the 
*spouse reports being in the same location and doing the same activity, code as spouses
*together (paid work and education are additional activities where many surveys do
*not collect who else is present information.
   if (a>50 and a<99) and b=a sppart=1.
End repeat.
freq sppart.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp.sav'
   /keep hldid persid epnum sppart.

compute sptime=0.
if sppart=1 sptime=time.
desc sptime.

aggregate 
   outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\sumsp.sav'
   /break hldid persid
   /sppart=SUM(sptime).

*put aggregate file together.

match files
      files='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\sum.sav'
      table='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa85ind.sav'
     /by hldid persid
     /drop sppart.
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 main10 main11 main47 main60 (0=-9).

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

match files files=* 
      files='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\sumsp.sav'
   /BY hldid persid.
recode sppart (sysmis=0).
if civstat=2 sppart=-7.
desc sppart.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\USA1985haf.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:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\USA1985haf.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.

recode disab (-7=-9).

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:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\USA1985haf.sav'
   /drop wtsum.

*final tidy of episode file.

get file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa1985hef.sav'
   /drop sppart.

match files files=* 
      files='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\withsp.sav'
   /BY hldid persid epnum.
execute.

rename variables (mainx secx=main sec).
recode sppart (sysmis=-7).

apply dictionary from 
   'C:\Users\Kimberly Fisher\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.

compute test=0.
if inout=2 and mtrav=1 test=1.
freq test.
*these are excursion cases, leave as participant originally coded these entries.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\usa1985hef.sav'
   /drop test.

GET file='C:\Users\Kimberly Fisher\Documents\Kimberly\time\ATUS\1985\USA85withrate.SAV'.

compute countrya=38.
compute survey=1985.
recode rating (99=-8).
var lab rating how much diarist enjoyed this activity.
val lab rating
   0 did not enjoy activity at all
  10 thoroughly enjoyed activity
  -8 missing.
freq rating.

save outfile='C:\Users\Kimberly Fisher\Documents\Kimberly\time\mtus\mtus data\W553\USA\1985\USA1985W6emotion.SAV'
   /keep countrya survey hldid persid epnum rating.