* Encoding: UTF-8.
* seas2021.sps.

*  Establishing drinking/abstainer status (F).

recode F_1b (10 THRU 11 = 0) (1 thru 9 = 1) into consumers.
EXECUTE.

VARIABLE LABELS consumers 'alcohol consumers'.

VALUE LABELS consumers 1 'alc consumers' 0 'abstiners'.



* Generic frequncy of drinking (F). 

recode  f_1b (0=0) (1=365) (2=286) (3=182) (4=78) (5=30) (6=12) (7=8.5) (8=3.5)
                   (9=1) (10=0) (11=0) (sysmis=0) INTO rf_1b.
EXECUTE .

Variable labels rf_1b 'generic frequency of drinking (number of days in year)'.


*******  Beverage Specific Quantity-Frequency (BSQF) ********.

recode bsqf_1b bsqf_3b bsqf_5b bsqf_7b (0=0) (1=365) (2=286) (3=182) (4=78) (5=30) (6=12) (7=8.5) (8=3.5)
                   (9=1) (10=0) (11=0) (sysmis=0) INTO rbsqf_1 rbsqf_3 rbsqf_5 rbsqf_7.
EXECUTE .

Variable labels rbsqf_1 'frequency of beer drinking (number of days in year)'
                        /rbsqf_3 'frequency of wine drinking (number of days in year)' 
                        /rbsqf_5 'frequency of spirits drinking (number of days in year)'
                        /rbsqf_7 'frequency of drinking of additional beverage (number of days in year)'. 

COMPUTE bsqf_beer = bsqf_2a *  rbsqf_1 .  
COMPUTE bsqf_wine = bsqf_4a *  rbsqf_3 .  
COMPUTE bsqf_spir = bsqf_6a *  rbsqf_5 .  
COMPUTE bsqf_adbev = bsqf_8a *  rbsqf_7 . 

COMPUTE bsqf_alc = bsqf_beer + bsqf_wine + bsqf_spir + bsqf_adbev.
EXECUTE .

Variable labels bsqf_beer 'beer annual consumption in cl of 100% alc'
                        /bsqf_wine 'wine annual consumption in cl of 100% alc' 
                        /bsqf_spir 'spirits annual consumption in cl of 100% alc'
                        /bsqf_adbev 'additional beverage annual consumption in cl of 100% alc'                  
                        /bsqf_alc 'alcohol annual consumption in cl of 100% alc'.  


IF (consumers eq 1 and bsqf_alc eq 0) bsqf_alc = -1. 
EXECUTE.
recode bsqf_alc (-1 = sysmis) (ELSE = copy). 
EXECUTE.

RECODE bsqf_alc (18250 THRU hi = 18250) (ELSE = COPY) into bsqf_alc1. 
EXECUTE.

VARIABLE LABELS bsqf_alc1 'alcohol annual consumption in cl of 100% alc - after capping (18250 cl. 100% alc)'. 


* Capping extreme values.

recode bsqf_2a bsqf_4a bsqf_6a bsqf_8a (sysmis=0). 
EXECUTE .


* Capping alcohol consumption per day (BSQF) 50-80 cl capped to 50 cl; 80 cl and highest put into sysmis. 

RECODE bsqf_2a (50 THRU 80 = 50) (80 thru Highest=sysmis)  (Lowest thru 50=copy) INTO cbsqf_2.
RECODE bsqf_4a  (50 THRU 80 = 50) (80 thru Highest=sysmis) (Lowest thru 50=copy) INTO cbsqf_4.
RECODE bsqf_6a (50 THRU 80 = 50) (80 thru Highest=sysmis)  (Lowest thru 50=copy) INTO cbsqf_6.
RECODE bsqf_8a  (50 THRU 80 = 50) (80 thru Highest=sysmis) (Lowest thru 50=copy) INTO cbsqf_8.
EXECUTE.

VARIABLE LABELS cbsqf_2 'beer consumption per day in 100% alc - capped'
    /cbsqf_4 'wine consumption per day in 100% alc - capped'
    /cbsqf_6 'spirits consumption per day in 100% alc - capped'
    /cbsqf_8 'additional beverage consumption per day in 100% alc - capped'.

