Yes that is a lot clearer now. Thank you very much for your help! ... View more

Thank you Rick that's really useful. I have managed to get the residuals out using outpm = and have now recreated the QQ plot using PROC UNIVARIATE which is very helpful ods graphics on; proc mixed data=data; class treatment; model x_level = treatment / residual outpm = resid; run; ods graphics off; proc univariate data = resid; var PearsonResid; histogram PearsonResid / normal(sigma=1); QQplot PearsonResid / normal(mu=0 sigma=1); run; I wanted to plot the residuals as a standard normal probability plot, so have used the code above. The line plotted using this matches up with the one plotted in PROC MIXED so I was just wondering does PROC MIXED always plot a line y = x for any normal data? ... View more