If you develop your SAS/IML programs in SAS/IML Studio, you can use high-level statistical graphics. For example, the following statements create three scatter plots that duplicate the low-resolution plots created in the previous section. Two of the plots are shown in Figure 4.11. The main steps in the program are indicated by numbered comments; these steps are explained in the list that follows the program.
x = {1 1 1, 1 2 4, 1 3 9, 1 4 16, 1 5 25}; /* 1 */ y = {1, 5, 9, 23, 36}; x1 = x[,2]; /* data = second column of X */ x = orpol(x1,2); /* generates orthogonal polynomials */ run Regress; /* runs the Regress module */ declare DataObject dobj; /* 2 */ dobj = DataObject.Create("Reg", /* 3 */ {"x" "y" "Residuals" "Predicted"}, x1 || y || resid || yhat); declare ScatterPlot p1, p2, p3; p1 = ScatterPlot.Create(dobj, "x", "Residuals"); /* 4 */ p1.SetTitleText("Plot of Residuals", true); p2 = ScatterPlot.Create(dobj, "x", "Predicted"); /* 5 */ p2.SetTitleText("Plot of Predicted Values", true); p3 = ScatterPlot.Create(dobj, "x", "y"); /* 6 */ p3.SetTitleText("Scatter Plot with Regression Line", true); p3.DrawUseDataCoordinates(); p3.DrawLine(x1,yhat); /* 7 */
To completely understand this program, you should read SAS/IML Studio for SAS/STAT Users. The following list describes the main steps of the program:
Use SAS/IML to create the data and run the Regress module.
Specify that the dobj
variable is an object of the DataObject class. SAS/IML Studio extends the SAS/IML language by adding object-oriented programming
techniques.
Create an object of the DataObject class from SAS/IML vectors.
Create a scatter plot of the residuals versus the values of the explanatory variable.
Create a scatter plot of the predicted values versus the values of the explanatory variable.
Create a scatter plot of the observed responses versus the values of the explanatory variable.
Overlay a line for the predicted values.