use V3 V51 ISCO88 using "C:\Users\User\work\data sets\issp - work orientations\work orientations iii (2005)\ZA4350_F1.dta", clear
// Generate variables
gen jobsatisf = 7 - V51
iskoisei isei, isko(ISCO88)
drop V51 ISCO88
preserve
// Calculating intercept and slope variation using OLS
statsby inter = _b[_cons] ///
slope = _b[isei] ///
, by(V3) ///
saving(ols, replace): ///
regress jobsatisf isei
merge m:1 V3 using ols
drop _merge
// Visualizing intercept and slope variation
gen yhat_ols = inter + slope*isei
separate jobsatisf, by(V3)
separate yhat_ols, by(V3)
twoway (line yhat_ols1-yhat_ols43 isei, sort(V3 isei)) ///
(lfit jobsatisf isei, clwidth(vvthick) clcolor(black)) ///
, legend(off) ytitle("Job satisfaction") xtitle("ISEI") ///
xlabel(16 25 50 75 90) ylabel(,format(%6.1f)) ///
caption("{it:Source:} ISSP 2005 (Work Orientations III), own calculations", span) ///
name(one, replace)
restore
// Calculating intercept and slope variation using a random effects model
center isei, inplace
mixed jobsatisf isei || V3: isei, var
predict u1 u0, reffects
// Visualizing intercept and slope variation
gen predRandomSlope = (_b[_cons] + u0) + ((_b[isei] + u1) * isei)
twoway (line predRandomSlope isei, connect(ascending) sort(V3 isei)), ///
ytitle("Job satisfaction") ///
xtitle("ISEI (centered)") ///
xlabel(-25 0 25 50) ///
ylabel(,format(%6.1f)) ///
caption("{it:Source:} ISSP 2005 (Work Orientations III), own calculations", span) ///
name(two, replace)
