data duncanprestige ; infile 'c:\temp\duncan.txt' firstobs = 2 ; length occupation $ 20 ; input occupation $ type $ income education prestige ; obs = _n_ ; title 'Prestige data'; run ; proc print data = duncanprestige (obs = 10) ; run ; proc means data = duncanprestige ; run ; * This code illustrates the use of partial regression plots with the PARTIAL option ; proc reg data = duncanprestige ; model prestige = education income / partial ; output out = rduncanprestige p = pred r = res student = stres rstudent = jkres h = hat cookd = cook covratio = covr dffits = dffits lclm=lcimu uclm=ucimu lcl=lpi ucl=upi ; run ; proc print data = rduncanprestige ; run ; proc plot data = rduncanprestige vpercent = 80 ; plot jkres*pred hat*obs jkres*hat ; proc capability noprint data = rduncanprestige lineprinter ; var jkres ; qqplot jkres /normal(mu = est sigma = est symbol='.') square ; run ;