COMPUTE cbsqf_beer = cbsqf_2 *  rbsqf_1 .  
COMPUTE cbsqf_wine = cbsqf_4 *  rbsqf_3 .  
COMPUTE cbsqf_spir = cbsqf_6 *  rbsqf_5 .  
COMPUTE cbsqf_adbev = cbsqf_8 *  rbsqf_7 . 
EXECUTE .

recode cbsqf_beer cbsqf_wine cbsqf_spir cbsqf_adbev (SYSMIS = 0). 
EXECUTE .

COMPUTE cbsqf_alc = cbsqf_beer + cbsqf_wine + cbsqf_spir + cbsqf_adbev.
EXECUTE .

Variable labels cbsqf_beer 'beer annual consumption in cl of 100% alc - capped'
                        /cbsqf_wine 'wine annual consumption in cl of 100% alc - capped'
                        /cbsqf_spir 'spirits annual consumption in cl of 100% alc - capped'
                        /cbsqf_alc 'alcohol annual consumption in cl of 100% alc - capped' 
                       /cbsqf_adbev 'additional beverages annual consumption in cl of 100% alc - capped'.

* Second capping BSQF - capping cbsqf on the level of 18250 cl of 100% alcohol. All above 18250 cl set to this figure. 

IF (consumers eq 1 and cbsqf_alc eq 0) cbsqf_alc = -1. 
EXECUTE.
recode cbsqf_alc (-1 = sysmis) (ELSE = copy). 
EXECUTE.

IF (sysmis(cbsqf_2) or sysmis(cbsqf_4) or sysmis(cbsqf_6) or sysmis(cbsqf_8)) cbsqf_alc = -1. 
EXECUTE.

MISSING VALUES cbsqf_alc (-1). 
EXECUTE.

RECODE cbsqf_alc (18250 THRU hi = 18250) (ELSE = COPY) into c1cbsqf_alc . 
EXECUTE.

VARIABLE LABELS c1cbsqf_alc 'alcohol annual consumption in cl of 100% alc with firts capping - after second capping (18250 cl. 100% alc)'. 


***** Risky Single Occasion Drinking (RSOD) *****.

* RSOD - conversion of frequency scale into number of days. 

recode rsod_2b rsod_7b (0=0) (1=365) (2=286) (3=182) (4=78) (5=30) (6=12) (7=8.5) (8=3.5)
                   (9=1) (10=0) (sysmis=0) INTO rrsod_2b rrsod_7b.
EXECUTE .

* Capping max vol on one occasion RSOD_1.

recode RSOD_1a (sysmis=sysmis) (lo THRU 400 = copy) (400 THRU hi = 400) 
            into cRSOD_1a. 
EXECUTE .

Variable labels cRSOD_1a 'Max vol of 100% alc in grams -capped'.


* Gender (SD_1_gender)

compute SD_1_gender = SD_1. 
if (SD_1 eq 3 and RSOD_1_Other eq 1) SD_1_gender = 1. 
if (SD_1 eq 3 and RSOD_1_Other eq 2) SD_1_gender = 2. 
EXECUTE.


*****  Social class status (SD_10 - SD_18) *********.

recode sd_11 sd_15 sd_17 (1 THRU 2=1) (4=2) (3=3) (5=3) (6=4) (7=5) (9=9) (0=0) into rsd_11 rsd_15 rsd_17.
EXECUTE.

recode sd_13 (0=1) (1=2) (2=3) (3=4) (4=5) (5=6) (9=9) into rsd_13. 

MISSING VALUES  rsd_13 (9).

COMPUTE social_class = rsd_11.
EXECUTE.

if (social_class EQ 0) social_class = rsd_15.
EXECUTE.

if (social_class EQ 0) social_class = rsd_17.
EXECUTE.

recode social_class (SYSMIS = 9) (0=9). 
EXECUTE.


**** Household income (SD_20, SD_20a). 

