** user=**** 
** password=**** 
** package=spss
** project=lis

set printback=listing .

* SPSS Slovakia LIS 1996 : Example CAMSIS matching .

get file= sk96h .

* Slovakia LIS 1996 : Example CAMSIS matching .
* (here to household file ).



compute isco88=-999 .
if (d14 > 999 & d14 < 10000) isco88=d14 .

* Current occ unit (isco88 one digit as 4 digit)=d14, 
* valid d14 in range > 999 and < 10000, 
* Employment status differences not applicable to Slovakia.


* Key indexing variable is isco88.
sort cases by isco88. 

temp.
select if (d14 > 999 & d14 < 10000).
means tables=isco88 by d3  .

sav out= cstemp .
get file="u:/camsis/csskis96.sav" . 
descriptives var=all.
sort cases by isco88.
match files table=* /in=csocc 
     /file=cstemp /in=source /by=isco88.
fre var=csocc source.
cro csocc by source.
select if (source=1).
* (this removes those CAMSIS isco88 by empst values which are not
*   represented in the LIS sample).
cro source by d3. 


* Current values are within gender. 
* also can make a cross-gender CAMSIS variable : .
compute cgentcs=-999.
if (d3=1) cgentcs=mcam .
if (d3=2) cgentcs=fcam .

descriptives var=isco88 cgentcs  .
temp.
select if (d14 > 999 & d14 < 10000).
descriptives var= isco88 cgentcs .


* Quick assessment: compare values by an education recode.
* some recode of education .
compute educ=d10. 

recode educ (1=1) (2,3,4=2) (5=3) (else=-9).
add value labels educ  
   1 "Elementary" 2 "Intermediate" 3 "University" . 
variable label educ "education categories".
fre var=educ .

temp.
select if (d14 > 999 & d14 < 10000).
descriptives var= educ cgentcs isei siops  . 

* Men .
temp. 
select if (d14 > 988  & d14 < 10000 & d3=1 & educ > 0).
means tables= cgentcs by educ  
             /statistics=anova.
temp. 
select if (d14 > 988  & d14 < 10000 & d3=1 & educ > 0 & siops > 0).
means tables= isei siops by educ  
             /statistics=anova.

* Women .
temp. 
select if (d14 > 999  & d14 < 10000 & d3=2 & educ > 0).
means tables= cgentcs  by educ 
             /statistics=anova.
temp. 
select if (d14 > 999  & d14 < 10000 & d3=2 & educ > 0 & siops > 0).
means tables= isei siops by educ  
             /statistics=anova.




/* For interest here are the Eta-squared statistics 
 with education from above:

Men: 
CAMSIS  : 0.408
ISEI  : 0.400
SIOPS : 0.408

Women: 
CAMSIS  :  0.345
ISEI : 0.367
SIOPS : 0.367
 
*/