The HPNLMOD Procedure

PROC HPNLMOD Statement

  • PROC HPNLMOD <options>;

The PROC HPNLMOD statement invokes the procedure. Table 1 summarizes important options in the PROC HPNLMOD statement by function. These and other options in the PROC HPNLMOD statement are then described fully in alphabetical order.

Table 1: PROC HPNLMOD Statement Options

Option Description
Basic Options
DATA= Specifies the input data set
OUT= Specifies the output data set
Output Options
CORR Specifies the correlation matrix
COV Specifies the covariance matrix
ECORR Specifies the correlation matrix of additional estimates
ECOV Specifies the covariance matrix of additional estimates
DF Specifies the default degrees of freedom
NOPRINT Suppresses ODS output
NOITPRINT Suppresses output about iterations within the optimization process
Optimization Options
ABSCONV= Tunes an absolute function convergence criterion
ABSFCONV= Tunes an absolute difference function convergence criterion
ABSGCONV= Tunes the absolute gradient convergence criterion
FCONV= Tunes the relative function convergence criterion
GCONV= Tunes the relative gradient convergence criterion
MAXITER= Chooses the maximum number of iterations in any optimization
MAXFUNC= Specifies the maximum number of function evaluations in any optimization
MAXTIME= Specifies the upper limit seconds of CPU time for any optimization
MINITER= Specifies the minimum number of iterations in any optimization
TECHNIQUE= Selects the optimization technique
Tolerance Options
SINGULAR= Tunes the general singularity criterion
User-Defined Format Options
FMTLIBXML= Specifies a file reference for a format stream
XMLFORMAT= Specifies a file name for a format stream


You can specify the following options in the PROC HPNLMOD statement.

ABSCONV=r
ABSTOL=r

specifies an absolute function convergence criterion. For minimization, termination requires f left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis less-than-or-equal-to sans-serif-italic r, where bold-italic psi is the vector of parameters in the optimization and f left-parenthesis dot right-parenthesis is the objective function. The default value of sans-serif-italic r is the negative square root of the largest double-precision value, which serves only as a protection against overflow.

ABSFCONV=r <n>
ABSFTOL=r<n>

specifies an absolute difference function convergence criterion. For all techniques except the Nelder-Mead simplex (NMSIMP) technique, termination requires a small change of the function value in successive iterations:

StartAbsoluteValue f left-parenthesis bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis Baseline right-parenthesis minus f left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis EndAbsoluteValue less-than-or-equal-to r

Here, bold-italic psi denotes the vector of parameters that participate in the optimization, and f left-parenthesis dot right-parenthesis is the objective function. The same formula is used for the NMSIMP technique, but bold-italic psi left-parenthesis k right-parenthesis is defined as the vertex that has the lowest function value, and bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis is defined as the vertex that has the highest function value in the simplex. The default value is sans-serif-italic r equals 0. The optional integer value sans-serif-italic n specifies the number of successive iterations for which the criterion must be satisfied before the process can be terminated.

ABSGCONV=r <n>
ABSGTOL=r<n>

specifies an absolute gradient convergence criterion. Termination requires the maximum absolute gradient element to be small:

max Underscript j Endscripts StartAbsoluteValue g Subscript j Baseline left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis EndAbsoluteValue less-than-or-equal-to r

Here, bold-italic psi denotes the vector of parameters that participate in the optimization, and g Subscript j Baseline left-parenthesis dot right-parenthesis is the gradient of the objective function with respect to the jth parameter. This criterion is not used by the NMSIMP technique. The default value is sans-serif-italic r equals1Eminus5. The optional integer value sans-serif-italic n specifies the number of successive iterations for which the criterion must be satisfied before the process can be terminated.

ALPHA=alpha

specifies the level of significance alpha that is used in the construction of 100 left-parenthesis 1 minus alpha right-parenthesis% confidence intervals. The value must be strictly between 0 and 1; the default value of alpha equals 0.05 results in 95% intervals. This value is used as the default confidence level for limits that are computed in the "Parameter Estimates" table and is used in the LOWER and UPPER options in the PREDICT statement.

CORR

requests the approximate correlation matrix for the parameter estimates.

COV

requests the approximate covariance matrix for the parameter estimates.

DATA=SAS-data-set

names the SAS data set to be used by PROC HPNLMOD. The default is the most recently created data set.

DF=n

specifies the default number of degrees of freedom to use in the calculation of p-values and confidence limits for additional parameter estimates.

ECORR

requests the approximate correlation matrix for all expressions that are specified in ESTIMATE statements.

ECOV

requests the approximate covariance matrix for all expressions that are specified in ESTIMATE statements.

FCONV=r<n>
FTOL=r<n>

specifies a relative function convergence criterion. For all techniques except NMSIMP, termination requires a small relative change of the function value in successive iterations:

StartFraction StartAbsoluteValue f left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis minus f left-parenthesis bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis Baseline right-parenthesis EndAbsoluteValue Over StartAbsoluteValue f left-parenthesis bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis Baseline right-parenthesis EndAbsoluteValue EndFraction less-than-or-equal-to r