RECODE SD_20A ( 0 = SYSMIS )
 ( 1 = 750 )
 ( 2 = 1250 )
 ( 3 = 1750 )
 ( 4 = 2250 )
 ( 5 = 2750 )
 ( 6 = 3250 )
 ( 7 = 3750 )
 ( 8 = 4250 )
 ( 9 = 4750 )
 ( 10 = 5250 )
 ( 11 = 5750 )
 ( 12 = 9000 )
 ( 13 = 375 )
 ( 14 = 750 )
 ( 15 = 1250 )
 ( 16 = 1750 )
 ( 17 = 2250 )
 ( 18 = 2750 )
 ( 19 = 3250 )
 ( 20 = 3750 )
 ( 21 = 4250 )
 ( 22 = 4750 )
 ( 23 = 5250 )
 ( 24 = 8250 )
 ( 25 = 750 )
 ( 26 = 1500 )
 ( 27 = 2500 )
 ( 28 = 3500 )
 ( 29 = 4500 )
 ( 30 = 7500 )
 ( 31 = 750 )
 ( 32 = 1500 )
 ( 33 = 2500 )
 ( 34 = 3500 )
 ( 35 = 4500 )
 ( 36 = 7500 )
 ( 47 = 187.5 )
 ( 48 = 300.5 )
 ( 49 = 425.5 )
 ( 50 = 625.5 )
 ( 51 = 875.5 )
 ( 52 = 1100.5 )
 ( 53 = 1350.5 )
 ( 54 = 1750.5 )
 ( 55 = 2250.5 )
 ( 56 = 2750.5 )
 ( 57 = 4500 )
 ( 63 = 11250 )
 ( 64 = 22500.5 )
 ( 65 = 37500.5 )
 ( 66 = 52500.5 )
 ( 67 = 67500.5 )
 ( 68 = 82500.5 )
 ( 69 = 105000.5 )
 ( 70 = 135000 )
 ( 71 = 225000 )
 ( 72 = 74999.25 )
 ( 73 = 149999.5 )
 ( 74 = 249999.5 )
 ( 75 = 349999.5 )
 ( 76 = 449999.5 )
 ( 77 = 549999.5 )
 ( 78 = 649999.5 )
 ( 79 = 749999.5 )
 ( 80 = 849999.5 )
 ( 81 = 949999.5 )
 ( 82 = 1500000 )
 ( 83 = 675 )
 ( 84 = 1250.5 )
 ( 85 = 2050.5 )
 ( 86 = 2750.5 )
 ( 87 = 4500 )
 ( 88 = 1500 )
 ( 89 = 3000.5 )
 ( 90 = 5000.5 )
 ( 91 = 7000.5 )
 ( 92 = 9000.5 )
 ( 93 = 11000.5 )
 ( 94 = 14000.5 )
 ( 95 = 18000.5 )
 ( 96 = 30001.5 )
 ( 97 = 6750 )
 ( 98 = 13500.5 )
 ( 99 = 21000.5 )
 ( 100 = 27000.5 )
 ( 101 = 45000 )
 ( 112 = 15000 )
 ( 113 = 27500.5 )
 ( 114 = 42500.5 )
 ( 115 = 67500.5 )
 ( 116 = 127500 )
 ( 117 = 1500 )
 ( 118 = 2250 )
 ( 119 = 2750 )
 ( 120 = 3750 )
 ( 121 = 5250 )
 ( 122 = 9000 )
 ( 123 = 7499.25 )
 ( 124 = 14999.5 )
 ( 125 = 24999.5 )
 ( 126 = 34999.5 )
 ( 127 = 44999.5 )
 ( 128 = 54999.5 )
 ( 129 = 64999.5 )
 ( 130 = 74999.5 )
 ( 131 = 84999.5 )
 ( 132 = 94999.5 )
 ( 133 = 150000 )
 ( 134 = 562.5 )
 ( 135 = 1125.5 )
 ( 136 = 1750.5 )
 ( 137 = 2250.5 )
 ( 138 = 3751.5 )
 ( 139 = 41250 )
 ( 140 = 67500 )
 ( 141 = 95000.5 )
 ( 142 = 125000.5 )
 ( 143 = 155000.5 )
 ( 144 = 185000.5 )
 ( 145 = 225000.5 )
 ( 146 = 275000.5 )
 ( 147 = 450000 )
 ( 154 = 7499.25 )
 ( 155 = 14999.5 )
 ( 156 = 24999.5 )
 ( 157 = 34999.5 )
 ( 158 = 44999.5 )
 ( 159 = 54999.5 )
 ( 160 = 64999.5 )
 ( 161 = 74999.5 )
 ( 162 = 84999.5 )
 ( 163 = 94999.5 )
 ( 164 = 150000 )
 ( 165 = 450 )
 ( 166 = 750.5 )
 ( 167 = 1050.5 )
 ( 168 = 1350.5 )
 ( 169 = 1650.5 )
 ( 170 = 1950.5 )
 ( 171 = 2250.5 )
 ( 172 = 2550.5 )
 ( 173 = 2850.5 )
 ( 174 = 3250.5 )
 ( 175 = 4000.5 )
 ( 176 = 5250.5 )
 ( 177 = 6750.5 )
 ( 178 = 11250 )
 ( 179 = 1500 )
 ( 180 = 3000.5 )
 ( 181 = 5000.5 )
 ( 182 = 7000.5 )
 ( 183 = 9000.5 )
 ( 184 = 11000.5 )
 ( 185 = 14000.5 )
 ( 186 = 18000.5 )
 ( 187 = 30001.5 )
 ( 188 = 1500 )
 ( 189 = 3000.5 )
 ( 190 = 5000.5 )
 ( 191 = 7000.5 )
 ( 192 = 9000.5 )
 ( 193 = 11000.5 )
 ( 194 = 14000.5 )
 ( 195 = 18000.5 )
 ( 196 = 30001.5 )
 ( 198 = 1125.5 )
 ( 199 = 1375.5 )
 ( 200 = 1750.5 )
 ( 201 = 2250.5 )
 ( 202 = 2750.5 )
 ( 203 = 3250.5 )
 ( 204 = 3750.5 )
 ( 205 = 4500.5 )
 ( 206 = 5500.5 )
 ( 207 = 7000.5 )
 ( 208 = 12000 )
 ( 209 = 750 )
 ( 212 = 1750.5 )
 ( 213 = 2250.5 )
 ( 214 = 2750.5 )
 ( 215 = 3250.5 )
 ( 216 = 3750.5 )
 ( 217 = 4500.5 )
 ( 218 = 5500.5 )
 ( 219 = 7000.5 )
 ( 250 = 150000 )
 ( 251 = 374500 )
 ( 252 = 524500 )
 ( 253 = 699500 )
 ( 254 = 899500 )
 ( 255 = 1500000 )
 ( 256 = 3450 )
 ( 257 = 5450 )
 ( 258 = 7150 )
 ( 259 = 8550 )
 ( 260 = 9950 )
 ( 261 = 11650 )
 ( 262 = 13400 )
 ( 263 = 14850 )
 ( 264 = 16250 )
 ( 265 = 18550 )
 ( 266 = 21700 )
 ( 267 = 24800 )
 ( 268 = 29350 )
 ( 269 = 35650 )
 ( 270 = 45050 )
 ( 271 = 58150 )
 ( 272 = 71250 )
 ( 273 = 90650 )
 ( 274 = 116600 )
 ( 275 = 142250 )
 ( 276 = 168200 )
 ( 277 = 193850 )
 ( 278 = 219500 )
 ( 279 = 245750 )
 ( 280 = 271700 )
 ( 281 = 297600 )
 ( 282 = 466050 )
 ( 283 = 1125 )
 ( 284 = 1700 )
 ( 285 = 2100.5 )
 ( 286 = 2525.5 )
 ( 287 = 2975.5 )
 ( 288 = 3400.5 )
 ( 289 = 3875.5 )
 ( 290 = 4575.5 )
 ( 291 = 5750.5 )
 ( 292 = 9750 )
 ( 293 = 5249.25 )
 ( 294 = 8499.5 )
 ( 295 = 11499.5 )
 ( 296 = 14499.5 )
 ( 297 = 17999.5 )
 ( 298 = 21499.5 )
 ( 299 = 24499.5 )
 ( 300 = 27999.5 )
 ( 301 = 32499.5 )
 ( 302 = 37499.5 )
 ( 303 = 42499.5 )
 ( 304 = 47500 )
 ( 305 = 75000 )
 ( 306 = 120000 )
 ( 307 = 70000.5 )
 ( 308 = 55000.5 )
 ( 309 = 47500.5 )
 ( 310 = 42500.5 )
 ( 311 = 37500.5 )
 ( 312 = 32500.5 )
 ( 313 = 29000.5 )
 ( 314 = 26500.5 )
 ( 315 = 24000.5 )
 ( 316 = 22000.5 )
 ( 317 = 20500.5 )
 ( 318 = 18500.5 )
 ( 319 = 16000.5 )
 ( 320 = 12500.5 )
 ( 321 = 7500 )
 ( 322 = 525 )
 ( 323 = 850 )
 ( 324 = 1250.5 )
 ( 325 = 2000.5 )
 ( 326 = 3000.5 )
 ( 327 = 5250 )
 ( 328 = 12000 )
 ( 329 = 20500.5 )
 ( 330 = 30000.5 )
 ( 331 = 40000.5 )
 ( 332 = 52500.5 )
 ( 333 = 75000.5 )
 ( 334 = 105000.5 )
 ( 335 = 135000.5 )
 ( 336 = 225000 )
 ( 337 = 6000 )
 ( 338 = 8999.5 )
 ( 339 = 10999.5 )
 ( 340 = 12999.5 )
 ( 341 = 15499.5 )
 ( 342 = 17999.5 )
 ( 343 = 20499.5 )
 ( 344 = 23499.5 )
 ( 345 = 27499.5 )
 ( 346 = 32499.5 )
 ( 347 = 38499.5 )
 ( 348 = 45999.5 )
 ( 349 = 57499.5 )
 ( 350 = 74499.5 )
 ( 351 = 126000 )
 ( 352 = 675 )
 ( 353 = 1250.5 )
 ( 354 = 2050.5 )
 ( 355 = 2750.5 )
 ( 356 = 4500 )
 ( 357 = 675 )
 ( 358 = 1250.5 )
 ( 359 = 2050.5 )
 ( 360 = 2750.5 )
 ( 361 = 4500 )
 ( 401 = 187.5 )
 ( 402 = 375 )
 ( 403 = 650.5 )
 ( 404 = 900.5 )
 ( 405 = 1125.5 )
 ( 406 = 1375.5 )
 ( 407 = 1625.5 )
 ( 408 = 1875.5 )
 ( 409 = 2250.5 )
 ( 410 = 3751.5 )
 ( 601 = 675 )
 ( 602 = 1250.5 )
 ( 603 = 2050.5 )
 ( 604 = 2750.5 )
 ( 605 = 4500 )
 ( 1001 = 210000 )
 ( 1002 = 350000 )
 ( 1003 = 485000 )
 ( 1004 = 645000 )
 ( 1005 = 825000 )
 ( 1006 = 1380000 )
 ( 2302 = 1125.5 )
 ( 2303 = 1375.5 )
 ( 2304 = 1750.5 )
 ( 2305 = 2250.5 )
 ( 2306 = 2750.5 )
 ( 2307 = 3250.5 )
 ( 2308 = 3750.5 )
 ( 2309 = 4500.5 )
 ( 2310 = 5500.5 )
 ( 2311 = 7000.5 )
 ( 2312 = 12000 )
 ( 2401 = 675 )
 ( 2402 = 1250.5 )
 ( 2403 = 2050.5 )
 ( 2404 = 2750.5 )
 ( 2405 = 4500 )
 ( 2501 = 225 )
 ( 2502 = 400.5 )
 ( 2503 = 750.5 )
 ( 2504 = 1250.5 )
 ( 2505 = 1750.5 )
 ( 2506 = 2500.5 )
 ( 2507 = 3500.5 )
 ( 2508 = 4500.5 )
 ( 2509 = 5500.5 )
 ( 2510 = 6500.5 )
 ( 2511 = 7500.5 )
 ( 2512 = 12001.5 )
 ( 9999 = SYSMIS )
 INTO SD_20B.
