data vocab1 ; infile 'c:\temp\vocabulary.txt' firstobs = 2 delimiter = '09'x ; input id year sex $ education vocabulary ; run ; proc print data = vocab1 (obs = 45) ; run ; proc freq data = vocab1 ; tables year sex education ; run ; * These two programs yield the information in Table 12.2 ; proc reg data = vocab1 ; model vocabulary = education ; run ; proc glm data = vocab1 ; class education ; model vocabulary = education ; run ; * The MEANS statement with the options below gives the Levene test from section 12.4.2 ; proc glm data = vocab1 ; class education ; model vocabulary = education ; means education / hovtest = levene (type = abs) ; run ;