Here, bold-italic psi denotes the vector of parameters that participate in the optimization, and f left-parenthesis dot right-parenthesis is the objective function. The same formula is used for the NMSIMP technique, but bold-italic psi Superscript left-parenthesis k right-parenthesis is defined as the vertex that has the lowest function value, and bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis is defined as the vertex that has the highest function value in the simplex. The default is sans-serif-italic r equals 10 Superscript minus normal upper F normal upper D normal upper I normal upper G normal upper I normal upper T normal upper S, where FDIGITS is by default minus log Subscript 10 Baseline left-brace epsilon right-brace and epsilon is the machine precision. The optional integer value sans-serif-italic n specifies the number of successive iterations for which the criterion must be satisfied before the process can terminate.

FMTLIBXML=file-ref

specifies the file reference for the XML stream that contains the user-defined format definitions. For information about how to generate a XML stream for your formats, see the section Working with Formats.

GCONV=r<n>
GTOL=r<n>

specifies a relative gradient convergence criterion. For all techniques except the conjugate gradient (CONGRA) and NMSIMP techniques, termination requires that the normalized predicted function reduction be small:

StartFraction bold g left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis prime left-bracket bold upper H Superscript left-parenthesis k right-parenthesis Baseline right-bracket Superscript negative 1 Baseline bold g left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis Over StartAbsoluteValue f left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis EndAbsoluteValue EndFraction less-than-or-equal-to r

Here, bold-italic psi denotes the vector of parameters that participate in the optimization, f left-parenthesis dot right-parenthesis is the objective function, and bold g left-parenthesis dot right-parenthesis is the gradient. For the CONGRA technique (where a reliable Hessian estimate bold upper H is not available), the following criterion is used:

StartFraction parallel-to bold g left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis parallel-to Subscript 2 Superscript 2 Baseline parallel-to bold s left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis parallel-to Over parallel-to bold g left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis minus bold g left-parenthesis bold-italic psi Superscript left-parenthesis k minus 1 right-parenthesis Baseline right-parenthesis parallel-to Subscript 2 Baseline StartAbsoluteValue f left-parenthesis bold-italic psi Superscript left-parenthesis k right-parenthesis Baseline right-parenthesis EndAbsoluteValue EndFraction less-than-or-equal-to r

This criterion is not used by the NMSIMP technique. The default value is sans-serif-italic r equals1Eminus8. The optional integer value sans-serif-italic n specifies the number of successive iterations for which the criterion must be satisfied before the process can terminate.

MAXFUNC=n
MAXFU=n

specifies the maximum number of function calls in the optimization process. The default values are as follows, depending on the optimization technique (which you specify in the TECHNIQUE= option):

  • TRUREG, NRRIDG, NEWRAP: sans-serif-italic n equals 125

  • QUANEW, DBLDOG: sans-serif-italic n equals 500

  • CONGRA: sans-serif-italic n equals 1,000

  • NMSIMP: sans-serif-italic n equals 3,000

Optimization can terminate only after completing a full iteration. Therefore, the number of function calls that are actually performed can exceed n.

MAXITER=n
MAXIT=n

specifies the maximum number of iterations in the optimization process. The default values are as follows, depending on the optimization technique (which you specify in the TECHNIQUE= option):

  • TRUREG, NRRIDG, NEWRAP: sans-serif-italic n equals 50

  • QUANEW, DBLDOG: sans-serif-italic n equals 200

  • CONGRA: sans-serif-italic n equals 400

  • NMSIMP: sans-serif-italic n equals 1,000

These default values also apply when sans-serif-italic n is specified as a missing value.

MAXTIME=r

specifies an upper limit of sans-serif-italic r seconds of CPU time for the optimization process. The default value is the largest floating-point double representation of your computer. This time that is specified by r is checked only once at the end of each iteration. Therefore, the actual running time can be longer than r.

MINITER=n
MINIT=n

specifies the minimum number of iterations. The default value is 0. If you request more iterations than are actually needed for convergence to a stationary point, the optimization algorithms can behave strangely. For example, the effect of rounding errors can prevent the algorithm from continuing for the required number of iterations.

NOITPRINT

suppresses the display of the "Iteration History" table.

NOPRINT

suppresses the generation of ODS output.

OUT=SAS-data-set

names the SAS data set to be created when one or more PREDICT statements are specified. A single OUT= data set is created to contain all predicted values when more than one PREDICT statement is specified. An error message is produced if a PREDICT statement is specified and an OUT= data set is not specified. For more information about output data sets in SAS high-performance analytical procedures, see the section Output Data Sets.

SINGULAR=number

tunes the general singularity criterion that is applied in sweeps and inversions. The default is 1 normal upper E Baseline 4 times the machine epsilon; this product is approximately 1Eminus12 on most computers.

TECHNIQUE=keyword
TECH=keyword

specifies the optimization technique for obtaining maximum likelihood estimates. You can choose from the following techniques by specifying the appropriate keyword:

CONGRA

performs a conjugate-gradient optimization.

DBLDOG

performs a version of double-dogleg optimization.

LEVMAR

performs a Levenberg-Marquardt optimization.

NEWRAP

performs a Newton-Raphson optimization that combines a line-search algorithm with ridging.

NMSIMP

performs a Nelder-Mead simplex optimization.

NONE

performs no optimization.

NRRIDG

performs a Newton-Raphson optimization with ridging.

QUANEW

performs a quasi-Newton optimization.

TRUREG

performs a trust-region optimization.

The default value is TECHNIQUE=LEVMAR for least squares regression models and TECHNIQUE=NRRIDG for models where the distribution is specified.

XMLFORMAT=filename

specifies the file name for the XML stream that contains the user-defined format definitions. For information about how to generate a XML stream for your formats, see the section Working with Formats.

Last updated: December 09, 2022