The MODEL Procedure

Iteration History

The options ITPRINT, ITDETAILS, XPX, I, and ITALL specify a detailed listing of each iteration of the minimization process.

ITPRINT Option

The ITPRINT information is selected whenever any iteration information is requested.

The following information is displayed for each iteration:

N

is the number of usable observations.

Objective

is the corrected objective function value.

Trace(S)

is the trace of the S matrix.

subit

is the number of subiterations required to find a lamda or a damping factor that reduces the objective function.

R

is the R convergence measure.

The estimates for the parameters at each iteration are also printed.

ITDETAILS Option

The additional values printed for the ITDETAILS option are as follows:

Theta

is the angle in degrees between normal upper Delta, the parameter change vector, and the negative gradient of the objective function.

Phi

is the directional derivative of the objective function in the normal upper Delta direction scaled by the objective function.

Stepsize

is the value of the damping factor used to reduce normal upper Delta if the Gauss-Newton method is used.

Lambda

is the value of lamda if the Marquardt method is used.

Rank(XPX)

is the rank of the bold upper X prime bold upper X matrix (output if the projected Jacobian crossproducts matrix is singular).

The definitions of PPC and R are explained in the section Convergence Criteria. When the values of PPC are large, the parameter associated with the criteria is displayed in parentheses after the value.

XPX and I Options

The XPX and the I options select the printing of the augmented bold upper X prime bold upper X matrix and the augmented bold upper X prime bold upper X matrix after a sweep operation (Goodnight 1979) has been performed on it. An example of the output from the following statements is shown in Figure 36:

proc model data=test2;
   y1 = a1 * x2 * x2 - exp( d1*x1);
   y2 = a2 * x1 * x1 + b2 * exp( d2*x2);
   fit y1 y2 / itall XPX I ;
run;

Figure 36: XPX and I Options Output

The MODEL Procedure
OLS Estimation

Cross Products for System At OLS Iteration 0
 a1d1a2b2d2Residual
a11839468-33818.350.00.000.0000003879959
d1-338181276.450.00.000.000000-76928
a200.0042925.01275.150.154739470686
b200.001275.250.010.00386716055
d200.000.20.000.0000642
Residual3879959-76928.14470686.316055.072.32971824576144

XPX Inverse for System At OLS Iteration 0
 a1d1a2b2d2Residual
a10.0000010.0000280.0000000.00000.002
d10.0000280.0015270.0000000.00000.00-9
a20.0000000.0000000.000097-0.0025-0.086
b20.0000000.000000-0.0024550.08250.95172
d20.0000000.000000-0.0849150.947615746.7111931
Residual1.952150-8.5468755.823969171.623411930.8910819902


The first matrix, labeled "Cross Products," for OLS estimation is

StartLayout 1st Row  Start 2 By 2 Matrix 1st Row 1st Column bold upper X prime bold upper X 2nd Column bold upper X prime bold r 2nd Row 1st Column bold r prime bold upper X 2nd Column bold r prime bold r EndMatrix EndLayout

The column labeled Residual in the output is the vector bold upper X prime bold r, which is the gradient of the objective function. The diagonal scalar value bold r prime bold r is the objective function uncorrected for degrees of freedom. The second matrix, labeled "XPX Inverse," is created through a sweep operation on the augmented bold upper X prime bold upper X matrix to get

StartLayout 1st Row  Start 2 By 2 Matrix 1st Row 1st Column left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript negative 1 2nd Column left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript negative 1 Baseline bold upper X prime bold r 2nd Row 1st Column left-parenthesis bold upper X prime bold r right-parenthesis prime left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript negative 1 2nd Column bold r prime bold r minus left-parenthesis bold upper X prime bold r right-parenthesis prime left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript negative 1 Baseline bold upper X prime bold r EndMatrix EndLayout

Note that the residual column is the change vector used to update the parameter estimates at each iteration. The corner scalar element is used to compute the R convergence criteria.

ITALL Option

The ITALL option, in addition to causing the output of all of the preceding options, outputs the S matrix, the inverse of the S matrix, the CROSS matrix, and the swept CROSS matrix. An example of a portion of the CROSS matrix for the preceding example is shown in Figure 37.

Figure 37: ITALL Option Crossproducts Matrix Output

The MODEL Procedure
OLS Estimation

Crossproducts Matrix At OLS Iteration 0
 1@PRED.y1/@a1@PRED.y1/@d1@PRED.y2/@a2@PRED.y2/@b2@PRED.y2/@d2RESID.y1RESID.y2
150.006409-239.161275.050.000.0038031470016053
@PRED.y1/@a16409.081839468-33818.35187766.16409.880.81393438799594065028
@PRED.y1/@d1-239.16-338181276.45-7253.0-239.19-0.026177-76928-85084
@PRED.y2/@a21275.00187766-7253.0042925.01275.150.154739420583470686
@PRED.y2/@b250.006410-239.191275.250.010.0038671470216055
@PRED.y2/@d20.001-0.030.20.000.00006422
RESID.y114699.973879959-76928.14420582.914701.771.8203561182710212234106
RESID.y216052.764065028-85083.68470686.316055.072.3297181223410612749042


Last updated: December 20, 2021