EXECUTE.


RECODE SD_20   (0=SYSMIS)  (9999 = SYSMIS) (99999 = SYSMIS) (999999 = SYSMIS) (999999999=SYSMIS).
EXECUTE.

* lower range.

if (COUNTRY eq 10 and SD_20 le 375)  SD_20 = 375. 
if (COUNTRY eq 11 and SD_20 le 750)  SD_20 = 750. 
if (COUNTRY eq 12 and SD_20 le 187.5)  SD_20 = 187.5. 
if (COUNTRY eq 13 and SD_20 le 187.5)  SD_20 = 187.5. 
if (COUNTRY eq 14 and SD_20 le 675)  SD_20 = 675. 
if (COUNTRY eq 15 and SD_20 le 11250)  SD_20 = 11250. 
if (COUNTRY eq 16 and SD_20 le 675)  SD_20 = 675. 
if (COUNTRY eq 17 and SD_20 le 7500)  SD_20 = 7500. 
if (COUNTRY eq 18 and SD_20 le 74999.25)  SD_20 = 74999.25. 
if (COUNTRY eq 19 and SD_20 le 1500)  SD_20 = 1500. 
if (COUNTRY eq 20 and SD_20 le 15000)  SD_20 = 15000. 

if (COUNTRY eq 21 and SD_20 le 1500)  SD_20 = 1500. 
if (COUNTRY eq 22 and SD_20 le 750)  SD_20 = 750. 
if (COUNTRY eq 23 and SD_20 le 562.5)  SD_20 = 562.5. 
if (COUNTRY eq 24 and SD_20 le 41250)  SD_20 = 41250. 
if (COUNTRY eq 25 and SD_20 le 2520000)  SD_20 = 2520000. 
if (COUNTRY eq 26 and SD_20 le 7499.25)  SD_20 = 7499.25. 
if (COUNTRY eq 27 and SD_20 le 450)  SD_20 = 450. 
if (COUNTRY eq 28 and SD_20 le 1500)  SD_20 = 1500. 
if (COUNTRY eq 29 and SD_20 le 1500)  SD_20 = 1500. 
if (COUNTRY eq 30 and SD_20 le 750)  SD_20 = 750. 

