** user=**** 
** password=**** 
** package=stata
** project=les

#delimit;

* US LES 1997 : STATA Example CAMSIS matching to US 97 data;
*(note that this uses different index var and educ from us90);
*(index var is soc80 and is only at title only level;


use $us97e;
* Current occ unit (soc80's)=fj04; 
* in this example only, can ignore status diffs; 
*[no status derivation];



* Derive key indexing variable soc80;
gen soc80=-999;
replace soc80=fj04 if (fj04 > 0 );
tabulate db02 if (fj04 > 0 ), 
       summarize(soc80) mean standard freq ;



* Match on CAMSIS scores from separate file;
sort soc80;
save cstemp, replace;
use soc80 htcs wtcs siops isei 
      using $mydata\camsis\csuses97;
summarize;
sort soc80;
merge soc80 using cstemp;
gen occinfo=_merge;
tab occinfo; 
drop if (occinfo==1);
* (this removes those iscobst values which are not represented in the LES sample);
tab occinfo db02; 




* Current values are within gender; 
* also can make a cross-gender CAMSIS variable : ;
gen cgentcs=-999;
replace cgentcs=htcs if (db02==1);
replace cgentcs=wtcs if (db02==2);
gen cgenscs=-999;
*replace cgenscs=hscs if (db02==1);
*replace cgenscs=wscs if (db02==2);



* Review scores distributions;
summarize soc80 cgentcs cgenscs siops isei ;
summarize soc80 cgentcs cgenscs siops isei if fj04 > 0;
summarize soc80 cgentcs cgenscs siops isei if (fj04 > 0 & siops > 0); 
* note that the last summary shows us that a handful of cases have occs 
*  but missing camsis and siopsisei scores by the CAMSIS / ISMF translation;




* Quick assessment: compare values by a US education recode;
gen educ=et01; 
recode educ 0 31/34=1 35/40=2 41/46=3 ;
label define edu  1 "Primary or less" 2 "Intermediate" 
 3 "Degree level"; 
label values educ edu;
label variable educ "education categories";
tab educ ;
sort db02;
by db02:
  summarize educ 
      cgentcs cgenscs siops isei if fj04 > 0; 
by db02: 
 oneway cgentcs educ  if (fj04 > 0 & siops > 0), tabulate;

by db02: 
 oneway siops educ  if (fj04 > 0 & siops > 0), tabulate;
by db02: 
 oneway isei educ  if (fj04 > 0 & isei > 0), tabulate;




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

Men: 
CAMSIS title only : 0.
CAMSIS title by status : 0.
SIOPS : 0.
ISEI: 0.
Women: 
CAMSIS title only : 0.
CAMSIS title by status : 0.
SIOPS : 0.
ISEI: 0.
 
*/