data table177 ; input block treatment $ y x @@ ; cards ; 1 Check 23.1 13 1 Full 30.1 7 1 N0 26.4 10 1 P0 26.2 8 2 Check 20.9 12 2 Full 31.8 5 2 N0 27.2 9 2 P0 25.3 9 3 Check 28.3 7 3 Full 32.4 6 3 N0 28.6 6 3 P0 29.7 7 4 Check 25.0 9 4 Full 30.6 7 4 N0 28.5 6 4 P0 26.0 7 5 Check 25.1 8 5 Full 27.5 9 5 N0 30.8 5 5 P0 24.9 9 ; run ; proc print data = table177 ; run ; proc plot data = table177 ; plot y*x = treatment ; run ; proc glm data = table177 ; class block treatment ; model y = x block treatment ; lsmeans treatment / stderr ; lsmeans treatment / adjust = dunnett pdiff = control ('Full') cl ; output out = r_table177 p = pred r = res rstudent = jkres ; run ; proc plot data = r_table177 ; plot res*pred jkres*pred ; proc capability noprint data = r_table177 lineprinter ; var res ; qqplot res /normal(mu = est sigma = est symbol='.') square ; run ;