if (COUNTRY eq 31 and SD_20 le 675)  SD_20 = 675. 
if (COUNTRY eq 32 and SD_20 le 225)  SD_20 = 225. 
if (COUNTRY eq 33 and SD_20 le 3450)  SD_20 = 3450. 
if (COUNTRY eq 34 and SD_20 le 150000)  SD_20 = 150000. 
if (COUNTRY eq 35 and SD_20 le 1125)  SD_20 = 1125. 
if (COUNTRY eq 36 and SD_20 le 5249.25)  SD_20 = 5249.25. 
if (COUNTRY eq 37 and SD_20 le 525)  SD_20 = 525. 
if (COUNTRY eq 38 and SD_20 le 12000)  SD_20 = 12000. 
if (COUNTRY eq 39 and SD_20 le 1500)  SD_20 = 1500. 
if (COUNTRY eq 40 and SD_20 le 6750)  SD_20 = 6750. 

if (COUNTRY eq 41 and SD_20 le 675)  SD_20 = 675. 
if (COUNTRY eq 42 and SD_20 le 6000)  SD_20 = 6000. 
if (COUNTRY eq 43 and SD_20 le 7499.25)  SD_20 = 7499.25. 
EXECUTE.

* higer range.

if (COUNTRY eq 10 and SD_20 ge 8250)  SD_20 = 8250. 
if (COUNTRY eq 11 and SD_20 ge 7500)  SD_20 = 7500. 
if (COUNTRY eq 12 and SD_20 ge 3751.5)  SD_20 = 3750.5. 
if (COUNTRY eq 13 and SD_20 ge 4500)  SD_20 = 4500. 
if (COUNTRY eq 14 and SD_20 ge 4500)  SD_20 = 4500. 
if (COUNTRY eq 15 and SD_20 ge 225000)  SD_20 = 225000. 
if (COUNTRY eq 16 and SD_20 ge 4500)  SD_20 = 4500. 
if (COUNTRY eq 17 and SD_20 ge 120000)  SD_20 = 120000. 
if (COUNTRY eq 18 and SD_20 ge 1500000)  SD_20 = 1500000. 
if (COUNTRY eq 19 and SD_20 ge 30001.5)  SD_20 = 30001.5. 
if (COUNTRY eq 20 and SD_20 ge 127500)  SD_20 = 127500. 

