The SEVSELECT Procedure
ODS Graphics
Statistical procedures use ODS Graphics to create graphs as part of their output. ODS Graphics is described in detail in the "Statistical Graphics Using ODS" chapter in SAS/STAT User's Guide.
Before you create graphs, ODS Graphics must be enabled (for example, with the ODS GRAPHICS ON statement). For more information about enabling and disabling ODS Graphics, see the section "Enabling and Disabling ODS Graphics" in that chapter.
The overall appearance of graphs is controlled by ODS styles. Styles and other aspects of using ODS Graphics are discussed in the section "A Primer on ODS Statistical Graphics" in that chapter.
This section describes how the SEVSELECT procedure uses ODS to create graphics.
ODS Graph Names
PROC SEVSELECT assigns a name to each graph that it creates by using ODS. You can use these names to selectively reference the graphs. The names are listed in Table 14.8.
Table 14.8: ODS Graphics Produced by PROC SEVSELECT
ODS Graph Name | Plot Description | PLOTS= Option |
|---|---|---|
CDFDistPlot | CDF plot per distribution | CDFPERDIST |
CDFPlot | Comparative CDF plot | CDF |
PDFDistPlot | PDF plot per distribution | PDFPERDIST |
PDFPlot | Comparative PDF plot | |
PPPlot | P-P plot of CDF and EDF | PP |
QQPlot | Q-Q plot |
Comparative CDF Plot
The comparative CDF plot helps you visually compare the cumulative distribution function (CDF) estimates of all the candidate distribution models to the empirical distribution function (EDF) estimate. The plot does not contain CDF estimates for models whose parameter estimation process does not converge. The horizontal axis represents the values of the response variable. The vertical axis represents the values of the CDF or EDF estimates.
If you specify truncation, then conditional CDF estimates are plotted. Otherwise, unconditional CDF estimates are plotted. The conditional estimates are computed by using the method that is described in the section Truncation and Conditional CDF Estimates.
If you specify regression effects, then the plotted CDF estimates are those of a mixture distribution. For more information, see the section CDF and PDF Estimates with Regression Effects.
CDF Plot per Distribution
The CDF plot per distribution shows the CDF estimates of each candidate distribution model unless that model’s parameter estimation process does not converge. The plot also contains estimates of the EDF. The horizontal axis represents the values of the response variable. The vertical axis represents the values of the CDF or EDF estimates.
This plot shows the lower and upper pointwise confidence limits for the EDF estimates. For an EDF estimate with standard error , they are computed as and , respectively, where is the pth quantile from the standard normal distribution and denotes the confidence level that you specify in the EDFALPHA= option (the default is ).
If you specify truncation, then conditional CDF estimates are plotted. Otherwise, unconditional CDF estimates are plotted. The conditional estimates are computed by using the method that is described in the section Truncation and Conditional CDF Estimates.
If you specify regression effects, then the plotted CDF estimates are those of a mixture distribution. For more information, see the section CDF and PDF Estimates with Regression Effects.
Comparative PDF Plot
The comparative PDF plot helps you visually compare the probability density function (PDF) estimates of all the candidate distribution models. The plot does not contain PDF estimates for models whose parameter estimation process does not converge. The horizontal axis represents the values of the response variable. The vertical axis represents the values of the PDF estimates.
If you specify the HISTOGRAM option, then the plot also contains the histogram of response variable values. If you specify the KERNEL option, then the plot also contains the kernel density estimate of the response variable values.
If you specify regression effects, then the plotted PDF estimates are those of a mixture distribution. For more information, see the section CDF and PDF Estimates with Regression Effects.
PDF Plot per Distribution
The PDF plot per distribution shows the PDF estimates of each candidate distribution model unless that model’s parameter estimation process does not converge. The horizontal axis represents the values of the response variable. The vertical axis represents the values of the PDF estimates.
If you specify the HISTOGRAM option, then the plot also contains the histogram of response variable values. If you specify the KERNEL option, then the plot also contains the kernel density estimate of the response variable values.
If you specify regression effects, then the plotted PDF estimates are those of a mixture distribution. For more information, see the section CDF and PDF Estimates with Regression Effects.
P-P Plot of CDF and EDF
The P-P plot of CDF and EDF is the probability-probability plot that compares the CDF estimates of a distribution to the EDF estimates. A plot is not prepared for models whose parameter estimation process does not converge. The horizontal axis represents the CDF estimates of a candidate distribution, and the vertical axis represents the EDF estimates.
This plot can be interpreted as displaying the data that are used for computing the EDF-based statistics of fit for the given candidate distribution. As described in the section EDF-Based Statistics, these statistics are computed by comparing the EDF, denoted by , to the CDF, denoted by , at each of the response variable values y. Using the probability inverse transform , this is equivalent to comparing the EDF of z, denoted by , to the CDF of z, denoted by (D’Agostino and Stephens 1986, Ch. 4). Because the CDF of z is a uniform distribution (), the EDF-based statistics can be computed by comparing the estimate of the EDF of z to the estimate of z itself, which is the estimate of the CDF. The horizontal axis of the plot represents . The vertical axis represents the estimated EDF of z, . The plot contains a scatter plot of (, ) points and a reference line, , that represents the expected uniform distribution of z. Points that are scattered closer to the reference line indicate a better fit than points that are scattered farther away from the reference line.
If you specify truncation, then the EDF estimates are conditional, as described in the section EDF Estimates and Truncation. So conditional estimates of CDF are displayed; these estimates are computed by using the method that is described in the section Truncation and Conditional CDF Estimates.
If you specify regression effects, then the displayed CDF estimates, both unconditional and conditional, are those of a mixture distribution. For more information, see the section CDF and PDF Estimates with Regression Effects.
Q-Q Plot
The Q-Q plot is a quantile-quantile scatter plot that compares the empirical quantiles to the quantiles from a candidate distribution. A plot is not prepared for models whose parameter estimation process does not converge. The horizontal axis represents the quantiles from a candidate distribution, and the vertical axis represents the empirical quantiles.
Each point in the plot corresponds to a specific value of the EDF estimate, . The Y coordinate is the value of the response variable for which is computed. The X coordinate is computed by using one of the two following methods for a candidate distribution named dist:
If you have defined the dist_QUANTILE function that satisfies the requirements listed in the section <phrase remap="Argument">dist</phrase>_QUANTILE, then that function is invoked by using and estimated distribution parameters as arguments. The QUANTILE function is defined in the
Sashelp.Svrtdistlibrary for all the predefined distributions.If the dist_QUANTILE function is not defined, then PROC SEVSELECT numerically inverts the dist_CDF function at the CDF value of for the estimated distribution parameters. If the dist_CDF function is not defined, then the exp(dist_LOGCDF) function is inverted. If the inversion fails, the corresponding point is not plotted in the Q-Q plot.
If you specify truncation, then the EDF estimates are conditional, as described in the section EDF Estimates and Truncation. The CDF inversion process, whether done numerically or by evaluating the dist_QUANTILE function, needs to accept an unconditional CDF value. So the value is first transformed to an unconditional estimate as
where and are as defined in the section Truncation and Conditional CDF Estimates.
If you specify regression effects, then the value of the first distribution parameter is determined by using the method that is equivalent to specifying the DFMIXTURE=MEAN option in the SCALEMODEL statement. This is described in the section CDF and PDF Estimates with Regression Effects.