One of the most important uses of PROC MODEL is to estimate unknown parameters in a nonlinear model. A simple nonlinear model has the form:

where x is a vector of exogenous variables. To estimate unknown parameters by using PROC MODEL, do the following:
Use the DATA= option in a PROC MODEL statement to specify the input SAS data set that contains and , the observed values of the variables.
Write the equation for the model by using SAS programming statements, including all parameters and arithmetic operators but leaving off the unobserved error component, .
Use a FIT statement to fit the model equation to the input data to determine the unknown parameters, .
The SASHELP library contains the data set CITIMON, which contains the variable LHUR, the monthly unemployment figures, and the variable IP, the monthly industrial production index. You suspect that the unemployment rates are inversely proportional to the industrial production index. Assume that these variables are related by the following nonlinear equation:

In this equation a, b, and c are unknown coefficients and is an unobserved random error.
The following statements illustrate how to use PROC MODEL to estimate values for a, b, and c from the data in SASHELP.CITIMON.
proc model data=sashelp.citimon; lhur = 1/(a * ip + b) + c; fit lhur; run;
Notice that the model equation is written as a SAS assignment statement. The variable LHUR is assumed to be the dependent variable because it is named in the FIT statement and is on the lefthand side of the assignment.
PROC MODEL determines that LHUR and IP are observed variables because they are in the input data set. A, B, and C are treated as unknown parameters to be estimated from the data because they are not in the input data set. If the data set contained a variable named A, B, or C, you would need to explicitly declare the parameters with a PARMS statement.
In response to the FIT statement, PROC MODEL estimates values for A, B, and C by using nonlinear least squares and prints the results. The first part of the output is a “Model Summary” table, shown in Figure 19.1.
Figure 19.1: Model Summary Report
Model Summary  

Model Variables  1 
Parameters  3 
Equations  1 
Number of Statements  1 
Model Variables  LHUR 

Parameters  a b c 
Equations  LHUR 
This table details the size of the model, including the number of programming statements that define the model, and lists the dependent variables (LHUR in this case), the unknown parameters (A, B, and C), and the model equations. In this case the equation is named for the dependent variable, LHUR.
PROC MODEL then prints a summary of the estimation problem, as shown in Figure 19.2.
Figure 19.2: Estimation Problem Report
The Equation to Estimate is  

LHUR =  F(a, b, c(1)) 
The notation used in the summary of the estimation problem indicates that LHUR is a function of A, B, and C, which are to be estimated by fitting the function to the data. If the partial derivative of the equation with respect to a parameter is a simple variable or constant, the derivative is shown in parentheses after the parameter name. In this case, the derivative with respect to the intercept C is 1. The derivatives with respect to A and B are complex expressions and so are not shown.
Next, PROC MODEL prints an estimation summary as shown in Figure 19.3.
Figure 19.3: Estimation Summary Report
Data Set Options  

DATA=  SASHELP.CITIMON 
Minimization Summary  

Parameters Estimated  3 
Method  Gauss 
Iterations  10 
Final Convergence Criteria  

R  0.000737 
PPC(b)  0.003943 
RPC(b)  0.00968 
Object  4.784E6 
Trace(S)  0.533325 
Objective Value  0.522214 
Observations Processed  

Read  145 
Solved  145 
Used  144 
Missing  1 
The estimation summary provides information on the iterative process used to compute the estimates. The heading “OLS Estimation Summary” indicates that the nonlinear ordinary least squares (OLS) estimation method is used. This table indicates that all three parameters were estimated successfully by using 144 nonmissing observations from the data set SASHELP.CITIMON. Calculating the estimates required 10 iterations of the GAUSS method. Various measures of how well the iterative process converged are also shown. For example, the “RPC(B)” value 0.00968 means that on the final iteration the largest relative change in any estimate was for parameter B, which changed by 0.968 percent. See the section Convergence Criteria for details.
PROC MODEL then prints the estimation results. The first part of this table is the summary of residual errors, shown in Figure 19.4.
Figure 19.4: Summary of Residual Errors Report
Nonlinear OLS Summary of Residual Errors  

Equation  DF Model  DF Error  SSE  MSE  Root MSE  RSquare  Adj RSq  Label 
LHUR  3  141  75.1989  0.5333  0.7303  0.7472  0.7436  UNEMPLOYMENT RATE: ALL WORKERS, 16 YEARS 
This table lists the sum of squared errors (SSE), the mean squared error (MSE), the root mean squared error (root MSE), and the R and adjusted R statistics. The R value of 0.7472 means that the estimated model explains approximately 75 percent more of the variability in LHUR than a mean model explains.
Following the summary of residual errors is the parameter estimates table, shown in Figure 19.5.
Figure 19.5: Parameter Estimates
Nonlinear OLS Parameter Estimates  

Parameter  Estimate  Approx Std Err  t Value  Approx Pr > t 
a  0.009046  0.00343  2.63  0.0094 
b  0.57059  0.2617  2.18  0.0309 
c  3.337151  0.7297  4.57  <.0001 
Because the model is nonlinear, the standard error of the estimate, the t value, and its significance level are only approximate. These values are computed using asymptotic formulas that are correct for large sample sizes but only approximately correct for smaller samples. Thus, you should use caution in interpreting these statistics for nonlinear models, especially for small sample sizes. For linear models, these results are exact and are the same as standard linear regression.
The last part of the output produced by the FIT statement is shown in Figure 19.6.
Figure 19.6: System Summary Statistics
Number of Observations  Statistics for System  

Used  144  Objective  0.5222 
Missing  1  Objective*N  75.1989 
This table lists the objective value for the estimation of the nonlinear system. Since there is only a single equation in this case, the objective value is the same as the residual MSE for LHUR except that the objective value does not include a degreesoffreedom correction. This can be seen in the fact that “Objective*N” equals the residual SSE, 75.1989. N is 144, the number of observations used.
Computing parameter estimates for nonlinear equations requires an iterative process. Starting with an initial guess for the parameter values, PROC MODEL tries different parameter values until the objective function of the estimation method is minimized. (The objective function of the estimation method is sometimes called the fitting function.) This process does not always succeed, and whether it does succeed depends greatly on the starting values used. By default, PROC MODEL uses the starting value 0.0001 for all parameters.
Consequently, in order to use PROC MODEL to achieve convergence of parameter estimates, you need to know two things: how to recognize convergence failure by interpreting diagnostic output, and how to specify reasonable starting values. The MODEL procedure includes alternate iterative techniques and grid search capabilities to aid in finding estimates. See the section Troubleshooting Convergence Problems for more details.