/* set the graphics environment */ goptions reset=global gunit=pct border ftext=swissb htitle=6 htext=3; data turbine; input x y ; label x='Gas Temperature' y='Heat, BTU per Kilowatt hours'; cards; 100 99.1 150 98.5 200 98.2 250 98.0 300 97.8 350 97.6 400 97.5 450 97.0 500 96.8 ; /* define symbol characteristics */ /*symbol interpol=none regression analysis with */ symbol1 interpol=none value=diamond height=3 cv=red ci=blue co=green width=2; symbol2 interpol=join value=dot height=3 cv=blue ci=blue co=green width=2; symbol3 interpol=none value=dot height=3 cv=green ci=blue co=green width=2; proc reg; model y = x /all; output out=hold p=pred r=resid; run; proc gplot data=hold; plot y*x pred*x / overlay; plot resid*pred/vref=0 cvref=blue; run; /*find the normal probability quantiles for the ratios*/ proc rank data=hold out=ranks; var resid ; ranks residrks; data n1; set ranks; r1=residrks/10; /*add one to sample size*/ z1=probit(r1); label z1='normal score for residuals'; proc gplot data=n1; title 'normal probability plot of residuals'; plot resid*z1; run;