* Data from "Analysis of Messy Data, Vol. 1" by Milliken and Johnson 2nd edition page 477 ; * see the diagram for this lecture to show the layout of the experiment ; data splitb ; input plot nitro irrig y @@ ; cards ; 1 1 1 55 1 1 2 71 1 3 1 69 1 3 2 78 1 2 1 62 1 2 2 77 2 3 2 81 2 3 1 77 2 1 2 77 2 1 1 63 2 2 2 79 2 2 1 66 3 2 1 70 3 2 2 78 3 3 1 79 3 3 2 80 3 1 1 63 3 1 2 77 4 3 1 76 4 3 2 79 4 2 1 66 4 2 2 76 4 1 1 65 4 1 2 75 ; proc print data = splitb ; run ; proc mixed data = splitb cl covtest method = type3 ; class plot nitro irrig ; model y = nitro irrig nitro*irrig / ddfm = kr ; random plot nitro*plot irrig*plot ; lsmeans nitro irrig nitro*irrig / diff ; run ; proc mixed data = splitb cl covtest method = reml ; class plot nitro irrig ; model y = nitro irrig nitro*irrig / ddfm = kr outp = rsplitb ; random plot nitro*plot irrig*plot ; lsmeans nitro irrig nitro*irrig / diff ; run ; proc print data = rsplitb ; run ; proc plot data = rsplitb vpercent = 80 ; plot resid*pred ; run ; proc capability noprint data = rsplitb lineprinter ; var resid ; qqplot resid /normal(mu = est sigma = est symbol='.') square ; run ;