data slid1 ; infile 'c:\temp\SLID-Ontario.txt' firstobs = 2 delimiter = '09'x ; input age sex $ wages yearsed ; male = (sex eq 'Male') ; log2wages = log2(wages) ; * 2^(mean(log2wages)) = 2^3.7798 = 13.74 ; g = wages*(log(wages/13.74) - 1) ; obs = _n_ ; title 'SLID data'; run ; proc print data = slid1 (obs = 10) ; run ; proc means data = slid1 ; var log2wages ; run ; ods html ; ods graphics on; proc transreg data = slid1 ss2 plots = all; model BoxCox(wages / lambda=-.4 to .4 by 0.02) = identity(male age yearsed)/ pboxcoxtable ; run; * Here the constructed variable approach is illustrated ; proc reg data = slid1 ; model wages = male age yearsed g / partial ; run ; ods graphics off ; ods html close ;