if (COUNTRY eq 21 and SD_20 ge 9000)  SD_20 = 9000. 
if (COUNTRY eq 22 and SD_20 ge 9000)  SD_20 = 9000. 
if (COUNTRY eq 23 and SD_20 ge 3751.5)  SD_20 = 3751.5. 
if (COUNTRY eq 24 and SD_20 ge 450000)  SD_20 = 450000. 
if (COUNTRY eq 25 and SD_20 ge 16560000)  SD_20 = 16560000. 

if (COUNTRY eq 26 and SD_20 ge 150000)  SD_20 = 150000. 
if (COUNTRY eq 27 and SD_20 ge 11250)  SD_20 = 11250. 
if (COUNTRY eq 28 and SD_20 ge 30001.5)  SD_20 = 30001.5. 
if (COUNTRY eq 29 and SD_20 ge 30001.5)  SD_20 = 30001.5. 
if (COUNTRY eq 30 and SD_20 ge 12000)  SD_20 = 12000. 

if (COUNTRY eq 31 and SD_20 ge 4500)  SD_20 = 4500. 
if (COUNTRY eq 32 and SD_20 ge 12001.5)  SD_20 = 12001.5. 
if (COUNTRY eq 33 and SD_20 ge 466050)  SD_20 = 466050. 
if (COUNTRY eq 34 and SD_20 ge 1500000)  SD_20 = 11500000. 
if (COUNTRY eq 35 and SD_20 ge 9750)  SD_20 = 9750. 
if (COUNTRY eq 36 and SD_20 ge 75000)  SD_20 = 75000. 
if (COUNTRY eq 37 and SD_20 ge 5250)  SD_20 = 5250. 
if (COUNTRY eq 38 and SD_20 ge 225000)  SD_20 = 225000. 
if (COUNTRY eq 39 and SD_20 ge 30001.5)  SD_20 = 30001.5. 
if (COUNTRY eq 40 and SD_20 ge 45000)  SD_20 = 45000. 

