ANALYZE Statement |
You use the ANALYZE statement to estimate the parameters of the probability distribution specified in the DISTRIBUTION statement without producing any graphical output. The ANALYZE statement performs the same analysis as the PROBPLOT statement, but it does not produce any plots. In addition, you can use the ANALYZE statement to analyze data with the binomial and Poisson distributions. The third format for the preceding ANALYZE statement applies only to Poisson and binomial data. You can use any number of ANALYZE statements after a PROC RELIABILITY statement; each ANALYZE statement produces a separate analysis. You can specify an optional label to distinguish between multiple ANALYZE statements in the output.
You must specify one variable. If your data are right censored, you must specify a censor-variable and, in parentheses, the values of the censor-variable that correspond to censored data values.
If you are using the binomial or Poisson distributions, you must specify variable1 to represent a binomial or Poisson count and variable2 to provide an exposure measure for the Poisson distribution or the binomial sample size for the binomial distribution.
You can optionally specify one or two group-variables. The ANALYZE statement produces an analysis for each level combination of the group-variable values. The observations in a given level are referred to as a cell.
The elements of the ANALYZE statement are described as follows.
represents the data for which an analysis is to be produced. A variable must be a numeric variable in the input data set.
indicates which observations in the input data set are right censored. You specify the values of censor-variable that represent censored observations by placing those values in parentheses after the variable name. If your data are not right censored, then you omit the specification of censor-variable; otherwise, censor-variable must be a numeric variable in the input data set.
is another method of specifying the data. You can use this syntax in a situation where uncensored, interval-censored, left-censored, and right-censored values occur in the same set of data. Table 14.27 shows how you use this syntax to specify different types of censoring by using combinations of missing and nonmissing values. See the section Lognormal Analysis with Arbitrary Censoring for an example of using this syntax to create a probability plot.
represents the count data for which a Poisson or binomial analysis is to be produced. A variable1 must be a numeric variable in the input data set.
provides either an exposure measure for a Poisson analysis or a binomial number of trials for a binomial analysis. A variable2 must be a numeric variable in the input data set.
are one or two group variables. If no group variables are specified, a single analysis is produced. The group-variables can be numeric or character variables in the input data set.
Note that the parentheses surrounding the textitgroup-variables are needed only if two group variables are specified.
control the features of the analysis. All options are specified after a slash (/) in the ANALYZE statement.
The following tables summarize the options available in the ANALYZE statement. You can specify one or more of these options to control the parameter estimation and provide optional analyses.
Option |
Option Description |
||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CONFIDENCE=number |
Specifies the confidence coefficient for all confidence intervals. Specify a number between 0 and 1. The default value is 0.95. |
||||||||||||||||
CONVERGE=number |
Specifies the convergence criterion for maximum likelihood fit. See the section Maximum Likelihood Estimation for details. |
||||||||||||||||
CONVH=number |
Specifies the convergence criterion for the relative Hessian convergence criterion. See the section Maximum Likelihood Estimation for details. |
||||||||||||||||
CORRB |
Requests the parameter correlation matrix. |
||||||||||||||||
COVB |
Requests the parameter covariance matrix. |
||||||||||||||||
FITTYPE | FIT= fit-specification |
Specifies the method of estimating distribution parameters. The available fit-specifications and their meanings are shown in the following table. |
||||||||||||||||
|
|||||||||||||||||
INEST | IN=SAS-data-set |
Specifies a SAS data set that can contain initial values, equality constraints, upper bounds, or lower bounds for the scale, shape, and threshold parameters in a three-parameter Weibull model for lifetime data. Applies only to three-parameter Weibull models. See the section INEST Data Set for the Three-Parameter Weibull for details. |
||||||||||||||||
ITPRINT |
Requests the iteration history for maximum likelihood fit. |
||||||||||||||||
ITPRINTEM |
Requests the iteration history for the Turnbull algorithm. |
||||||||||||||||
LRCL |
Requests likelihood ratio confidence intervals for distribution parameters. |
||||||||||||||||
LRCLPER |
Requests likelihood ratio confidence intervals for distribution percentiles. |
||||||||||||||||
LOCATION=number <LINIT> |
Specifies fixed or initial value of location parameter. |
||||||||||||||||
MAXIT=number |
Specifies the maximum number of iterations allowed for maximum likelihood fit. |
||||||||||||||||
MAXITEREM | |
|||||||||||||||||
MAXITEM=number1 <,number2> |
number1 specifies the maximum number of iterations allowed for Turnbull algorithm. Iteration history is printed in increments of number2 if requested with ITPRINTEM. See the section Interval-Censored Data for details. |
||||||||||||||||
NOPCTILES |
Suppresses computation of percentiles. |
||||||||||||||||
NOPOLISH |
Suppresses the setting of small interval probabilities to 0 in the Turnbull algorithm. See the section Interval-Censored Data for details. |
||||||||||||||||
PCTLIST=number-list |
Specifies a list of percentages for which to compute percentile estimates. number-list must be a list of numbers separated by blanks or commas. Each number in the list must be between 0 and 100. If this option is not specified, percentiles are computed for a standard list of percentages. |
||||||||||||||||
PPOS=plotting-position |
Specifies the plotting-position type used to compute nonparametric estimates of the probability distribution function. See the section Probability Plotting for details. The available plotting-position types are shown in the following table. |
||||||||||||||||
|
|||||||||||||||||
PPOUT |
Requests a table of cumulative probabilities. |
||||||||||||||||
PRINTPROBS |
Print intervals and associated probabilities for the Turnbull algorithm. |
||||||||||||||||
PROBLIST=number-list |
Specifies a list of initial values for Turnbull algorithm. See the section Interval-Censored Data for details. |
||||||||||||||||
PSTABLE=number |
Specifies stable parameterization. The number must be between 0 and 1. See the section Stable Parameters for further information. |
||||||||||||||||
READOUT |
Analyzes readout data. |
||||||||||||||||
SCALE=number <SCINIT> |
Specifies the fixed or initial value of scale parameter. |
||||||||||||||||
SHAPE=number <SHINIT> |
Specifies the fixed or initial value of shape parameter. |
||||||||||||||||
SINGULAR=number |
Specifies the singularity criterion for matrix inversion. |
||||||||||||||||
Requests that the survival function, cumulative distribution function, and confidence limits be computed for values in number-list. See the section Reliability Function for details. |
|||||||||||||||||
THRESHOLD=number |
Specifies a fixed threshold parameter. See Table 14.51 for the distributions with a threshold parameter. |
||||||||||||||||
TOLLIKE=number |
Specifies the criterion for convergence in the Turnbull algorithm. The default is . See the section Interval-Censored Data for details. |
||||||||||||||||
TOLPROB=number |
Specifies the criterion for setting interval probability to 0 in the Turnbull algorithm. Default is . See the section Interval-Censored Data for details. |
Option |
Option Description |
---|---|
CONFIDENCE=number |
Specifies the confidence coefficient for all confidence intervals. Specify a number between 0 and 1. The default value is 0.95. |
PREDICT(number) |
Requests predicted counts for exposure number for Poisson or sample size number for binomial. |
TOLERANCE(number) |
Requests exposure for Poisson or sample size for binomial to estimate Poisson rate or binomial probability within number with probability given by the CONFIDENCE= option. |