***************************************************** *SYNTAX TO CREATE World 5.53 FILE (11 November 2008)* *Kimberly Fisher, Centre for Time Use Research * *University of Oxford, United Kingdom * *earlier version constructed by Tingting Lui at the * *University of Calgary, Canada, July 2006 * *****************************************************. Get file='C:\time-use-library\Germany\1992\daten\zbe_hh.sav'. compute countreg = 13. compute survey = 1991. compute swave = 0. compute hldid = idhh. freq m71 m72. temp. select if m71=2. freq hldid. temp. select if m72=2. freq hldid. cro e1x by typ6. *e1x has more categories. compute hhldsize = e1x. freq e301x e302x e303. compute nmotor= e301x+e302x. freq nmotor. Do if nmotor>1. compute vehicle=4. Else if nmotor=1. compute vehicle=3. Else if e303>0. compute vehicle=2. Else. compute vehicle=0. End if. freq vehicle. cro vehicle by e303 nmotor. recode e3014x (0=0) (1,2=1) into computer. cro computer by e3014x. recode typ5x (1,2=1) (3 thru 10=4) (11 thru 13=2) (14 thru 30=3) (32 thru 34=4) into hhtype. recode typ5x (1,2,11,12,13,34=-7) (3,4,7,8,14,15,20,21,26,27,32=1) (5,6,9,10,17,18,19,23,24,25,29,30,33=3) into agekidx. cro typ5x by hhtype agekidx. recode typ5x (3 thru 10=4) into singphh. freq typ9x. compute incorig=typ9x. recode typ9x (50=-8) (1 thru 5=1) (10 thru 15=3) (else=2) into income. cro income by typ9x. freq income incorig. if r12x=1 msamp=10. if r12x=2 msamp=11. cro msamp by r12x. compute nchild= typ7x. freq s2 typ1 typ2. compute ownhome=-9. compute urban=-9. sort cases by hldid. SAVE OUTFILE='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\hh.sav' /KEEP countreg survey swave msamp hldid hhldsize vehicle computer hhtype agekidx incorig income nchild ownhome urban singphh s2 typ1 typ2. get file ='C:\time-use-library\Germany\1992\daten\zbe_ps.sav'. compute hldid = idhh. desc id1 e3x. *id1 is a numbered list from first to last case in the dataset, e3x is *like persid if most MTUS surveys. Use e3x for persid, but id1 needed for *matching to the diary file. Can be dropped after this match. compute persid=e3x. compute sex=e7. recode persid (sysmis=-8). freq typpart e8 e10. cro typpart by e8. Do if e8=1. compute civstat=1. compute cohab=0. Else if e10=1. compute civstat=1. compute cohab=1. Else if e8>1. compute civstat=2. compute cohab=-7. End if. freq civstat cohab. cro civstat cohab by e8 typpart e10. compute spid=0. if civstat=1 spid=persid. cro civstat by spid. cro e51x by e52x. *e52x distinguishes 15-17 from 18-20, use this variable for age. recode e52x (1=1) (2=4) (3=8) (4=13) (5=16) (6=19) (7=22) (8=27) (9=32) (10=37) (11=42) (12=47) (13=52) (14=57) (15=62) (16=67) (17=72) (18=75) into age. AGGREGATE /OUTFILE='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav' /BREAK=hldid persid /spid2 = MAX(spid) / idpart1= max(idpart1) /id1x=max(id1) /sex=max(sex) /age=max(e52x). SAVE OUTFILE='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav'. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav'. temp. select if idpart1=99997. freq spid2. *idpart maps to the reference person's spouse. There are the other *135 persons in couples who are not the reference persons or the spouse of the *reference persons. select if spid2>0. *this step recodes the id of the partner to = the id of the respondent. *this will allow the matching of the person id to the partner. rename variables (idpart1=id1). execute. compute m1=0. compute m2=999999. if id1=99997 m1=id1x. if id1=99997 m2=id1x. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /newpid1= MAX(m1) /newpid2= min(m2) /ncase=n. freq ncase. temp. select if ncase=1. freq id1. *all 20 cases of only person in couple in household are partner of *someone other than the reference person and no other peron in a *couple completed a diary. if ncase=1 spid2=-8. temp. select if ncase=3. freq id1. *5 instances of one party in the second couple kept a diary and 1 did not. if ncase=3 and id1=99997 spid2=-8. temp. select if ncase=6. freq id1. *1 household with 3 couples who completed diaries. If hldid=2798 and persid=3 id1=8041. If hldid=2798 and persid=4 id1=8040. If hldid=2798 and persid=5 id1=8043. If hldid=2798 and persid=6 id1=8042. temp. select if id1=99997. freq ncase spid2. *some households where the only couple does not include the *household reference person. Also some households where there *are two couples. if id1=99997 and (ncase=2 or ncase=4) and id1x=newpid1 id1=newpid2. if id1=99997 and (ncase=2 or ncase=4) and id1x=newpid2 id1=newpid1. execute. sort cases by hldid id1. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav' /drop persid id1x sex age newpid1 newpid2 ncase. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav'. match files files=* files='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav' /by hldid id1. cro spid by spid2. rename variables (spid2=partid). recode partid idhh (sysmis=-7). cro partid by civstat. *25 false lines added by the partners for whom there is no information. *remove these lines. select if idhh>-7. cro partid by civstat. freq e4x. cro e4x by persid. recode e4x (1=1) (2=2) (3=3) (5=4) (4,7=9) (8=10) into relrefp. compute pin=0. if e4x=5 pin=1. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /pinhh= MAX(pin). freq pinhh. compute cphome=0. Do if civstat=2. if e4x=3 or e4x=4 cphome=1. if (e4x=1 or e4x=2) and pinhh=1 cphome=1. End if. freq cphome. cro cphome by e4x civstat. compute pnt1=0. compute pnt2=99. compute p1prt=0. if e4x=5 pnt1=persid. if e4x=5 pnt2=persid. if persid=1 p1prt=partid. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /parent1= MAX(pnt1) /parent2= MIN(pnt2) /p1prtid= MAX(p1prt). freq parent1 parent2 p1prtid. compute parntid1=-7. compute parntid2=-7. if e4x=3 parntid1=1. if e4x=3 parntid2=p1prtid. execute. Do if parntid1=-7. if persid=1 or (persid=2 and p1prtid>0) parntid1=parent1. if e4x=4 parntid1=parent1. End if. Do if parntid2=-7. if persid=1 or (persid=2 and p1prtid>0) parntid2=parent2. if e4x=4 parntid2=parent2. End if. recode parntid1 parntid2 (0=-7) (99=-7). freq parntid1 parntid2. cro parntid1 parntid2 by cphome e4x. *only grand-children not mapped to parents - 34 such cases. temp. select if cphome=1 and parntid1=-7. freq e4x. *these 34 grandchild cases visually examined. *one case of siblings and one inappropriately coded as a grandchild. Do if hldid=60 and persid=4. compute parntid1=1. compute parntid2=2. compute relrefp=3. End if. *3 generations, 2 couples. Do if (hldid=179 and persid=5) or (hldid=587 and persid=5) or (hldid=798 and persid=5) or (hldid=2193 and persid>4). compute parntid1=3. compute parntid2=4. End if. *single parents living with parents. Do if (hldid=1361 or hldid=1406 or hldid=1702 or hldid=3236 or hldid=3740 or hldid=4378) and persid>3. compute parntid1=3. End if. *adult living with great grandparent or grandparent. Do if (hldid=1212 and persid=2) or (hldid=1679 and persid=3) or (hldid=2488 and persid=2). compute cphome=0. End if. *3 generations, 2 couples. Do if (hldid=3192 and persid>5). compute parntid1=4. compute parntid2=5. End if. *single parents living with parents. Do if hldid=3195 and persid=5. compute parntid1=4. End if. *single parents living with parents. Do if (hldid=3339 and persid=5) or (hldid=5006 and persid>4) or (hldid=5082 and persid=5). compute parntid1=3. End if. *3 generation 2 couple household with 2nd couple cohabitting. Do if hldid=3405 and persid=4. compute parntid1=3. compute parntid2=5. End if. *4 cases of grandchild or great grandchild being looked after by grandparent. *hldid=3421 hldid=4811 hldid=5046 hldid=5517. *4 generation household with grandchild child of single parent and parent id miscoded. Do if hldid=3901 and persid=3. compute parntid1=2. compute parntid2=-7. End if. *3 generation household. Do if hldid=3943 and persid>3. compute parntid1=2. compute parntid2=3. End if. *3 generation household with half of couple with young child not in household. Do if hldid=5494 and persid=5. compute parntid1=5. End if. *3 generations, 2 couples. Do if (hldid=5622 and persid=5). compute parntid1=3. compute parntid2=4. End if. cro cphome by parntid1 parntid2. temp. select if cphome=1 and parntid1=-7. freq e4x. temp. select if cphome=1 and parntid1>-7. freq civstat. temp. select if cphome=0 and parntid1>-7. freq civstat. compute test=0. if civstat=2 and cphome=0 and parntid1>-7 test=1. freq test. *correct errors. Do if hldid=2844 and persid=1. compute parntid1=2. compute parntid2=-7. Else if hldid=2844 and persid=3. compute parntid1=2. compute parntid2=-7. compute relrefp=2. Else if hldid=2844. compute parntid1=-7. compute parntid2=-7. end if. execute. cro parntid1 by parntid2. compute test=0. if parntid1>1 and parntid1=parntid2 test=1. freq test. *237 problems to fix. if parntid1>1 and parntid1=parntid2 parntid2=-7. cro parntid1 by parntid2. compute citizen=-9. compute rushed=-9. compute health=-9. *carer. cro s723ax by s723b2x. *most cases on both variables are people who do unpaid care of *adults, a handful unclear if they get pay, but not clear if *those paid get paid a minor incentive or this is their *main job. As the performance of care likely to impact daily *routines all these cases recorded as carers. recode s723ax (1=1) (2=0) (else=-8) into carer. cro s723ax by carer. freq carer. recode s3 (1 thru 3=1) (else=0) into disab. freq disab. temp. select if disab=1. freq age. freq e21x e26x. cro e21x by e26x. recode e21x (1=3) (2=7) (3=12) (4=17) (5=22) (6=27) (7=32) (8=37) (9=42) (10=47) (11=52) (12=57) (13=62) (14=67) (15=72) (16=75) (99997=-7) into workhrs. If e26x=1 workhrs=workhrs+1. If e26x=2 workhrs=workhrs+2. If e26x=3 workhrs=workhrs+3. If e26x=4 workhrs=workhrs+4. If e26x=5 workhrs=workhrs+7. If e26x=6 workhrs=workhrs+12. If e26x=7 workhrs=workhrs+17. If e26x=8 workhrs=workhrs+22. If e26x=9 workhrs=workhrs+25. freq workhrs. temp. select if workhrs>0. freq e28. cro e29x by e28. temp. select if workhrs>0. cro e29x by e28. Do if e28=1. compute empinclm=e29x. Else if workhrs>0. compute empinclm=e29x*100. Else. compute empinclm=-7. End if. recode empinclm (5000=50) (9900=99). freq empinclm. Recode e12 (0=0) (1 thru 4=1) into student. freq student. cro student by e12 e14 e281 sic1712 sic1722. cro e13 by e14. compute educa = (e13*100)+e14. Do if e13=5 or e14>4. compute edtry=3. Else if e13=3 or e13=4. compute edtry=2. Else if e14=4. compute edtry=2. Else if e13<3 or e13=6. compute edtry=1. End if. freq edtry. cro edtry by e13 e14. means workhrs by e20 /cells=min max mean. Do if e20=1 or workhrs>34. compute empstat=1. Else if (e20>1 and e20<7). compute empstat=2. Else if workhrs>0. compute empstat=-8. Else. compute empstat=4. End if. freq empstat. cro empstat by e20 e281. recode empstat (1,2,3 =1) (4=0) (else=-8) into emp. compute unemp=0. if e281=2 unemp=1. cro unemp by e281. compute retired=0. if e281=1 or e28=2 or age>70 retired=1. cro retired by e281 e28. Do if e17=3 or e17=8 or e17=9. compute sector=1. Else if e17<3 or (e17>3 and e17<8) or e17=10. compute sector=2. Else. compute sector=-7. End if. freq sector. cro sector by e17 empstat. cro e29x by e18x. temp. select if e17=3. freq e18x. Do if e17=3 and (e18x=27 or e18x=28 or e18x=24). compute occup=4. Else if e18x=27. compute occup=1. Else if e18x=30 or e17=8. compute occup=10. Else if e18x=28. compute occup=2. Else if e18x=24. compute occup=3. Else if e18x=34. compute occup=5. Else if e18x=32. compute occup=6. Else if e18x=31. compute occup=7. Else if (e18x=14 or e18x=33) and (e17>2 and e17<10). compute occup=8. Else if e18x=29 and (e17>2 and e17<10). compute occup=9. Else if (e18x=25 or (e18x>34 and e18x<38)) and (e17>2 and e17<10). compute occup=11. Else if e18x<3. compute occup=12. Else if ((e18x>2 and e18x<14) or (e18x>14 and e18x<24) or e18x=26) and (e17>2 and e17<10). compute occup=13. Else if e17<3 or e17=10 or e18x=38. compute occup=14. Else. compute occup=-7. end if. freq occup. cro occup by e17. temp. select if e18x<13. cro occup by e18x. temp. select if e18x>12 and e18x<25. cro occup by e18x. temp. select if e18x>24. cro occup by e18x. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav'. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav' /keep hldid partid empstat. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav'. select if partid>-7. rename variables (empstat partid=empsp persid). sort cases by hldid persid. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav'. match files files='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\sp.sav' files='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav' /by hldid persid. recode empsp (sysmis=-7). cro empsp by partid civstat empstat. if civstat=1 and empsp=-7 empsp=-8. if civstat=1 and partid=-7 partid=-8. freq empsp. match files files=* table='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\hh.sav' /by hldid. execute. freq age agekidx. compute ak=99. if parntid1>0 ak=age. freq ak. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /agekid2= MIN(ak). freq agekid2. cro agekid2 by agekidx. *a number of inconsitencies corrected by household grid information, *also full codes for agekidx now possible. if agekid2<5 and agekidx~=1 agekidx=1. if agekid2=8 and agekidx=3 agekidx=2. if (agekid2>17 and agekid2<99) and agekidx=-7 agekidx=4. if agekid2=99 and agekidx=1 agekid2=3. if agekid2=99 and agekidx=3 agekid2=15. Do if agekid2=8 and agekidx=1. compute agekidx=2. Else if agekid2=13 and agekidx=1. compute agekidx=2. compute agekid2=8. Else if agekid2>13 and agekidx=1. compute agekid2=3. Else if agekid2>16 and agekidx=3. compute agekid2=15. End if. execute. recode agekid2 (99=-7). freq agekidx agekid2. cro agekid2 by agekidx. *check nchild and hhsize. compute sz=1. compute nc=0. if age<19 nc=1. means age by nc /cells = min max. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /nc2= SUM(nc) /hhsz=SUM(sz). cro nc2 by nchild. cro nc2 nchild by agekidx. *the household level nchild largely overlaps, but there are *a number of inconsistencies. Is the variable from the individual *level information. rename variables (nchild=nchildhh). rename variables (nc2=nchild). cro nchild by agekidx. cro hhsz by hhldsize. *these match perfectly. Do if (age>18 and age<40) and (agekidx=-7 or agekidx=4). compute famstat=0. Else if age>39 and (agekidx=-7 or agekidx=4). compute famstat=3. Else if age>18 and agekidx=1. compute famstat=1. Else if age>18 and (agekidx=2 or agekidx=3). compute famstat=2. Else if age<19 and cphome=1. compute famstat=4. Else if age<19 and cphome=0. compute famstat=5. End if. freq famstat. cro famstat by cphome agekidx. means age by famstat /cells=min max. recode singphh (sysmis=0). cro singphh by nc hhldsize hhtype nchild. temp. select if singphh=4. cro hhldsize by nchild. compute singpar=0. if singphh=4 and nc=0 singpar=1. freq singpar. cro singpar by singphh nc hhtype. *pull in formation from coding of parent id above to mark *11 single parents who live in complex households. Do if (hldid=1361 or hldid=1406 or hldid=1702 or hldid=3236 or hldid=3740 or hldid=4378) and persid=3. compute singpar=1. Else if hldid=3195 and persid=4. compute singpar=1. Else if (hldid=3339 or hldid=5006 or hldid=5082) and persid=3. compute singpar=1. Else if hldid=3901 and persid=2. compute singpar=1. End if. freq singpar. *is a household level variable marking if any member needs care. *see if this variable might identify where missing cases for the *variable disab. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=hldid /ndis= SUM(disab). freq ndis. cro ndis by s2. *disab has a particularly low value for this country, but all cases *which could have been identified have been identified. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav' /keep countreg survey swave msamp hldid persid partid parntid1 parntid2 hhldsize vehicle computer hhtype agekidx incorig income nchild ownhome urban agekid2 empsp sex civstat cohab age relrefp cphome citizen rushed health carer disab workhrs empinclm student educa edtry empstat emp unemp retired sector occup famstat singpar sic151 sic152 t11 t12 t21 t22 id1. get file='C:\time-use-library\Germany\1992\daten\zbe_tak.sav'. compute id=idtag. rename variables (idtag=diary). freq id diary. compute hldid = idhh. freq zha1 zga1 zfw1 zwo1 zmk1 zmh1 zmv1 zmc1 zms1 zma1. *note location only at home, not at home. * *note with whom appears as 1 variable with *0=not reported, 1=relatives/neighbours, 2=friends *+ 5 0=no 1=yes present for: *children, other household members, colleagues, others, alone. *the is a problem for coding time with spouse, as time *with other household members will be spouse only for *couple alone, while with with other household members may *include spouse or may not include spouse but only other *household members. * *mode of transport coded within activity codes, also some *location codes in activity codes. * *missing activity code 920. *diary processing. compute miss=0. Do repeat a= zha1 to zha288. if a=920 miss=miss+5. End repeat. freq miss. *.3% missing time, 9 diaries with >91 missing minutes. compute miss=0. Do repeat a= zha1 to zha288 / b=zga1 to zga288. Do if a=920 and ((b>0 and b<920) or b>920). compute a=b. compute b=0. End if. End repeat. execute. Do repeat a= zha1 to zha288. if a=920 miss=miss+5. End repeat. freq miss. *.3% missing time, 9 diaries with >91 missing minutes, but some *minor improvement in the number of diaries with no missing time. compute nep=1. Vector zha=zha1 to zha288 / zga=zga1 to zga288 /zwo=zwo1 to zwo288. Loop i=1 to 277. compute j=i+1. If zha(i)~=zha(j) or zga(i)~=zga(j) or zwo(i)~=zwo(j) nep=nep+1. End loop. freq nep. *26 low-episode diaries. recode miss (91 thru high=1) (else=0) into miss2. recode nep (1 thru 6=1) (else=0) into nep2. cro nep2 by miss2. *Only 1 diary bad on both these cases. compute residual=0. compute carflag=0. compute eatdrk=0. compute trav=0. compute rst=0. compute pcare=0. compute sppart=0. vector av(41) main(69). do repeat a=av1 to av41. compute a=0. End repeat. Do repeat a=main1 to main69. compute a=0. End repeat. freq carflag sppart av35 main28. Do repeat a=zmh1 to zmh288. If a=1 sppart=sppart+5. End repeat. freq sppart. *test for how to identify walking the dog. compute test=0. Do repeat a= zha1 to zha288 /b= zga1 to zga288 /c= zwo1 to zwo288. if ((a>43 and a<49) and (b=81 or b=91 or b=181 or b=191 or b=281 or b=291 or b=381 or b=481 or b=581 or b=681 or b=781 or b=881 or b=891 or b=981)) or ((a=721 or a=722 or a=81 or a=91 or a=181 or a=191 or a=281 or a=291 or a=381 or a=481 or a=581 or a=681 or a=781 or a=881 or a=891 or a=981) and (b>43 and b<49)) test=test+5. End repeat. freq test. *not much of this activity, but record separately anyway. *look at cases with missing diary time (80 have some missing time). *no cases in the first episode, only 5 missing the last episode. compute test=0. if zha288=920 test=1. freq test. *person who goes out in evening in first diary, comes home, goes to bed; *in second diary, goes out, walks home, then 90 missing minutes - recode *part of imputed personal and household care and part as sleep. Do repeat a=zha271 to zha279. if idhh=1157 and id1=3277 and diary=2 a=1. End repeat. Do repeat a=zha280 to zha288. if idhh=1157 and id1=3277 and diary=2 a=3. End repeat. *the other 4 are people who kept one good diary then gave up half-way through the *second diary. No change made to huge gaps. sort cases by miss(d). *asleep then outside cleaning. Do repeat a=zha210 to zha251. if idhh=5443 and id1=15621 and diary=1 a=3. End repeat. *many activities all recorded at home, other people in and out of house. Do repeat a=zha107 to zha112 /b=zha175 to zha180. if idhh=1791 and id1=5139 and diary=1 a=1. if idhh=1791 and id1=5139 and diary=1 b=1. end repeat. Do repeat a=zha127 to zha134. if idhh=1791 and id1=5139 and diary=1 a=1. end repeat. Do repeat a=zha226 to zha228. if idhh=1791 and id1=5139 and diary=1 a=1. end repeat. Do repeat a=zha235 to zha239 /b=zha271 to zha275. if idhh=1791 and id1=5139 and diary=1 a=1. if idhh=1791 and id1=5139 and diary=1 b=1. end repeat. *housework and gardening followed by gap at home then leisure at home. Do repeat a=zha217 to zha240. if idhh=5548 and id1=15914 and diary=1 a=1. end repeat. *household activities, reading at home, then sleep at home. Do repeat a=zha237 to zha264. if idhh=5292 and id1=15202 and diary=1 a=1. end repeat. sort cases by idhh id1 diary. vector zha=zha1 to zha288 / zwo=zwo1 to zwo288. Loop i=1 to 285. compute j=i+1. compute k=i+2. compute l=i+3. Do if zha(i)=511 and zha(j)=920 and zha(k)=920 and zha(l)=920. compute zha(j)=1. compute zha(k)=1. compute zha(l)=1. Else if zha(i)=920 and zha(j)=920 and zha(k)=920 and zha(l)=511. compute zha(i)=1. compute zha(j)=1. compute zha(k)=1. Else if zha(i)=920 and zha(j)=920 and zha(k)=920 and zwo(i)=1 and zwo(j)=1 and zwo(k)=1 and zwo(l)=2. compute zha(i)=1. compute zha(j)=1. compute zha(k)=1. Else if zha(j)=920 and zha(k)=920 and zha(l)=920 and zwo(i)=2 and zwo(j)=1. compute zha(j)=1. compute zha(k)=1. compute zha(l)=1. End if. End loop. execute. vector zha=zha1 to zha288 / zwo=zwo1 to zwo288. Loop i=1 to 286. compute j=i+1. compute k=i+2. Do if zha(i)=511 and zha(j)=920 and zha(k)=920. compute zha(j)=1. compute zha(k)=1. Else if zha(i)=920 and zha(j)=920 and zha(k)=511. compute zha(i)=1. compute zha(j)=1. Else if zha(i)=920 and zha(j)=920 and zwo(i)=1 and zwo(j)=1 and zwo(k)=2. compute zha(i)=1. compute zha(j)=1. Else if zha(j)=920 and zha(k)=920 and zwo(i)=2 and zwo(j)=1. compute zha(j)=1. compute zha(k)=1. End if. End loop. execute. vector zha=zha1 to zha288 / zwo=zwo1 to zwo288. Loop i=1 to 287. compute j=i+1. Do if zha(i)=511 and zha(j)=920. compute zha(j)=1. Else if zha(i)=920 and zha(j)=511. compute zha(i)=1. Else if zha(i)=920 and zwo(i)=1 and zwo(j)=2. compute zha(i)=1. Else if zha(j)=920 and zwo(i)=2 and zwo(j)=1. compute zha(j)=1. End if. End loop. execute. compute miss=0. Do repeat a= zha1 to zha288. if a=920 miss=miss+5. End repeat. freq miss. *.3% missing time, 5 diaries with >91 missing minutes, but some *minor improvements in the number of diaries with no missing time *and reductions in missing time for other diaries. Do repeat a= zha1 to zha288 /b= zga1 to zga288 /c= zwo1 to zwo288. Do if a=11 or (a>13 and a<20). compute av6=av6+5. compute main18=main18+5. compute eatdrk=1. Else if a=12 or a=13. compute av6=av6+5. compute main19=main19+5. compute eatdrk=1. Else if a>20 and a<30. compute av8=av8+5. compute main21=main21+5. Else if a>30 and a<40. compute av7=av7+5. compute main20=main20+5. Else if (a>40 and a<44) or a=49. compute av9 = av9+5. compute main46 = main46+5. Else if a>43 and a<49. compute av8 = av8+5. compute main27 = main27+5. compute carflag=1. Else if ((a>43 and a<49) and (b=81 or b=91 or b=181 or b=191 or b=281 or b=291 or b=381 or b=481 or b=581 or b=681 or b=781 or b=881 or b=891 or b=981)) or ((a=721 or a=722 or a=81 or a=91 or a=181 or a=191 or a=281 or a=291 or a=381 or a=481 or a=581 or a=681 or a=781 or a=881 or a=891 or a=981) and (b>43 and b<49)). compute av8= av8+5. compute main47 = main47+5. compute pcare=1. compute trav=1. Else if a>43 and a<49. compute av8= av8+5. compute main27 = main27+5. compute pcare=1. Else if a=50 or a=150.. compute av10 = av10+5. compute main24 = main24+5. Else if a=63. compute av8 = av8+5. compute main23 = main23+5. compute pcare=1. Else if a=61 or (a>63 and a<71) or a=99. compute av8 = av8+5. compute main23 = main23+5. Else if a=62 or a=162. compute av10 = av10+5. compute main26 = main26+5. Else if (a>70 and a<99) or (a>170 and a<199). compute av12 = av12+5. compute main67 = main67+5. compute trav=1. Else if (a>99 and a<124) or (a>124 and a<150) or a=161 or a=169 or a=170 or a=199. compute av8 = av8+5. compute main22 = main22+5. Else if a=124. compute av40 = av40+5. compute main54 = main54+5. Else if (a>170 and a<199). compute av12 = av12+5. compute main = main +5. Else if a=220. compute av3 = av3+5. compute main9 = main9+5. Else if a=234 or a=235. compute av3 = av3+5. compute main10 = main10+5. Else if (a>200 and a<220) and (b>280 and b<297). compute av5 = av5+5. compute main11 = main11+5. Else if ((a>200 and a<220) and c=1) or a=233. compute av2 = av2+5. compute main8 = main8+5. Else if (a>200 and a<220) or a=270 or a=299. compute av1 = av1+5. compute main7 = main7+5. Else if a=250. compute av1 = av1+5. compute main12 = main12+5. compute eatdrk=1. compute rst=1. Else if a=231 or a=232 or a=239. compute av1 = av1+5. compute main13 = main13+5. Else if a=240. compute av2 = av2+5. compute main14 = main14+5. Else if (a>280 and a<299). compute av5 = av5+5. compute main63 = main63+5. compute trav=1. Else if (a>300 and a<380)or a=399. compute av23 = av23+5. compute main33 = main33+5. Else if a>380 and a<399. compute av17 = av17+5. compute main65 = main65+5. compute trav=1. Else if (a>400 and a<450) and c=1. compute av33 = av33+5. compute main16 = main16+5. Else if (a>400 and a<440) or a=470 or a=499. compute av4 = av4+5. compute main15 = main15+5. Else if a=440. compute av4 = av4+5. compute main17 = main17+5. Else if a=450. compute av4 = av4+5. compute main5 = main5+5. compute eatdrk=1. Else if a>480 and a<499. compute av5 = av5+5. compute main64 = main64+5. compute trav=1. Else if a=511 or a=519. compute av16 = av16+5. compute main2 = main2+5. compute rst=1. Else if a=512. compute av36 = av36+5. compute main55 = main55+5. compute rst=1. Else if (a>520 and a<530) or a=570 or a=599. compute av13 = av13+5. compute main4 = main4+5. compute pcare=1. Else if a=531 or a=539. compute av15 = av15+5. compute main6 = main6+5. compute eatdrk=1. Else if a=532. compute av28 = av28+5. compute main39 = main39+5. compute eatdrk=1. Else if a=533. compute av28 = av28+5. compute main5 = main5+5. compute eatdrk=1. Else if a>580 and a<599. compute av12 = av12+5. compute main67 = main67+5. compute trav=1. Else if a>600 and a<630. compute av37 = av37+5. compute main49 = main49 +5. Else if a=631 and c=1. compute av38 = av38+5. compute main48 = main48+5. Else if a=631. compute av29 = av29+5. compute main48 = main48+5. Else if (a>631 and a<640) or a=670 or a=699. compute av25 = av25+5. compute main40 = main40+5. Else if a=640. compute av18 = av18+5. compute main38 = main38+5. compute trav=1. Else if a>680 and a<699. compute av17 = av17+5. compute main68 = main68+5. compute trav=1. Else if a=711. compute av31 = av31+5. compute main59 = main59+5. Else if a=712. compute av30 = av30+5. compute main58 = main58+5. Else if a=713. compute av32 = av32+5. compute main57 = main57+5. Else if a=714. compute av34 = av34+5. compute main56 = main56+5. Else if a>714 and a<718. compute av35 = av35+5. compute main56 = main56+5. Else if a=718. compute av33 = av33+5. compute main61 = main61+5. Else if a=719 or a=770. compute av40 = av40+5. compute main51 = main51+5. Else if a=721. compute av19 = av19+5. compute main42 = main42+5. compute trav=1. Else if a=722 or ((a=729 or (a>733 and a<770) or a=799) and b=781). compute av21 = av21+5. compute main43 = main43+5. compute trav=1. Else if (a=729 or (a>733 and a<770) or a=799) and b=782. compute av19 = av19+5. compute main44 = main44+5. compute trav=1. Else if a=729 or (a>733 and a<770) or a=799. compute av40 = av40+5. compute main35 = main35+5. Else if a=723. compute av40 = av40+5. compute main50 = main50+5. Else if a=724. compute av11 = av11+5. compute main31 = main31+5. compute carflag=1. Else if a=731. compute av40 = av40+5. compute main52 = main52+5. Else if a=732. compute av20 = av20+5. compute main36 = main36+5. Else if a=733. compute av24 = av24+5. compute main37 = main37+5. Else if a>780 and a<799. compute av17 = av17+5. compute main68 = main68+5. compute trav=1. Else if a=811. compute av11 = av11+5. compute main29 = main29+5. compute carflag=1. Else if a=812. compute av11 = av11+5. compute main30 = main30+5. compute carflag=1. Else if a=813 or a=824. compute av11 = av11+5. compute main28 = main28+5. compute carflag=1. Else if a=819 or a=870 or a=899. compute av11 = av11+5. compute main31 = main31+5. compute carflag=1. Else if (a>819 and a<824) or a=829 or a=875. compute av8 = av8+5. compute main32 = main32+5. compute carflag=1. Else if a>875 and a<899. compute av12 = av12+5. compute main66 = main66+5. compute carflag=1. compute trav=1. Else if a=910. compute av23 = av23+5. compute main33 = main33+5. Else if a=920. compute av41 = av41+5. compute main69 = main69+5. Else if a=970. compute av13 = av13+5. compute main4 = main4+5. compute pcare=1. Else if a>970 and a<999. compute av17 = av17+5. compute main68 = main68+5. compute trav=1. Else if a=999 and c=1. compute av40 = av40+5. compute main51 = main51+5. Else if a=999 and c=2. compute av40 = av40+5. compute main35 = main35+5. Else if a=1. compute av41 = av41+5. compute main1 = main1+5. compute pcare=1. Else if a=3. compute av16 = av16+5. compute main3 = main3+5. compute rst=1. Else. compute residual=a. End if. End repeat. freq residual. compute misbasic=0. If eatdrk=0 misbasic=misbasic+1. If trav=0 misbasic=misbasic+1. If rst=0 misbasic=misbasic+1. If pcare=0 misbasic=misbasic+1. freq misbasic. *90.4% of diaries not missing any basic information; *72 diaries missing 2+ categories of basic information. cro carflag by misbasic. *24 cases carers missing 2 activities, 1 carer missing *all 4 basic activities. temp. select if misbasic>1 and carflag=1. freq eatdrk trav rst pcare miss nep. temp. select if misbasic=4 and carflag=1. freq nep. *all carers missing 2+ activities recorded no travel, only 1 * recorded personal care. Only 1 recorded no eating or drinking, *and only the 1 carer diarist missing all 4 basic activities *recorded no rest. None of these diaries is missing too much *activity time to count as good diaries. 24 have good episode *counts, the diarist missing all 4 basic activities has only *1 episode so would be excluded anyway. compute baddiary=0. if miss>90 or nep<7 or misbasic>1 baddiary=1. if misbasic=2 and carflag=1 baddiary=0. freq baddiary. *73 - 0.3% bad diaries. recode nep (0 thru 6=1) (7 thru high=0). recode miss (91 thru high=1) (else=0). if carflag=1 and misbasic=2 misbasic=0. recode misbasic (0,1=0) (2 thru 4=1). cro miss by misbasic by nep. *no diaries bad on all 3 criterion. *6 diaries bad on 2 criterion - 1 low episode and missing 91+ *minutes, 5 low episode and missing basic activities. *4 diaries only missing 91+ minutes *43 diaries missing only basic activities *20 diaries only low episode. compute tot1=sum (av1 to av41). compute tot2=sum (main1 to main69). freq tot1 tot2. *happy news, both total 1440. sort cases by idhh id1 diary. match files files=* table='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ind.sav' /by hldid id1. freq age. *youngest people no longer included. cro t12 by t22. *mostly both diaries completed in the same month, but a few cases where *the second diary completed in the second month. Do if diary=1. compute month=t12. Else if diary=2. compute month=t22. End if. freq month. cro month by t12 t22. compute day=0. Do if diary=1 and t11=1 or sic151=1. compute day=2. Else if diary=1 and t11=2 or sic151=2. compute day=3. Else if diary=1 and t11=3 or sic151=3. compute day=4. Else if diary=1 and t11=4 or sic151=4. compute day=5. Else if diary=1 and t11=5 or sic151=5. compute day=6. Else if diary=1 and t11=6 or sic151=6. compute day=7. Else if diary=1 and t11=7 or sic151=7. compute day=1. Else if diary=2 and t21=1. compute day=2. Else if diary=2 and t21=2. compute day=3. Else if diary=2 and t21=3. compute day=4. Else if diary=2 and t21=4. compute day=5. Else if diary=2 and t21=5. compute day=6. Else if diary=2 and t21=6. compute day=7. Else if diary=2 and t21=7. compute day=1. End if. freq day. temp. select if diary=1. cro day by sic151 t11. temp. select if diary=2. cro day by sic152 t21. compute test1=0. compute test2=0. compute test3=0. if diary=2 and day=t21 test1=1. if hldid=56 or hldid=316 or hldid=667 or hldid=1011 or hldid=1174 or hldid=1705 or hldid=2262 or hldid=2693 or hldid=2836 or hldid=3303 or hldid=3679 or hldid=3816 or hldid=4036 or hldid=4096 or hldid=4169 or hldid=4470 or hldid=4559 or hldid=4731 or hldid=4822 or hldid=5088 or hldid=5292 or hldid=5722 test2=1. if persid=lag(persid) and day=lag(day) test3=1. cro test1 by test2 test3. *mystery as to why day does not code correctly for 11 cases - fix these, *which all are for the second day. if test=1 day=lag(day)+1. recode day (8=1). temp. select if diary=2. cro day by sic152 t21. compute year=1992. if month>9 year=1991. cro month by year. if civstat=2 sppart=-7. freq sppart. temp. select if civstat=2. freq sppart. Do If baddiary=0 and sex>-8 and age>-8 and day>-8. compute badcase=0. Else if baddiary=0 and (sex=-8 or age=-8) and day>-8. compute badcase=1. Else if baddiary=0 and sex>-8 and age>-8 and day=-8. compute badcase=2. Else if sex>-8 and age>-8 and day>-8 and nep=0 and miss=1 and misbasic=0. compute badcase=3. Else if sex>-8 and age>-8 and day>-8 and nep=1 and miss=0 and misbasic=0. compute badcase=4. Else if sex>-8 and age>-8 and day>-8 and nep=0 and miss=0 and misbasic=1. compute badcase=5. Else if (sex=-8 or age=-8) and day=-8 and nep=0 and miss=0 and misbasic=0. compute badcase=6. Else if (sex=-8 or age=-8) and day>-8 and nep=0 and miss=1 and misbasic=0. compute badcase=7. Else if (sex=-8 or age=-8) and day>-8 and nep=1 and miss=0 and misbasic=0. compute badcase=8. Else if (sex=-8 or age=-8) and day>-8 and nep=0 and miss=0 and misbasic=1. compute badcase=9. Else if sex>-8 and age>-8 and day=-8 and nep=0 and miss=1 and misbasic=0. compute badcase=10. Else if sex>-8 and age>-8 and day=-8 and nep=1 and miss=0 and misbasic=0. compute badcase=11. Else if sex>-8 and age>-8 and day=-8 and nep=0 and miss=0 and misbasic=1. compute badcase=12. Else if sex>-8 and age>-8 and day>-8 and nep=1 and miss=1 and misbasic=0. compute badcase=13. Else if sex>-8 and age>-8 and day>-8 and nep=0 and miss=1 and misbasic=1. compute badcase=14. Else if sex>-8 and age>-8 and day>-8 and nep=1 and miss=0 and misbasic=1. compute badcase=15. Else if (sex=-8 or age=-8) and day=-8 and nep=0 and miss=1 and misbasic=0. compute badcase=16. Else if (sex=-8 or age=-8) and day=-8 and nep=1 and miss=0 and misbasic=0. compute badcase=17. Else if (sex=-8 or age=-8) and day=-8 and nep=0 and miss=0 and misbasic=1. compute badcase=18. Else if (sex=-8 or age=-8) and day>-8 and nep=1 and miss=1 and misbasic=0. compute badcase=19. Else if (sex=-8 or age=-8) and day>-8 and nep=0 and miss=1 and misbasic=1. compute badcase=20. Else if (sex=-8 or age=-8) and day>-8 and nep=1 and miss=0 and misbasic=1. compute badcase=21. Else if sex>-8 and age>-8 and day=-8 and nep=1 and miss=1 and misbasic=0. compute badcase=22. Else if sex>-8 and age>-8 and day=-8 and nep=0 and miss=1 and misbasic=1. compute badcase=23. Else if sex>-8 and age>-8 and day=-8 and nep=1 and miss=0 and misbasic=1. compute badcase=24. Else if sex>-8 and age>-8 and day>-8 and nep=1 and miss=1 and misbasic=1. compute badcase=25. Else if (sex=-8 or age=-8) and day=-8 and nep=1 and miss=1 and misbasic=0. compute badcase=26. Else if (sex=-8 or age=-8) and day=-8 and nep=0 and miss=1 and misbasic=1. compute badcase=27. Else if (sex=-8 or age=-8) and day=-8 and nep=1 and miss=0 and misbasic=1. compute badcase=28. Else if (sex=-8 or age=-8) and day>-8 and nep=1 and miss=1 and misbasic=1. compute badcase=29. Else if sex>-8 and age>-8 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. *check that all cases have a code. freq badcase. *double check that all correctly coded. cro badcase by sex day nep miss misbasic. freq day. weight by m11. freq day. weight by m12. freq day. weight off. *m11, the inflation factor, does not correct for day distribution, *though the uninflated m11 does correct for day distribution. compute ocombwt = m12. compute agesexgp=sex*100+age. temp. select if badcase=0. freq agesexgp. weight by ocombwt. freq agesexgp. weight off. ******************************** *Distributions by age and sex * * weighted no * * by m12 weight * *men * *113 448 889 * *116 420 761 * *119 325 426 * *122 1129 883 * *127 1018 681 * *132 1395 1121 * *137 957 1396 * *142 1049 1364 * *147 896 1070 * *152 1048 1191 * *157 1001 946 * *162 905 678 * *167 700 383 * *172 424 210 * *175 512 270 * * * *women * *213 445 889 * *216 395 708 * *219 366 381 * *222 1069 743 * *227 1051 975 * *232 1300 1543 * *237 1000 1701 * *242 963 1452 * *247 894 1026 * *252 1128 1255 * *257 995 879 * *262 1008 612 * *267 1359 584 * *272 794 341 * *275 818 381 * ********************************. *constructed population weight =expected/actual. if agesexgp=113 agesexwt =( 448/ 889). if agesexgp=116 agesexwt =( 420/ 761). if agesexgp=119 agesexwt =( 325/ 426). if agesexgp=122 agesexwt =(1129/ 883). if agesexgp=127 agesexwt =(1018/ 681). if agesexgp=132 agesexwt =(1395/1121). if agesexgp=137 agesexwt =( 957/1396). if agesexgp=142 agesexwt =(1049/1364). if agesexgp=147 agesexwt =( 896/1070). if agesexgp=152 agesexwt =(1048/1191). if agesexgp=157 agesexwt =(1001/ 946). if agesexgp=162 agesexwt =( 905/ 678). if agesexgp=167 agesexwt =( 700/ 383). if agesexgp=172 agesexwt =( 424/ 210). if agesexgp=175 agesexwt =( 512/ 270). if agesexgp=213 agesexwt =( 445/ 889). if agesexgp=216 agesexwt =( 395/ 708). if agesexgp=219 agesexwt =( 366/ 381). if agesexgp=222 agesexwt =(1069/ 743). if agesexgp=227 agesexwt =(1051/ 975). if agesexgp=232 agesexwt =(1300/1543). if agesexgp=237 agesexwt =(1000/1701). if agesexgp=242 agesexwt =( 963/1452). if agesexgp=247 agesexwt =( 894/1026). if agesexgp=252 agesexwt =(1128/1255). if agesexgp=257 agesexwt =( 995/ 879). if agesexgp=262 agesexwt =(1008/ 612). if agesexgp=267 agesexwt =(1359/ 584). if agesexgp=272 agesexwt =( 794/ 341). if agesexgp=275 agesexwt =( 818/ 381). desc agesexwt. sort cases by agesexgp day. aggregate out='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\group.sav' /break=agesexgp /grouptot=sum(agesexwt)/ groupn=nu. execute. aggregate out='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\day.sav' /break= agesexgp day /daytot=sum(agesexwt) /dayn=nu. execute. match files file=* table='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\group.sav' /by agesexgp . execute. match files file=* table='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\day.sav' /by agesexgp day. execute. compute exptot=(grouptot/7). do if badcase=0. compute propwt=(exptot/daytot)/(dayn/daytot). else if badcase>0. compute propwt=0. end if. weight off. desc propwt ocombwt. weight by propwt. cro day by sex /cel=col. freq agesexgp day. weight off. *no need to adjust weight further as producing the correct number of cases. compute pwt=propwt*(25812/25739). desc pwt. weight by pwt. freq day. weight off. *attempting further adjustment to the weight makes the weight overestimate the *number of diaries. apply dictionary from 'C:\Kimberly\time\mtus\MTUS database\MTUSW58schema.sav'. val lab empinclm 1 under 300 2 300-600 3 600-1,000 4 1,000-1,400 5 1,400-1,800 6 1,800-2,200 7 2,200-2,500 8 2,500-3,000 9 3,000-3,500 10 3,500-4,000 11 4,000-4,500 12 4,500-5,000 13 5,000-5,500 14 5,500-6,000 15 6,000-6,500 16 6,500-7,000 17 7,000 and over 50 no data 99 no income. val lab incorig 1 'under 800' 2 '800-1,000' 3 '1,000-1,400' 4 '1,400-1,800' 5 '1,800-2,200' 6 '2,200-3,000' 7 '3,000-3,500' 8 '3,500-4,000' 9 '4,000-4,500' 10 '4,500-5,000' 11 '5,000-5,500' 12 '5,500-6,000' 13 '6,000-8,000' 14 '8,000-10,000' 15 '10,000 and over' 50 'no data'. val lab educa 101 'lower 2ndary leaving cert, no vocational training' 102 'lower 2ndary leaving cert & apprenticeship' 103 'lower 2ndary leaving cert & traineeship' 104 'lower 2ndary leaving cert & higher vocational diploma' 105 'lwr 2ndary leaving cert&fachschule, ddr profess diplma' 201 'intermediate 2ndary leaving cert & no vocational' 202 'intermediate 2ndary leaving cert & apprenticeship' 203 'intermediate 2ndary leaving cert & traineeship' 204 'intermediate 2ndary leaving cert & higher voc diplom' 205 'intermed 2ndary leaving cert&fachschule,ddrprofess dip' 301 'allgemeinbildendepolytechnischeoberschule,ddr,no voc' 302 'allgemeinbildendepolytechnischeoberschule,ddr,apprent' 303 'allgemeinbildendepolytechnischeoberschule,ddr,trnshp' 304 'allgemeinbildendepolytechnischeoberschule,ddr,hgr voc' 305 'allgemeinbildendepolytechnischeoberschule,ddr,f ddr' 401 'tech college entry-level leaving cert,no voc' 402 'tech college entry-level leaving cert &apprentship' 403 'tech college entry-level leaving cert & traineeship' 404 'tech college entry-level leaving cert & higher voc diploma' 405 'tech college entry-level leaving cert&f, ddr professdip' 406 'techcollege entry-level leaving cert&tech college dgree' 501 'university entry-level leaving cert, no vocational' 502 'university entry-level leaving cert & apprenticeship' 503 'university entry-level leaving cert & traineeship' 504 'university entry-level leaving cert & hghr voc diplma' 505 'university entry-level leaving cert & f, ddrprofess dip' 506 'university entry-level leaving cert&tec college degree' 507 'university entry-level leaving cert&university degree' 601 'still without school leaving cert, no voc' 602 'still without school leaving cert & apprenticeship' 604 'still without school leaving cert & higher voc diplma'. *recoding time use variables that cannot be created to -9. recode av14 av22 av26 av27 av39 main25 main34 main41 main45 main53 main60 main62 (0=-9). execute. sort cases by survey hldid persid id. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w553.sav' /keep countreg 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 edtry 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 ocombwt propwt. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58.sav' /keep countreg 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 edtry rushed health carer disab 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 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 sppart ocombwt propwt. *final quality checks. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w553.sav'. desc hldid persid 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. AGGREGATE /OUTFILE=* MODE=ADDVARIABLES /BREAK=countreg /wtsum = SUM(propwt). freq wtsum. freq countreg survey swave msamp persid id 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 edtry rushed health carer disab. cro empsp by civstat. 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. cro civstat by hhtype relrefp partid cohab. *note no relrefp=2 should appear if civstat=2. cro hhtype by hhldsize. cro parntid1 by parntid2. weight by propwt. freq day. cro day by sex /cells= count column. weight off. aggregate /outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\aggr.sav' /break=hldid /urbanmin = min(urban) /urbanmax = max(urban) /incmin = min(income) /incmax = max(income) /akidmin = min(agekidx) /akidmax = max(agekidx) /akid2min = min(agekid2) /akid2max = max(agekid2) /nchimin = min(nchild) /nchimax = max(nchild) /hsizemin = min(hhldsize) /hsizemax = max(hhldsize) /htypmin = min(hhtype) /htypmax = max(hhtype) /ownhomin = min(ownhome) /ownhomax = max(ownhome) /compmin = min(computer) /compmax = max(computer) /vehicmin = min(vehicle) /vehicmax = max(vehicle). get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\aggr.sav'. cro urbanmin by urbanmax. cro incmin by incmax. cro akidmin by akidmax. cro akid2min by akid2max. cro nchimin by nchimax. cro hsizemin by hsizemax. cro htypmin by htypmax. cro ownhomin by ownhomax. cro compmin by compmax. cro vehicmin by vehicmax. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58.sav'. desc main1 to main69 sppart. *check all variables have variable and value labels. Quality of other variables *already checked in the review of the W53 version. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58kid.sav'. select if age>17. freq age badcase. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58.sav'. get file='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58kid.sav'. select if age<18. freq age badcase. save outfile='C:\Kimberly\time\mtus\mtus data\W553\Germany\1992\ger1991w58kid.sav'.