if (COUNTRY eq 41 and SD_20 ge 4500)  SD_20 = 4500. 
if (COUNTRY eq 42 and SD_20 ge 126000)  SD_20 = 126000. 
if (COUNTRY eq 43 and SD_20 ge 150000)  SD_20 = 150000. 
EXECUTE.

COMPUTE SD_20C = SD_20.
EXECUTE.
if (sysmis(SD_20)) SD_20C = SD_20B.
EXECUTE.

variable labels SD_20B 'Household income calculated on the basis of categories'
  /SD_20C 'Integrated household income (SD_20 and SD_20B)'.


********* Income per capita *******************.

* creationg variable with number of haushold members. 

recode SD_6A (sysmis = 1) (else = copy) into SD_6Anumber. 
EXECUTE.

VARIABLE LABELS SD_6Anumber 'Number of haushold members'.

* recalculation of yearly income into mothly income. 

COMPUTE sd_20month = SD_20C.
IF (SD_20_income_type2 = 1) sd_20month = SD_20C/12. 
EXECUTE.

* correction for Iceland (SD_20 - yearly, SD_20B monthly) . 

DO IF (COUNTRY = 25 and sysmis(SD_20)).
COMPUTE sd_20month = SD_20B.
END IF.
EXECUTE.


* calculating monthly income per capita of haushold mambers. 
COMPUTE inc_mth = sd_20month/SD_6Anumber.  
EXECUTE.

VARIABLE LABELS sd_20month 'Monthly household income final'
    /inc_mth 'Monthly income per capita of household members final'.


* correction for Slovenia (everything to sysmis) . 

DO IF (COUNTRY = 40).
RECODE SD_20 SD_20A SD_20B sd_20month inc_mth SD_20C (SYSMIS=SYSMIS) (ELSE=SYSMIS).
END IF.
EXECUTE.


* Variable filter. 

compute SD_20_incom_type = SD_20_income_type2.
recode SD_20_incom_type (1=1) (else = 0). 
EXECUTE.
VARIABLE LABELS SD_20_incom_type 'Yearly - monthly'.
VALUE LABELS SD_20_incom_type 1 'yearly' 0 'monthly'. 


* correction for Iceland.

DO IF (COUNTRY = 25 and sysmis(SD_20)).
COMPUTE SD_20_incom_type = 0.
END IF.
EXECUTE.