goptions device=lj4sips gaccess=gsasfile autofeed reset=all gunit=pct border ftext=swissb htitle=6 htext=3; options ls=64 ps=45 nodate nonumber; options linesize = 90 nodate;data bolts; input diameter break; lackofit=diameter; cards; 0.1 1.62 0.1 1.73 0.1 1.7 0.1 1.66 0.1 1.74 0.1 1.72 0.2 1.71 0.2 1.78 0.2 1.79 0.2 1.86 0.2 1.7 0.2 1.84 0.3 1.86 0.3 1.86 0.3 1.9 0.3 1.95 0.3 1.96 0.3 2 0.4 2.14 0.4 2.07 0.4 2.11 0.4 2.18 0.4 2.17 0.4 2.07 0.5 2.45 0.5 2.42 0.5 2.33 0.5 2.36 0.5 2.38 0.5 2.31 ; TITLE1 ' The Breaking Strengths'; Title2 'Breaking Strength vs. Diameters'; proc reg; model break = diameter/XPX I P CLI CLM; OUTPUT OUT = HOLD R=RESID P=PRED; run; /* define symbol characteristics */ symbol interpol=none /*regression analysis with */ /*symbol interpol=rlclm95 regression analysis with */ /* confidence limits */ value=dot /* plot symbol */ height=1 /* plot symbol height */ cv=red /* plot symbol color */ ci=blue /* regression line color */ co=green /* confidence limits color */ width=2; /* line width */ /* produce plot */ PROC GPLOT DATA=HOLD; TITLE3 'SCATTERPLOTS'; PLOT break*diameter; PLOT RESID*diameter RESID*break; run; %include 'nqplot.sas'; TITLE3 'Normal Probability Plots'; %nqplot(data=hold,var=resid) run; proc glm; classes lackofit; model break = diameter lackofit /ss1; title1 'Testing for lack of fit'; run;