Hello,

I'm trying to run an estimation for my honours thesis detailed in this paper (open-source version): https://www.lifecoursecentre.org.au/...Siminski-2.pdf, though with a longer dataset.

I am attempting to replicate equation 3 (page 13), however am unsure on how to save the scalars for occupation from equation 2 (page 12).

My initial regression is as follows:

reg income i.occupation age age^2

I have estimated the coefficients of different occupations by creating tabulated dummies with the data, however where creating a scalar for the constant, age and age^2 is simple enough, the function does not work with the variable i.occupation.

Does anyone know how I can save the coefficients of these dummies? Many thanks in advance!